Exploring spatio-temporal patterns of mortality using mixed effects models.

A linear mixed effects (LME) model previously used for a spatial analysis of mortality data for a single time period is extended to include time trends and spatio-temporal interactions. This model includes functions of age and time period that can account for increasing and decreasing death rates over time and age, and a change-point of rates at a predetermined age. A geographic hierarchy is included that provides both regional and small area age-specific rate estimates, stabilizing rates based on small numbers of deaths by sharing information within a region. The proposed log-linear analysis of rates allows the use of commercially available software for parameter estimation, and provides an estimator of overdispersion directly as the residual variance. Because of concerns about the accuracy of small area rate estimates when there are many instances of no observed deaths, we consider potential sources of error, focusing particularly on the similarity of likelihood inferences using the LME model for rates as compared to an exact Poisson-normal mixed effects model for counts. The proposed LME model is applied to breast cancer deaths which occurred among white women during 1979-1996. For this example, application of diagnostics for multiparameter likelihood comparisons suggests a restriction of age to a minimum of either 25 or 35, depending on whether small area rate estimates are required. Investigation into a convergence problem led to the discovery that the changes in breast cancer geographic patterns over time are related more to urbanization than to region, as previously thought. Published in 2000 by John Wiley & Sons, Ltd.

[1]  D. A. Sprott.  Normal likelihoods and their relation to large sample theory of estimation , 1973 .

[2]  Stephen E. Fienberg,et al.  Discrete Multivariate Analysis: Theory and Practice , 1976 .

[3]  Salomon Minkin.  Assessing the quadratic approximation to the log likelihood function in nonnormal linear models , 1983 .

[4]  Dennis E. Jennings.  Judging Inference Adequacy in Logistic Regression , 1986 .

[5]  David S. Stoffer,et al.  Estimation and Identification of Space-Time ARMAX Models in the Presence of Missing Data , 1986 .

[6]  Linda W. Pickle,et al.  Atlas of U.S. cancer mortality among whites: 1950-1980. , 1987 .

[7]  M. Woodbury,et al.  Empirical Bayes procedures for stabilizing maps of U.S. cancer mortality rates. , 1989, Journal of the American Statistical Association.

[8]  Rural Economy Division.,et al.  Rural-Urban Continuum Codes for Metro and Nonmetro Counties, 1993 , 1990 .

[9]  R. Cook,et al.  Diagnostics for Assessing the Accuracy of Normal Approximations in Exponential Family Nonlinear Models , 1990 .

[10]  J. R. Wallis,et al.  An Approach to Statistical Spatial-Temporal Modeling of Meteorological Fields , 1994 .

[11]  Spatiotemporal Analysis of Area Mortality , 1994 .

[12]  A van der Linde,et al.  Spatial-temporal analysis of mortality using splines. , 1995, Biometrics.

[13]  D. Clayton,et al.  Bayesian analysis of space-time variation in disease risk. , 1995, Statistics in medicine.

[14]  R. Littell.  SAS System for Mixed Models , 1996 .

[15]  Knut Sølna,et al.  Time Trend Estimation for a Geographic Region , 1996 .

[16]  Bradley P. Carlin,et al.  Hierarchical Spatio-Temporal Mapping of Disease Rates , 1997 .

[17]  B P Carlin,et al.  Spatio-temporal models with errors in covariates: mapping Ohio lung cancer mortality. , 1998, Statistics in medicine.