STREAM-SMR: Sequential Bayesian state-space monitoring of standardized mortality ratios - a simulation comparison with risk-adjusted CUSUM and EWMA
Ohno, K.
Show abstract
Background. Sequential monitoring of risk-adjusted mortality in intensive care typically relies on alarm-generating control charts - the risk-adjusted CUSUM, EWMA, or VLAD. These charts signal deterioration but do not return what clinicians and registry stewards ultimately need to interpret: a calibrated, continuously updated estimate of the standardized mortality ratio (SMR) itself. Methods. STREAM-SMR is a conjugate gamma-Poisson dynamic generalized linear model in which the latent log-SMR evolves through a discount factor delta and the alarm statistic is the posterior exceedance probability P(SMR > 1). The construction is closed-form, exact for zero-death months, and computationally trivial at registry scale. Under a protocol frozen before any evaluation runs and calibrated to published national ICU registry aggregates, we compared STREAM-SMR (delta in {0.90, 0.95, 0.97}) with the risk-adjusted CUSUM and risk-adjusted EWMA across three facility-volume strata (50, 200, and 800 annual admissions) and five change scenarios (sustained steps, gradual drift, transient deterioration, and improvement). All methods were Monte-Carlo-calibrated to a common 5% false-alarm probability over a 60-month horizon, with 1,000 replications per cell. Results. At delta = 0.90, STREAM-SMR matched the detection frontier of the risk-adjusted CUSUM to within one to two months across sustained-shift scenarios - median delay for an SMR step to 1.5 of 6 versus 5 months in large facilities, 14 versus 14 in medium, and 20 versus 22 in small - while returning filtered SMR estimates whose 95% credible intervals held at least 91% empirical coverage in every scenario-stratum cell. Empirical false-alarm probabilities were close to the 5% nominal target for all methods (range 0.035-0.066). For a three-month transient deterioration, detection was faster with STREAM-SMR conditional on occurring, but overall detection probability favored the CUSUM in large facilities. Conclusions. STREAM-SMR unifies monitoring and estimation in a single Bayesian object: for a detection-delay premium of at most one to two months against the theoretically optimal CUSUM, it returns an interpretable, uncertainty-quantified SMR trajectory at every time point. The discount factor is an explicit dial between estimate smoothness and detection speed. Simulation code and the frozen protocol are publicly archived.
Matching journals
The top 6 journals account for 50% of the predicted probability mass.
Similar papers in this journal
Similar papers in this journal
- Joint Modeling of Longitudinal Biomarker and Survival Outcomes with the Presence of Competing Risk in Nested Case-Control Studies with Application to the TEDDY Microbiome Dataset 91%
- Robustifying genomic classifiers to batch effects via ensemble learning 90%
- High-dimensional Biomarker Identification for Scalable and Interpretable Disease Prediction via Machine Learning Models 90%
Similar papers in this journal
- A Bayesian Susceptible-Infectious-Hospitalized-Ventilated-Recovered Model to Predict Demand for COVID-19 Inpatient Care in a Large Healthcare System 93%
- A scaling approach to estimate the COVID-19 infection fatality ratio from incomplete data 92%
- Was R < 1 before the English lockdowns? On modelling mechanistic detail, causality and inference about Covid-19 91%
Similar papers in this journal
- Probabilistic Cause-of-disease Assignment using Case-control Diagnostic Tests: A Latent Variable Regression Approach 91%
- Penalized reduced rank regression for multi-outcome survival data supports a common metabolic risk score for age-related diseases 91%
- A Stability-Enhanced Lasso Approach for Covariate Selection in Non-Linear Mixed Effect Model 91%
Similar papers in this journal
- Machine Learning Generalizability Across Healthcare Settings: Insights from multi-site COVID-19 screening 92%
- Continuous-Time and Dynamic Suicide Attempt Risk Prediction with Neural Ordinary Differential Equations 92%
- Generation of realistic synthetic data using multimodal neural ordinary differential equations 91%
"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.