31:T5e1,function phaseChartInc(incdata, opts_param) { const options = { hist_lookback: 0, hist_len: 7, ref_offset: 0, smooth: 2, smooth_growth: 2, postproc: (data => data), filter: (a => true), ...opts_param }; const data = incdata.sort((a,b) => b.population - a.population).map(reg => { let srs = reg.series.map(a => a.c7/7); if (options.hist_lookback > 0) { srs = srs.slice(0, -options.hist_lookback); } for (let ii = 1; ii <= options.smooth; ii++) { srs = srs.map((a, idx, arr) => idx idx { const aPr = srsGrowth[idx-1]; const aN = idx < srsGrowth.length-1 ? srsGrowth[idx+1] : srsGrowth[idx]*srsGrowth[idx]/aPr; return { x: a * pop_norm, y: index_from_ratio(Math.pow(aN / aPr, 7/2)) }; }).slice(-options.hist_len-1); const refIdx = series.length - 1 - options.ref_offset; return { name: reg.name, region: reg.region, x: series[refIdx].x, y: series[refIdx].y, series: series, r: /*reg.population //*/series[refIdx].x * reg.population / pop_norm }; }); return phaseChart(options.postproc(data).filter(options.filter), options); }32:T17f4,function phaseChart(data, opts_param) { const minX = d3.min(data.map(a => d3.min(a.series.map(s => s.x).filter(a => a > 0)))); const maxX = d3.max(data.map(a => d3.max(a.series.map(s => s.x).filter(a => a >= 0)))); const minY = d3.min(data.map(a => d3.min(a.series.map(s => s.y).filter(a => a>=-1)))); const maxY = d3.max(data.map(a => d3.max(a.series.map(s => s.y).filter(a => a<=100)))); const maxR = d3.max(data.map(a => a.r)); //return [minX, maxX, minY, maxY]; const options = { bndsX: [Math.max(0.01, minX * 0.95), Math.min(1000, maxX * 1.2)], bndsY: [Math.max(0, (minY+1)*0.95)-1, Math.min(6, (maxY+1)*1.1)-1], w: width, h: width*9/16, scaleR: 1, clampX: true, clampY: true, ...opts_param }; const margin = ({top: 10 + (options.title ? 10 : 0), right: 5, bottom: 35, left: 45}); const nTicksX = options.w / 40; const smoothW = Math.round(7 + Math.sqrt(options.smooth)); const scaleX = d3.scaleLog(options.bndsX, [margin.left, options.w - margin.right]); scaleX.clamp(options.clampX); const scaleXAxis = d3.axisBottom(scaleX).ticks(nTicksX, ",").tickFormat(d3.format('.2s')); const xAxis = g => g .attr("font-family", "'IBM Plex Sans',sans-serif") .attr("transform", `translate(0,${options.h - margin.bottom})`) .call(scaleXAxis) .call(g => g.select(".domain").remove()) .call(g => g.append("text") .attr("x", options.w) .attr("y", margin.bottom - 4) .attr("fill", "currentColor") .attr("text-anchor", "end") .text("Primeri na " + (pop_norm/1000) + "k preb. (" + smoothW + "d glajenje)→")); const yTicks = indexTicks.filter(a => inRange(a, options.bndsY)); const scaleY = d3.scaleLinear(options.bndsY, [options.h - margin.bottom, margin.top]).domain(options.bndsY); scaleY.clamp(options.clampY); const yAxis = g => g .attr("transform", `translate(${margin.left},0)`) .call(d3.axisLeft(scaleY).tickValues(yTicks).tickFormat(idx => d3.format('.00%')( ratio_from_index(idx) - 1 ))) .call(g => g.select(".domain").remove()) .call(g => g.append("text") .attr("x", -margin.left) .attr("y", 7) .attr("fill", "currentColor") .attr("text-anchor", "start") .text("↑ Tedenska rast (" + smoothW + "d glajenje)")); const radius = d3.scaleSqrt([0, 314], [0, options.w / 20 * options.scaleR]); const grid = g => g .attr("stroke", "black") .attr("stroke-width", 0.05) .call(g => g.append("g") .selectAll("line") .data(scaleX.ticks(nTicksX)) .join("line") .attr("x1", d => 0.5 + scaleX(d)) .attr("x2", d => 0.5 + scaleX(d)) .attr("y1", margin.top) .attr("y2", options.h - margin.bottom)) .call(g => g.append("g") .selectAll("line") .data(yTicks) .join("line") .attr("y1", d => 0.5 + scaleY(d)) .attr("y2", d => 0.5 + scaleY(d)) .attr("x1", margin.left) .attr("x2", options.w - margin.right)); const svg = d3.create("svg") .attr("viewBox", [0, 0, options.w, options.h]) .attr("font-family", "'IBM Plex Sans',sans-serif"); svg.append("g").append("line") // attach a line .style("stroke", "black") // colour the line .style("stroke-width", 0.2) .attr("x1", scaleX(0)) // x position of the first end of the line .attr("y1", scaleY(0)) // y position of the first end of the line .attr("x2", scaleX(100000)) // x position of the second end of the line .attr("y2", scaleY(0)); svg.append("g") .selectAll("rect") .data([ {x1: 0.01, x2: 14.5, c: 'green'}, {x1: 14.5, x2: 29, c: 'yellow'}, {x1: 29, x2: 48, c: 'orange'}, {x1: 48, x2: 65, c: 'red'}, {x1: 65, x2: 1000, c: 'darkgray'}]) .join("rect") .attr("opacity",0.2) .attr("x", d => scaleX(d.x1)) .attr("y", options.h - margin.bottom) .attr("width", d => scaleX(d.x2) - scaleX(d.x1)) .attr("height", 20) .attr("fill", d => d.c); svg.append("g") .call(xAxis); svg.append("g") .call(yAxis); svg.append("g") .call(grid); const line = d3.line() .defined(d => d && (d.x || isFinite(d.x)) && (d.y || isFinite(d.y))) .x((d,i) => scaleX(d.x || 0) || 0) .y((d,i) => scaleY(d.y || 1) || 0); svg.append("g") .selectAll("path") .data(data) .join("path") .attr("fill", "none") .attr("stroke", d => reg_colors[d.region]) .attr("stroke-width", d => radius(d.r)/8) .attr("stroke-linecap", "round") .attr("stroke-linejoin", "round") .attr("opacity", '0.4') //.style("mix-blend-mode", "multiply") .attr("d", d => line(d.series)); svg.append("g") .attr("stroke", "none") .selectAll("circle") .data(data) .join("circle") .attr("fill", d => reg_colors[d.region]) .attr("opacity", "0.65") .attr("cx", d => scaleX(d.x)) .attr("cy", d => scaleY(d.y)) .attr("r", d => radius(d.r)) .append("title") .text(d => d.name); if (options.labels) { svg.append("g") .attr("font-size", 10) .selectAll("text") .data(data) .join("text") .attr("fill", "black") .attr("dy", "0.35em") .attr("x", d => scaleX(d.x) + 7) .attr("y", d => scaleY(d.y) - radius(d.r)) .text(d => d.name); } svg.append("g") .append("text") .attr("font-size", 10) .attr("text-anchor", "end") .attr("x", options.w) .attr("y", margin.top) .text(strFromDate(d3.timeDay.offset(last_data_update, -options.hist_lookback))); const title_text = (typeof options.title === 'string' || options.title instanceof String) ? options.title : options.title(data); if (title_text) { svg.append("g") .attr("font-size", 20) .attr("font-weight", "bold") .attr("text-anchor", "middle") .append("text") .attr("x", 0.5 * (options.w + margin.left - margin.right)) .attr("y", margin.top-2) .text(title_text); } return svg.node(); }33:T89d,regions_table = function(data_arr, title, sortFun = (() => 0)) { const form = html`

${data_arr .sort(sortFun) .map(m => m) .map((m, index) => '' ).join('')}

#	Ime	1d	Aktivni primeri	Inc. [na 100k]	Pr.7d [7d smooth]	Pr.7d [${short_offset}d smooth]
' + (index+1) + '	' + m.name + '	' + m.cases_1d + ' [' + Math.round(m.incidence1d * 1400000) + ']' + '	' + inc_eqgrow_str(m.cases_prevWk, m.cases_lastWk) + '	' + inc_eqgrow_str(m.incidencePrevWk * 100000, m.incidenceLastWk * 100000) + '	' + inc_eqgrow_str(m.incidenceLastWk * 100000, m.incidenceNextWk * 100000) + '	' + inc_eqgrow_str(m.incidenceLastWk * 100000, m.incidenceWkNextS*100000) + '

`; const toggleContent = function(event) { event.target.classList.toggle("active"); var content = event.target.nextElementSibling; if (content.style.display === "block") { content.style.display = "none"; } else { content.style.display = "block"; } }; form.button.onclick = toggleContent; return form; }34:T8e2,function produce_series(cases_todate, population, dates) { const cases_1d = cases_todate.map((c, idx) => c - (cases_todate[idx-1] || 0)); const cases_7dg = cases_1d.map((c, idx, arr) => 7 * geomavg(cases_1d, idx-6, idx)); const cases_7d = cases_1d.map((c, idx, arr) => 7 * average(cases_1d, idx-6, idx)); const cases_14d = cases_7d.map((c, idx, arr) => 14 * average(cases_1d, idx-13, idx)); const cases_21d = cases_14d.map((c, idx, arr) => 21 * average(cases_1d, idx-20, idx)); const growth_series_1d = cases_21d.map((c, idx) => Math.pow(c / cases_21d[idx-1], 7) || 1); const growth_series_7d = cases_7d.map((c, idx) => (c / cases_7d[idx-7]) || 1); const growth_series_14d = cases_14d.map((c, idx) => Math.pow(c / cases_14d[idx-1], 7) || 1); const growth_series_s = cases_7d.map((c, idx) => Math.pow(c / cases_7d[idx - short_offset], 7/short_offset) || 1); const lastIdx = cases_todate.length - 1; const casesLastWk = cases_7d[lastIdx]; const casesLastWkS = cases_7d[lastIdx - short_offset] || 0; const casesPrevWk = cases_7d[lastIdx - 7] || 0; const growth = (casesLastWk / casesPrevWk) || 0; const growthS = (casesLastWk / casesLastWkS) || 0; const growthWkS = Math.pow(growthS, 7/short_offset) || 0; const casesNextWk = casesLastWk * growth; const casesNextWkS = casesLastWk * growthWkS; return { cases_active: cases_14d[lastIdx] || 0, cases_total: cases_todate[lastIdx] || 0, cases_total7: cases_todate[lastIdx - 7] || 0, cases_total14: cases_todate[lastIdx - 14] || 0, cases_lastWk: casesLastWk, cases_lastWkS: casesLastWkS, cases_prevWk: casesPrevWk, cases_1d: cases_1d[lastIdx], incidence1d: cases_1d[lastIdx] / population, incidence14d: cases_14d[lastIdx] / population, incidenceLastWk: casesLastWk / population, incidencePrevWk: casesPrevWk / population, growth: growth, incidenceNextWk: casesNextWk / population, incidenceLastWkS: casesLastWkS / population, growthS: growthS, growthWkS: growthWkS, incidenceWkNextS: casesNextWkS / population, series: cases_7d.map((c,idx) => ({ c7: cases_7d[idx]/population, date: dates[idx] })) }; }35:[] 2f:{"title":"Municipalities - sledilnik.org","cells":[{"id":0,"value":"md`# Municipalities - sledilnik.org`","pinned":false,"hidden":false,"mode":"ojs"},{"id":197,"value":"md`## Tabele`","pinned":false,"hidden":false,"mode":"ojs"},{"id":80,"value":"munis_table('Urejeno po incidenci 14d', (a,b) => b.incidence14d - a.incidence14d)","pinned":false,"hidden":false,"mode":"ojs"},{"id":89,"value":"munis_table('Urejeno po incidenci 7d', (a,b) => b.incidenceLastWk - a.incidenceLastWk)","pinned":false,"hidden":false,"mode":"ojs"},{"id":91,"value":"munis_table('Urejeno po hitrosti naraščanja 7d incidence', (a,b) => b.growth - a.growth)","pinned":false,"hidden":false,"mode":"ojs"},{"id":113,"value":"munis_table('Urejeno po projekciji 7d-inc. čez 1 teden, s 7d razmerjem',(a,b) => b.incidenceNextWk - a.incidenceNextWk)","pinned":false,"hidden":false,"mode":"ojs"},{"id":318,"value":"munis_table('Urejeno po projekciji čez 1 teden s ' + short_offset + 'd razmerjem',(a,b) => b.incidenceWkNextS - a.incidenceWkNextS)","pinned":false,"hidden":false,"mode":"ojs"},{"id":375,"value":"munis_table('Urejeno po 1d incidenci',(a,b) => b.incidence1d - a.incidence1d)","pinned":false,"hidden":false,"mode":"ojs"},{"id":267,"value":"regions_table(regions_incidence, 'REGIJE (urejene po projekciji 14d inc.)', (a,b) => \n b.incidenceLastWk + b.incidenceWkNextS - a.incidenceLastWk - a.incidenceWkNextS)","pinned":true,"hidden":false,"mode":"ojs"},{"id":409,"value":"regions_incidence","pinned":true,"hidden":false,"mode":"ojs"},{"id":413,"value":"d3.csvFormat(regions_incidence)","pinned":true,"hidden":false,"mode":"ojs"},{"id":392,"value":"{\n const dt = d3.timeFormat('%Y-%m-%d')(last(muni_data.timeSeries).date);\n \n return html`DOWNLOAD REGIONS CSV ${dt}`\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":231,"value":"md`## Diagram`","pinned":false,"hidden":false,"mode":"ojs"},{"id":661,"value":"pop_norm = 100000;","pinned":true,"hidden":false,"mode":"ojs"},{"id":847,"value":"$31","pinned":true,"hidden":false,"mode":"ojs"},{"id":923,"value":"function index_from_ratio(ratio) {\n const idx = (ratio - 1)/(ratio + 1);\n return idx;//Math.pow(Math.abs(idx), 1/2) * Math.sign(idx);\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":934,"value":"function ratio_from_index(index) {\n const idx = index;// Math.pow(index, 2);\n return (1 + idx) / (1 - idx);\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":931,"value":"indexTicks=[0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1,1.1,1.2,1.3,1.4,1.5,1.6,1.75,2.0,2.5,3, 3.5, 6, 11,61,Number.POSITIVE_INFINITY].map(index_from_ratio);","pinned":true,"hidden":false,"mode":"ojs"},{"id":471,"value":"$32","pinned":true,"hidden":false,"mode":"ojs"},{"id":213,"value":"d3.timeDay.range(d3.isoParse('2020-03-04'), new Date())","pinned":true,"hidden":false,"mode":"ojs"},{"id":1449,"value":"regions_diagram = phaseChartInc(\n regions_incidence,\n {title: 'Regije', labels: true, h: width*9/16, hist_len: 14, smooth: 21, hist_lookback: date_back, ref_offset: 8, scaleR: 2, bndsX: [0.3, 250], bndsY: [-1, 1]})","pinned":true,"hidden":false,"mode":"ojs"},{"id":1691,"value":"regions = phaseChartInc(regions_incidence/*.filter(r => 'mb,ce,sg,ms'.includes(r.id))*/,\n {title: \"Statistične regije\", labels: true, hist_len: 21, scaleR: 0.5, smooth: 100, ref_offset: 7, smooth_growth: 0, hist_lookback: date_back, bndsX: [10, 160], bndsY: [-0.5,0.5]})","pinned":false,"hidden":false,"mode":"ojs"},{"id":191,"value":"viewof date = Scrubber(d3.timeDay.range(d3.isoParse('2020-09-01'), d3.isoParse('2021-03-12'/*last_data_update*/)), {\n autoplay: false,\n loop: false,\n delay: 100,\n initial: 0,\n format: date => strFromDate(date)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1581,"value":"munis_diagram = phaseChartInc(\n muni_incidence.filter(m => true/*5000 <= m.population*/),\n {\n title: 'Občine', hist_len: 7, smooth: 3, ref_offset: 3,\n scaleR: 4, bndsX:[0.1, 100], bndsY: [-1, 1]\n });","pinned":true,"hidden":false,"mode":"ojs"},{"id":1468,"value":"d3.isoParse(last_data_update)","pinned":true,"hidden":false,"mode":"ojs"},{"id":1451,"value":"date_back=d3.timeDay.count(date, last_data_update);","pinned":true,"hidden":false,"mode":"ojs"},{"id":1264,"value":"reg_colors = ({\n lj: 'rgb( 69, 120, 68)',\n mb: 'rgb(249, 93, 106)',\n ce: 'rgb(102, 81, 145)',\n kr: 'rgb(255, 166, 0)',\n nm: 'rgb(175, 165, 63)',\n ng: 'rgb(219, 165, 29)',\n kp: 'rgb(112, 164, 113)',\n ms: 'rgb( 2, 74, 102)',\n kk: 'rgb(212, 80, 135)',\n za: 'rgb( 16, 130, 154)',\n sg: 'rgb(119, 124, 41)',\n po: 'rgb(160, 81, 149)'\n});","pinned":true,"hidden":false,"mode":"ojs"},{"id":104,"value":"md`## CODE`","pinned":false,"hidden":false,"mode":"ojs"},{"id":453,"value":"short_offset = 3;","pinned":false,"hidden":false,"mode":"ojs"},{"id":297,"value":"regions_table(muni_incidence.slice(1,5),\"AA\", ((a,b) => b.incidenceWkNext4 - a.incidenceWkNext4))","pinned":true,"hidden":false,"mode":"ojs"},{"id":247,"value":"munis_table = function(title, sortFun) {\n return regions_table(muni_incidence, title, sortFun);\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":43,"value":"$33","pinned":false,"hidden":false,"mode":"ojs"},{"id":306,"value":"function inc_eq_str(prev, last) {\n const frmt = d3.format('.2d');\n return frmt(+prev) + ' + ' + frmt(+last) + ' = ' + frmt((last + prev) || 0)\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":309,"value":"function inc_eqgrow_str(prev, last) {\n return inc_eq_str(prev, last) + ' ' + grow_str2(prev, last);\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":333,"value":"function grow_str2(prev, last) {\n return grow_str(last/prev);\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":349,"value":"function grow_str(ratio) {\n return '(' + (+(ratio > 1.0) ? '+' : '') + d3.format('.0%')(+ratio - 1) + ')';\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":269,"value":"d3.sum(regions_incidence.map(ri => ri.population))","pinned":true,"hidden":false,"mode":"ojs"},{"id":243,"value":"regions_incidence = d3.groups(muni_incidence, m => m.region).map(reg => {\n const obj = {\n ...regions_dict[reg[0]],\n population_munis: 0,\n cases_todate: []\n };\n reg[1].forEach(ob => {\n ob.cases_todate.forEach((c, idx) => {\n if (obj.cases_todate.length <= idx) {\n obj.cases_todate.push(0);\n }\n obj.cases_todate[idx] += (c || 0);\n })\n obj.population_munis += ob.population;\n });\n const series_data = produce_series(obj.cases_todate, obj.population, reg[1][0].series.map(s => s.date));\n return {...obj, ...series_data, region: obj.id};\n});","pinned":true,"hidden":false,"mode":"ojs"},{"id":72,"value":"muni_incidence.filter(a => a.incidence14d>0.02)","pinned":true,"hidden":false,"mode":"ojs"},{"id":25,"value":"muni_incidence = {\n const muni_objs = muni_data.dict.map(d => {\n const cases_todate = muni_data.timeSeries.map(dt => dt['region.'+d.region+'.'+d.id+'.cases.confirmed.todate']);\n const series_data = produce_series(cases_todate, d.population, muni_data.timeSeries.map(a => a.date));\n return {\n ...d,\n ...series_data,\n cases_todate: cases_todate\n };\n });\n return muni_objs.filter(m => m.id != 'slovenija');\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":628,"value":"$34","pinned":false,"hidden":false,"mode":"ojs"},{"id":717,"value":"function geomavg(arr, startIdx, endIdx) {\n let cum = 1;\n let cnt = 0;\n for (let idx = startIdx; idx <= endIdx; idx++) {\n if (arr[idx]) {\n cum = cum * arr[idx];\n cnt++;\n }\n }\n return Math.pow(cum, 1/cnt);\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":981,"value":"function average(arr, startIdx, endIdx) {\n let cum = 0;\n for (let idx = startIdx; idx <= endIdx; idx++) {\n if (arr[idx]) {\n cum = cum + arr[idx];\n }\n }\n return cum/(endIdx - startIdx + 1);\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":832,"value":"function geomgauss3(arr, idx) {\n const weights = [1, 2, 1];\n const wOff = 1;\n let cum = 1;\n let cnt = 0;\n weights.forEach((w, ii) => {\n const aIdx = idx + ii - wOff;\n if (arr[aIdx]) {\n cum = cum * Math.pow(arr[aIdx], w);\n cnt += w;\n }\n });\n return Math.pow(cum, 1/cnt); \n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":396,"value":"function last(arr) {\n return arr[arr.length-1];\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":600,"value":"function strFromDate(d) {\n return d3.timeFormat('%Y\\-%m\\-%d')(d);\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1433,"value":"function inRange(value, range) {\n return range[0] <= value && value <= range[1];\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":398,"value":"last_data_update = last(muni_data.timeSeries).date","pinned":false,"hidden":false,"mode":"ojs"},{"id":446,"value":"regions_dict = d3.csv(\"https://raw.githubusercontent.com/sledilnik/data/master/csv/dict-region.csv\", d3.autoType).then(a => {\n const ret = {};\n a.forEach(reg => {\n ret[reg.id] = reg;\n });\n return ret;\n});","pinned":false,"hidden":false,"mode":"ojs"},{"id":8,"value":"muni_data = Promise.all([\n d3.csv(\"https://raw.githubusercontent.com/sledilnik/data/master/csv/municipality-cases.csv\", d3.autoType),\n d3.csv(\"https://raw.githubusercontent.com/sledilnik/data/master/csv/dict-municipality.csv\", d3.autoType)\n]).then(([a,b]) => ({timeSeries: a, dict: b}));","pinned":false,"hidden":false,"mode":"ojs"},{"id":18,"value":"si_data = d3.csv(\"https://raw.githubusercontent.com/sledilnik/data/master/csv/stats.csv\", d3.autoType)","pinned":false,"hidden":false,"mode":"ojs"},{"id":22,"value":"regression = require('https://bundle.run/regression@2.0.1');","pinned":false,"hidden":false,"mode":"ojs"},{"id":193,"value":"import {Scrubber} from \"@mbostock/scrubber\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":1393,"value":"style = html`\n\n`","pinned":false,"hidden":false,"mode":"ojs"},{"id":5,"value":"d3 = require('d3','d3-time')","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W35","stdlib":"1","titling":"auto"}