Triple-goal estimates for disease mapping.

Maps of regional morbidity and mortality rates play an important role in assessing environmental equity. They provide effective tools for identifying areas with potentially elevated risk, determining spatial trend, and formulating and validating aetiological hypotheses about disease. Bayes and empirical Bayes methods produce stable small-area rate estimates that retain geographic and demographic resolution. The beauty of the Bayesian approach lies in its ability to structure complicated models, inferential goals and analyses. Three inferential goals are relevant to disease mapping and risk assessment: (i) computing accurate estimates of disease rates in small geographic areas; (ii) estimating the distribution of disease rates over the region; (iii) ranking the disease rates so that environmental investigation can be prioritized. No single set of estimates can simultaneously optimize these three goals, and Shen and Louis propose a set of estimates that perform well on all three goals. These are optimal for estimating the distribution of rates and for ranking, and maintain a high accuracy in estimating area-specific rates. However, the Shen/Louis method is sensitive to choice of priors. To address this issue we introduce a robustified version of the method based on a smoothed non-parametric estimate of the prior. We evaluate the performance of this method through a simulation study, and illustrate it using a data set of county-specific lung cancer rates in Ohio.

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

[2]  Wei Shen,et al.  Empirical Bayes Estimation via the Smoothing by Roughening Approach , 1999 .

[3]  Noel A Cressie,et al.  Statistics for Spatial Data. , 1992 .

[4]  D. Clayton,et al.  Empirical Bayes estimates of age-standardized relative risks for use in disease mapping. , 1987, Biometrics.

[5]  T. Louis,et al.  Empirical Bayes methods for stabilizing incidence rates before mapping. , 1994, Epidemiology.

[6]  Thomas A. Louis,et al.  Empirical Bayes Ranking Methods , 1989 .

[7]  Bradley P. Carlin,et al.  Hierarchical models for mapping Ohio lung cancer rates , 1997 .

[8]  T. Louis,et al.  Empirical bayes estimators for spatially correlated incidence rates , 1994 .

[9]  L. Bernardinelli,et al.  Bayesian methods for mapping disease risk , 1996 .

[10]  T. Louis.  Estimating a population of parameter values using Bayes and empirical Bayes methods , 1984 .

[11]  S. Zeger,et al.  A Smooth Nonparametric Estimate of a Mixing Distribution Using Mixtures of Gaussians , 1996 .

[12]  L Knorr-Held,et al.  Modelling risk from a disease in time and space. , 1998, Statistics in medicine.

[13]  D. Berry,et al.  Empirical Bayes Estimation of a Binomial Parameter Via Mixtures of Dirichlet Processes , 1979 .

[14]  Two semi-parametric empirical Bayes estimators , 1999 .

[15]  N. Laird.  Nonparametric Maximum Likelihood Estimation of a Mixing Distribution , 1978 .

[16]  L Bernardinelli,et al.  Empirical Bayes versus fully Bayesian analysis of geographical variation in disease risk. , 1992, Statistics in medicine.

[17]  J. Besag,et al.  Bayesian image restoration, with two applications in spatial statistics , 1991 .

[18]  Nan M. Laird,et al.  Empirical bayes estimates using the nonparametric maximum likelihood estimate for the priort , 1982 .

[19]  T. Louis,et al.  Triple‐goal estimates in two‐stage hierarchical models , 1998 .

[20]  T. Louis,et al.  Smoothing the non-parametric estimate of a prior distribution by roughening: A computational study , 1991 .

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

[22]  C. Morris.  Parametric Empirical Bayes Inference: Theory and Applications , 1983 .

[23]  M. Ghosh.  Constrained Bayes Estimation with Applications , 1992 .

[24]  M. Escobar.  Estimating Normal Means with a Dirichlet Process Prior , 1994 .