Few things are more anticipated in 2020 than the arrival of an effective and safe COVID vaccine. On that note, I wanted to compare the deaths from COVID to peak deaths from commonly known vaccine preventable diseases prior to vaccines existing for those diseases.
I found a paper in JAMA that looked at the peak deaths from diseases for which we now have vaccines that are recommended as part of routine immunization programs. These diseases include:
*Part of DTaP vaccine
**Part of MMR vaccine
From that JAMA paper I took the peak death counts and converted them to the 2020 US population using data from the US Census Bureau. I did not include Hemophilus influenzae B because no peak death statistics were included.
#COVID deaths from CDC as of July 6, 2020 disease<-c("Diphtheria","Measles","Mumps","Pertussis","Polio, Acute","Polio,Paralytic","Rubella","Cong Rubella","Smallpox","Tetanus","Hep A","Hep B","Pneumo","Varicella","COVID") deaths<-c(3065,552,50,7518,2720,3145,24,2160,2510,511,298,267,7300,138,129811) #populations for 2020 from https://www.census.gov/content/dam/Census/library/publications/2020/demo/p25-1144.pdf and for all other years from https://www.census.gov/population/estimates/nation/popclockest.txt and are all in 100,000s years<-c(1936,1958,1964,1934,1949,1952,1968,1964,1902,1947,1971,1985,1999,1973,2020) pops<-c(1280.5318,1748.81904,1918.88791,1263.73773,1491.8813,1575.5274,2007.06052,1918.88791,791.63,1441.26071,2076.60677,2379.23795,2726.90813,2119.08788,3326.39) deaths_per_pop<-deaths/pops df<-data.frame(disease,deaths,years,pops,deaths_per_pop) #now I lumped a few diseases together to make fewer categories for better visualization; diphtheria, tetanus and pertussis are given together as the DTaP vaccine, measles, mumps and rubella as MMR, and then the two types of polio into one diseases<-c("DTaP","MMR","Polio","Smallpox","Hep A","Hep B","Pneumo","Varicella","COVID") deathrates<-c(df$deaths_per_pop+df$deaths_per_pop+df$deaths_per_pop,df$deaths_per_pop+df$deaths_per_pop+df$deaths_per_pop+df$deaths_per_pop,df$deaths_per_pop+df$deaths_per_pop,df$deaths_per_pop,df$deaths_per_pop,df$deaths_per_pop,df$deaths_per_pop,df$deaths_per_pop,df$deaths_per_pop) #convert to 2020 pop deaths<-ceiling(deathrates*3326.390) df<-data.frame(diseases=diseases,deaths=deaths) df$group<-c(rep("Vaccine-Preventable Illnesses",8),"COVID")
I will arrange the data frame by number of deaths so that when the bar graph is generated, everything looks in order.
From the graph above you can see that were all of these vaccine-preventable illnesses killing Americans at the same rate as they were in their peak pre-vaccine year, they still would only kill roughly half of what COVID has killed so far only HALFWAY through this year.