Biostatistics Advance Access originally published online on February 16, 2006
Biostatistics 2006 7(4):530-550; doi:10.1093/biostatistics/kxj024
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
A Bayesian forecasting model: predicting U.S. male mortality
Division of Biostatistics, University of Texas School of Public Health at Houston, 1200 Herman Pressler Drive, E831, Houston, TX 77030, USA claudia.pedroza{at}uth.tmc.edu
This article presents a Bayesian approach to forecast mortality rates. This approach formalizes the LeeCarter method as a statistical model accounting for all sources of variability. Markov chain Monte Carlo methods are used to fit the model and to sample from the posterior predictive distribution. This paper also shows how multiple imputations can be readily incorporated into the model to handle missing data and presents some possible extensions to the model. The methodology is applied to U.S. male mortality data. Mortality rate forecasts are formed for the period 19901999 based on data from 19591989. These forecasts are compared to the actual observed values. Results from the forecasts show the Bayesian prediction intervals to be appropriately wider than those obtained from the LeeCarter method, correctly incorporating all known sources of variability. An extension to the model is also presented and the resulting forecast variability appears better suited to the observed data.
Keywords: Bayesian prediction; LeeCarter method; Missing data; Mortality forecasting; State-space model