31:T7d6,caseCountMap = map({ title: `${d3.timeFormat('%B %d, %Y')(new Date())}: COVID-19 reports in Montana by county (${totalCases} cases / ${deaths} death${(deaths === 1) ? '' : 's'})`, height: 400, projection: d3.geoAlbers().rotate([110, 0, 0]), layers: [ new FeatureLayer({ geodata: counties, featureStyle: { fill: d => { // const countiesWithCases = countyCounts.map(c => c.county) const hasCases = (countiesWithCases.includes(d.NAMELABEL)) if (hasCases) return '#ffd92f' else return '#ddd' } } }), new FeatureLayer({ geodata: lakes, featureStyle: { fill: d => '#bbb', opacity: d => 1, stroke: d => 'none', } }), new FeatureLayer({ geodata: highways, featureStyle: { fill: d => 'none', stroke: d => '#bbb', } }), new PointLayer({ geodata: cities, labelStyle: { text: d => d.NAME, dx: d => d.offset ? d.offset[0] : '0.5em', dy: d => d.offset ? d.offset[1] : '0.5em', }, }), new PointLayer({ geodata: markCounties, featureStyle: { radius: d => rScale(d.cases), fill: d => '#CA532B', fillOpacity: d => 0.6, stroke: d => (d.deaths > 0) ? '#000' : '#a84320', strokeWidth: d => (d.deaths > 0) ? 2 : 1, }, labelStyle: { text: d => (d.deaths > 0) ? `${d.cases}/${d.deaths}` : `${d.cases}`, dx: d => 0, dy: d => '0.3em', textAnchor: d => 'middle', fontSize: d => '12px', fill: d => '#fff', }, }), new AnnotationLayer({ annotations: [ {x: 670, y: 330, text: `Data: covid19.mt.gov`}, {x: 670, y: 345, text: `Graphic: Eric Dietrich / MTFP`}, ], annotationStyle: { textSize: '12px', textStyle: 'italic', textAnchor: 'end', } }) ], // title: '2020 proposed grants by county' })32:Tb42,activeCaseCountMap = () => map({ // title: (width > 500) ? `Known active COVID-19 cases in Montana, ${d3.timeFormat('%B %d, %Y')(new Date())}` : `Known active COVID-19 cases`, width: largeMapWidth(), height: largeMapWidth() * 0.6, margin: {top: 10, left: 10, right: 10, bottom: 10}, projection: d3.geoAlbers().rotate([110, 0, 0]), layers: [ new FeatureLayer({ geodata: counties, featureStyle: { fill: d => { // const masksRequired = (d.activeCases >= 4) // const masksRequired = true // as of Nov. 20, 2020 const masksRequired = false // as of Jan. 12, 2021 const hasCases = (countiesWithCases.includes(d.NAMELABEL)) if (masksRequired) return '#e8a96b' if (hasCases) return '#f2e091' else return '#ddd' } } }), new FeatureLayer({ geodata: lakes, featureStyle: { fill: d => '#aaa', opacity: d => 1, stroke: d => 'none', } }), new FeatureLayer({ geodata: highways, featureStyle: { fill: d => 'none', stroke: d => '#aaa', } }), new PointLayer({ geodata: cities, labelStyle: { text: d => (width > 500) ? d.NAME : '', dx: d => d.offset ? d.offset[0] : '0.5em', dy: d => d.offset ? d.offset[1] : '0.5em', }, }), new PointLayer({ geodata: markCounties, featureStyle: { radius: d => rScale(d.activeCases), fill: d => '#1f78b4', fillOpacity: d => 0.6, // stroke: d => (d.deaths > 0) ? '#000' : '#a84320', stroke: d => '#1f78b4', // strokeWidth: d => (d.deaths > 0) ? 2 : 1, strokeWidth: d => 1, }, labelStyle: { text: d => `${d.activeCases}`, dx: d => ((d.activeCases > 9) || (d.activeCases <= 0)) ? 0 : '0.7em', dy: d => '0.3em', textAnchor: d => 'middle', fontSize: d => '12px', fill: d => (d.activeCases > 9) ? '#fff' : '#444', }, }), // new AnnotationLayer({ // annotations: [ // {x: Math.min(width, 800) - 50, y: (Math.min(width, 800) * 0.6) - 60, text: `Data: covid19.mt.gov`}, // {x: Math.min(width, 800) - 50, y: (Math.min(width, 800) * 0.6) - 45, text: `Graphic: Eric Dietrich / MTFP`}, // ], // annotationStyle: { // textSize: '12px', // textStyle: 'italic', // textAnchor: 'end', // } // }), // new AnnotationLayer({ // annotations: [ // {x: 40, y: 315, text: `${activeCases} active cases`}, // {x: 40, y: 330, text: `${totalCases} total cases`}, // {x: 40, y: 345, text: `${deaths} deaths`}, // ], // annotationStyle: { // textSize: '12px', // // textStyle: 'italic', // textAnchor: 'start', // } // }) ], })33:T5eb,totalDeathMap = () => map({ title: ``, width: largeMapWidth(), height: largeMapWidth() * 0.6, margin: {top: 10, left: 10, right: 10, bottom: 10}, projection: d3.geoAlbers().rotate([110, 0, 0]), layers: [ new FeatureLayer({ geodata: counties, featureStyle: { fill: d => { if (d.deaths > 0) return '#b5915e' else return '#ddd' }, opacity: d => 0.7, } }), new FeatureLayer({ geodata: lakes, featureStyle: { fill: d => '#bbb', opacity: d => 1, stroke: d => 'none', } }), new FeatureLayer({ geodata: highways, featureStyle: { fill: d => 'none', stroke: d => '#bbb', } }), new PointLayer({ geodata: cities, customMarker: squareMarker, labelStyle: { text: d => d.NAME, dx: d => d.offset ? d.offset[0] : '0.5em', dy: d => d.offset ? d.offset[1] : '0.5em', }, }), new PointLayer({ geodata: markCounties, featureStyle: { radius: d => deathSizeScale(d.deaths), fill: d => '#636363', fillOpacity: d => 0.7, // stroke: d => '#CA532B', stroke: d => '#222', strokeWidth: d => 2, }, labelStyle: { text: d => (d.deaths > 0) ? `${d.deaths}` : ``, dx: d => 0, dy: d => '0.3em', textAnchor: d => 'middle', fontSize: d => '12px', fill: d => '#fff', }, }), ], })34:T6ce,vaccinationMap = () => map({ // title: (width > 500) ? `Known active COVID-19 cases in Montana, ${d3.timeFormat('%B %d, %Y')(new Date())}` : `Known active COVID-19 cases`, width: largeMapWidth(), height: largeMapWidth() * 0.6, margin: {top: 10, left: 10, right: 10, bottom: 10}, projection: d3.geoAlbers().rotate([110, 0, 0]), layers: [ new FeatureLayer({ geodata: counties, featureStyle: { fill: d => { return '#ccc' } } }), new FeatureLayer({ geodata: lakes, featureStyle: { fill: d => '#aaa', opacity: d => 1, stroke: d => 'none', } }), new FeatureLayer({ geodata: highways, featureStyle: { fill: d => 'none', stroke: d => '#aaa', } }), new PointLayer({ geodata: cities, labelStyle: { text: d => (width > 500) ? d.NAME : '', dx: d => d.offset ? d.offset[0] : '0.5em', dy: d => d.offset ? d.offset[1] : '0.5em', }, }), new PointLayer({ geodata: vaccinationCounties, customMarker: vaccinationPieMarker, featureStyle: { radius: d => rScale(d.activeCases), fill: d => '#1f78b4', fillOpacity: d => 0.6, // stroke: d => (d.deaths > 0) ? '#000' : '#a84320', stroke: d => '#1f78b4', // strokeWidth: d => (d.deaths > 0) ? 2 : 1, strokeWidth: d => 1, }, labelStyle: { // text: d => `${d.totalDosesAdministered}`, dx: d => 0, dy: d => '0.3em', textAnchor: d => 'middle', fontSize: d => '12px', // fill: d => (d.activeCases > 0) ? '#fff' : '#222', fill: d => '#222' }, }), ], })35:T681,communityTransmissionMap = () => map({ // title: (width > 500) ? `Known active COVID-19 cases in Montana, ${d3.timeFormat('%B %d, %Y')(new Date())}` : `Known active COVID-19 cases`, width: largeMapWidth(), height: largeMapWidth() * 0.6, margin: {top: 10, left: 10, right: 10, bottom: 10}, projection: d3.geoAlbers().rotate([110, 0, 0]), layers: [ new FeatureLayer({ geodata: cdcMatchedCounties, featureStyle: { fill: d => { const color = { 'low': '#9ecae1', 'moderate': '#fee391', 'substantial': '#fd8d3c', 'high': '#d7301f' }[d.cdcTransmissionCategory] || '#ccc' return color // const masksRequired = (d.activeCases >= 4) // const masksRequired = true // as of Nov. 20, 2020 // const masksRequired = false // as of Jan. 12, 2021 // const hasCases = (countiesWithCases.includes(d.NAMELABEL)) // if (masksRequired) return '#e8a96b' // if (hasCases) return '#f2e091' // else return '#ddd' } } }), // new FeatureLayer({ // geodata: lakes, // featureStyle: { // fill: d => '#aaa', // opacity: d => 1, // stroke: d => 'none', // } // }), new FeatureLayer({ geodata: highways, featureStyle: { fill: d => 'none', stroke: d => '#aaa', } }), new PointLayer({ geodata: cities, labelStyle: { text: d => (width > 500) ? d.NAME : '', dx: d => d.offset ? d.offset[0] : '0.5em', dy: d => d.offset ? d.offset[1] : '0.5em', }, }), ], })36:T615,totalDeathMapStatic20210921 = () => map({ title: ``, width: largeMapWidth(), height: largeMapWidth() * 0.6, margin: {top: 10, left: 10, right: 10, bottom: 10}, projection: d3.geoAlbers().rotate([110, 0, 0]), layers: [ new FeatureLayer({ geodata: counties, featureStyle: { fill: d => { if (d.deaths > 0) return '#b5915e' else return '#ddd' }, opacity: d => 0.7, } }), new FeatureLayer({ geodata: lakes, featureStyle: { fill: d => '#888', opacity: d => 1, stroke: d => 'none', } }), new FeatureLayer({ geodata: highways, featureStyle: { fill: d => 'none', stroke: d => '#888', } }), new PointLayer({ geodata: cities, customMarker: squareMarker, labelStyle: { text: d => d.NAME, dx: d => d.offset ? d.offset[0] : '0.5em', dy: d => d.offset ? d.offset[1] : '0.5em', }, }), new PointLayer({ geodata: merge(clean(snapshot20210929), counties), featureStyle: { radius: d => deathSizeScale(d.deaths), fill: d => '#636363', fillOpacity: d => 0.7, // stroke: d => '#CA532B', stroke: d => '#222', strokeWidth: d => 2, }, labelStyle: { text: d => (d.deaths > 0) ? `${d.deaths}` : ``, dx: d => 0, dy: d => '0.3em', textAnchor: d => 'middle', fontSize: d => '12px', fill: d => '#fff', }, }), ], })37:Ta2f,cities = ({ "type": "FeatureCollection", "features": [ { "type": "Feature", "geometry": { "type": "Point", "coordinates": [ -114.00941481367617, 46.86845609815254 ] }, "properties": { "GEOID": "3050200", "NAME": "Missoula", offset: [-15,12], } }, { "type": "Feature", "geometry": { "type": "Point", "coordinates": [ -105.83850296885869, 46.40588180031621 ] }, "properties": { "GEOID": "3049525", "NAME": "Miles City" } }, { "type": "Feature", "geometry": { "type": "Point", "coordinates": [ -114.32618037830662, 48.21553771988646 ] }, "properties": { "GEOID": "3040075", "NAME": "Kalispell", offset: [-50,5] } }, { "type": "Feature", "geometry": { "type": "Point", "coordinates": [ -109.68049816861198, 48.54275568290601 ] }, "properties": { "GEOID": "3035050", "NAME": "Havre" } }, { "type": "Feature", "geometry": { "type": "Point", "coordinates": [ -112.54225090137798, 45.99630474669668 ] }, "properties": { "GEOID": "3011397", "NAME": "Butte" } }, { "type": "Feature", "geometry": { "type": "Point", "coordinates": [ -111.05476616601253, 45.682775182946244 ] }, "properties": { "GEOID": "3008950", "NAME": "Bozeman" } }, { "type": "Feature", "geometry": { "type": "Point", "coordinates": [ -112.0196972521252, 46.59642057169749 ] }, "properties": { "GEOID": "3035600", "NAME": "Helena", offset: [-15,-5], } }, { "type": "Feature", "geometry": { "type": "Point", "coordinates": [ -108.55027669152761, 45.788885426373696 ] }, "properties": { "GEOID": "3006550", "NAME": "Billings", offset: [-20,15], } }, { "type": "Feature", "geometry": { "type": "Point", "coordinates": [ -111.29948415503986, 47.50216159847621 ] }, "properties": { "GEOID": "3032800", "NAME": "Great Falls" } } ] })38:T566,map = (opts) => { const layers = opts.layers if (layers.length === 0) return 'No layers to map' const title = opts.title const height = opts.height || 600 const width = opts.width || 800 // width comes predefined const topMargin = title ? 40 : 10 const margin = opts.margin || {top: topMargin, left: 40, right: 40, bottom: 10} const projection = opts.projection || d3.geoMercator() // Dimensional calcs const plotWidth = width - margin.left - margin.right const plotHeight = height - margin.top - margin.bottom const svg = d3.select(DOM.svg(width, height)) .style('background-color', '#fff') // .style("width", "100%") .style("height", "auto") // .style('border', '1px solid #ddd') .attr('font-family','sans-serif') const plot = svg.append('g') .attr('transform', `translate(${margin.left},${margin.top})`) .attr('font-size',10) .attr('text-anchor','middle') // fit map based on first layer (i.e. base layer) projection.fitSize([plotWidth, plotHeight], layers[0].geodata) layers.forEach(layer => layer.plot(plot, projection)) if (title) { svg.append('text') .attr('y', margin.top - 14) .attr('x', margin.left + plotWidth / 2) .text(title) .attr('text-anchor','middle') .attr('font-size', 14) .attr('font-weight', 'bold') } return svg.node(); }39:T47c,class FeatureLayer { constructor(opts) { // required opts this.geodata = opts.geodata this.features = opts.geodata.features // optional opts this.featureStyle = opts.featureStyle || {} // calculated this.style = { opacity: d => 1, fill: d => '#ddd', fillOpacity: d => 1, stroke: d => '#fff', strokeWidth: d => 1, strokeOpacity: d => 1, ...this.featureStyle } } plot(parent, projection){ const path = d3.geoPath().projection(projection); const style = this.style const shapes = parent.append('g').attr('class','feature-layer') shapes.selectAll('.feature') .data(this.features).enter() .append('path') .attr('d', path) .attr('class', 'feature') .attr('opacity', (f => style.opacity(f.properties))) .attr('fill', (f => style.fill(f.properties))) .attr('fillOpacity', (f => style.fillOpacity(f.properties))) .attr('stroke', (f => style.stroke(f.properties))) .attr('strokeWidth', (f => style.strokeWidth(f.properties))) .attr('strokeOpacity', (f => style.strokeOpacity(f.properties))) } }3a:T94e,class PointLayer { constructor(opts) { // required opts this.geodata = opts.geodata this.features = opts.geodata.features // optional opts this.featureStyle = opts.featureStyle || {} this.labelStyle = opts.labelStyle this.customMarker = opts.customMarker || null // calculated this.style = { radius: d => 2, opacity: d => 1, fill: d => '#222', fillOpacity: d => 1, stroke: d => null, strokeWidth: d => 1, strokeOpacity: d => 1, ...this.featureStyle } this.label = { text: d => null, fill: d => '#000', fontSize: d=> '12px', dx: d => '0.5em', dy: d => '0.5em', textAnchor: d => 'start', ...this.labelStyle } } plot(parent, projection){ const path = d3.geoPath().projection(projection); const style = this.style const label = this.label // const shapes = parent.append('g').attr('class','feature-layer') const customMarker = this.customMarker const markerGroups = parent.append('g').attr('class','feature-layer') .selectAll('.feature') .data(this.features).enter() .append('g') .attr('transform', d => { const point = path.centroid(d) return `translate(${point[0]},${point[1]})` }) if (customMarker) { markerGroups.each(customMarker, label) } else { markerGroups.append('circle') .attr('r', (f => style.radius(f.properties))) .attr('cx', 0) .attr('cy', 0) .attr('class', 'feature') .attr('opacity', (f => style.opacity(f.properties))) .attr('fill', (f => style.fill(f.properties))) .attr('fill-opacity', (f => style.fillOpacity(f.properties))) .attr('stroke', (f => style.stroke(f.properties))) .attr('stroke-width', (f => style.strokeWidth(f.properties))) .attr('stroke-opacity', (f => style.strokeOpacity(f.properties))) } if (label) { markerGroups.append('text') .attr('x',0) .attr('y',0) .attr('font-size', f => label.fontSize(f.properties)) .attr('text-anchor', f => label.textAnchor(f.properties)) .attr('dy', f => label.dy(f.properties)) .attr('dx', f => label.dx(f.properties)) .attr('fill', f => label.fill(f.properties)) .text(f => label.text(f.properties)) } } }3b:T493,class AnnotationLayer { constructor(opts) { // required opts this.annotations = opts.annotations // optional opts this.annotationStyle = opts.annotationStyle this.style = { fill: d => '#000', dx: d => '0.5em', dy: d => '0.5em', textAnchor: 'start', ...this.annotationStyle } } plot(parent, projection){ // const path = d3.geoPath().projection(projection); const style = this.style // const shapes = parent.append('g').attr('class','feature-layer') const annotationGroups = parent.append('g').attr('class','feature-layer') .selectAll('.feature') .data(this.annotations).enter() .append('g') .attr('transform', d => { return `translate(${d.x},${d.y})` }) annotationGroups.append('text') .attr('x',0) .attr('y',0) .attr('text-anchor', style.textAnchor) .attr('font-size', style.textSize) .attr('font-style', style.textStyle) // .attr('dy', f => style.dy(f.properties)) // .attr('dx', f => style.dx(f.properties)) // .attr('fill', f => style.fill(f.properties)) .text(f => f.text) } }3c:Tc15,class GroupedPointLayer { constructor(opts) { // required opts this.geodata = opts.geodata this.features = opts.geodata.features this.joinData = opts.joinData this.markerStyle = { // defaults radius: d => 5, opacity: d => 1, fill: d => '#444', fillOpacity: d => 1, stroke: d => '#222', strokeWidth: d => 1, strokeOpacity: d => 1, ...opts.markerStyle } // optional opts this.groupDx = opts.groupDx || 0 this.groupDy = opts.groupDy || 0 } plot(parent, projection){ this.path = d3.geoPath().projection(projection); const sortedFeatures = this.features.sort((a,b) => this.path.centroid(a)[0] - this.path.centroid(b)[0] ) // Awkward internal function b/c .each() screws with this. calls const markerStyle = this.markerStyle const sizer = markerStyle.radius const countyNameKey = item => item.county.toUpperCase() const joinData = this.joinData const groupDx = this.groupDx const groupDy = this.groupDy function buildMapGroup(d,i,nodes) { // console.log('group', joinData) const thisGroupData = joinData .filter(row => countyNameKey(row) === d.properties.NAME) if (thisGroupData.length === 0) return const labelRadiusThreshold = 20 const pack = (data, width, height) => d3.pack() .size([width - 2, height - 2]) .radius(d => sizer(d.data)) .padding(2) (d3.hierarchy({children: data}).sum(sizer)) const root = pack(thisGroupData, 1, 1) const leaf = d3.select(this) .append('g') .selectAll("g") .data(root.leaves()) .enter().append("g") // assumes pre-translated centroids // with jitter // .attr("transform", (d,i) => `translate(${d.x + 1 + Math.random()*2},${d.y + 1 + Math.random()*2})`) // without jitter .attr("transform", (d,i) => `translate(${d.x + 1 + groupDx},${d.y + 1 + groupDy})`) leaf.append("circle") .attr("id", d => (d.leafUid = DOM.uid("leaf")).id) .attr("r", d => d.r) .attr('opacity', f => markerStyle.opacity(f.data)) .attr('fill', f => markerStyle.fill(f.data)) .attr("fill-opacity", f => markerStyle.fillOpacity(f.data)) .attr('stroke', f => markerStyle.stroke(f.data)) .attr('stroke-width', f => markerStyle.strokeWidth(f.data)) .attr('stroke-opacity', f => markerStyle.strokeOpacity(f.data)) leaf.append("clipPath") .attr("id", d => (d.clipUid = DOM.uid("clip")).id) .append("use") .attr("xlink:href", d => d.leafUid.href); // leaf.append("title") // .text(makeTitleText) } const groups = parent.append('g').attr('class','grouped-point-layer') .selectAll('.point-group') .data(this.features).enter() .append('g').attr('class','point-group') .attr('transform', d => `translate(${this.path.centroid(d)[0]},${this.path.centroid(d)[1]})`) .each(buildMapGroup) } }3d:[{"id":"c1ba9ac597918c7eff95486da091d82b0bda248131f42ac4c220b8dded5b4fa8babd5bfbecdcdb11a44023dd0446d7cf0b499cf62961cadc5798053ab9656fb4","href":"https://static.observableusercontent.com/files/c1ba9ac597918c7eff95486da091d82b0bda248131f42ac4c220b8dded5b4fa8babd5bfbecdcdb11a44023dd0446d7cf0b499cf62961cadc5798053ab9656fb4","name":"state-covid-dashboard-6-30-2021.csv","type":"text/csv","lastModified":1625090149162.35,"size":18599},{"id":"4e55d25df775569d47439f9bd71995fbbd449b1cd7a945e68186885401992e6d36923f38d708ea06a17518ae0eae0b71fbeb2ea56635492b91793d786e63fcdb","href":"https://static.observableusercontent.com/files/4e55d25df775569d47439f9bd71995fbbd449b1cd7a945e68186885401992e6d36923f38d708ea06a17518ae0eae0b71fbeb2ea56635492b91793d786e63fcdb","name":"snapshot-20210929-cleaned.csv","type":"text/csv","lastModified":1632939568965.043,"size":2910}] 2f:{"title":"Montana COVID-19 cases","cells":[{"id":481,"value":"md`# Montana COVID-19 cases\n\nUsing case data from state government, accesed via the API behind [the state's official map](https://montana.maps.arcgis.com/apps/MapSeries/index.html?appid=7c34f3412536439491adcc2103421d4b) and cached via a [Github scraper repository](https://github.com/eidietrich/mt-covid-dashboard-history).\n`","pinned":false,"hidden":false,"mode":"ojs"},{"id":443,"value":"fullData = {\n\n // Pull data direct from state dashboard \n // const raw = await d3.json('https://services.arcgis.com/qnjIrwR8z5Izc0ij/ArcGIS/rest/services/COVID_Cases_Production_View/FeatureServer/0/query?f=json&where=Total%20%3C%3E%200&returnGeometry=false&spatialRel=esriSpatialRelIntersects&outFields=*&outSR=102100&resultOffset=0&resultRecordCount=56&cacheHint=true')\n\n // Pull data as cached in github\n const raw = await d3.json('https://raw.githubusercontent.com/eidietrich/mt-covid-dashboard-history/main/mt-county-case-snapshot.json')\n \n const features = raw.features\n return features.map(d => d.attributes)\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1056,"value":"vaccinationData = {\n // Data direct from state dashboard\n // const raw = await d3.json('https://services.arcgis.com/qnjIrwR8z5Izc0ij/arcgis/rest/services/COVID_Vaccination_PRD_View/FeatureServer/0/query?f=json&where=1%3D1&returnGeometry=false&spatialRel=esriSpatialRelIntersects&outFields=*&orderByFields=NAMELABEL%20asc&outSR=102100&resultOffset=0&resultRecordCount=57&resultType=standard&cacheHint=true')\n // Data as cached by git scraper\n const raw = await d3.json('https://raw.githubusercontent.com/eidietrich/mt-covid-dashboard-history/main/mt-county-vaccination-snapshot.json')\n const features = raw.features\n return features.map(d => ({\n ...d.attributes,\n dateReported: new Date(d.attributes.Date_Reported),\n }))\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1091,"value":"vaccinationCounts = vaccinationData.map(d => {\n // const deathCountMatch = deathsManual.find(c => c.county === d.NAMELABEL)\n return {\n county: d.NAMELABEL,\n eligiblePop: d.Eligible_Pop,\n firstDosesAdministered: d.Dose_1,\n secondDosesAdministered: d.Dose_2,\n totalDosesAdministered: d.Total_Doses_Admin,\n fractionFirstDose: (d.Dose_1 - d.Dose_2) / d.Eligible_Pop,\n fractionSecondDose: d.Dose_2 / d.Eligible_Pop,\n }\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1593,"value":"vaccinationCounts.sort((a,b) => a.fractionSecondDose - b.fractionSecondDose)","pinned":true,"hidden":false,"mode":"ojs"},{"id":577,"value":"deathsManual = [\n {county: 'Lincoln', deaths: 1},\n {county: 'Madison', deaths: 1},\n {county: 'Toole', deaths: 6},\n {county: 'Missoula', deaths: 1},\n {county: 'Flathead', deaths: 2},\n {county: 'Cascade', deaths: 2},\n {county: 'Yellowstone', deaths: 1},\n]","pinned":false,"hidden":false,"mode":"ojs"},{"id":633,"value":"itemizedDeaths = {\n const countiesWithDeaths = Array.from(deathsManual.map(d => d.county))\n let deaths = []\n countiesWithDeaths.forEach(county => {\n deathsManual.filter(d => d.county === county).forEach(d => {\n for (let i=0;i html`

\n \n Active cases in county\n \n\n \n\n

`","pinned":false,"hidden":false,"mode":"ojs"},{"id":922,"value":"activeCaseCountEmbed = html`\n

\n\n

${(width > 500) ? `Known active COVID-19 cases in Montana, ${d3.timeFormat('%B %d, %Y')(new Date())}` : `Known active COVID-19 cases`}

\n${activeCaseMapLegend()}\n

${activeCaseCountMap()}

Data: Montana COVID-19 dashboard. Map: Eric Dietrich / MTFP

\n`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1294,"value":"md`### Cumulative death counts map`","pinned":false,"hidden":false,"mode":"ojs"},{"id":626,"value":"$33","pinned":false,"hidden":false,"mode":"ojs"},{"id":977,"value":"totalDeathLegend = () => html`

\n \n COVID-19 deaths in county\n\n \n Counties with at least one death\n

`","pinned":false,"hidden":false,"mode":"ojs"},{"id":974,"value":"totalDeathEmbed = html`\n

\n\n

COVID-19 deaths reported in Montana by county as of ${d3.timeFormat('%B %d, %Y')(new Date())}

${deaths} deaths total

\n${totalDeathLegend()}\n

${totalDeathMap()}

Data: Montana COVID-19 dashboard. Map: Eric Dietrich / MTFP

\n`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1297,"value":"md`### Vaccination progress map`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1183,"value":"vaccinationColors = ['#08519c', '#78c679', '#fee391']","pinned":false,"hidden":false,"mode":"ojs"},{"id":1180,"value":"vaccinationLegend = () => html`

\n\n \n Vaccine-eligible population\n\n \n Fraction with 1st dose\n\n \n Fraction fully immunized\n\n

`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1082,"value":"vaccinationPieMarker = (d, i, nodes) => {\n // console.log(d, i, nodes)\n const data = d.properties\n const pied = d3.pie()([\n data.fractionSecondDose,\n data.fractionFirstDose,\n 1 - data.fractionFirstDose - data.fractionSecondDose\n ])\n const w = 3\n const radius = popScale(data.eligiblePop)\n d3.select(nodes[i]).append('g')\n .selectAll('arcs')\n .data(pied).enter()\n .append('path')\n .attr('d', d3.arc().innerRadius(0).outerRadius(radius))\n .attr('fill', (d, i) => vaccinationColors[i])\n .attr('opacity', 0.8)\n .attr('class','feature')\n .attr('transform', 'scaleX(-1)')\n \n d3.select(nodes[i]).append('circle')\n .attr('r', radius)\n .attr('cx', 0)\n .attr('cy', 0)\n .attr('class', 'feature')\n .attr('stroke', '#222')\n .attr('stroke-width', 0.75)\n .attr('fill', 'none')\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1074,"value":"vaccinationEmbed = html`\n

\n\n

${(width > 500) ? `COVID-19 vaccination progress by county, ${d3.timeFormat('%B %d, %Y')(new Date())}` : `COVID-19 vaccination progress`}

\n${vaccinationLegend()}\n

${await vaccinationMap()}

Data: Montana COVID-19 dashboard. Note that the source data doesn't include vaccinations issued by the U.S. Dept. of Veterans Affairs or Indian Health Service. Map: Eric Dietrich / MTFP

\n`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1072,"value":"$34","pinned":false,"hidden":false,"mode":"ojs"},{"id":1313,"value":"md`### CDC community transmission map`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1326,"value":"cdcData = d3.json('https://raw.githubusercontent.com/eidietrich/cdc-mask-map-history/main/montana-cleaned.json')","pinned":false,"hidden":false,"mode":"ojs"},{"id":1330,"value":"cdcMatchedCounties = ({\n ...counties,\n features: counties.features.map(feature => {\n const match = cdcData.find(d => +d.fips_code === +feature.properties.ALLFIPS)\n if (!match) throw `Missing cdc Data match ${feature.properties.county}/${feature.properties.ALLFIPS}`\n return {\n ...feature,\n properties: {\n ...feature.properties,\n cdcTransmissionCategory: match.community_transmission_level,\n }\n }\n })\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1316,"value":"communityTransmissionMapEmbed = html`\n

\n\n

Current CDC COVID-19 community transmission rating

\n${communityTransmissionLegend()}\n

${communityTransmissionMap()}

Source: Live data from the CDC's community transmission map

\n`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1357,"value":"communityTransmissionLegend = () => html`

\n\n\n \n Low\n \n\n\n \n Moderate\n \n\n\n \n Substantial*\n \n\n\n \n High*\n \n\n

*Masks recommended in indoor public spaces as of July 27

\n\n

`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1318,"value":"$35","pinned":false,"hidden":false,"mode":"ojs"},{"id":1393,"value":"communityTransmissionCounts = {\n const buckets = ['high','substantial','moderate','low']\n return buckets.map(bucket => ({\n bucket,\n count: cdcData.filter(d => d.community_transmission_level === bucket).length,\n }))\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1422,"value":"md`## Static graphics for 9/29/21 story`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1465,"value":"snapshot20210929 = FileAttachment(\"snapshot-20210929-cleaned.csv\").csv()","pinned":false,"hidden":false,"mode":"ojs"},{"id":1427,"value":"totalDeathEmbedStatic20210921 = html`\n

\n\n

COVID-19 deaths reported in Montana by county as of September 29, 2021

\n${totalDeathLegend()}\n

${totalDeathMapStatic20210921()}

Data: Montana COVID-19 dashboard. Map: Eric Dietrich / MTFP

\n`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1527,"value":"deathDistributionEmbedStatic20210921 = html`\n

\n\n

Montana's 2020 COVID-19 deaths

${deathAgeDistribution2020()}

${deathRaceDistribution2020()}

${deathPlaceDistribution2020()}

Data: Montana Dept. Public Health & Human Services. Counts are based on death certificate data. In contrast to some reported figures, they include both resident and non-resident deaths.

\n`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1432,"value":"$36","pinned":false,"hidden":false,"mode":"ojs"},{"id":1475,"value":"deathsByAgeBracket2020 = {\n // data from https://dphhs.mt.gov/assets/publichealth/Epidemiology/VSU/VSUDeathCertificateICDDx2020.pdf\n // Represents residents and non-resident deaths\n const data = [\n {bracket: '0-18', count: 1},\n {bracket: '18–29', count: 1},\n {bracket: '30–39', count: 17},\n {bracket: '40–49', count: 30},\n {bracket: '50–64', count: 167},\n {bracket: '65–74', count: 276},\n {bracket: '75–84', count: 372},\n {bracket: '85+', count: 415},\n ]\n const total = d3.sum(data, d => d.count)\n data.forEach((d,i) => {\n d.index = i\n d.fraction = d.count / total\n })\n return data\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1544,"value":"deathsByRace2020 = {\n // data from https://dphhs.mt.gov/assets/publichealth/Epidemiology/VSU/VSUDeathCertificateICDDx2020.pdf\n // Represents residents and non-resident deaths\n const data = [\n {bracket: 'White', count: 1002},\n {bracket: 'American Indian', count: 263},\n {bracket: 'Other', count: 14},\n ]\n const total = d3.sum(data, d => d.count)\n data.forEach((d,i) => {\n d.index = i\n d.fraction = d.count / total\n })\n return data\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1568,"value":"deathsByPlace2020 = {\n // data from https://dphhs.mt.gov/assets/publichealth/Epidemiology/VSU/VSUDeathCertificateICDDx2020.pdf\n // Represents residents and non-resident deaths\n const data = [\n {bracket: 'Hospital', count: 765},\n {bracket: 'Nursing home', count: 355},\n {bracket: 'Personal home', count: 103},\n {bracket: 'Other', count: 56},\n ]\n const total = d3.sum(data, d => d.count)\n data.forEach((d,i) => {\n d.index = i\n d.fraction = d.count / total\n })\n return data\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1479,"value":"deathAgeDistribution2020 = () => Plot.plot({\n width: Math.min(width - 10, 600), \n marginLeft: 90,\n x: {\n label: 'Recorded deaths',\n domain: [0,1000]\n },\n y: {\n label: 'Age bracket',\n labelAnchor: 'top',\n },\n marks: [\n Plot.barX(deathsByAgeBracket2020, {x: \"count\", y: \"bracket\", fill: '#665' }),\n Plot.ruleX([0])\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1547,"value":"deathRaceDistribution2020 = () => Plot.plot({\n width: Math.min(width - 10, 600),\n marginLeft: 90,\n x: {\n label: 'Recorded deaths',\n domain: [0,1000]\n },\n y: {\n label: 'Race',\n labelAnchor: 'top',\n domain: deathsByRace2020.map(d => d.bracket)\n },\n marks: [\n Plot.barX(deathsByRace2020, {x: \"count\", y: \"bracket\", fill: '#665' }),\n Plot.ruleX([0])\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1570,"value":"deathPlaceDistribution2020 = () => Plot.plot({\n width: Math.min(width - 10, 600),\n marginLeft: 90,\n x: {\n label: 'Recorded deaths',\n domain: [0,1000]\n },\n y: {\n label: 'Place of death',\n labelAnchor: 'top',\n domain: deathsByPlace2020.map(d => d.bracket)\n },\n marks: [\n Plot.barX(deathsByPlace2020, {x: \"count\", y: \"bracket\", fill: '#665' }),\n Plot.ruleX([0])\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1308,"value":"md`### Case age distribution`","pinned":false,"hidden":false,"mode":"ojs"},{"id":521,"value":"embed({\n \"$schema\": \"https://vega.github.io/schema/vega-lite/v4.json\",\n \"description\": \"A simple bar chart with embedded data.\",\n width: 500,\n \n background: \"#fff\",\n \"data\": {\n \"values\": ageBracketCounts\n },\n title: `Montana COVID-19 cases by age bracket ${d3.timeFormat('%B %d, %Y')(new Date())} (${totalCases} total)`,\n layer: [\n {\n \"mark\": { type:\"bar\", color: \"#CA532B\"},\n \"encoding\": {\n \"x\": {\"field\": \"ageBracket\", \"type\": \"ordinal\", title: \"Age bracket\"},\n \"y\": {\"field\": \"count\", \"type\": \"quantitative\", title: \"COVID-19 case count\",\n // scale: {domain:[0,50]}\n }\n }\n },\n {\n \"mark\": {type: \"text\", dy: -5},\n \"encoding\": {\n \"x\": {\"field\": \"ageBracket\", \"type\": \"ordinal\", title: \"Age bracket\"},\n \"y\": {\"field\": \"count\", \"type\": \"quantitative\"},\n \"text\": {\"field\": \"count\", \"type\": \"quantitative\", },\n }\n },\n ]\n \n}\n\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":1304,"value":"md`### Miscellanea`","pinned":false,"hidden":false,"mode":"ojs"},{"id":914,"value":"largeMapWidth = () => {\n if (width <= 600) {\n return width - 20\n } else if (width >= 1200 ){\n return 600\n } else {\n return width * 0.9 * 0.66 - 10\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1012,"value":"secondaryMapWidth = () => {\n if (width <= 600) {\n return width - 20\n } else if (width >= 1400 ){\n return 1400 * 0.9 * 0.5\n } else {\n return width * 0.9 * 0.5 - 10\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":16,"value":"rawDataOld = [\n {date: '3/13/20', county: 'Gallatin', count: 1},\n {date: '3/13/20', county: 'Yellowstone', count: 1},\n {date: '3/13/20', county: 'Silver Bow', count: 1},\n {date: '3/13/20', county: 'Broadwater', count: 1},\n {date: '3/14/20', county: 'Missoula', count: 2},\n {date: '3/16/20', county: 'Missoula', count: 1},\n {date: '3/16/20', county: 'Yellowstone', count: 1},\n {date: '3/17/20', county: 'Gallatin', count: 2}, // includes one case reclassified as New Hampshire resident\n {date: '3/18/20', county: 'Missoula', count: 1},\n {date: '3/18/20', county: 'Gallatin', count: 1},\n {date: '3/19/20', county: 'Yellowstone', count: 3},\n {date: '3/19/20', county: 'Roosevelt', count: 1},\n]","pinned":false,"hidden":false,"mode":"ojs"},{"id":1445,"value":"clean = (fullData) => {\n return fullData.map(d => {\n // const deathCountMatch = deathsManual.find(c => c.county === d.NAMELABEL)\n return {\n county: d.NAMELABEL,\n cases: d.Total,\n newCases: d.NewCases,\n // deaths: deathCountMatch ? deathCountMatch.deaths : 0,\n deaths: d.TotalDeaths,\n hospitalized: d.HospitalizationCount,\n recoveredCases: d.TotalRecovered,\n activeCases: d.TotalActive,\n }\n }).filter(d => d.cases > 0)\n} ","pinned":false,"hidden":false,"mode":"ojs"},{"id":488,"value":"countyCounts = clean(fullData)\n// countyCounts = fullData.map(d => {\n// // const deathCountMatch = deathsManual.find(c => c.county === d.NAMELABEL)\n// return {\n// county: d.NAMELABEL,\n// cases: d.Total,\n// newCases: d.NewCases,\n// // deaths: deathCountMatch ? deathCountMatch.deaths : 0,\n// deaths: d.TotalDeaths,\n// hospitalized: d.HospitalizationCount,\n// recoveredCases: d.TotalRecovered,\n// activeCases: d.TotalActive,\n// }\n// }).filter(d => d.cases > 0)","pinned":false,"hidden":false,"mode":"ojs"},{"id":499,"value":"ageBracketCounts = {\n const ageBrackets = Object.keys(fullData[0]).filter(d => d.match(\"^T_\"))\n const ageBracketCounts = ageBrackets.map(bracketKey => {\n return {\n key: bracketKey,\n ageBracket: bracketKey.replace(\"T_\",\"\").replace('_','–'),\n count: d3.sum(fullData, d => d[bracketKey]),\n countMale: d3.sum(fullData, d => d[`M_${bracketKey.replace(\"T_\",\"\")}`]),\n countFemale: d3.sum(fullData, d => d[`F_${bracketKey.replace(\"T_\",\"\")}`]),\n }\n })\n return ageBracketCounts\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":125,"value":"rScale = {\n if (width < 800) {\n return d3.scaleSqrt().domain([0,20000]).range([0,75])\n }\n else {\n return d3.scaleSqrt().domain([0,20000]).range([0,200])\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1111,"value":"popScale = {\n // return d3.scaleSqrt().domain([0,100000]).range([4,20])\n if (width < 800) {\n return d3.scaleSqrt().domain([0,100000]).range([4,20])\n }\n else {\n return d3.scaleSqrt().domain([0,100000]).range([4,25])\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":715,"value":"deathSizeScale = d3.scaleSqrt().domain([0,200]).range([2,25])","pinned":false,"hidden":false,"mode":"ojs"},{"id":280,"value":"// filterDate = new Date('3/21/20')\nfilterDate = new Date()","pinned":false,"hidden":false,"mode":"ojs"},{"id":364,"value":"totalCases = d3.sum(countyCounts, d => d.cases)","pinned":false,"hidden":false,"mode":"ojs"},{"id":775,"value":"activeCases = d3.sum(countyCounts, d => d.activeCases)","pinned":false,"hidden":false,"mode":"ojs"},{"id":584,"value":"deaths = {\n const countyDeaths = d3.sum(countyCounts, d => d.deaths)\n const unlocatedDeathsObj = deathsManual.find(d => d.county === 'Unknown')\n if (unlocatedDeathsObj) {\n return countyDeaths + unlocatedDeathsObj.deaths\n } else {\n return countyDeaths\n }\n \n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":881,"value":"countyCounts.map(d => ({county: d.county, deaths: d.deaths})).filter(d => d.deaths > 0).sort((a,b) => (b.deaths - a.deaths))","pinned":false,"hidden":false,"mode":"ojs"},{"id":493,"value":"countiesWithCases = countyCounts.map(c => c.county)","pinned":false,"hidden":false,"mode":"ojs"},{"id":53,"value":"md`## Geodata`","pinned":false,"hidden":false,"mode":"ojs"},{"id":47,"value":"counties = {\n const mtCounties = await d3.json('https://gist.githubusercontent.com/eidietrich/1f4727a8fb1a32bda14c921b29028bfd/raw/2479b3ae8ca3aa09c2b9f2cb9c9ca7f892250a96/mt-counties.geojson')\n \n return mtCounties\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1443,"value":"merge = (dataCleaned, geojson) => {\n const features = geojson.features\n features.forEach(feature => {\n const match = dataCleaned.find(d => d.county === feature.properties.NAMELABEL)\n feature.properties = {\n ...feature.properties,\n ...match,\n }\n })\n\n return {\n type: geojson.type,\n crs: geojson.crs,\n features: features.sort((a,b) => b.properties.count - a.properties.count),\n }\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":110,"value":"markCounties = merge(clean(fullData), counties)\n// markCounties = {\n// counties.features.forEach(county => {\n// const match = countyCounts.find(d => d.county === county.properties.NAMELABEL)\n// county.properties = {\n// ...county.properties,\n// ...match,\n// }\n// })\n\n// return {\n// type: counties.type,\n// crs: counties.crs,\n// features: counties.features.sort((a,b) => b.properties.count - a.properties.count),\n// }\n// }","pinned":true,"hidden":false,"mode":"ojs"},{"id":1086,"value":"vaccinationCounties = {\n const mergedFeatures = counties.features.map(county => {\n const match = vaccinationCounts.find(d => d.county === county.properties.NAMELABEL)\n return {\n ...county,\n properties: {\n ...county.properties,\n ...match,\n }\n }\n })\n\n return {\n type: counties.type,\n crs: counties.crs,\n features: mergedFeatures.sort((a,b) => b.properties.eligiblePop - a.properties.eligiblePop),\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":50,"value":"$37","pinned":false,"hidden":false,"mode":"ojs"},{"id":57,"value":"highways = {\n const highways = await d3.json('https://gist.githubusercontent.com/eidietrich/67cbe6196b6b87698fc42f50e96076ab/raw/4d51ceb627c76360497837e1ac5338698f7df752/ne_mt_highways_state_resolution.topojson.json')\n \n const topo = topojson.feature(highways, highways.objects.highways)\n topo.features = topo.features.filter(d => d.properties.type === 'Major Highway')\n return topo\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":61,"value":"lakes = {\n const topo = await d3.json('https://gist.githubusercontent.com/eidietrich/753ea92617c64c938ad36fddad3b0b7f/raw/63273fccb218ffdcf1d14b0c0d31cbe658588661/mt-major-lakes.json')\n \n const geo = topojson.feature(topo, topo.objects.lakes)\n return geo\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":29,"value":"md`## Libraries`","pinned":false,"hidden":false,"mode":"ojs"},{"id":33,"value":"d3 = require(\"d3@5\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":37,"value":"topojson = require(\"topojson-client@3\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":523,"value":"embed = require('vega-embed@4')","pinned":false,"hidden":false,"mode":"ojs"},{"id":39,"value":"$38","pinned":false,"hidden":false,"mode":"ojs"},{"id":41,"value":"$39","pinned":false,"hidden":false,"mode":"ojs"},{"id":261,"value":"// Note: PointLayere here has been updated vs earlier version, now handles non-point geometries","pinned":true,"hidden":false,"mode":"ojs"},{"id":94,"value":"$3a","pinned":false,"hidden":false,"mode":"ojs"},{"id":309,"value":"$3b","pinned":false,"hidden":false,"mode":"ojs"},{"id":628,"value":"$3c","pinned":false,"hidden":false,"mode":"ojs"},{"id":668,"value":"squareMarker = (d, i, nodes) => {\n // console.log(d, i, nodes)\n const w = 3\n d3.select(nodes[i]).append('rect')\n .attr('x', -w/2)\n .attr('y', -w/2)\n .attr('width', w)\n .attr('height', w)\n .attr('class', 'feature')\n .attr('stroke', '#666')\n .attr('fill', 'none')\n}","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W3d","stdlib":"1","titling":"auto"}