31:T1080,{ let ageGroupIndices = [selectedAgeRange]; let geography = stateEDdata.geography; let regionName = geographyToName(geography); let ageRangeNames = stateEDdata.age_groups.map(d=>d.age_range); function dataForAgeGroup(i) { return stateEDdata.ed_trends_by_Geography[0].age_groups[i].data; } function meanForAgeGroup(i) { return stateEDdata.ed_trends_by_Geography[0].age_groups[i].mean_percent; } let dataAll = dataForAgeGroup(9); let dataBaseline = dataForAgeGroup(5); let ageRangeString = ageGroupIndices.map(a=>ageRangeNames[a]); let key = "7day_Avg_Percent"; let ageMeanPercentage = meanForAgeGroup(selectedAgeRange); let allMeanPercentage = meanForAgeGroup(9); let ageAdjust = allMeanPercentage / ageMeanPercentage; let ratioData = dataAll.map((d,i)=>{ let all = d[key]; let baseline = dataBaseline[i][key]; let s = 0; ageGroupIndices.forEach(ageI=>s+=dataForAgeGroup(ageI)[i][key]||0); return { Date: d.Date, ageGroup:ageRangeString, ratio: s/all/ageGroupIndices.length, baselineRatio: s/baseline/ageGroupIndices.length, percentage: s/ageGroupIndices.length, ageAdjustedPercentage: s*ageAdjust/ageGroupIndices.length, allAgesPercentage: all, baselinePercentage: baseline, } }); console.log(ratioData); let x = d=>new Date(d.Date); let y = d=>d.ratio; let y2 = d=>d.baselineRatio; let channels = { Percentage: d=>d.percentage, "All ages %": d=>d.allAgesPercentage, "Ratio to all ages": y, }; let plot1 = Plot.plot({ title: html`${regionName} Percentage of ED visits with diagnosed COVID-19`, width: width, height: Math.min(window.screen.availHeight/3,width), color: { legend: true, }, x: { axis: 'bottom', grid: true, label: "Date", //domain: [0, 100], }, y: { grid: 15, label: `Percentage of ED visits with COVID diagnoses for ${ageRangeString} year olds and all ages.` }, marks: [ Plot.ruleY([0]), Plot.lineY(ratioData, { x, y: d=>d.percentage, stroke:d=>d.ageGroup, channels, tip: true, }), Plot.lineY(ratioData, { x, y: d=>d.allAgesPercentage, stroke:d=>"All ages", tip: false, }), //Plot.lineY(ratioData, { // x, // y: d=>d.baselinePercentage, // stroke:d=>"Baseline (40-49)" //}), ] }); let plot2 = Plot.plot({ title: html`${regionName} Age-adjusted percentage of ED visits with diagnosed COVID-19`, width: width, height: Math.min(window.screen.availHeight/3,width), color: { legend: true, }, x: { axis: 'bottom', grid: true, label: "Date", //domain: [0, 100], }, y: { grid: 15, label: `Age-adjusted percentage of ED visits with COVID diagnoses for ${ageRangeString} year olds and all ages.` }, marks: [ Plot.ruleY([0]), Plot.lineY(ratioData, { x, y: d=>d.ageAdjustedPercentage, stroke:d=>d.ageGroup, channels, tip: true, }), Plot.lineY(ratioData, { x, y: d=>d.allAgesPercentage, stroke:d=>"All ages", tip: false, }), ] }); let plot3 = Plot.plot({ title: html`${regionName} Ratio of percentage of ED visits with diagnosed COVID-19, for ${ageRangeString} year olds versus all.`, width: width, height: Math.min(window.screen.availHeight/3,width), x: { axis: 'bottom', grid: true, label: "Date", //domain: [0, 100], }, y: { grid: 15, label: `Ratio of the ED COVID % for ${ageRangeString} year olds compared to all.` }, marks: [ Plot.ruleY([0]), Plot.lineY(ratioData, {x, y, stroke:d=>"Ratio versus all",channels,tip:true}), // Comparing to 40-49 as a baseline did not result in a better plot, probably because of noise //Plot.lineY(ratioData, {x, y:y2, stroke:d=>"Ratio versus baseline"}), ] }); //console.log("Age adjustment",ageAdjust); return html`${plot1}
${plot2}
${plot3}` }32:[] 2f:{"title":"CDC COVID Emergency Department visits by age","cells":[{"id":0,"value":"# CDC COVID Emergency Department visits by age","pinned":false,"hidden":false,"mode":"md"},{"id":101,"value":"* Data from https://covid.cdc.gov/covid-data-tracker/#ed-visits_separated_by_age_group\n
(Note the coverage information in the footnotes)\n* X thread: https://twitter.com/JPWeiland/status/1712958387770536288","pinned":false,"hidden":false,"mode":"md"},{"id":128,"value":"**Ratio for age group relative to all ages** of percentage of ED visits with diagnosed COVID-19, for ${regionToName(selectedState)}.","pinned":false,"hidden":false,"mode":"md"},{"id":35,"value":"viewof selectedState = Inputs.select(regionList, {format:regionToName,label: \"Select state or region\",disabled:regionsDisabled})","pinned":false,"hidden":false,"mode":"ojs"},{"id":81,"value":"viewof selectedAgeRange = Inputs.select([...ageRangeNames.keys()].slice(0,-1), {format:d=>ageRangeNames[d],label: \"Select age range\",value:0})","pinned":false,"hidden":false,"mode":"ojs"},{"id":237,"value":"{\n if (loadingData) {\n return md`Loading data for ${regionToName(selectedState)}...`\n }\n if (regionToName(selectedState) !== geographyToName(geography)) {\n return md`Loading data for ${regionToName(selectedState)}......`\n }\n return md`${regionToName(selectedState)} loaded. Last date: ${lastData.Date}`;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":26,"value":"$31","pinned":false,"hidden":false,"mode":"ojs"},{"id":56,"value":"________\n### Code Appendix","pinned":false,"hidden":false,"mode":"md"},{"id":234,"value":"mutable loadingData = 1","pinned":true,"hidden":false,"mode":"ojs"},{"id":84,"value":"ageRangeNames = [\"00-11\", \"12-15\", \"16-17\", \"18-24\", \"25-39\", \"40-49\", \"50-64\", \"65-74\", \"75+\", \"all_ages\"]\n// Will reset the input choice if you get the names from the data:\n//stateEDdata.ed_trends_by_Geography[0].age_groups.map(d=>d.age_range);\n","pinned":false,"hidden":false,"mode":"ojs"},{"id":250,"value":"mutable geography = \"\"","pinned":true,"hidden":false,"mode":"ojs"},{"id":254,"value":"meanRateByAgeGroup = {\n stateEDdata.ed_trends_by_Geography[0].age_groups.forEach(d=>{\n d.mean_percent = d3.mean(d.data,dd=>dd.Daily_Percent || 0);\n });\n return stateEDdata;\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":320,"value":"lastDataOld = stateEDdata?.ed_trends_by_Geography?.[0].age_groups?.slice(-1)[0]?.data?.slice(-1)[0]","pinned":true,"hidden":false,"mode":"ojs"},{"id":415,"value":"lastData = stateEDdata.filter(d=>Math.abs(d.week_end - dataDateExtent[1]) < 60*1000)","pinned":true,"hidden":false,"mode":"ojs"},{"id":400,"value":"dataDateExtent = d3.extent(stateEDdata, d=>d.week_end)","pinned":true,"hidden":false,"mode":"ojs"},{"id":408,"value":"stateEDdata.slice(10)","pinned":true,"hidden":false,"mode":"ojs"},{"id":18,"value":"stateEDdata= {\n let data = stateNSSPdataDemo;\n //mutable geography = data[0].geography;\n return data;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":417,"value":"Plot.plot({\n width: width,\n x: {\n grid: true\n },\n y: {\n grid: true\n },\n color: {\n legend: true\n },\n marks: [\n Plot.auto(dataForPlot, {x: \"week_end\", y: \"Ratio 18-64\", color: \"demographics_values\"})\n ]\n})","pinned":true,"hidden":false,"mode":"ojs"},{"id":477,"value":"Plot.plot({\n width: width,\n x: {\n grid: true\n },\n y: {\n grid: true\n },\n color: {\n legend: true\n },\n marks: [\n Plot.auto(dataForPlot, {x: \"week_end\", y: \"Ratio White\", color: \"demographics_values\"})\n ]\n})","pinned":true,"hidden":false,"mode":"ojs"},{"id":484,"value":"Plot.plot({\n width: width,\n x: {\n grid: true\n },\n y: {\n grid: true\n },\n color: {\n legend: true\n },\n marks: [\n Plot.auto(dataForPlot, {x: \"week_end\", y: \"Ratio Male\", color: \"demographics_values\"})\n ]\n})","pinned":true,"hidden":false,"mode":"ojs"},{"id":426,"value":"dataForPlot = ratios.filter(d=>d.pathogen==\"COVID-19\" && \n d.demographics_values !== \"Other\" && \n d.demographics_values !== \"Invalid\" &&\n {Female:1, \"5-17 years\":1}[d.demographics_values]\n )","pinned":true,"hidden":false,"mode":"ojs"},{"id":524,"value":"maleFemaleDataAllPathogens = ratios.filter(d=>d.pathogen!=\"Combined\" && \n {Female:1, Male:1}[d.demographics_values]\n )","pinned":true,"hidden":false,"mode":"ojs"},{"id":527,"value":"Plot.plot({\n width: width,\n x: {\n grid: true\n },\n y: {\n grid: true\n },\n color: {\n legend: true\n },\n marks: [\n Plot.auto(maleFemaleDataAllPathogens, {x: \"week_end\", y: \"percent_visits\", fx: \"demographics_values\", color: \"pathogen\"})\n ]\n})","pinned":true,"hidden":false,"mode":"ojs"},{"id":506,"value":"maleFemaleData = ratios.filter(d=>d.pathogen==\"COVID-19\" && \n {Female:1, Male:1}[d.demographics_values]\n )","pinned":true,"hidden":false,"mode":"ojs"},{"id":508,"value":"Plot.plot({\n width: width,\n x: {\n grid: true\n },\n y: {\n grid: true\n },\n color: {\n legend: true\n },\n marks: [\n Plot.auto(maleFemaleData, {x: \"week_end\", y: \"percent_visits\", color: \"demographics_values\"})\n ]\n})","pinned":true,"hidden":false,"mode":"ojs"},{"id":522,"value":"Plot.plot({\n width: width,\n x: {\n grid: true\n },\n y: {\n grid: true\n },\n color: {\n legend: true\n },\n marks: [\n Plot.auto(maleFemaleData, {x: \"week_end\", y: \"Ratio Male\", color: \"demographics_values\"})\n ]\n})","pinned":true,"hidden":false,"mode":"ojs"},{"id":436,"value":"dataForPlot","pinned":true,"hidden":false,"mode":"ojs"},{"id":437,"value":"dataByPathogenDemographicsDate = d3.group(stateEDdata, d=>d.pathogen, d=>d.demographics_values, d=>d.week_end)","pinned":true,"hidden":false,"mode":"ojs"},{"id":439,"value":"dates = d3.flatRollup(stateEDdata, v=>v.length, d=>new Date(Math.round(d.week_end.valueOf()/1000)*1000))","pinned":true,"hidden":false,"mode":"ojs"},{"id":456,"value":"dataByPathogenDemographicsDate.get(\"COVID-19\").get(\"18-64 years\").get(new Date(\"2025-09-06T07:00Z\"))[0].percent_visits","pinned":true,"hidden":false,"mode":"ojs"},{"id":450,"value":"ratios = stateEDdata.map(d=> ({\n ...d,\n \"Ratio 18-64\": d.percent_visits / dataByPathogenDemographicsDate.get(d.pathogen)?.get(\"18-64 years\")?.get(d.week_end)?.[0]?.percent_visits,\n \"Ratio Male\": d.percent_visits / dataByPathogenDemographicsDate.get(d.pathogen)?.get(\"Male\")?.get(d.week_end)?.[0]?.percent_visits,\n \"Ratio White\": d.percent_visits / dataByPathogenDemographicsDate.get(d.pathogen)?.get(\"White, NH\")?.get(d.week_end)?.[0]?.percent_visits,\n}))","pinned":true,"hidden":false,"mode":"ojs"},{"id":364,"value":"stateEDdataOld = {\n mutable loadingData = 1;\n let url = urlForState(selectedState) ;\n console.log(\"Fetching\",url);\n let data = await d3.json(url);\n mutable loadingData = 0;\n mutable geography = data.ed_trends_by_Geography[0].Geography\n return data;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":5,"value":"function urlForState(state) {\n return corsProxy+encodeURIComponent(\"https://covid.cdc.gov/covid-data-tracker/COVIDData/getAjaxData?id=ed_trends_by_\" + state.toLowerCase());\n}\n","pinned":false,"hidden":false,"mode":"ojs"},{"id":10,"value":"urlForState(selectedState)","pinned":false,"hidden":false,"mode":"ojs"},{"id":20,"value":"corsProxy = 'https://cors.lvab.workers.dev/?u='","pinned":false,"hidden":false,"mode":"ojs"},{"id":32,"value":"regionList = ['united_states', //'allStates', \n 'alabama', 'alaska', \n //'american_samoa', \n 'arizona', 'arkansas', 'california', 'colorado', 'connecticut', 'delaware', 'florida', 'georgia', \n //'guam', \n 'hawaii', 'idaho', 'illinois', 'indiana', 'iowa', 'kansas', 'kentucky', 'louisiana', 'maine', 'maryland', 'massachusetts', 'michigan', 'minnesota', 'mississippi', 'missouri', 'montana', 'nebraska', 'nevada', 'new_hampshire', 'new_jersey', 'new_mexico', 'new_york', 'north_carolina', 'north_dakota',\n //'northern_mariana islands',\n 'ohio', 'oklahoma', 'oregon', 'pennsylvania', 'puerto_rico', 'rhode_island', 'south_carolina', 'south_dakota', 'tennessee', 'texas', 'utah', 'vermont', 'virginia', \n //'virgin_islands', \n 'washington', 'west_virginia', 'wisconsin', 'wyoming', 'region_1', 'region_2', 'region_3', 'region_4', 'region_5', 'region_6', 'region_7', 'region_8', 'region_9', 'region_10'];","pinned":false,"hidden":false,"mode":"ojs"},{"id":195,"value":"regionName = regionList.map(regionToName)","pinned":true,"hidden":false,"mode":"ojs"},{"id":197,"value":"function regionToName(s) {\n return s.charAt(0).toUpperCase() + s.slice(1).replace(/_./,function(match) {\n return match.toUpperCase();\n }).replace(\"_\",\" \");\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":304,"value":"function geographyToName(s) {\n return s.charAt(0).toUpperCase() + s.slice(1).replace(/ ./,function(match) {\n return match.toUpperCase();\n }).replace(\"_\",\" \");\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":105,"value":"regionsDisabled=[\"minnesota\",\"missouri\",\"oklahoma\",\"puerto_rico\"];","pinned":false,"hidden":false,"mode":"ojs"},{"id":270,"value":"cdcWastewaterURL = \"https://www.cdc.gov/wcms/vizdata/NCEZID_DIDRI/NWSSRegionalLevel.json\"","pinned":true,"hidden":false,"mode":"ojs"},{"id":334,"value":"urlNSSP = \"https://data.cdc.gov/api/v3/views/rdmq-nq56/query.csv?app_token=\"+appToken","pinned":true,"hidden":false,"mode":"ojs"},{"id":376,"value":"appToken = \"PsDHDtxnq8boCzAm6E5n8jYoc\"","pinned":true,"hidden":false,"mode":"ojs"},{"id":374,"value":"urlNSSPdemo = \"https://data.cdc.gov/api/v3/views/7xva-uux8/query.csv?app_token=\"+appToken","pinned":true,"hidden":false,"mode":"ojs"},{"id":379,"value":"nsspDataDemo = d3.csv(getProxyURL(urlNSSPdemo), d3.autoType)","pinned":true,"hidden":false,"mode":"ojs"},{"id":336,"value":"nsspData = d3.csv(getProxyURL(urlNSSP), d3.autoType)","pinned":true,"hidden":false,"mode":"ojs"},{"id":353,"value":"d3.extent(nsspData, d=>d.week_end)","pinned":true,"hidden":false,"mode":"ojs"},{"id":355,"value":"counties = d3.rollup(stateNSSPdata, v=>d3.extent(v,d=>d.week_end), d=>d.county)","pinned":true,"hidden":false,"mode":"ojs"},{"id":361,"value":"states = d3.rollup(nsspData, v=>d3.extent(v,d=>d.week_end), d=>d.geography)","pinned":true,"hidden":false,"mode":"ojs"},{"id":387,"value":"statesDemo = d3.rollup(nsspDataDemo, v=>d3.extent(v,d=>d.week_end), d=>d.geography)","pinned":true,"hidden":false,"mode":"ojs"},{"id":346,"value":"stateNSSPdata = nsspData.filter(d=>d.geography == regionToName(selectedState))","pinned":true,"hidden":false,"mode":"ojs"},{"id":382,"value":"stateNSSPdataDemo = nsspDataDemo.filter(d=>d.geography == regionToName(selectedState))","pinned":true,"hidden":false,"mode":"ojs"},{"id":349,"value":"selectedState","pinned":true,"hidden":false,"mode":"ojs"},{"id":338,"value":"import {getProxyURL} from \"@labvegas/cors\"","pinned":true,"hidden":false,"mode":"ojs"},{"id":341,"value":"apiKeyID = \"87cztuikrrkgzymmik2avc5dy\"","pinned":true,"hidden":false,"mode":"ojs"}],"files":"$W32","stdlib":"1","titling":"auto"}