2b:null 2c:["$","title",null,{"children":"Covid-simulator | Bart Mesuere | Observable"}] 31:T561,{ const line = vl.markLine() .encode( vl.x().fieldT("DATE").title(null), vl.y().fieldQ("AVG_CASES").title("Gemiddeld aantal gevallen per dag").axis({grid: true}), vl.tooltip([ {field: "DATE", type: "temporal", title: 'Datum'}, {field: "AVG_CASES", type: "quantitative", title: 'Cases', format:".1f"}, ]) ); const projectedLine = vl.markLine().data(projected_national_cases) .encode( vl.x().fieldT("DATE").title(null), vl.y().fieldQ("AVG_CASES"), vl.strokeDash().value([3, 2]), vl.tooltip([ {field: "DATE", type: "temporal", title: 'Datum'}, {field: "AVG_CASES", type: "quantitative", title: 'Cases', format:".1f"}, ]) ); const dots = vl.markBar({opacity: 0.4, width: "2"}) .encode( vl.x().fieldT("DATE").title(null), vl.y().fieldQ("CASES"), vl.tooltip([ {field: "DATE", type: "temporal", title: 'Datum'}, {field: "CASES", type: "quantitative", title: 'Cases', format:".1f"}, ]) ); const rule = vl.markRule({stroke: "red", strokeWidth: 1}) .encode( vl.y().datum(800), ); return vl.layer(rule, dots, line, projectedLine) .data(national_cases.filter(d => new Date(d.DATE) >= new Date("2020-09-01"))) .title("Projectie aantal bevestigde gevallen per dag") .width(800) .render(); }32:T6fa,{ const bars = vl.markBar({width: 2}) .encode( vl.x().fieldT("DATE").title(null), vl.y().fieldQ("INCIDENCE").title("Incidentie").axis({labelPadding: 7}), vl.opacity().fieldN("PREDICTED").scale({range: [1, 0.4]}).legend(false), vl.color().fieldN("LEVEL").scale({domain: ["< 50", "< 200", "< 500", "500+"], range: ["#016351", "#ffbf00", "#c71e1d", "#66023C"]}).title("Incidentie (14 dagen)").legend({orient: "bottom"}), vl.tooltip([ {field: "DATE", type: "temporal", title: 'Datum'}, {field: "INCIDENCE", type: "quantitative", title: 'Incidentie', format:".1f"}, ]) ); const xmas = vl .markText({ text: "🎄", fontSize: 35, baseline: "bottom", align: "center", dx: 0.5, dy: 3 }).data([{DATE: "2020-12-25", VALUE: projected_incidence.find(d => d.DATE === "2020-12-25").INCIDENCE}]) .encode( vl.y().fieldQ("VALUE"), vl.x().fieldT("DATE") ); const val = vl .markText({ text: "❤️", fontSize: 35, baseline: "bottom", align: "center", dx: 0.5, dy: 3 }).data([{DATE: "2021-02-14", VALUE: projected_incidence.find(d => d.DATE === "2021-02-14").INCIDENCE}]) .encode( vl.y().fieldQ("VALUE"), vl.x().fieldT("DATE") ); const easter = vl .markText({ text: "🐣", fontSize: 35, baseline: "bottom", align: "center", dx: 0.5, dy: 3 }).data([{DATE: "2021-04-04", VALUE: projected_incidence.find(d => d.DATE === "2021-04-04").INCIDENCE}]) .encode( vl.y().fieldQ("VALUE"), vl.x().fieldT("DATE") ); return vl.layer(bars, xmas, val, easter) .data(projected_incidence.filter(d => new Date(d.DATE) >= new Date("2020-09-01"))) .title("Projectie incidentie") .width(800) .render(); }33:[] 2f:{"id":"9c205a1fa1868d50","slug":"corona-8-mei-simulator","title":"Covid-simulator","thumbnail":"19086eb25d62b69e86362e9c794a2c2ca6dea576a4481902ceec9696300169b3","access_level":"public","create_time":"$D2020-10-26T18:23:05.754Z","update_time":"$D2021-11-18T10:22:38.947Z","view_time":null,"first_public_time":"$D2020-10-26T19:57:29.361Z","fork":{"id":"dc34afa05963afb8","slug":"projectie-hospitalisaties","title":"Projectie hospitalisaties","owner":{"id":"1f20537af0cde79f","login":"bmesuere","name":"Bart Mesuere"}},"stars":0,"starred":false,"owner":{"type":"team","id":"1f20537af0cde79f","login":"bmesuere","name":"Bart Mesuere","avatar_url":"https://avatars-next.observableusercontent.com/eb2388e453b68d683131930d7e0b914d9c3d26dc21802c4e3c852e91ae6e1799","tier":"starter_2024"},"version":1982,"publish_version":1982,"role":null,"authored":false,"authors":[{"type":"individual","id":"d42baac23e3ffb5f","login":"bmesuere","name":"Bart Mesuere","avatar_url":"https://avatars-next.observableusercontent.com/eb2388e453b68d683131930d7e0b914d9c3d26dc21802c4e3c852e91ae6e1799","tier":null}],"members":[{"type":"individual","id":"d42baac23e3ffb5f","login":"bmesuere","name":"Bart Mesuere","avatar_url":"https://avatars-next.observableusercontent.com/eb2388e453b68d683131930d7e0b914d9c3d26dc21802c4e3c852e91ae6e1799","tier":null,"role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"subscribed":null,"forks":0,"collection":null,"type":"notebook","body":{"title":"Covid-simulator","cells":[{"id":0,"value":"md`# Covid-simulator\n\n**Dit is geen wetenschappelijk model of voorspelling van de toekomst. Het is enkel een tool om snel het effect van verschillende groeipercentages in te schatten.**\n\nDeze pagina simuleert het gemiddeld aantal bevestigde gevallen tot het einde van het jaar aan de hand van enkele zelf gekozen parameters:\n* met de eerste parameter kan je de uiteindelijke groei- of afnamesnelheid kiezen\n* met de tweede slider kan je kiezen binnen hoeveel dagen we die groei- of afnamesnelheid behalen`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1711,"value":"md`Gekozen groeisnelheid: ${getRateInWords(final_rate)}`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1686,"value":"viewof final_rate = html`\n`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1708,"value":"md`Gekozen vertraging: **${days_till_rate}** dagen`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1688,"value":"viewof days_till_rate = html``","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1315,"value":"$31","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1742,"value":"$32","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1188,"value":"md`**Berekening van de onderliggende waarden**`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":7,"value":"raw_cases = (await fetch('https://epistat.sciensano.be/Data/COVID19BE_CASES_AGESEX.json?date=${today}')).json()","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1309,"value":"national_cases = {\n let result = [];\n const grouped = d3a.rollup(raw_cases, v => {\n return {\n CASES : d3.sum(v, d => d.CASES)\n };\n }, d => d.DATE);\n dates.forEach(date => {\n const v = grouped.get(date);\n result.push({DATE: date, CASES: v.CASES});\n });\n for (let i = 6; i < result.length; i++) {\n result[i].AVG_CASES = (result[i - 6].CASES + result[i - 5].CASES + result[i - 4].CASES + result[i - 3].CASES + result[i - 2].CASES + result[i - 1].CASES + result[i].CASES) / 7;\n }\n for (let i = 13; i < result.length; i++) {\n result[i].INCIDENCE = 100000 * (7*result[i - 7].AVG_CASES + 7*result[i].AVG_CASES) / 11492641;\n }\n result = result.filter(d => new Date(d.DATE) <= new Date(d_4))\n return result;\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1329,"value":"days_to_project = Math.floor((new Date(\"2021-11-01\") - Date.now())/(60*60*24*1000)) + 5 + 7 + 2","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1323,"value":"avg_cases_rate = national_cases.slice(-1)[0].AVG_CASES / national_cases.slice(-8)[0].AVG_CASES","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1332,"value":"projected_national_cases = {\n const result = [];\n let previous = national_cases.slice(-1)[0];\n result.push(previous);\n for (let i = 0; i < days_to_project; i++) {\n const date = d3.timeFormat(\"%Y-%m-%d\")(d3.timeDay.offset(Date.parse(d_3), i));\n const avg_cases = previous.AVG_CASES * Math.pow(getRate(i), 1/7);\n previous = {DATE: date, AVG_CASES: avg_cases, PREDICTED : true};\n result.push(previous);\n }\n return result;\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1746,"value":"projected_incidence = {\n const result = national_cases.concat(projected_national_cases.slice(1));\n for (let i = 0; i < 13; i++) {\n result[i].INCIDENCE = 0;\n }\n for (let i = 13; i < result.length; i++) {\n result[i].INCIDENCE = 100000 * (7*result[i - 7].AVG_CASES + 7*result[i].AVG_CASES) / 11492641;\n if (result[i].INCIDENCE < 50) result[i].LEVEL = \"< 50\";\n else if (result[i].INCIDENCE < 200) result[i].LEVEL = \"< 200\";\n else if (result[i].INCIDENCE < 500) result[i].LEVEL = \"< 500\";\n else result[i].LEVEL = \"500+\";\n }\n return result;\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1693,"value":"function getRate(day) {\n if (day >= days_till_rate) return final_rate;\n return avg_cases_rate - (day * (avg_cases_rate - final_rate) / days_till_rate);\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1714,"value":"function getRateInWords(rate) {\n if (rate < 1) {\n return `**afname** van **${Math.round((1-rate)*100)}%** per week`;\n }\n if (rate > 1) {\n return `**toename** van **${Math.round((rate-1)*100)}%** per week`;\n }\n return \"constante snelheid\";\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":928,"value":"md`**Other stuff**`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":607,"value":"today = d3.timeFormat(\"%Y-%m-%d\")(d3.timeHour.offset(Date.now(), -4))","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1280,"value":"d_1 = d3.timeFormat(\"%Y-%m-%d\")(d3.timeDay.offset(Date.parse(today), -1))","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1668,"value":"d_3 = d3.timeFormat(\"%Y-%m-%d\")(d3.timeDay.offset(Date.parse(today), -3))","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1647,"value":"d_4 = d3.timeFormat(\"%Y-%m-%d\")(d3.timeDay.offset(Date.parse(today), -4))","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":75,"value":"dates = Array.from(new Set(raw_cases.map(c => c.DATE))).filter(d => d)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":78,"value":"provinces = Array.from(new Set(raw_cases.map(c => c.PROVINCE)))","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":380,"value":"extendedProvinces = new Array(...provinces, \"Vlaanderen\", \"Wallonië\", \"België\")","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":20,"value":"d3 = require(\"d3@5\")","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":82,"value":"d3a = require(\"d3-array@^2.2\")","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":679,"value":"import {vl} from '@vega/vega-lite-api'","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"}],"files":"$W33","stdlib":"1","titling":"auto"}} 30:"$2f:body"