The following data are from the CDC WONDER Multiple Cause of Death Database. This is a searchable public-use database of aggregated data for all deaths in the United States. Each death record is extracted from the a death certificate and given a primary uncderlying cause of death (UCOD) and up to 19 additional multiple causes of death (MCOD). These are originally assigned by the physician on the death certificate but are recategorized following a multisystems algorithms by the National Center for Health Statistics.
Data were extracted for 1999-2014 for deaths with MCOD diagnoses for septicemia and influenza & pneumonia codes.
- Below we will adhere to the following definitions:
- Septicemia = MCOD codes for A40-A41
- Influenza and Pneumonia = MCOD codes for J09-J18
- Sepsis = MCOD codes for either of the above codes
Citation for Data: Centers for Disease Control and Prevention, National Center for Health Statistics. Multiple Cause of Death 1999-2014 on CDC WONDER Online Database, released 2015. Data are from the Multiple Cause of Death Files, 1999-2014, as compiled from data provided by the 57 vital statistics jurisdictions through the Vital Statistics Cooperative Program. Accessed at http://wonder.cdc.gov/mcd-icd10.html on May 10, 2016 10:17:01 AM
library(dplyr)
library(ggplot2)
library(knitr)
kable(septicemia, type="markdown")
age.group | Race | MCOD | deaths | population | crude.rate | L95 | U95 | SE |
---|---|---|---|---|---|---|---|---|
15-24 years | American Indian or Alaska Native | Septicemia (A40-A41) | 185 | 10787577 | 1.7 | 1.5 | 2.0 | 0.1 |
15-24 years | Asian or Pacific Islander | Septicemia (A40-A41) | 427 | 36154211 | 1.2 | 1.1 | 1.3 | 0.1 |
15-24 years | Black or African American | Septicemia (A40-A41) | 3020 | 107822451 | 2.8 | 2.7 | 2.9 | 0.1 |
15-24 years | White | Septicemia (A40-A41) | 8156 | 522206951 | 1.6 | 1.5 | 1.6 | 0.0 |
25-34 years | American Indian or Alaska Native | Septicemia (A40-A41) | 517 | 9053351 | 5.7 | 5.2 | 6.2 | 0.3 |
25-34 years | Asian or Pacific Islander | Septicemia (A40-A41) | 849 | 43006670 | 2.0 | 1.8 | 2.1 | 0.1 |
25-34 years | Black or African American | Septicemia (A40-A41) | 7242 | 91289236 | 7.9 | 7.8 | 8.1 | 0.1 |
25-34 years | White | Septicemia (A40-A41) | 17173 | 504875354 | 3.4 | 3.4 | 3.5 | 0.0 |
35-44 years | American Indian or Alaska Native | Septicemia (A40-A41) | 1303 | 8480501 | 15.4 | 14.5 | 16.2 | 0.4 |
35-44 years | Asian or Pacific Islander | Septicemia (A40-A41) | 1731 | 39700779 | 4.4 | 4.2 | 4.6 | 0.1 |
35-44 years | Black or African American | Septicemia (A40-A41) | 18454 | 91137461 | 20.2 | 20.0 | 20.5 | 0.1 |
35-44 years | White | Septicemia (A40-A41) | 46627 | 544960014 | 8.6 | 8.5 | 8.6 | 0.0 |
45-54 years | American Indian or Alaska Native | Septicemia (A40-A41) | 2669 | 7356880 | 36.3 | 34.9 | 37.7 | 0.7 |
45-54 years | Asian or Pacific Islander | Septicemia (A40-A41) | 3880 | 32492694 | 11.9 | 11.6 | 12.3 | 0.2 |
45-54 years | Black or African American | Septicemia (A40-A41) | 44571 | 83283474 | 53.5 | 53.0 | 54.0 | 0.3 |
45-54 years | White | Septicemia (A40-A41) | 129274 | 553220646 | 23.4 | 23.2 | 23.5 | 0.1 |
55-64 years | American Indian or Alaska Native | Septicemia (A40-A41) | 3459 | 4611773 | 75.0 | 72.5 | 77.5 | 1.3 |
55-64 years | Asian or Pacific Islander | Septicemia (A40-A41) | 7382 | 22253840 | 33.2 | 32.4 | 33.9 | 0.4 |
55-64 years | Black or African American | Septicemia (A40-A41) | 69558 | 55389764 | 125.6 | 124.6 | 126.5 | 0.5 |
55-64 years | White | Septicemia (A40-A41) | 246795 | 432707999 | 57.0 | 56.8 | 57.3 | 0.1 |
65-74 years | American Indian or Alaska Native | Septicemia (A40-A41) | 3693 | 2281495 | 161.9 | 156.6 | 167.1 | 2.7 |
65-74 years | Asian or Pacific Islander | Septicemia (A40-A41) | 11699 | 12631311 | 92.6 | 90.9 | 94.3 | 0.9 |
65-74 years | Black or African American | Septicemia (A40-A41) | 84655 | 31286246 | 270.6 | 268.8 | 272.4 | 0.9 |
65-74 years | White | Septicemia (A40-A41) | 385175 | 283869849 | 135.7 | 135.3 | 136.1 | 0.2 |
75-84 years | American Indian or Alaska Native | Septicemia (A40-A41) | 3410 | 998396 | 341.5 | 330.1 | 353.0 | 5.8 |
75-84 years | Asian or Pacific Islander | Septicemia (A40-A41) | 16437 | 6346595 | 259.0 | 255.0 | 262.9 | 2.0 |
75-84 years | Black or African American | Septicemia (A40-A41) | 95561 | 16533737 | 578.0 | 574.3 | 581.6 | 1.9 |
75-84 years | White | Septicemia (A40-A41) | 601555 | 183946653 | 327.0 | 326.2 | 327.9 | 0.4 |
85+ years | American Indian or Alaska Native | Septicemia (A40-A41) | 1952 | 301745 | 646.9 | 618.2 | 675.6 | 14.6 |
85+ years | Asian or Pacific Islander | Septicemia (A40-A41) | 13600 | 1982010 | 686.2 | 674.6 | 697.7 | 5.9 |
85+ years | Black or African American | Septicemia (A40-A41) | 72733 | 5796853 | 1254.7 | 1245.6 | 1263.8 | 4.7 |
85+ years | White | Septicemia (A40-A41) | 546878 | 72489228 | 754.4 | 752.4 | 756.4 | 1.0 |
g0<-
ggplot(septicemia, aes(age.group, crude.rate, group=Race, color=Race))+
scale_y_continuous(breaks=seq(0, 1500, 100))+
geom_line(size=.5) +
geom_point()+
geom_ribbon(aes(ymin=L95, ymax=U95), alpha=0.2)+
theme_classic()+
ggtitle("Septicemia Mortality (A40-A41) by Age and Race, USA 1999-2014")+
labs(x="Age Group", y = "Crude Mortality per 100,000")
print(g0)
g01<-
ggplot(septicemia2, aes(age.group, RR, group=Race, color=Race))+
geom_line(size=.5) +
geom_point()+
geom_line(aes(y=1), color='black', linetype="dashed")+
theme_classic()+
ggtitle("Septicemia Mortality (A40-A41) by Age and Race, USA 1999-2014")+
labs(x="Age Group", y = "Risk Ratio to each Age Group's Average")
print(g01)
g1<-
ggplot(pneumonia, aes(age.group, crude.rate, group=Race, color=Race))+
scale_y_continuous(breaks=seq(0, 1600, 100))+
geom_line(size=.5) +
geom_point()+
geom_ribbon(aes(ymin=L95, ymax=U95), alpha=0.2)+
theme_classic()+
ggtitle("Influenza & Pneumonia Mortality (J09-J18) by Age and Race, USA 1999-2014")+
labs(x="Age Group", y = "Crude Mortality per 100,000")
print(g1)
g2 <-
ggplot(pneumonia2, aes(age.group, RR, group=Race, color=Race))+
geom_line(size=.5) +
geom_point()+
geom_line(aes(y=1), color='black', linetype="dashed")+
theme_classic()+
ggtitle("Influenza & Pneumonia Mortality (J09-J18) by Age and Race, USA 1999-2014")+
labs(x="Age Group", y = "Risk Ratio to each Age Group's Average")
print(g2)
kable(s2, type="markdown")
age.group | Race | deaths | population | crude.rate | L95 | U95 | SE | RR1 | RR2 |
---|---|---|---|---|---|---|---|---|---|
15-24 years | American Indian or Alaska Native | 334 | 10787577 | 3.1 | 2.8 | 3.4 | 0.2 | 1.00 | 1.63 |
15-24 years | Asian or Pacific Islander | 677 | 36154211 | 1.9 | 1.7 | 2.0 | 0.1 | 0.61 | 1.00 |
15-24 years | Black or African American | 4749 | 107822451 | 4.4 | 4.3 | 4.5 | 0.1 | 1.42 | 2.32 |
15-24 years | White | 15674 | 522206951 | 3.0 | 3.0 | 3.0 | 0.0 | 0.97 | 1.58 |
25-34 years | American Indian or Alaska Native | 800 | 9053351 | 8.8 | 8.2 | 9.4 | 0.3 | 1.19 | 3.03 |
25-34 years | Asian or Pacific Islander | 1257 | 43006670 | 2.9 | 2.8 | 3.1 | 0.1 | 0.39 | 1.00 |
25-34 years | Black or African American | 10954 | 91289236 | 12.0 | 11.8 | 12.2 | 0.1 | 1.62 | 4.14 |
25-34 years | White | 29761 | 504875354 | 5.9 | 5.8 | 6.0 | 0.0 | 0.80 | 2.03 |
35-44 years | American Indian or Alaska Native | 1969 | 8480501 | 23.2 | 22.2 | 24.2 | 0.5 | 1.27 | 3.68 |
35-44 years | Asian or Pacific Islander | 2505 | 39700779 | 6.3 | 6.1 | 6.6 | 0.1 | 0.35 | 1.00 |
35-44 years | Black or African American | 27020 | 91137461 | 29.6 | 29.3 | 30.0 | 0.2 | 1.62 | 4.70 |
35-44 years | White | 75420 | 544960014 | 13.8 | 13.7 | 13.9 | 0.1 | 0.76 | 2.19 |
45-54 years | American Indian or Alaska Native | 3795 | 7356880 | 51.6 | 49.9 | 53.2 | 0.8 | 1.14 | 2.93 |
45-54 years | Asian or Pacific Islander | 5723 | 32492694 | 17.6 | 17.2 | 18.1 | 0.2 | 0.39 | 1.00 |
45-54 years | Black or African American | 62706 | 83283474 | 75.3 | 74.7 | 75.9 | 0.3 | 1.67 | 4.28 |
45-54 years | White | 201500 | 553220646 | 36.4 | 36.3 | 36.6 | 0.1 | 0.80 | 2.07 |
55-64 years | American Indian or Alaska Native | 5150 | 4611773 | 111.7 | 108.6 | 114.7 | 1.6 | 1.05 | 2.24 |
55-64 years | Asian or Pacific Islander | 11102 | 22253840 | 49.9 | 49.0 | 50.8 | 0.5 | 0.47 | 1.00 |
55-64 years | Black or African American | 96046 | 55389764 | 173.4 | 172.3 | 174.5 | 0.6 | 1.63 | 3.47 |
55-64 years | White | 388578 | 432707999 | 89.8 | 89.5 | 90.1 | 0.1 | 0.85 | 1.80 |
65-74 years | American Indian or Alaska Native | 6081 | 2281495 | 266.5 | 259.8 | 273.2 | 3.4 | 1.02 | 1.74 |
65-74 years | Asian or Pacific Islander | 19355 | 12631311 | 153.2 | 151.1 | 155.4 | 1.1 | 0.58 | 1.00 |
65-74 years | Black or African American | 122361 | 31286246 | 391.1 | 388.9 | 393.3 | 1.1 | 1.49 | 2.55 |
65-74 years | White | 676174 | 283869849 | 238.2 | 237.6 | 238.8 | 0.3 | 0.91 | 1.55 |
75-84 years | American Indian or Alaska Native | 6725 | 998396 | 673.6 | 657.5 | 689.7 | 8.2 | 0.96 | 1.28 |
75-84 years | Asian or Pacific Islander | 33370 | 6346595 | 525.8 | 520.2 | 531.4 | 2.9 | 0.75 | 1.00 |
75-84 years | Black or African American | 150522 | 16533737 | 910.4 | 905.8 | 915.0 | 2.3 | 1.30 | 1.73 |
75-84 years | White | 1286803 | 183946653 | 699.6 | 698.3 | 700.8 | 0.6 | 1.00 | 1.33 |
85+ years | American Indian or Alaska Native | 5214 | 301745 | 1727.9 | 1681.0 | 1774.9 | 23.9 | 0.87 | 1.00 |
85+ years | Asian or Pacific Islander | 36424 | 1982010 | 1837.7 | 1818.9 | 1856.6 | 9.6 | 0.93 | 1.06 |
85+ years | Black or African American | 130921 | 5796853 | 2258.5 | 2246.3 | 2270.7 | 6.2 | 1.14 | 1.31 |
85+ years | White | 1538209 | 72489228 | 2122.0 | 2118.6 | 2125.3 | 1.7 | 1.07 | 1.23 |
We will plot crude rate by age-group and race.
g3 <-
ggplot(s2, aes(age.group, crude.rate, group=Race, color=Race))+
scale_y_continuous(breaks=seq(0, 2500, 200))+
geom_line(size=.5) +
geom_point()+
geom_ribbon(aes(ymin=L95, ymax=U95), alpha=0.2)+
theme_classic()+
ggtitle("Sepsis Mortality by Age and Race, USA 1999-2014")+
labs(x="Age Group", y = "Crude Mortality per 100,000")
print(g3)
Now we will plot the risk ratio of each race within each age-group. Here we will make the referent within each age-group the average crude rate of the the age-group.
g4 <-
ggplot(s2, aes(age.group, RR1, group=Race, color=Race))+
geom_line(size=.5) +
geom_point()+
geom_line(aes(y=1), color='black', linetype="dashed")+
theme_classic()+
ggtitle("Sepsis by Age and Race, USA 1999-2014")+
labs(x="Age Group", y = "Risk Ratio to each Age Group's Average")
print(g4)
Now we will plot the risk ratio of each race within each age-group. Here we will make the referent within each age-group the minimum crude rate of the the age-group. This is the Asian/Pacific Islander group for every age-group except the >/= 85 years-old age-group.
g5<-
ggplot(s2, aes(age.group, RR2, group=Race, color=Race))+
geom_line(size=.5) +
geom_point()+
theme_classic()+
theme(legend.justification = c(1,1), legend.position=c(1,1),
legend.title = element_blank(), text = element_text(size = 18))+
ggtitle("Sepsis Mortality Risk by Age and Race, USA 1999-2014")+
labs(x="Age Group", y ="Risk Ratio")
print(g5)