Bayesian Multiresolution Modeling Of Georeferenced Data.

Current implementations of multiresolution methods are limited in terms of possible types of responses and approaches to inference. We provide a multiresolution approach for spatial analysis of non-Gaussian responses using latent Gaussian models and Bayesian inference via integrated nested Laplace approximation (INLA). The approach builds on `LatticeKrig', but uses a reparameterization of the model parameters that is intuitive and interpretable so that modeling and prior selection can be guided by expert knowledge about the different spatial scales at which dependence acts. The priors can be used to make inference robust and integration over model parameters allows for more accurate posterior estimates of uncertainty. The extended LatticeKrig (ELK) model is compared to a standard implementation of LatticeKrig (LK), and a standard Matern model, and we find modest improvement in spatial oversmoothing and prediction for the ELK model for counts of secondary education completion for women in Kenya collected in the 2014 Kenya demographic health survey. Through a simulation study with Gaussian responses and a realistic mix of short and long scale dependencies, we demonstrate that the differences between the three approaches for prediction increases with distance to nearest observation.

[1]  H. Künsch.  Gaussian Markov random fields , 1979 .

[2]  J. Hyman.  Accurate Monotonicity Preserving Cubic Interpolation , 1983 .

[3]  Henry C. Tuckwell.  Elementary applications of probability theory , 1988 .

[4]  M. Stein,et al.  A Bayesian analysis of kriging , 1993 .

[5]  L Knorr-Held,et al.  Bayesian Detection of Clusters and Discontinuities in Disease Maps , 2000, Biometrics.

[6]  Turid Follestad,et al.  GMRFLib: a C-library for fast and exact simulation of Gaussian Markov random elds , 2002 .

[7]  Stephen E. Fienberg,et al.  Testing Statistical Hypotheses , 2005 .

[8]  Glen D Meeden,et al.  Fuzzy and randomized confidence intervals and P-values , 2005 .

[9]  S. Martino.  Approximate Bayesian Inference for Latent Gaussian Models , 2007 .

[10]  A. Gelfand,et al.  Gaussian predictive process models for large spatial data sets , 2008, Journal of the Royal Statistical Society. Series B, Statistical methodology.

[11]  N. Cressie,et al.  Fixed rank kriging for very large spatial data sets , 2008 .

[12]  H. Rue,et al.  Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations , 2009 .

[13]  Toshihiro Hirano,et al.  Covariance tapering for prediction of large spatial data sets in transformed random fields , 2013 .

[14]  Jianhua Z. Huang,et al.  A full scale approximation of covariance functions for large spatial data sets , 2012 .

[15]  Daniel W. Apley,et al.  Local Gaussian Process Approximation for Large Computer Experiments , 2013, 1303.0383.

[16]  Alan E Gelfand,et al.  A multivariate spatial mixture model for areal data: examining regional differences in standardized test scores , 2014, Journal of the Royal Statistical Society. Series C, Applied statistics.

[17]  Bani K. Mallick,et al.  Adaptive Bayesian Nonstationary Modeling for Large Spatial Datasets Using Covariance Approximations , 2014 .

[18]  D. Nychka,et al.  A Multiresolution Gaussian Process Model for the Analysis of Large Spatial Datasets , 2015 .

[19]  Matthias Katzfuss,et al.  A Multi-Resolution Approximation for Massive Spatial Datasets , 2015, 1507.04789.

[20]  Finn Lindgren,et al.  Bayesian Spatial Modelling with R-INLA , 2015 .

[21]  Sudipto Banerjee,et al.  Hierarchical Nearest-Neighbor Gaussian Process Models for Large Geostatistical Datasets , 2014, Journal of the American Statistical Association.

[22]  Haavard Rue,et al.  Bayesian Computing with INLA: A Review , 2016, 1604.00860.

[23]  Dorit Hammerling,et al.  Parallel inference for massive distributed spatial data using low-rank models , 2017, Stat. Comput..

[24]  Matthias Katzfuss,et al.  A General Framework for Vecchia Approximations of Gaussian Processes , 2017, Statistical Science.

[25]  M. Fuentes,et al.  Circulant Embedding of Approximate Covariances for Inference From Gaussian Data on Large Lattices , 2017 .

[26]  Sudipto Banerjee,et al.  Web Appendix: Meta-Kriging: Scalable Bayesian Modeling and Inference for Massive Spatial Datasets , 2018 .

[27]  Joseph Guinness,et al.  Spectral density estimation for random fields via periodic embeddings. , 2019, Biometrika.

[28]  Jon Wakefield,et al.  Design- and Model-Based Approaches to Small-Area Estimation in A Low- and Middle-Income Country Context: Comparisons and Recommendations , 2019, Journal of Survey Statistics and Methodology.

[29]  Dorit Hammerling,et al.  A Case Study Competition Among Methods for Analyzing Large Spatial Data , 2017, Journal of Agricultural, Biological and Environmental Statistics.

[30]  Jon Wakefield,et al.  Changes in the spatial distribution of the under-five mortality rate: Small-area analysis of 122 DHS surveys in 262 subregions of 35 countries in Africa , 2019, PloS one.