31:Tb72,chart = { const height = 400; const svg = d3.create("svg") .attr("viewBox", [0, 0, width, height]); const data = runningFraction; const samples = pastSamples; const indices = pastIndices; const pastData = d3.group(samples.map((sample,i) => { return sample; }), (d,i) => indices[i]); const max_index = d3.max(pastData.keys()); const x = d3.scaleLinear() .domain(d3.extent(data, (d, i) => i)) .range([margin.left, width - margin.right]) const meanRate = d3.mean(rates); const y = d3.scaleLog() .domain(d3.extent([...samples, ...data.map((d,i) => Math.pow(Math.E, meanRate*i))], (d, i) => d)).nice() .range([height - margin.bottom, margin.top]); const xAxis = g => g .attr("transform", `translate(0,${height - margin.bottom})`) .call(d3.axisBottom(x).ticks(width / 80).tickSizeOuter(0)); const yAxis = g => g .attr("transform", `translate(${margin.left},0)`) .call(d3.axisLeft(y)) .call(g => g.select(".domain").remove()) .call(g => g.select(".tick:last-of-type text").clone() .attr("x", 3) .attr("text-anchor", "start") .attr("font-weight", "bold") .text("Individual wealth trajectory (log scale)")); const line = d3.line() .defined(d => !isNaN(d)) .x((d,i) => x(i)) .y(d => y(d)) const ensembleAvg = d3.line() .defined(d => !isNaN(d)) .x((d,i) => x(i)) .y((d,i) => y(Math.pow(Math.E, i*meanRate))) const baseLine = d3.line() .defined(d => !isNaN(d)) .x((d,i) => x(i)) .y((d,i) => y(1)) svg.append("g") .call(xAxis); svg.append("g") .call(yAxis); // svg.append("path") // .datum(data) // .attr("fill", "none") // .attr("stroke", "red") // .attr("stroke-width", 5) // .attr("stroke-linejoin", "round") // .attr("stroke-linecap", "round") // .attr("d", line); svg.append("path") .datum(data) .attr("fill", "none") .attr("stroke", "grey") .attr("stroke-width", 1) .attr("stroke-linejoin", "round") .attr("stroke-linecap", "round") .attr("d", baseLine); svg.append("path") .datum(data) .attr("fill", "none") .attr("stroke", "blue") .attr("stroke-width", 1) .attr("stroke-linejoin", "round") .attr("stroke-linecap", "round") .attr("d", ensembleAvg); //console.log(pastData); Array.from(pastData.entries()).forEach(entry => { const v = entry[1]; const k = entry[0]; //console.log(k,v); const line_color = (k == max_index) ? "red": "lightgrey"; const stroke_width = (k == max_index) ? 5: 1; svg.append("path") .datum(v) .attr("fill", "none") .attr("stroke", line_color) .attr("stroke-width", stroke_width) .attr("stroke-linejoin", "round") .attr("stroke-linecap", "round") .attr("d", line); }); return svg.node(); }32:T47a,runningFraction = { //let samples = []; resample; let wealth = 1.0; let temporalWealths = [wealth]; let temporalAverages = [wealth]; const chunk_size = 10; const rng = d3.randomInt(rates.length); while(true){ // For improved drawing add 100 iterations at a time before yield for(let i=0; i < chunk_size; i++){ const randIdx = rng(); const rate = rates[randIdx]; wealth = wealth * (1.0 + rate); //samples.push(sample); temporalWealths.push(wealth); const newAvg = d3.mean(temporalWealths); temporalAverages.push(newAvg); //console.log(randIdx, rate, wealth, newAvg); //yield newAvg; mutable pastSamples.push(wealth); mutable pastIndices.push(parseInt(resample)); } mutable pastSamples = mutable pastSamples; mutable pastIndices = mutable pastIndices; // mutable pastSamples = [...pastSamples, fraction]; // mutable pastIndices = [...pastIndices, parseInt(resample)]; if(noLoop || temporalWealths.length > 1000) { yield Promises.delay(1, temporalWealths); break; }; yield Promises.delay(1, temporalWealths); } }33:[] 2f:{"title":"Multiplicative coin toss: Ergodicity breaking","cells":[{"id":322,"value":"md`# Multiplicative coin toss: Ergodicity breaking`","pinned":true,"hidden":false,"mode":"ojs"},{"id":286,"value":"viewof form_data = form(html``)","pinned":false,"hidden":false,"mode":"ojs"},{"id":85,"value":"viewof resample = {\n const btn = html`