This is reproducible code and data for our PNAS Brief Report, Dynamics of racial disparities in all-cause mortality during the COVID-19 pandemic, which uses the CDC WONDER provisional data to estimate monthly age-standardized mortality rates and rate ratios. We use the NCI Joinpoint Regression Program to fit segmented regressions (joinpoint models) to summarize trends in the monthly rate ratios. The full citation is:
Aschmann, Hélène E., Alicia R. Riley, Ruijia Chen, Yea-Hung Chen, Kirsten Bibbins-Domingo, Andrew C. Stokes, M. Maria Glymour, and Mathew V. Kiang. 2022. “Dynamics of Racial Disparities in All-Cause Mortality during the COVID-19 Pandemic.” Proceedings of the National Academy of Sciences 119(40):e2210941119. doi: 10.1073/pnas.2210941119.
As research documenting disparate impacts of COVID-19 by race and ethnicity grows, little attention has been given to dynamics in mortality disparities during the pandemic and whether changes in disparities persist. We estimate age-standardized monthly all-cause mortality in the United States from January 2018 through February 2022 for seven racial/ethnic populations. Using joinpoint regression, we quantify trends in race-specific rate ratios relative to non-Hispanic White mortality to examine the magnitude of pandemic-related shifts in mortality disparities. Prepandemic disparities were stable from January 2018 through February 2020. With the start of the pandemic, relative mortality disadvantages increased for American Indian or Alaska Native (AIAN), Native Hawaiian or other Pacific Islander (NHOPI), and Black individuals, and relative mortality advantages decreased for Asian and Hispanic groups. Rate ratios generally increased during COVID-19 surges, with different patterns in the summer 2021 and winter 2021/2022 surges, when disparities approached prepandemic levels for Asian and Black individuals. However, two populations below age 65 fared worse than White individuals during these surges. For AIAN people, the observed rate ratio reached 2.25 (95% CI = 2.14, 2.37) in October 2021 vs. a prepandemic mean of 1.74 (95% CI = 1.62, 1.86), and for NHOPI people, the observed rate ratio reached 2.12 (95% CI = 1.92, 2.33) in August 2021 vs. a prepandemic mean of 1.31 (95% CI = 1.13, 1.49). Our results highlight the dynamic nature of racial/ethnic disparities in mortality and raise alarm about the exacerbation of mortality inequities for Indigenous groups due to the pandemic.
All analyses are conducted using R
, which can be downloaded via
CRAN. We also recommend the use of
RStudio when
running R
, which will allow users to take advantage of
renv
for dependency
management. If not using renv
, we provide session information at the
end of this file to ensure packages are the correct version. Note that
for the joinpoint analyses, you must use an external program, Joinpoint
Regression Program which
can be downloaded from the National Cancer Institute website.
We use publicly available data — the README.md
file in data_raw
has
links to our CDC WONDER query and each query file has the dates of data
access.
The code files (in code
) are designed to be run in order. Files 01
to 04
manipulate the raw data into files that can be run in the
Joinpoint regression program. After running these files, one must run
joinpoint analyses in Joinpoint Regression Program. For convenience, we
have provided our saved analysis settings (files ending in .jps
) as
well as our model results (files ending in .jpo
). Once files have been
exported out of the Joinpoint Regression Program, one can run code files
05
and higher to gather the results and create plots.
The numeric representations of the figures in our paper are available in
output
.
- Hélène Aschmann
(
: heleneaschmann | (
: @LnAschmann)
- Kirsten
Bibbins-Domingo
(
: @KBibbinsDomingo)
- Ruijia Chen
(
: 2018Ruijia |
: @RuijiaChen19)
- Yea-Hung Chen
(
: yea-hung |
: @ch272n)
- Maria Glymour
(
: @MariaGlymour)
- Mathew Kiang
(
: mkiang |
: @mathewkiang)
- Alicia Riley
(
: aliciarr |
: @aliciacita)
- Andrew Stokes
(
: @astokespop)
> sessioninfo::session_info()
─ Session info ─────────────────────────────────────────────────────────
setting value
version R version 4.1.0 (2021-05-18)
os macOS 12.4
system x86_64, darwin17.0
ui RStudio
language (EN)
collate en_US.UTF-8
ctype en_US.UTF-8
tz America/Los_Angeles
date 2022-07-18
rstudio 2022.02.1+461 Prairie Trillium (desktop)
pandoc 2.17.1.1 @ /Applications/RStudio.app/Contents/MacOS/quarto/bin/ (via rmarkdown)
─ Packages ─────────────────────────────────────────────────────────────
package * version date (UTC) lib source
assertthat 0.2.1 2019-03-21 [2] CRAN (R 4.1.0)
backports 1.4.1 2021-12-13 [2] CRAN (R 4.1.0)
broom 0.8.0 2022-04-13 [2] CRAN (R 4.1.2)
cellranger 1.1.0 2016-07-27 [2] CRAN (R 4.1.0)
cli 3.3.0 2022-04-25 [2] CRAN (R 4.1.2)
colorspace 2.0-3 2022-02-21 [2] CRAN (R 4.1.2)
crayon 1.5.1 2022-03-26 [2] CRAN (R 4.1.2)
DBI 1.1.2 2021-12-20 [2] CRAN (R 4.1.0)
dbplyr 2.1.1 2021-04-06 [2] CRAN (R 4.1.0)
digest 0.6.29 2021-12-01 [2] CRAN (R 4.1.0)
dplyr * 1.0.9 2022-04-28 [2] CRAN (R 4.1.2)
ellipsis 0.3.2 2021-04-29 [2] CRAN (R 4.1.0)
evaluate 0.15 2022-02-18 [2] CRAN (R 4.1.2)
fansi 1.0.3 2022-03-24 [2] CRAN (R 4.1.2)
fastmap 1.1.0 2021-01-25 [2] CRAN (R 4.1.0)
forcats * 0.5.1 2021-01-27 [2] CRAN (R 4.1.0)
fs 1.5.2 2021-12-08 [2] CRAN (R 4.1.0)
generics 0.1.2 2022-01-31 [2] CRAN (R 4.1.2)
ggplot2 * 3.3.6 2022-05-03 [2] CRAN (R 4.1.2)
glue 1.6.2 2022-02-24 [2] CRAN (R 4.1.2)
gtable 0.3.0 2019-03-25 [2] CRAN (R 4.1.0)
haven 2.5.0 2022-04-15 [2] CRAN (R 4.1.2)
here * 1.0.1 2020-12-13 [2] CRAN (R 4.1.0)
hms 1.1.1 2021-09-26 [2] CRAN (R 4.1.0)
htmltools 0.5.2 2021-08-25 [2] CRAN (R 4.1.0)
httr 1.4.3 2022-05-04 [2] CRAN (R 4.1.2)
janitor * 2.1.0 2021-01-05 [2] CRAN (R 4.1.0)
jsonlite 1.8.0 2022-02-22 [2] CRAN (R 4.1.2)
knitr 1.39 2022-04-26 [2] CRAN (R 4.1.2)
lattice 0.20-45 2021-09-22 [2] CRAN (R 4.1.0)
lifecycle 1.0.1 2021-09-24 [2] CRAN (R 4.1.0)
lubridate 1.8.0 2021-10-07 [2] CRAN (R 4.1.0)
magrittr 2.0.3 2022-03-30 [2] CRAN (R 4.1.2)
modelr 0.1.8 2020-05-19 [2] CRAN (R 4.1.0)
munsell 0.5.0 2018-06-12 [2] CRAN (R 4.1.0)
narcan * 0.1.1 2022-03-28 [2] Github (mkiang/narcan@a9c3d53)
pillar 1.7.0 2022-02-01 [2] CRAN (R 4.1.2)
pkgconfig 2.0.3 2019-09-22 [2] CRAN (R 4.1.0)
purrr * 0.3.4 2020-04-17 [2] CRAN (R 4.1.0)
R6 2.5.1 2021-08-19 [2] CRAN (R 4.1.0)
readr * 2.1.2 2022-01-30 [2] CRAN (R 4.1.2)
readxl 1.4.0 2022-03-28 [2] CRAN (R 4.1.2)
renv 0.15.5 2022-05-26 [1] CRAN (R 4.1.2)
reprex 2.0.1 2021-08-05 [2] CRAN (R 4.1.0)
rlang 1.0.2 2022-03-04 [2] CRAN (R 4.1.2)
rmarkdown 2.14 2022-04-25 [2] CRAN (R 4.1.2)
rprojroot 2.0.3 2022-04-02 [2] CRAN (R 4.1.2)
rsconnect 0.8.25 2021-11-19 [2] CRAN (R 4.1.0)
rstudioapi 0.13 2020-11-12 [2] CRAN (R 4.1.0)
rvest 1.0.2 2021-10-16 [2] CRAN (R 4.1.0)
scales 1.2.0 2022-04-13 [2] CRAN (R 4.1.2)
sessioninfo 1.2.2 2021-12-06 [2] CRAN (R 4.1.0)
snakecase 0.11.0 2019-05-25 [2] CRAN (R 4.1.0)
stringi 1.7.6 2021-11-29 [2] CRAN (R 4.1.0)
stringr * 1.4.0 2019-02-10 [2] CRAN (R 4.1.0)
tibble * 3.1.7 2022-05-03 [2] CRAN (R 4.1.2)
tidyr * 1.2.0 2022-02-01 [2] CRAN (R 4.1.2)
tidyselect 1.1.2 2022-02-21 [2] CRAN (R 4.1.2)
tidyverse * 1.3.1 2021-04-15 [2] CRAN (R 4.1.0)
tzdb 0.3.0 2022-03-28 [2] CRAN (R 4.1.0)
utf8 1.2.2 2021-07-24 [2] CRAN (R 4.1.0)
vctrs 0.4.1 2022-04-13 [2] CRAN (R 4.1.2)
withr 2.5.0 2022-03-03 [2] CRAN (R 4.1.2)
xfun 0.31 2022-05-10 [2] CRAN (R 4.1.2)
xml2 1.3.3 2021-11-30 [2] CRAN (R 4.1.0)
yaml 2.3.5 2022-02-21 [2] CRAN (R 4.1.2)
zoo 1.8-10 2022-04-15 [2] CRAN (R 4.1.2)
[1] /Users/mvk/Library/CloudStorage/Dropbox/Projects_Active/dynamic_inequality/renv/library/R-4.1/x86_64-apple-darwin17.0
[2] /Library/Frameworks/R.framework/Versions/4.1/Resources/library
────────────────────────────────────────────────────────────────────────