REVIEW 3 major objections 4 minor 35 references
A fast approach for analyzing spatio-temporal patterns in ischemic heart disease mortality across US counties (1999-2021)
T0 review · 3 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read Using 23 years of county-level data for every contiguous US county, this paper claims the long decline in ischemic heart disease mortality slowed after 2009, flattened after 2014, and turned into a slight rise after 2019.
desk verdict Solid, reproducible bigDM application; honest about limits, but the headline change-point claim rests on single imputation that could create the pattern. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central mechanism is a hierarchical Bayesian spatio-temporal disease mapping model. Death counts follow $Y_{it} \sim \text{Poisson}(e_{it} r_{it})$, and the log relative risk is decomposed as $\log r_{it} = \alpha_0 + \xi_i + \gamma_t + \delta_{it}$, with $\xi_i$ a spatially structured county effect, $\gamma_t$ a first-order random walk over years, and $\delta_{it}$ a Type II space-time interaction (structured in time, unstructured in space). To handle 3105 counties over 23 years, the paper replaces one global fit with a 'divide and conquer' scheme: partition the US into state-based subregions, add first-order neighboring counties across borders, fit the model separately in each partition, and allow different smoothing parameters per partition. Before fitting, suppressed counts (values below 10) are imputed with rounded, truncated predictions from a per-year Bayesian spatial model, so the smoothing machinery operates on a complete data grid.
What would settle it
Fit the same model to the complete, un-suppressed death records for a subset of counties; if the post-2014 flattening or the post-2019 uptick disappears, the reported trajectory is an artifact of imputation.
Extended reading notes
Core claim
The paper's central claim is that the risk of dying from ischemic heart disease in US counties followed a four-phase national trajectory over 1999-2021: a clear decline to about 2009, a slowing decline from 2009 to 2014, a flat trend from 2014 to 2019, and an increase after 2019 with a peak in 2020. Beneath that national curve, no single trend describes the country. The West remains below the national risk, the Northeast above it; rural counties carry higher risk than urban counties in the West, Midwest, and South; and counties such as Valley (Montana) and Menominee (Michigan) show rising trends even while the national risk falls. The authors identify 2009 and 2019 as candidate change points while cautioning that the two post-2019 years cannot separate a persistent increase from a transient pandemic effect.
Load-bearing premise
The load-bearing premise is that the roughly 12.5% of county-year death counts suppressed because they are below 10 can be replaced by spatially predicted, rounded, and truncated values without systematically distorting the smoothed trends; the paper itself concedes the results may be sensitive to this imputation.
Editorial extensions
If this is right
- If the trajectory is right, the decades-long US decline in heart disease mortality has levelled off, so public health planning should target a stalled epidemic rather than assume continuation of historical improvements.
- The 2009 and 2019 change points give surveillance systems concrete moments to monitor, especially whether the post-2019 increase persists outside the pandemic years.
- Rural counties in the West, Midwest, and South emerge as the highest-risk areas, reinforcing that national campaigns need county-level and rural-specific strategies.
- The partitioned fitting strategy makes it practical to screen spatio-temporal disease patterns for all US counties in minutes on a desktop, not just for heart disease but for any condition with similar data.
Reading between the lines
- A consequence the authors leave implicit is that a real post-2014 flattening strengthens the case that cardiovascular mortality is the main brake on US life-expectancy gains, since ischemic heart disease is the leading cause of death.
- The imputation step could be audited by fitting the same pipeline to complete, un-suppressed death records for a subset of states; if the flattening or the post-2019 uptick disappears, the finding would be an artifact of filling in small counts.
- Extending the series through 2023 and later would settle whether the post-2019 increase is persistent or a one-off COVID-19 spike, something the paper's data ending in 2021 cannot resolve.
- The same partitioned workflow transfers to other suppressed-count chronic diseases such as stroke or heart failure, where county-level screening could reuse the imputation and smoothing machinery.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript analyzes county-level ischaemic heart disease mortality in the 48 contiguous US states over 1999-2021 using CDC WONDER data, in which counts below 10 are suppressed. Missing counts are imputed by fitting a year-specific BYM2 spatial model and replacing each censored count with the rounded posterior median predicted count truncated at 9. The imputed data are then modeled with a partition-based 'divide and conquer' spatio-temporal model (bigDM) with ICAR spatial prior, RW1 temporal prior, and Type II interaction. The main finding is that the US decline in IHD mortality slowed after 2014 and there was a slight increase after 2019, with regional and urban/rural disparities. The paper also reports model selection favoring the partitioned ICAR/Type II model over alternatives and provides reproducible code and data.
Significance. If correct, the paper offers a computationally practical pipeline for Bayesian spatio-temporal disease mapping on a large domain with censored counts, and its descriptive finding of a flattening IHD mortality trend is consistent with prior literature (Mehta et al.; Vaughan et al.). The full reproducibility (code and data on GitHub) and the comparison with global models are strengths. The main risk is that the central temporal claim is identified only after a single imputation of 12.5% of the observations, and the paper's own discussion acknowledges that the results 'may be sensitive to this method' and that the post-2019 rise cannot be separated from COVID-19 with only two years of data. Given the time-varying missingness, this is a load-bearing limitation rather than a presentation issue.
major comments (3)
- [Section 2.2 and Eq. (2)] The imputation treats synthetic counts as observed, and no imputation uncertainty is propagated into the posterior credible intervals in Figures 5-7 or into the DIC/WAIC values in Table 2. Because suppression is deterministic (count < 10) and the fraction of suppressed cells increases over time as IHD mortality falls, the imputation error is time-varying and can distort the very slopes used to identify the post-2014 deceleration and post-2019 uptick. The rounded posterior medians truncated at 9 are not draws from the posterior predictive distribution, so the headline change-point claim could be an artifact of the imputation scheme. A censoring-appropriate reanalysis (e.g., multiple imputation, a censored likelihood as in Valeriano et al. 2021, or a sensitivity analysis excluding counties with any imputed counts) is needed before the central claim can be accepted.
- [Section 3, Figure 7] The paper's sanity check for imputation effects, based on comparing model-based rural trends with observed-only SMRs, only checks aggregate levels; it does not target the post-2014 slope or the post-2019 increase. At the aggregation level used, a time-varying bias in imputed counts can cancel or be masked, so this check does not address whether the change-point pattern is an artifact. The check should be repeated for the slope estimates, for example by comparing the model's year-to-year trend with an observed-only SMR trend restricted to counties without any suppression.
- [Abstract, Section 3, Discussion] The claimed change-point dates are stated inconsistently. The abstract says 'slowdown in the decrease of IHD mortality in the US after 2014'; the text near Figure 6 says 'a slowdown in the decrease is observed after 2009 and a flattening in the trend is observed after 2014'; and the Discussion says 'our analysis suggests 2009 and 2019 as potential change points.' Because these dates are the central result, the paper needs to state a single, precisely defined claim and either justify the date or formally test for a change point rather than relying on 'mere visual inspection.'
minor comments (4)
- [Section 3, Figure 4] The text 'from 1991 to 2021' should read 'from 1999 to 2021'.
- [Section 2.2] The phrase 'As suggested by the reviewers' is an artifact of the review process and should be removed or rephrased.
- [Abstract] The phrase 'differences exists' should be 'differences exist'.
- [Section 3] The expression 'the starting of the COVID-19 pandemic' should be 'the start of the COVID-19 pandemic', and the sentence 'From this year onwards, the average risks in rural areas is the highest' needs grammatical correction.
Circularity Check
No circular derivation; the central trends are posterior summaries of a standard spatio-temporal model, not fitted inputs relabeled as predictions.
full rationale
The paper's central claim (deceleration after 2014 and increase after 2019) is a descriptive summary of the posterior temporal trend from model (2), a standard Poisson ICAR/BYM2 + RW1 + Type II interaction model. No parameter of that model is defined in terms of the target trend, and the trend is not a fitted quantity that is then relabeled as a prediction; the change points are identified by visual inspection of the fitted curves, and previous literature is cited for corroboration, not as the source of the estimate. The imputation step (Section 2.2) does replace suppressed counts with spatial-model predictions, and the paper explicitly warns that the results may be sensitive to this method; however, the imputation model contains no temporal trend component, so the post-2014 flattening and post-2019 uptick are not imposed by construction. Any bias from truncation at 9 is a robustness or identifiability issue, not circularity. Self-citations to the bigDM package and the divide-and-conquer papers provide computational machinery and are not used to justify the epidemiological conclusion, so they are not load-bearing. Overall, no step in the derivation reduces to its own input.
Assumptions & free parameters
free parameters (2)
- truncation threshold for imputed counts =
9
- neighbor order k for border buffers =
1
assumptions (4)
- standard math Poisson likelihood for county-year death counts
- domain assumption ICAR/RW1/Type II interaction priors adequately capture residual spatio-temporal structure
- domain assumption Missing counts can be imputed using spatial borrowing of strength from neighboring counties
- domain assumption Expected counts computed with a single national rate make risks comparable across counties and years
Cite this review
Pith. "Pith review of A fast approach for analyzing spatio-temporal patterns in ischemic heart disease mortality across US counties (1999-2021)." pith.science (2026). https://pith.science/paper/I7AEM4FO
@misc{pith2026241114849,
author = {Pith},
title = {Pith review of: A fast approach for analyzing spatio-temporal patterns in ischemic heart disease mortality across US counties (1999-2021)},
year = {2026},
howpublished = {\url{https://pith.science/paper/I7AEM4FO}},
note = {Machine review of arXiv:2411.14849}
}
read the original abstract
Ischaemic heart disease (IHD) remains the primary cause of mortality in the US. This study focuses on using spatio-temporal disease mapping models to explore the temporal trends of IHD at the county level from 1999 to 2021. To manage the computational burden arising from the high-dimensional data, we employ scalable Bayesian models using a "divide and conquer" strategy. This approach allows for fast model fitting and serves as an efficient procedure for screening spatio-temporal patterns. Additionally, we analyze trends in four regional subdivisions, West, Midwest, South and Northeast, and in urban and rural areas. The dataset on IHD contains missing data, and we propose a procedure to impute the omitted information. The results show a slowdown in the decrease of IHD mortality in the US after 2014 with a slight increase noted after 2019. However, differences exists among the counties, the four big geographical regions, and rural and urban areas.
Figures
Figures from the paper (4 more)
Reference graph
Works this paper leans on
-
[1]
Berrett, C., Gurney, B., Arthur, D., Moon, T., & Williams, G. P. (2021). A bayesian change point model for spatio-temporal data
work page 2021
-
[2]
Besag, J. (1974). Spatial interaction and the statistical analysis of lattice systems (with discussion) . Journal of the Royal Statistical Society: Series B (Statistical Methodology) , 36, 192--236. https://doi.org/10.1111/j.2517-6161.1974.tb00999.x
arXiv 1974
-
[3]
Cooper, R., Cutler, J., Desvigne-Nickens, P., Fortmann, S. P., Friedman, L., Havlik, R., Hogelin, G., Marler, J., McGovern, P., Morosco, G., Mosca, L., Pearson, T., Stamler, J., Stryer, D., & Thom, T. (2000). Trends and disparities in coronary heart disease, stroke, and other cardiovascular diseases in the united states. Circulation , 102(25), 3137--3147....
-
[4]
Ford, E. S., Ajani, U. A., Croft, J. B., Critchley, J. A., Labarthe, D. R., Kottke, T. E., Giles, W. H., & Capewell, S. (2007). Explaining the decrease in u.s. deaths from coronary disease, 1980–2000. New England Journal of Medicine , 356(23), 2388--2398. https://doi.org/10.1056/NEJMsa053935
-
[5]
D., Etxeberria, J., & Militino, A
Goicoa, T., Ugarte, M. D., Etxeberria, J., & Militino, A. F. (2012). Comparing CAR and P-spline models in spatial disease mapping . Environmental and Ecological Statistics , 19, 573--599
work page 2012
-
[6]
Hall, E. W., Vaughan, A. S., Ritchey, M. D., Schieb, L., & Casper, M. (2019). Stagnating national declines in stroke mortality mask widespread county-level increases, 2010–2016. Stroke , 50(12), 3355--3359. https://doi.org/10.1161/STROKEAHA.119.026695
-
[7]
Heaton, M. J., Datta, A., Finley, A. O., Furrer, R., Guinness, J., Guhaniyogi, R., Gerber, F., Gramacy, R. B., Hammerling, D., Katzfuss, M., Lindgren, F., Nychka, D. W., Sun, F., & Zammit-Mangion, A. (2019). A case study competition among methods for analyzing large spatial data. Journal of Agricultural, Biological and Environmental Statistics , 24, 398--...
-
[8]
Knorr-Held, L. (2000). Bayesian modelling of inseparable space-time variation in disease risk. Statistics in Medicine , 19(17-18), 2555--2567
2000
Show all 35 references
-
[9]
Kulshreshtha, A., Goyal, A., Dabhadkar, K., Veledar, E., & Vaccarino, V. (2014). Urban-rural differences in coronary heart disease mortality in the united states: 1999–2009. Public Health Reports , 129(1), 19--29. https://doi.org/10.1177/003335491412900105
2014 doi
-
[10]
& Rue, H
Lindgren, F. & Rue, H. (2015). Bayesian Spatial Modelling with R - INLA . Journal of Statistical Software , 63(19), 1--25
2015
-
[11]
K., Abrams, L
Mehta, N. K., Abrams, L. R., & Myrskylä, M. (2020). Us life expectancy stalls due to cardiovascular disease, not drug deaths. Proc Natl Acad Sci U S A. , 31(117(13)), 6998--7000. https://doi.org/10.1073/pnas.1920391117
2020 doi
-
[12]
A., Wei, G
Mensah, G. A., Wei, G. S., Sorlie, P. D., Fine, L. J., Rosenberg, Y., Kaufmann, P. G., Mussolino, M. E., Hsu, L. L., Addou, E., Engelgau, M. M., & Gordon, D. (2017). Decline in cardiovascular mortality: Possible causes and implications. Circulation research , 120(2), 366–380. ...
2017 doi
-
[13]
Murphy, M. J. & Grundy, E. M. D. (2022). Slowdown in mortality improvement in the past decade: A us/uk comparison. The Journals of Gerontology: Series B , 77(Sup2), S138--S147. https://doi.org/10.1093/geronb/gbab220
2022 doi
-
[14]
Niekerk, J. V. & Rue, H. (2023). Low-rank variational bayes correction to the laplace method. https://doi.org/10.48550/arXiv.2111.12945
2023 doi
-
[15]
Orozco-Acosta, E., Adin, A., & Ugarte, M. D. (2021). Scalable bayesian modelling for smoothing disease risks in large spatial data sets using inla. Spatial Statistics , 41, 100496. https://doi.org/10.1016/j.spasta.2021.100496
2021
-
[16]
Orozco-Acosta, E., Adin, A., & Ugarte, M. D. (2023). Big problems in spatio-temporal disease mapping: Methods and software. Computer Methods and Programs in Biomedicine , 231, 107403. https://doi.org/10.1016/j.cmpb.2023.107403
2023
-
[17]
Quick, H. (2019). Estimating County- Level Mortality Rates Using Highly Censored Data From CDC WONDER. Preventing Chronic Disease , 16:180441. https://doi.org/10.5888/pcd16.180441
2019 doi
-
[18]
A., & Casper, M
Quick, H., Waller, L. A., & Casper, M. (2017). Multivariate spatio-temporal modeling of age-specific stroke mortality . The Annals of Applied Statistics , 11(4), 2165--2177. https://doi.org/10.1214/17-AOAS1068
2017 doi
-
[19]
A., & Casper, M
Quick, H., Waller, L. A., & Casper, M. (2018). A Multivariate Space–Time Model for Analysing County Level Heart Disease Death Rates by Race and Sex . Journal of the Royal Statistical Society Series C: Applied Statistics , 67(1), 291--304. https://doi.org/10.1111/rssc.12215
2018 doi
-
[20]
H., Simpson, D., & Rue, H
Riebler, A., Sørbye, S. H., Simpson, D., & Rue, H. (2016). An intuitive bayesian spatial model for disease mapping that accounts for scaling. Statistical Methods in Medical Research , 25(4), 1145--1165. https://doi.org/10.1177/0962280216660421
2016 doi
-
[21]
A., Vaduganathan, M., & Mensah, G
Roth, G. A., Vaduganathan, M., & Mensah, G. A. (2022). Impact of the covid-19 pandemic on cardiovascular health in 2020: Jacc state-of-the-art review. Journal of the American College of Cardiology , 80(6), 631--640. https://doi.org/10.1016/j.jacc.2022.06.008
2022 doi
-
[22]
& Held, L
Rue, H. & Held, L. (2005). Gaussian Markov Random Fields: Theory and Applications . CRC press
2005
-
[23]
D., Goicoa, T., Etxeberria, J., & Militino, A
Ugarte, M. D., Goicoa, T., Etxeberria, J., & Militino, A. F. (2012). A P-spline ANOVA type model in space-time disease mapping . Stochastic Environmental Research and Risk Assesment , 26(6), 835--845
2012
-
[24]
D., Goicoa, T., & Militino, A
Ugarte, M. D., Goicoa, T., & Militino, A. F. (2010). Spatio-temporal modelling of mortality risks using penalized splines . Environmetrics , 21(3--4), 270--289
2010
-
[25]
Valeriano, K. A. L., Lachos, V. H., Prates, M. O., & Matos, L. A. (2021). Likelihood-based inference for spatiotemporal data with censored and missing responses. Environmetrics , 32(3), e2663. https://doi.org/10.1002/env.2663
2021 doi
-
[26]
Van Niekerk , J., Krainski, E., Rustand, D., & Rue, H. (2023). A new avenue for bayesian inference with inla. Computational Statistics & Data Analysis , 181, 107692. https://doi.org/10.1016/j.csda.2023.107692
2023
-
[27]
S., Coronado, F., Casper, M., Loustalot, F., & Wright, J
Vaughan, A. S., Coronado, F., Casper, M., Loustalot, F., & Wright, J. S. (2022a). County‐level trends in hypertension‐related cardiovascular disease mortality—united states, 2000 to 2019. Journal of the American Heart Association , 11(7), e024785. https://doi.org/10.1161/JAHA....
2022 doi
-
[28]
S., Flynn, A., & Casper, M
Vaughan, A. S., Flynn, A., & Casper, M. (2022b). The where of when: Geographic variation in the timing of recent increases in us county-level heart disease death rates. Annals of Epidemiology , 72, 18--24. https://doi.org/10.1016/j.annepidem.2022.05.002
2022 doi
-
[29]
S., George, M
Vaughan, A. S., George, M. G., Jackson, S. L., Schieb, L., & Casper, M. (2021). Changing spatiotemporal trends in county‐level heart failure death rates in the united states, 1999 to 2018. Journal of the American Heart Association , 10(4), e018125. https://doi.org/10.1161/JAHA...
2021 doi
-
[30]
S., Quick, H., Schieb, L., Kramer, M
Vaughan, A. S., Quick, H., Schieb, L., Kramer, M. R., Taylor, H. A., & Casper, M. (2019). Changing rate orders of race-gender heart disease death rates: An exploration of county-level race-gender disparities. SSM - Population Health , 7, 100334. https://doi.org/10.1016/j.ssmph...
2019
-
[31]
S., Ritchey, M
Vaughan, A. S., Ritchey, M. D., Hannan, J., Kramer, M. R., & Casper, M. (2017). Widespread recent increases in county-level heart disease mortality across age groups. Annals of Epidemiology , 27(12), 796--800. https://doi.org/10.1016/j.annepidem.2017.10.012
2017 doi
-
[32]
S., Schieb, L., & Casper, M
Vaughan, A. S., Schieb, L., & Casper, M. (2020). Historic and recent trends in county-level coronary heart disease death rates by race, gender, and age group, united states, 1979-2017. PLoS One , 15(7). https://doi.org/10.1371/journal.pone.0235839
2020 doi
-
[33]
S., Schieb, L., Quick, H., Kramer, M
Vaughan, A. S., Schieb, L., Quick, H., Kramer, M. R., & Casper, M. (2018). Before the here and now: What we can learn from variation in spatiotemporal patterns of changing heart disease mortality by age group, time period, and birth cohort. Social Science & Medicine , 217, 97-...
2018 doi
-
[34]
Vicente, G., Adin, A., Goicoa, T., & Ugarte, M. D. (2023). High-dimensional order-free multivariate spatial disease mapping. Statistics and Computing , 33(5), 104
2023
-
[35]
A., O’Flaherty, M., Capewell, S., Ford, E
Wilmot, K. A., O’Flaherty, M., Capewell, S., Ford, E. S., & Vaccarino, V. (2015). Coronary heart disease mortality declines in the united states from 1979 through 2011: Evidence for stagnation in young adults, especially women. Circulation , 132(11), 997–1002. https://doi.org/...
2015 doi
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.