2a:null 2b:["$","title",null,{"children":"Boston-Area SARS-CoV-2 Wastewater Levels | Harris Lapiroff | Observable"}] 30:T1722,chart = { // Filter data to only be within the window selected in the brush chart const data_ = data.filter(d => ( +d.date > +brush.startDate && +d.date < +brush.endDate && systems.includes(d.system) )) const milestones_ = ( accompanyingData.includes('Annotations') ? milestones.filter(d => +d.date > +brush.startDate && +d.date < +brush.endDate) : [] ) // Returns a function to update a tooltip using data from the specified system const updateTooltip = (event, datum) => { const date = d3.utcDay.floor(datum.x) const fetchedData = data_.filter(d => +d.date === +date) event.currentTarget._tippy.setContent( d3.utcFormat('%-m/%d/%Y')(date) + '
' + fetchedData.map(d => `${capitalize(d.system)}: ${d3.format(',')(d.concentrationAvg)} cp/mL `).join('
') ) } // Render a plot const plot = Plot.plot({ y: { grid: true, type: scaleType, label: 'Copies/mL ↑', tickFormat: 's', domain: staticYAxis ? [ scaleType === 'log' ? d3.min(data, d => d.concentration) : 0, d3.max(data, d => d.ciHigh), ] : [ scaleType === 'log' ? d3.min(data_, d => d.concentration) : 0, d3.max(data_, d => d.ciHigh), ] }, color: { legend: true, domain: ['northern', 'southern'], tickFormat: t => `${capitalize(t)} System`, }, width: chartWidth, height: chartWidth / 1.6, marginRight: accompanyingData.includes('Annotations') ? 100 : 10, marks: [ // Milestones Marks (the ticks go up here to go behind everything else, but the text is at the bottom to go above) Plot.tickX(milestones_, { x: 'date', stroke: '#999', strokeDasharray: '4,2' }), // Daily Measurements accompanyingData.includes('Daily Measurements') ? Plot.dot(data_, { x: 'date', y: 'concentration', fill: 'system', opacity: 0.33, r: 2, }) : null, // Intervals accompanyingData.includes('Intervals') ? Plot.link(data_, { x1: 'date', x2: 'date', y1: 'ciLow', y2: 'ciHigh', strokeWidth: 3, stroke: 'system', opacity: 0.1, strokeLinecap: 'round', }) : null, // Average Intervals accompanyingData.includes('Average Intervals') ? Plot.areaY(data_,{ x: 'date', y1: 'avgHigh', y2: 'avgLow', fill: 'system', opacity: 0.15, }) : null, // 7-Day Average Line on( Plot.line(data_, { x: 'date', y: 'concentrationAvg', stroke: 'system', strokeWidth: 4 }), { // This updates the tooltip (added below) with inspector data for the point being hovered over pointermove: updateTooltip, pointerenter: updateTooltip, pointerdown: updateTooltip, } ), accompanyingData.includes('Annotations') ? [ Plot.dot( data_, Plot.selectMaxY({ x: 'date', y: 'concentrationAvg', fill: 'system', r: 5, stroke: '#FFF' }) ), Plot.text( data_, Plot.selectMaxY({ x: 'date', y: 'concentrationAvg', z: 'system', fill: '#000', strokeWidth: 3, stroke: '#FFF', text: (d, i, data) => { const sameSysData = data.filter(d_ => d.system === d_.system) return d3.format(',d')(d.concentrationAvg) + ' cp/Ml\n' + formatRatio(d.concentrationAvg / d3.min(sameSysData, d => d.concentrationAvg)) + ' min' }, dy: -12, textAnchor: 'start' }) ), // Min Plot.dot(data_, Plot.selectMinY({ x: 'date', y: 'concentrationAvg', fill: 'system', r: 5, stroke: '#FFF' })), Plot.text( data_, Plot.selectMinY({ x: 'date', y: 'concentrationAvg', z: 'system', fill: '#000', strokeWidth: 3, stroke: '#FFF', text: d => d3.format(',d')(d.concentrationAvg) + ' cp/Ml\n', dy: 12, textAnchor: 'start' }) ), // Curr Plot.dot(data_, Plot.selectMaxX({ x: 'date', y: 'concentrationAvg', fill: 'system', r: 5, stroke: '#FFF' })), Plot.text( data_, Plot.selectMaxX({ x: 'date', y: 'concentrationAvg', z: 'system', fill: '#000', strokeWidth: 3, stroke: '#FFF', text: d => { const sameSysData = data_.filter(d_ => d.system === d_.system) return `${formatRatio(d.concentrationAvg / d3.max(sameSysData, d => d.concentrationAvg))} peak` + '\n' + `${formatRatio(d.concentrationAvg / d3.min(sameSysData, d => d.concentrationAvg))} min` }, dx: 8, textAnchor: 'start' }) ), ] : null, // Milestone Text Plot.text(milestones_, { x: 'date', textAnchor: 'start', text: 'name', dx: -4, dy: 0, rotate: 90, frameAnchor: 'top', stroke: '#FFF', fill: '#333', strokeWidth: 4, }), ] }) // Change the paint order on text so the stroke shows up below the fill d3.select(plot).selectAll('text').attr('paint-order', 'stroke') // Make lines have crips edges d3.select(plot).selectAll('line').attr('shape-rendering', 'crispEdges') // Add tooltips showing inspector to the plot const tooltip = tippy(plot.querySelectorAll('[aria-label="line"] path'), { followCursor: true, allowHTML: true, theme: 'light', animation: 'fade', delay: [0, 200], }) return plot }31:T49c,viewof brush = { const data_ = data.filter(d => d.concentration) const plot = Plot.plot({ marginLeft: 0, marginRight: 0, marginTop: 0, y: { type: 'log', ticks: false }, x: { grid: true, ticks: 10 }, height: 80, width: width, marks: [ Plot.dot(data_, { x: 'date', y: 'concentration', fill: 'system', opacity: 0.4, r: 1.5 }), Plot.line(data_, { x: 'date', y: 'concentrationAvg', stroke: 'system', strokeWidth: 1.5 }), Plot.frame({ stroke: '#CCC' }) ] }) const xScale = plot.scale('x') const yScale = plot.scale('y') const brushed = ({ selection }) => { plot.value = { startDate: xScale.invert(selection[0]), endDate: xScale.invert(selection[1]) } plot.dispatchEvent(new Event('input'), { bubbles: true }) } const brush = d3.brushX() .extent([[xScale.range[0], yScale.range[1]], [xScale.range[1] + 1, yScale.range[0]]]) .on("start brush end", brushed) d3.select(plot) .append('g') .call(brush) .call( brush.move, [d3.utcDay.offset(d3.max(data_, d => d.date), -269), d3.utcDay.offset(d3.max(data_, d => d.date), 1)].map(xScale.apply) ) return plot }32:T463,mapAndLegend = htl.html`

${map.node()}

${Array.from(SYSTEMS_BY_TOWN).map(([system, towns]) => htl.html.fragment`

${capitalize(system)}

${t}

`)}

`33:T433,{ const data_ = data.filter(d => d.concentration !== null).sort((a, b) => a.date - b.date) return Plot.plot({ height: d3.utcYear.count( d3.utcYear.floor(d3.min(data_, d => d.date)), d3.utcYear.floor(d3.max(data_, d => d.date)) ) * 150, width: width, y: { type: scaleType, tickFormat: 's', grid: true, ticks: 3, domain: [100, 10000] }, x: { tickFormat: t => d3.timeFormat('%b %-d')(dateFromPct(t)) }, fy: { tickFormat: d => d.toFixed() }, facet: {data: data_, y: d => d3.utcDay.offset(d.date, 1).getFullYear() }, marks: [ Plot.link(data_, { x1: d => pctThroughYear(d.date), x2: d => pctThroughYear(d.date), y1: 'ciLow', y2: 'ciHigh', strokeWidth: 3, stroke: 'system', opacity: 0.25, strokeLinecap: 'round', }), Plot.line( data_, { x: d => pctThroughYear(d.date), y: 'concentrationAvg', stroke: 'system', strokeWidth: 3, } ), Plot.frame({ stroke: '#CCC' }), ], }) }34:T790,data = (await d3.csv(DATA_URL)).map(d => ({ northernVariant: parseVariantResults(d['Northern\rVariant\rResults']), southernVariant: parseVariantResults(d['Southern\rVariant Results']), ...d, date: d3.utcParse('%m/%d/%Y')(d['Sample Date']), southern: toIntOrNull(d['Southern\r(copies/mL)']), northern: toIntOrNull(d['Northern\r(copies/mL)']), southernAvg: toIntOrNull(d['Southern\r7 day avg']), northernAvg: toIntOrNull(d['Northern\r7 day avg']), })) .map(d => ({ ...d, southernLow: d.southern - toIntOrNull(d['Southern\rLow Confidence\rInterval']), northernLow: d.northern - toIntOrNull(d['Northern\rLow Confidence\rInterval']), southernHigh: d.southern + toIntOrNull(d['Southern\rHigh Confidence\rInterval']), northernHigh: d.northern + toIntOrNull(d['Northern\rHigh Confidence\rInterval']), })) .map((d, i, arr) => { const windowStart = Math.max(0, i - 7) const windowEnd = i const window = arr.slice(windowStart, windowEnd) return { ...d, southernLowAvg: d3.mean(window.filter(d => d.southernLow), d => d.southernLow), northernLowAvg: d3.mean(window.filter(d => d.northernLow), d => d.northernLow), southernHighAvg: d3.mean(window.filter(d => d.southernHigh), d => d.southernHigh), northernHighAvg: d3.mean(window.filter(d => d.northernHigh), d => d.northernHigh), } }) // Unfold the data so each system gets its own datapoints .map((d, i, arr) => ( [ { system: 'northern', date: d.date, concentration: d.northern, concentrationAvg: d.northernAvg, ciHigh: d.northernHigh, ciLow: d.northernLow, avgHigh: d.northernHighAvg, avgLow: d.northernLowAvg, }, { system: 'southern', date: d.date, concentration: d.southern, concentrationAvg: d.southernAvg, ciHigh: d.southernHigh, ciLow: d.southernLow, avgHigh: d.southernHighAvg, avgLow: d.southernLowAvg, } ] )) .flat()35:[{"id":"92d05ffcec5489aee036f3f348863d7af491525f94ca0044908c41b7da89d1ab766fc2ab1ad9c6082035d633a33bbe2cde1259209af166270d8e557729f25c38","href":"https://static.observableusercontent.com/files/92d05ffcec5489aee036f3f348863d7af491525f94ca0044908c41b7da89d1ab766fc2ab1ad9c6082035d633a33bbe2cde1259209af166270d8e557729f25c38","name":"mwra-data-20220219.csv","type":"text/csv","lastModified":1645545752881.248,"size":32069},{"id":"e2abb2fc4d10235a0b037a5fc7544b0db35fbc59d9970176fbb2cc1a809cda04d0c8821c1202961017d3a4cb639fc5d0c80f4dc98341916f373a5196b56f75c1","href":"https://static.observableusercontent.com/files/e2abb2fc4d10235a0b037a5fc7544b0db35fbc59d9970176fbb2cc1a809cda04d0c8821c1202961017d3a4cb639fc5d0c80f4dc98341916f373a5196b56f75c1","name":"mwra-data-20220217.csv","type":"text/csv","lastModified":1645164124246.283,"size":32031},{"id":"bbe7c93899b9c3ba14929fb55b89869e2f8c54972d2e6fbcd98d6998d06d4203c04f00406df9cfed73f93f4f2d04a5a9ef30b413d824a38aabad4fb9faa52dc4","href":"https://static.observableusercontent.com/files/bbe7c93899b9c3ba14929fb55b89869e2f8c54972d2e6fbcd98d6998d06d4203c04f00406df9cfed73f93f4f2d04a5a9ef30b413d824a38aabad4fb9faa52dc4","name":"ma-towns-state-plane-topo-1.json","type":"application/json","lastModified":1645685010969.58,"size":9168288},{"id":"92238e15c3f972fb181c0b9c59377f74ed21c302af222802db158fa410b2b86d3d775fccc41115c03db267ca24f76344df0f911e0ab51fbaa4948e470ad8df99","href":"https://static.observableusercontent.com/files/92238e15c3f972fb181c0b9c59377f74ed21c302af222802db158fa410b2b86d3d775fccc41115c03db267ca24f76344df0f911e0ab51fbaa4948e470ad8df99","name":"mwra-data-20220406.csv","type":"text/csv","lastModified":1649303307228.267,"size":35090},{"id":"97b0c8a2d5ad1cb464020b31bf5b9579f617f28aec6bdca1b88044d64175650aaf00cbbc18de68b3e64310365d4c952c8d5653944ef26f49150e35daf59189b9","href":"https://static.observableusercontent.com/files/97b0c8a2d5ad1cb464020b31bf5b9579f617f28aec6bdca1b88044d64175650aaf00cbbc18de68b3e64310365d4c952c8d5653944ef26f49150e35daf59189b9","name":"mwra-data-20220409.csv","type":"text/csv","lastModified":1649551000482.724,"size":35162},{"id":"f81dfc9a842398620d29eb9b77ba1732480e9bd1e467306dce829dca3ffe7390f1498fa28edf2498c97d6a9ba64c991adef4661f9fbd444d75764ff22b3731bd","href":"https://static.observableusercontent.com/files/f81dfc9a842398620d29eb9b77ba1732480e9bd1e467306dce829dca3ffe7390f1498fa28edf2498c97d6a9ba64c991adef4661f9fbd444d75764ff22b3731bd","name":"mwra-data-20220412.csv","type":"text/csv","lastModified":1649792663136.412,"size":35264},{"id":"84d0b70c118217af0ec6cb2f5d2595a739ca72676a5297520d5ed8f51d06abdc90ae47fc53bf0d73f12d9242e220c1a10bb24f0e9de94ace3f5321cc6d22e4dc","href":"https://static.observableusercontent.com/files/84d0b70c118217af0ec6cb2f5d2595a739ca72676a5297520d5ed8f51d06abdc90ae47fc53bf0d73f12d9242e220c1a10bb24f0e9de94ace3f5321cc6d22e4dc","name":"mwra-data-20220423.csv","type":"text/csv","lastModified":1650767714175.356,"size":35499},{"id":"b684a16af949b04764c1314aea6a6b470ccaeb254f0b71b2c3a69d76d87d8e76851e86d85d4b416a4147687a453b57556bc0d3b2b7339e91db2bf6c428fb0898","href":"https://static.observableusercontent.com/files/b684a16af949b04764c1314aea6a6b470ccaeb254f0b71b2c3a69d76d87d8e76851e86d85d4b416a4147687a453b57556bc0d3b2b7339e91db2bf6c428fb0898","name":"mwra-data-20220428.csv","type":"text/csv","lastModified":1651194910129.995,"size":35620},{"id":"b54fe595538ec7ac4cf4ecd64fb916e02adc81e4fb476717ef7518f2c7e7013668dd80d766d38edcbf00ec6e808ff4d67f1c19b73af2462f81e3aa52992af007","href":"https://static.observableusercontent.com/files/b54fe595538ec7ac4cf4ecd64fb916e02adc81e4fb476717ef7518f2c7e7013668dd80d766d38edcbf00ec6e808ff4d67f1c19b73af2462f81e3aa52992af007","name":"mwra-data-20220429.csv","type":"text/csv","lastModified":1651263356384.741,"size":35674},{"id":"c8fff05837f9f53fffd04b07ffab2da8f13629c722ae5915feaab3824748c7bc8e8b2499299f32cab426e734e2759442f58a1414fe9876e002d7fcbcc89fba67","href":"https://static.observableusercontent.com/files/c8fff05837f9f53fffd04b07ffab2da8f13629c722ae5915feaab3824748c7bc8e8b2499299f32cab426e734e2759442f58a1414fe9876e002d7fcbcc89fba67","name":"mwra-data-20220503.csv","type":"text/csv","lastModified":1651611264372.26,"size":36418},{"id":"6bddeb02eefb450a040750dd8b34e3f2ad226253fc6b4ca9b0b8aca18b5710d9b8714fc820eecc8631ebb7aa395d034166018901d2e1ffabc6e04388410b2842","href":"https://static.observableusercontent.com/files/6bddeb02eefb450a040750dd8b34e3f2ad226253fc6b4ca9b0b8aca18b5710d9b8714fc820eecc8631ebb7aa395d034166018901d2e1ffabc6e04388410b2842","name":"mwra-data-20220504.csv","type":"text/csv","lastModified":1651721998550.249,"size":36442},{"id":"a60de0ae01822f475e917bef6859b594543212f0aca6b4aa7c31bfdd52e85236dac54686107f2c14c0e9ac6b81e6ffeee986a32f6db87839100c2402b9d23664","href":"https://static.observableusercontent.com/files/a60de0ae01822f475e917bef6859b594543212f0aca6b4aa7c31bfdd52e85236dac54686107f2c14c0e9ac6b81e6ffeee986a32f6db87839100c2402b9d23664","name":"mwra-data-20220505.csv","type":"text/csv","lastModified":1651854220741.358,"size":36466},{"id":"e8379d7b1016dbd463e8563684c4789811db67bc9a0e36c93c345d03bdc59d4cc05f226d1a65e3af6c1682d5741aa9cf3d6c5aff415348d76c3b16dc216df0a5","href":"https://static.observableusercontent.com/files/e8379d7b1016dbd463e8563684c4789811db67bc9a0e36c93c345d03bdc59d4cc05f226d1a65e3af6c1682d5741aa9cf3d6c5aff415348d76c3b16dc216df0a5","name":"mwra-data-20220506.csv","type":"text/csv","lastModified":1651862190200.68,"size":36490},{"id":"e0ee8914fb55ec0c4d29aae68f0c77c62ddac45ccb4690f21910c683b86427675223c4ed383a9aa3f279706d0e73251315d52faad449037c6256dea33be2b340","href":"https://static.observableusercontent.com/files/e0ee8914fb55ec0c4d29aae68f0c77c62ddac45ccb4690f21910c683b86427675223c4ed383a9aa3f279706d0e73251315d52faad449037c6256dea33be2b340","name":"tabula-MWRAData20220215-data.csv","type":"text/csv","lastModified":1644953930626.159,"size":31539},{"id":"52f0958fd06a974bf4fa28f991084ec98bfa8d245af687b40797ed6b9723acd006ab068c53ba50385ce8481869cd81cdbd41db3ed3eda2890b7eca99be2aeebb","href":"https://static.observableusercontent.com/files/52f0958fd06a974bf4fa28f991084ec98bfa8d245af687b40797ed6b9723acd006ab068c53ba50385ce8481869cd81cdbd41db3ed3eda2890b7eca99be2aeebb","name":"mwra-data-20220224.csv","type":"text/csv","lastModified":1645733644930.919,"size":32261},{"id":"2186d712d789c9df86209493d663eac484417e620f968cbd67267fbca18aba3fd2b0de4a0b066f75d1cf33a3fed4964539a9cde4fb98c55998113e167f4ecf6e","href":"https://static.observableusercontent.com/files/2186d712d789c9df86209493d663eac484417e620f968cbd67267fbca18aba3fd2b0de4a0b066f75d1cf33a3fed4964539a9cde4fb98c55998113e167f4ecf6e","name":"mwra-data-20220308.csv","type":"text/csv","lastModified":1646758727634.238,"size":33983},{"id":"50ecbc6a7b5a0306bed5291ddefbe72ca90982f355e03675f531f5351cc9288d7848151c3dccc6f0cb184a548c505672a63cb5745121381b9341e4493385567b","href":"https://static.observableusercontent.com/files/50ecbc6a7b5a0306bed5291ddefbe72ca90982f355e03675f531f5351cc9288d7848151c3dccc6f0cb184a548c505672a63cb5745121381b9341e4493385567b","name":"mwra-data-20220309.csv","type":"text/csv","lastModified":1646880520688.07,"size":34003},{"id":"a8688ff2cd880827b0e74a154d6d05c7f6266d16b1a6be7c64e84b6146b026153e396604550b618b5b902a71ed4def50f762bf927e6d666da55eea94eb5e79e5","href":"https://static.observableusercontent.com/files/a8688ff2cd880827b0e74a154d6d05c7f6266d16b1a6be7c64e84b6146b026153e396604550b618b5b902a71ed4def50f762bf927e6d666da55eea94eb5e79e5","name":"mwra-data-20220305.csv","type":"text/csv","lastModified":1646666968635.175,"size":32432},{"id":"9ed6c73539725df2b10e8b9392f393e19f65152e8b9cab3f171fb4ac869f43970971afe228a2acfb658a6193578a578f547fede3374507943adb740f815537ce","href":"https://static.observableusercontent.com/files/9ed6c73539725df2b10e8b9392f393e19f65152e8b9cab3f171fb4ac869f43970971afe228a2acfb658a6193578a578f547fede3374507943adb740f815537ce","name":"mwra-data-20220218.csv","type":"text/csv","lastModified":1645247003026.441,"size":32050},{"id":"beb32b39d0514545911e6444c33101b21d07d104b958956be4fa65f0ce42b55fee891ed7cb4b810538e43bf307f9ffb97dc4c08dbb7be895dcc2f4a298671d6e","href":"https://static.observableusercontent.com/files/beb32b39d0514545911e6444c33101b21d07d104b958956be4fa65f0ce42b55fee891ed7cb4b810538e43bf307f9ffb97dc4c08dbb7be895dcc2f4a298671d6e","name":"mwra-data-20220301.csv","type":"text/csv","lastModified":1646167656175.737,"size":32357},{"id":"b12e92b06a3d8d565a2f4249b0148e29c4a9986966b68671c272f780dfd8a8342fde2fac00af94ba26c6a25db2520c3230d495b87ca2becaec141c17251825e7","href":"https://static.observableusercontent.com/files/b12e92b06a3d8d565a2f4249b0148e29c4a9986966b68671c272f780dfd8a8342fde2fac00af94ba26c6a25db2520c3230d495b87ca2becaec141c17251825e7","name":"mwra-data-20220223.csv","type":"text/csv","lastModified":1645658901158.232,"size":32242},{"id":"c0702062525aa49c8a0604f0ea0ce9614e918bf8888ff9d2a06c3f821f15437898047e58497d184c082fbcbe3535ef487ca8f019edc78ab1e38efc2751e1acee","href":"https://static.observableusercontent.com/files/c0702062525aa49c8a0604f0ea0ce9614e918bf8888ff9d2a06c3f821f15437898047e58497d184c082fbcbe3535ef487ca8f019edc78ab1e38efc2751e1acee","name":"mwra-data-20220303.csv","type":"text/csv","lastModified":1646333893766.801,"size":32393},{"id":"e97e405d4263f2218e032e8cdd487b3c41821a7baf079733868e10bd803552bc2a618539c15da26a3bf5c6004a36674e75a330aa5c6f796f47d09d50445030c1","href":"https://static.observableusercontent.com/files/e97e405d4263f2218e032e8cdd487b3c41821a7baf079733868e10bd803552bc2a618539c15da26a3bf5c6004a36674e75a330aa5c6f796f47d09d50445030c1","name":"tabula-MWRAData20220216-data.csv","type":"text/csv","lastModified":1645045515210.063,"size":31558},{"id":"4e37e6eb857b5c7f3d290c8f0011e616ed3d44d0699b64877d8dcd8bd41655502957609aa2ddd6a77ee93eda7cb46caf9c5f7788a848f8b77f9b3fbfa8bf3745","href":"https://static.observableusercontent.com/files/4e37e6eb857b5c7f3d290c8f0011e616ed3d44d0699b64877d8dcd8bd41655502957609aa2ddd6a77ee93eda7cb46caf9c5f7788a848f8b77f9b3fbfa8bf3745","name":"mwra-data-20220225.csv","type":"text/csv","lastModified":1645827231020.692,"size":32280},{"id":"5b0c6d3fc1d265c2d743791e1752b1bb8421cf8115d9424123f7bb8c2b11709b60c519deff16796792168266d3048839a045c6f5de616dd6d16ea4ca36c8a384","href":"https://static.observableusercontent.com/files/5b0c6d3fc1d265c2d743791e1752b1bb8421cf8115d9424123f7bb8c2b11709b60c519deff16796792168266d3048839a045c6f5de616dd6d16ea4ca36c8a384","name":"mwra-data-20220310.csv","type":"text/csv","lastModified":1646976632654.355,"size":33993},{"id":"c97994602352c96be424e03fccc9f932622fc7cf05ef55df5279c105e00826fa1570e4a47fe528418460b214c573d872b9c9e4ab04020ad19840df6761cf3390","href":"https://static.observableusercontent.com/files/c97994602352c96be424e03fccc9f932622fc7cf05ef55df5279c105e00826fa1570e4a47fe528418460b214c573d872b9c9e4ab04020ad19840df6761cf3390","name":"mwra-data-20220311.csv","type":"text/csv","lastModified":1647027485918.796,"size":34042},{"id":"a87ea5181927784b70c48c8a7bc220ae87be7cc6547bc265d8471f4e6b96d50d6cb773f64e0e12164ad05fdd493ce55435a66c79ec260c430b76a2b11a8cd310","href":"https://static.observableusercontent.com/files/a87ea5181927784b70c48c8a7bc220ae87be7cc6547bc265d8471f4e6b96d50d6cb773f64e0e12164ad05fdd493ce55435a66c79ec260c430b76a2b11a8cd310","name":"mwra-data-20220313.csv","type":"text/csv","lastModified":1647262338801.172,"size":34032},{"id":"0eac11a0ee586ae483b8deff090d5749b18b733c8a8baf30a0320b98636efa1bfd0fb8e96b2ee9e99d808e9d2b4860cb0ac2eb5edf71c58eef3d2a12153d9dc2","href":"https://static.observableusercontent.com/files/0eac11a0ee586ae483b8deff090d5749b18b733c8a8baf30a0320b98636efa1bfd0fb8e96b2ee9e99d808e9d2b4860cb0ac2eb5edf71c58eef3d2a12153d9dc2","name":"mwra-data-20220315.csv","type":"text/csv","lastModified":1647440793019.589,"size":34122},{"id":"4b03113a8152ee68b51656542fa28e5b4617ae9d868ba163723e35cc27715e636db958f4fd6c246f7aad3497f3640b3f02c07c5c21eeeefa53eba5f15f0ecf2c","href":"https://static.observableusercontent.com/files/4b03113a8152ee68b51656542fa28e5b4617ae9d868ba163723e35cc27715e636db958f4fd6c246f7aad3497f3640b3f02c07c5c21eeeefa53eba5f15f0ecf2c","name":"mwra-data-20220316.csv","type":"text/csv","lastModified":1647456967122.495,"size":34142},{"id":"7df2dd6ee5d14cbc95dda608863d58e69ff4fe1c0c3a8f57df95bec894adb4bee5c3001332b25dd56d950b7d20170915db1486eca7d8b6f8f14c84952fce0b25","href":"https://static.observableusercontent.com/files/7df2dd6ee5d14cbc95dda608863d58e69ff4fe1c0c3a8f57df95bec894adb4bee5c3001332b25dd56d950b7d20170915db1486eca7d8b6f8f14c84952fce0b25","name":"mwra-data-20220317.csv","type":"text/csv","lastModified":1647578556939.451,"size":34162},{"id":"824421c6fb8db73e7c70cb943e6c8b8e0d80968285fb8de900a51e68579b484830fca5f97934ca789cbf87cd401769ef5f6d375b401942ce235b6478a476fac6","href":"https://static.observableusercontent.com/files/824421c6fb8db73e7c70cb943e6c8b8e0d80968285fb8de900a51e68579b484830fca5f97934ca789cbf87cd401769ef5f6d375b401942ce235b6478a476fac6","name":"mwra-data-20220319.csv","type":"text/csv","lastModified":1647733659117.483,"size":34734},{"id":"bc58ab0432599fba3f130dbfd1ba02634598d06bab545819d1634b8faaa97727baf13a8c1b78e3d016eaf6958554858e0ffd3e893cd969cc7742cd1ddcfd0e43","href":"https://static.observableusercontent.com/files/bc58ab0432599fba3f130dbfd1ba02634598d06bab545819d1634b8faaa97727baf13a8c1b78e3d016eaf6958554858e0ffd3e893cd969cc7742cd1ddcfd0e43","name":"mwra-data-20220322.csv","type":"text/csv","lastModified":1647966465339.52,"size":34763},{"id":"192c5085799069d311381c801a330dbe1594e3d8a5f56f64dfeb491d175ecb1c0a6fa22a2921075c31374a7edc99cad33def134fac25ff699ff6cc0805a47120","href":"https://static.observableusercontent.com/files/192c5085799069d311381c801a330dbe1594e3d8a5f56f64dfeb491d175ecb1c0a6fa22a2921075c31374a7edc99cad33def134fac25ff699ff6cc0805a47120","name":"tabula-MWRAData20220323-data.csv","type":"text/csv","lastModified":1648094841784.903,"size":33423},{"id":"db49ff347364956620da9960a034e4f184c6f5031eded69c8e2e1493aaa7dd54c54432bda4402216eab4e2e863f56c84d72cd1eb6d1de6fbdb5e5b8728cd6a43","href":"https://static.observableusercontent.com/files/db49ff347364956620da9960a034e4f184c6f5031eded69c8e2e1493aaa7dd54c54432bda4402216eab4e2e863f56c84d72cd1eb6d1de6fbdb5e5b8728cd6a43","name":"tabula-MWRAData20220324-data.csv","type":"text/csv","lastModified":1648170247008.694,"size":34346},{"id":"8f89d3e9575b6d5a56f608e50fafb5c1bb387d1ad9a36637913a9142bfded9ee76da9fcc1c8a55089b0bbc66b41f916286fa1cdbdd2b8796bd405acafd3f0ef1","href":"https://static.observableusercontent.com/files/8f89d3e9575b6d5a56f608e50fafb5c1bb387d1ad9a36637913a9142bfded9ee76da9fcc1c8a55089b0bbc66b41f916286fa1cdbdd2b8796bd405acafd3f0ef1","name":"tabula-MWRAData20220325-data.csv","type":"text/csv","lastModified":1648261953001.017,"size":34367},{"id":"a9e33c7626a28af193f8bf57478dacb949c363c52621d6ffae099a0da2e2ca1ea61921f843016ea0d9c9612c5dcfc99298c3e3aceee88f7d94b83888f463db67","href":"https://static.observableusercontent.com/files/a9e33c7626a28af193f8bf57478dacb949c363c52621d6ffae099a0da2e2ca1ea61921f843016ea0d9c9612c5dcfc99298c3e3aceee88f7d94b83888f463db67","name":"mwra-data-20220328.csv","type":"text/csv","lastModified":1648527749746.564,"size":34877},{"id":"d348ee07da8f73f2f106a42fec294d1d9c0b1084b5aa76969dcaf03f0c6b8ef650dbf85113efb737c426d19530809687c823c6a042b6688f09f89143f067ef79","href":"https://static.observableusercontent.com/files/d348ee07da8f73f2f106a42fec294d1d9c0b1084b5aa76969dcaf03f0c6b8ef650dbf85113efb737c426d19530809687c823c6a042b6688f09f89143f067ef79","name":"mwra-data-20220329.csv","type":"text/csv","lastModified":1648569422113.86,"size":34898},{"id":"2f7e6d7c27cf80be7b1e84e7ae8370c9d2ef06bc840cf4e94cdf7db566ffa36cf93345ec0c6cb4bac0c0e1f2fcbe650cc3f186f22600b1b15ed6e3c84af094c5","href":"https://static.observableusercontent.com/files/2f7e6d7c27cf80be7b1e84e7ae8370c9d2ef06bc840cf4e94cdf7db566ffa36cf93345ec0c6cb4bac0c0e1f2fcbe650cc3f186f22600b1b15ed6e3c84af094c5","name":"mwra-data-20220329a.csv","type":"text/csv","lastModified":1648650725241.43,"size":34909},{"id":"7d15862316554d7f81a96f05922883d07f6460358e1ff767422c0ad156e3139c372790c9526d6a2907a1ed1d829dd4697776555b265b780a6cd1c95372916bc7","href":"https://static.observableusercontent.com/files/7d15862316554d7f81a96f05922883d07f6460358e1ff767422c0ad156e3139c372790c9526d6a2907a1ed1d829dd4697776555b265b780a6cd1c95372916bc7","name":"mwra-data-20220330.csv","type":"text/csv","lastModified":1648679495233.964,"size":34930},{"id":"d9d591288e8a6eb177273190125b42259575373962c508d3af7672d22a8162596c5adc94163d8ac052fdb39781b025cf247f5bddbfead157f2a8dfa1478a3c7f","href":"https://static.observableusercontent.com/files/d9d591288e8a6eb177273190125b42259575373962c508d3af7672d22a8162596c5adc94163d8ac052fdb39781b025cf247f5bddbfead157f2a8dfa1478a3c7f","name":"mwra-data-20220331.csv","type":"text/csv","lastModified":1648756026245.627,"size":34982},{"id":"718c218da8b2ddf2d9fc527b4f698dbc70fdb2c5f3ca4c31cff0cab33f6693518436ae2ffeface75f0cf8edbfce0dcb691c7c8233d293bc7e6276bc4a21520d4","href":"https://static.observableusercontent.com/files/718c218da8b2ddf2d9fc527b4f698dbc70fdb2c5f3ca4c31cff0cab33f6693518436ae2ffeface75f0cf8edbfce0dcb691c7c8233d293bc7e6276bc4a21520d4","name":"mwra-data-20220402.csv","type":"text/csv","lastModified":1648930474506.57,"size":34998},{"id":"d29faa2b97a855668fa35140595ca73131fc87ddfa89e913dede7b850dc0f3d09fe879c52d21b72411701fa062c0454be0c2900408c96d31a3f276bfdcce01c5","href":"https://static.observableusercontent.com/files/d29faa2b97a855668fa35140595ca73131fc87ddfa89e913dede7b850dc0f3d09fe879c52d21b72411701fa062c0454be0c2900408c96d31a3f276bfdcce01c5","name":"mwra-data-20220405.csv","type":"text/csv","lastModified":1649215007901.695,"size":35097},{"id":"0624a438c7c26d32b04e5dd7be4413c9ab282730d9256821c94ebabb560f5a9521b4a9643cace905c1b68d668731aa652b647dff4ac574a96cd9f2458ea6117a","href":"https://static.observableusercontent.com/files/0624a438c7c26d32b04e5dd7be4413c9ab282730d9256821c94ebabb560f5a9521b4a9643cace905c1b68d668731aa652b647dff4ac574a96cd9f2458ea6117a","name":"mwra-data-20220413.csv","type":"text/csv","lastModified":1649947240714.031,"size":35288},{"id":"6635e16fd99dbe1b90ee6627a023d18b9042ed92c291f85159726f58c37b11a8e07b068b8814b3080ec5dbfcbace92889b90c27c6f9bf1bb5116e36b75006f7f","href":"https://static.observableusercontent.com/files/6635e16fd99dbe1b90ee6627a023d18b9042ed92c291f85159726f58c37b11a8e07b068b8814b3080ec5dbfcbace92889b90c27c6f9bf1bb5116e36b75006f7f","name":"mwra-data-20220415.csv","type":"text/csv","lastModified":1650068877702.393,"size":35336},{"id":"e6342ecf15da5dc531f62be24852bf16aa2c7c275a888eee2ecafdfb5fac44019cb23e4e2d6c8ea21b57b292e0e20a297b5335fbd40eb89591f4e220e206cee0","href":"https://static.observableusercontent.com/files/e6342ecf15da5dc531f62be24852bf16aa2c7c275a888eee2ecafdfb5fac44019cb23e4e2d6c8ea21b57b292e0e20a297b5335fbd40eb89591f4e220e206cee0","name":"mwra-data-20220420.csv","type":"text/csv","lastModified":1650498521667.054,"size":35457},{"id":"b5a9fd3224431aacc558b072a17b0db457080d8267827c2b210d839a946fa1636510946817f60a33e5b6bdce651c5158aebb08d69d95fa88c604e2d104a6dc17","href":"https://static.observableusercontent.com/files/b5a9fd3224431aacc558b072a17b0db457080d8267827c2b210d839a946fa1636510946817f60a33e5b6bdce651c5158aebb08d69d95fa88c604e2d104a6dc17","name":"mwra-data-20220422.csv","type":"text/csv","lastModified":1650680388660.341,"size":35505},{"id":"49091f8b8fe1b9b295b279782fe749d7028b71acc5476853c37af86d2756515f9ed45b0cd25295f973b90aa371ee567d1e9dd3f575c603c4088abb58167822df","href":"https://static.observableusercontent.com/files/49091f8b8fe1b9b295b279782fe749d7028b71acc5476853c37af86d2756515f9ed45b0cd25295f973b90aa371ee567d1e9dd3f575c603c4088abb58167822df","name":"mwra-data-20220426.csv","type":"text/csv","lastModified":1650991160491.393,"size":35602},{"id":"f62d5460d27f4146255e3d945b7afdc81a30f7cea460d6f9b11c0adb82a5dbcd93f2aed2c4e6fb96b3cfd0b9795ac1b53830006b3413a94e95005416847cd228","href":"https://static.observableusercontent.com/files/f62d5460d27f4146255e3d945b7afdc81a30f7cea460d6f9b11c0adb82a5dbcd93f2aed2c4e6fb96b3cfd0b9795ac1b53830006b3413a94e95005416847cd228","name":"mwra-data-20220427.csv","type":"text/csv","lastModified":1651100011264.608,"size":35626},{"id":"7c6ea253c7b39aba8217c70adfceb32ea8ab83e42a04b064b5f1789641368e5f2e3aa48ef48db7003ed68dd2628e7e7a1ecaf169e56991f642e9b9e4a1a76c99","href":"https://static.observableusercontent.com/files/7c6ea253c7b39aba8217c70adfceb32ea8ab83e42a04b064b5f1789641368e5f2e3aa48ef48db7003ed68dd2628e7e7a1ecaf169e56991f642e9b9e4a1a76c99","name":"mwra-data-20220430.csv","type":"text/csv","lastModified":1651336574200.293,"size":36316}] 2e:{"id":"d213529454eff018","slug":"mwra-wastewater-data","title":"Boston-Area SARS-CoV-2 Wastewater Levels","thumbnail":"6e2a49cb05422fb5c938cc9ba70f6eab37a925ab0f66ff7a66c27d5c9e658074","access_level":"public","create_time":"$D2022-02-15T19:38:38.302Z","update_time":"$D2023-12-04T19:29:38.656Z","view_time":null,"first_public_time":"$D2022-02-15T20:27:47.364Z","has_unpublished_edits":null,"fork":false,"stars":13,"starred":false,"owner":{"type":"team","id":"3587ecfe3e28cf64","login":"harrislapiroff","name":"Harris Lapiroff","avatar_url":"https://avatars-next.observableusercontent.com/9777c782bd50f7d54c128b6ac9313f8731fb712666730d9e662539eb3f5a5f07","tier":"starter_2024"},"version":1857,"publish_version":null,"role":null,"authored":false,"authors":[{"type":"individual","id":"666e3ebbff9b41a2","login":"harrislapiroff","name":"Harris Lapiroff","avatar_url":"https://avatars-next.observableusercontent.com/9777c782bd50f7d54c128b6ac9313f8731fb712666730d9e662539eb3f5a5f07","tier":null}],"members":[{"type":"individual","id":"666e3ebbff9b41a2","login":"harrislapiroff","name":"Harris Lapiroff","avatar_url":"https://avatars-next.observableusercontent.com/9777c782bd50f7d54c128b6ac9313f8731fb712666730d9e662539eb3f5a5f07","tier":null,"role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"lapsed":null,"subscribed":null,"forks":0,"collection":{"id":"2f3f46dbe149ed9b","slug":"coronavirus","title":"Coronavirus","owner":{"id":"3587ecfe3e28cf64","login":"harrislapiroff","name":"Harris Lapiroff"}},"type":"notebook","body":{"titling":"auto","title":"Boston-Area SARS-CoV-2 Wastewater Levels","cells":[{"id":0,"value":"# Boston-Area SARS-CoV-2 Wastewater Levels","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":425,"value":"**Latest reading:**\n${d3.utcDay.count(latestDate(data), d3.utcDay.floor(new Date()))}\n${pluralize('day', d3.utcDay.count(latestDate(data), d3.utcDay.floor(new Date())))} ago,\n${d3.utcFormat('%-m/%-d/%Y')(latestDate(data))}\n${\n latestSiteUpdate ?\n htl.html.fragment`\n (${+latestDate(data) >= +latestSiteUpdate ?\n 'Up to date' :\n `Outdated, latest website data is from ${d3.utcFormat('%-m/%-d/%Y')(latestSiteUpdate)}`\n })\n ` : \n ''\n}","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":415,"value":"> Samples are currently taken 3 - 7 times a week and analyzed by [Biobot Analytics](https://biobot.io), a wastewater epidemiology company. Data are updated as received by Biobot following an internal review process.\n>\n> — [Massachusetts Water Resources Authority](https://www.mwra.com/biobot/biobotdata.htm)","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":93,"value":"viewof scaleType = Inputs.radio(['log', 'linear'], { label: 'Scale Type', value: 'log' })","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":189,"value":"viewof systems = Inputs.checkbox([\"northern\", \"southern\"], { label: \"Systems\", value: ['northern', 'southern'] })","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":1178,"value":"viewof accompanyingData = Inputs.checkbox(\n ['Daily Measurements', 'Intervals', 'Average Intervals', 'Annotations'],\n { label: 'Data', value: ['Daily Measurements', 'Intervals', 'Annotations'] }\n)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":1210,"value":"Lines represent a 7-day average and are always shown","pinned":false,"hidden":false,"mode":"html","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":640,"value":"viewof staticYAxis = Inputs.toggle({label: \"Preserve y-axis over all time\"})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":10,"value":"$30","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":656,"value":"Drag selection window to change date range (currently selected: ${d3.utcDay.count(brush.startDate, brush.endDate)} days):","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":435,"value":"$31","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":867,"value":"$32","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":1712,"value":"## By Year","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":1714,"value":"$33","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":200,"value":"***\n\n## Implementation\n\nThe data comes from this [GitHub repo](https://github.com/harrislapiroff/mwra-wastewater-scraper), which I maintain. The code in the repo scrapes PDF from the MWRA website and converts the latest data to CSV. I run it whenever I notice the data is outdated.","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":3,"value":"$34","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":1701,"value":"data","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":1319,"value":"milestones = holidays.inRange(...d3.extent(data, d => d.date)).map(h => ({\n ...h,\n date: d3.utcDay.floor(h.date), // Convert dates to UTC\n}))","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":326,"value":"### Check most recent update\n\nThis bit checks the MWRA website to see if there's a more recent PDF than the data we have. If there is, [**@harrislapiroff**](https://observablehq.com/@harrislapiroff/) needs to convert it to CSV and upload to this notebook.","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":339,"value":"RECENT_DATE_RE = /Biobot Data - samples submitted through (\\d{1,2}\\/\\d{1,2}\\/\\d{4})/","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":330,"value":"mwra_page = (await fetch(CORS_PROXY + 'https://www.mwra.com/biobot/biobotdata.htm')).text()","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":342,"value":"latestSiteUpdate = {\n try {\n return d3.utcParse('%m/%d/%Y')(mwra_page.match(RECENT_DATE_RE)[1])\n } catch {\n return null\n }\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":711,"value":"### Map","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":725,"value":"TOWNS = FileAttachment(\"ma-towns-state-plane-topo-1.json\").json()","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":732,"value":"SYSTEMS_BY_TOWN = new Map([\n ['northern', ['Wilmington', 'Reading', 'Bedford', 'Burlington', 'Woburn', 'Stoneham', 'Lexington', 'Winchester', 'Melrose', 'Arlington', 'Medford', 'Malden', 'Waltham', 'Belmont', 'Somerville', 'Everett', 'Revere', 'Chelsea', 'Watertown', 'Cambridge', 'Winthrop']],\n ['both', ['Newton', 'Boston', 'Brookline', 'Milton']],\n ['southern', ['Ashland', 'Framingham', 'Natick', 'Wellesley', 'Needham', 'Dedham', 'Quincy', 'Westwood', 'Canton', 'Randolph', 'Braintree', 'Weymouth', 'Hingham', 'Norwod', 'Canton', 'Walpole', 'Stoughton', 'Holbrook', 'Norwood']]\n])","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":823,"value":"systemForTown = t => {\n const system = Array.from(SYSTEMS_BY_TOWN.entries())\n .find(([system, towns]) => towns.includes(t))\n\n if (!system) return null\n\n return system[0]\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":785,"value":"includedTowns = topojson.merge(\n TOWNS,\n TOWNS.objects.towns.geometries.filter(t => !!systemForTown(t.properties.TOWN2))\n)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":839,"value":"mapFillScale = (system, both, none) => {\n const colorScale = chart.scale('color')\n const color = colorScale.apply(system)\n if (color) return color\n\n if (system === 'both') return both\n \n return none\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":903,"value":"colorScale = chart.scale('color')","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":900,"value":"bothTexture = textures.lines()\n .size(8)\n .strokeWidth(2.5)\n .stroke(colorScale.range[0])\n .background(colorScale.range[1])","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":718,"value":"map = {\n const mapWidth = Math.min(600, width)\n const map = d3.select(DOM.svg(mapWidth, mapWidth))\n map.attr('style', 'max-width: 100%')\n\n const reprojection = d3.geoIdentity().fitExtent([[20, 20], [mapWidth - 20, mapWidth - 20]],includedTowns)\n const path = d3.geoPath(reprojection)\n\n map.call(bothTexture)\n\n const townShapes = map.append('g')\n .selectAll()\n .data(topojson.feature(TOWNS, TOWNS.objects.towns).features)\n .join('path')\n .attr('d', path)\n .attr('fill', d => mapFillScale(systemForTown(d.properties.TOWN2), bothTexture.url(), '#EEE'))\n .attr('stroke', '#FFF')\n \n return map\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":369,"value":"### Utils","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":703,"value":"Given the entire dataset, find the most recent date for which data exists:","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":374,"value":"latestDate = data => data.filter(d => !!d.concentrationAvg).slice(-1)[0].date","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":932,"value":"Capitalize the first char of a string:","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":207,"value":"capitalize = s => s.substring(0, 1).toUpperCase() + s.substring(1)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":936,"value":"Convert a string to a number (or null, if it's an empty string):","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":34,"value":"toIntOrNull = d => d === \"\" ? null : +d","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":948,"value":"Pluralize a word conditionally, based on a provided count:","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":940,"value":"pluralize = (word, count, plural) => count === 1 ? word : (plural || word + 's')","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":1050,"value":"parseVariantResults = str => {\n if (str === '') return\n const [variant, pctStr] = str.split('\\r')\n return { [variant]: +pctStr.substring(0, pctStr.length - 1) / 100 }\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":1416,"value":"Format a ratio as \"N× higher/lower\":","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":1419,"value":"formatRatio = (n, formats = {}) => {\n const { numberFormat, lowerFormat, higherFormat, exactString } = Object.assign({}, {\n numberFormat: d3.format('d'),\n higherFormat: n => `↑${n}× over`,\n lowerFormat: n => `↓${n}× below`,\n exactString: 'At'\n }, formats)\n if (n < 0) throw('Ratio cannot be smaller than 0')\n if (n < 1) return lowerFormat(numberFormat(1 / n))\n if (n > 1) return higherFormat(numberFormat(n))\n return exactString\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":1487,"value":"chartWidth = Math.max(width, 600)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":241,"value":"***\n\n## Appendix","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":245,"value":"import {on} from \"@fil/plot-onclick-experimental-plugin\"","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":846,"value":"textures = require('textures@1.2.3')","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":254,"value":"tippy = require('tippy.js@6.3.7')","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":1839,"value":"holidays = import('https://cdn.skypack.dev/@18f/us-federal-holidays@4.0.0?min')","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":306,"value":"","pinned":false,"hidden":false,"mode":"html","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":332,"value":"CORS_PROXY = 'https://corsproxy.harrislapiroff.com/'","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":1307,"value":"DATA_URL = `https://raw.githubusercontent.com/harrislapiroff/mwra-wastewater-scraper/main/output/latest.csv?cachebreaker=${d3.timeFormat('%Y-%m-%d')(new Date())}`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":1726,"value":"pctThroughYear = (date) => {\n const year = date.getFullYear()\n const isLeapYear = (0 == year % 4) && (0 != year % 100) || (0 == year % 400)\n return d3.utcDay.count(d3.utcYear.floor(date), d3.utcDay.ceil(date)) / (isLeapYear ? 366 :365)\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":1779,"value":"dateFromPct = (pct) => {\n const dayOfYear = pct * 365 // will be slightly inaccurate for leap years, fine for our purposes\n return d3.utcDay.offset(d3.utcYear.floor(new Date()), dayOfYear)\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"}],"files":"$W35","stdlib":"1"}} 2f:"$2e:body"