Bayesian spatial modeling of disease risk in relation to multivariate environmental risk fields

The relationship between exposure to environmental chemicals during pregnancy and early childhood development is an important issue that has a spatial risk component. In this context, we have examined mental retardation and developmental delay (MRDD) outcome measures for children in a Medicaid population in South Carolina and sampled measures of soil chemistry (e.g. As, Hg, etc.) on a network of sites that are misaligned to the outcome residential addresses during pregnancy. The true chemical concentration at the residential addresses is not observed directly and must be interpolated from soil samples. In this study, we have developed a Bayesian joint model that interpolates soil chemical fields and estimates the associated MRDD risk simultaneously. Having multiple spatial fields to interpolate, we have considered a low-rank Kriging method for the interpolation that requires less computation than the Bayesian Kriging. We performed a sensitivity analysis for a bivariate smoothing, changing the number of knots and the smoothing parameter. These analyses show that a low-rank Kriging method can be used as an alternative to a full-rank Kriging, reducing the computational burden. However, the number of knots for the low-rank Kriging model needs to be selected with caution as a bivariate surface estimation can be sensitive to the choice of the number of knots.

[1]  M. R. Osborne,et al.  Estimation of covariance parameters in kriging via restricted maximum likelihood , 1991 .

[2]  D. Ruppert.  Selecting the Number of Knots for Penalized Splines , 2002 .

[3]  David B. Dunson,et al.  Bayesian Data Analysis , 2010 .

[4]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[5]  Yuedong Wang,et al.  Semiparametric Nonlinear Mixed-Effects Models and Their Applications , 2001 .

[6]  M. Wand,et al.  Generalized additive models for cancer mapping with incomplete covariates. , 2004, Biostatistics.

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

[8]  Sw. Banerjee,et al.  Hierarchical Modeling and Analysis for Spatial Data , 2003 .

[9]  S C Darby,et al.  Some aspects of measurement error in explanatory variables for continuous and binary regression models. , 1998, Statistics in medicine.

[10]  A. Gelman.  Prior distributions for variance parameters in hierarchical models (comment on article by Browne and Draper) , 2004 .

[11]  Daniel Krewski,et al.  On measurement error adjustment methods in Poisson regression , 1999 .

[12]  Ying C. MacNab,et al.  Spline smoothing in Bayesian disease mapping , 2007 .

[13]  I. Heid,et al.  On the potential of measurement error to induce differential bias on odds ratio estimates: an example from radon epidemiology , 2002, Statistics in medicine.

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

[15]  Andrew B Lawson,et al.  A spatial analysis of mental retardation of unknown cause and maternal residence during pregnancy. , 2008, Geospatial health.

[16]  S. Wood.  Generalized Additive Models: An Introduction with R , 2006 .

[17]  M. Wand,et al.  Semiparametric Regression: Parametric Regression , 2003 .

[18]  V. Zadnik,et al.  Effects of Residual Smoothing on the Posterior of the Fixed Effects in Disease‐Mapping Models , 2006, Biometrics.

[19]  J. Schafer,et al.  Computational Strategies for Multivariate Linear Mixed-Effects Models With Missing Values , 2002 .

[20]  Brian J Smith,et al.  Iowa radon leukaemia study: a hierarchical population risk model for spatially correlated exposure measured with error , 2007, Statistics in medicine.

[21]  N. Schupf,et al.  Epidemiology and Etiology of Mental Retardation , 2007, Handbook of Intellectual and Developmental Disabilities.

[22]  Refik Soyer,et al.  Bayesian Methods for Nonlinear Classification and Regression , 2004, Technometrics.

[23]  S. Chib,et al.  Understanding the Metropolis-Hastings Algorithm , 1995 .

[24]  Geoffrey M. Laslett,et al.  Kriging and Splines: An Empirical Comparison of their Predictive Performance in Some Applications , 1994 .

[25]  Yuan Liu,et al.  Evaluation of Bayesian models for focused clustering in health data , 2007 .

[26]  Jennifer A Hoeting,et al.  Model selection for geostatistical models. , 2006, Ecological applications : a publication of the Ecological Society of America.

[27]  Noel A Cressie,et al.  Statistics for Spatial Data, Revised Edition. , 1994 .

[28]  Montserrat Fuentes,et al.  Spatial Association between Speciated Fine Particles and Mortality , 2006, Biometrics.

[29]  Andrew B Lawson,et al.  Metal concentrations in rural topsoil in South Carolina: potential for human health impact. , 2008, The Science of the total environment.