Abstract
In geostatistics, both kriging and smoothing splines are commonly used to generate an interpolated map of a quantity of interest. The geoadditive model proposed by Kammann and Wand (J R Stat Soc: Ser C (Appl Stat) 52(1):1–18, 2003) represents a fusion of kriging and penalized spline additive models. Complex data issues, including non-linear covariate trends, multiple measurements at a location and clustered observations are easily handled using the geoadditive model. We propose a likelihood based estimation procedure that enables the estimation of the range (spatial decay) parameter associated with the penalized splines of the spatial component in the geoadditive model. We present how the spatial covariance structure (covariogram) can be derived from the geoadditive model. In a simulation study, we show that the underlying spatial process and prediction of the spatial map are estimated well using the proposed likelihood based estimation procedure. We present several applications of the proposed methods on real-life data examples.
Similar content being viewed by others
References
Akaike H (1973) Maximum likelihood identification of Gaussian autoregressive moving average models. Biometrika 60:255–265
Altman N (2000) Krige, smooth, both or neither. Aust N Z J Stat 42(4):441–461
Banerjee S, Gelfand AE, Finley AO, Sang H (2008) Gaussian predictive process models for large spatial data sets. J R Stat Soc Ser B 70(4):825–848
Brumback BA, Ruppert D, Wand MP (1999) Comment on Shively, Kohn and Wood. J Am Stat Assoc 94:794–797
Crainiceanu CM, Ruppert D (2004) Likelihood ratio tests in linear mixed models with one variance component. J R Stat Soc Ser B 66:165–185
Cressie NAC (1985) Fitting variogram models by weighted least squares. J Int Assoc Math Geol 5:563–586
Cressie NAC (1993) Statistics for spatial data. Wiley, New York
Cressie NAC, Johannesson G (2008) Fixed rank kriging for very large spatial data sets. J R Stat Soc Ser B 70(1):209–226
Diggle PJ, Ribeiro P (2002) Bayesian inference in Gaussian model-based geostatistics. Geogr Environ Model 6:129–146
Dubrule O (1983) Two methods with different objectives: splines and kriging. J Int Assoc Math Geol 15:245–257
Dubrule O (1984) Comparing splines and kriging. Comput Geosci 10:327–338
French JL, Kammann EE, Wand MP (2001) Comment on paper by Ke and Wang. J Am Stat Assoc 96:1285–1288
Gelfand AE, Diggle PJ, Fuentes M, Guttorp P (eds) (2010) Handbook of spatial statistics. Chapman & Hall/CRC, Boca Raton
Green PJ, Silverman BW (1994) Nonparametric regression and generalized linear models: a roughness approach. Chapman & Hall, London
Handcock MS, Meier K, Nychka D (1994) Comment on paper by (Laslett). J Am Stat Assoc 89:401–403
Hutchinson MF, Gessler PE (1994) Splines—more than just a smooth interpolator. Geoderma 62:45–67
Johnson ME, Moore LM, Ylvisaker D (1990) Minimax and maximin distance designs. J Stat Plan Inference 26:131–148
Kammann EE, Wand MP (2003) Geoadditive models. J R Stat Soc Ser C (Appl Stat) 52(1):1–18
Kaufman CG, Schervish MJ, Nychka DW (2008) Covariance tapering for likelihoodbased estimation in large spatial data sets. J Am Stat Assoc 103:1545–1555
Kleiber W, Nychka DW (2015) Equivalent kriging. Spatial. Statistics 12:31–49
Laslett GM (1994) Kriging and splines: an empirical comparison of their predictive performance in some applications. J Am Stat Assoc 89:392–409 (with discussion)
Lee D-J, Durbán M (2009) Smooth-CAR mixed models for spatial count data. Comput Stat Data Anal 53:2968–2979
Lee D-J, Durbán M (2011) P-spline ANOVA-type interaction models for spatio-temporal smoothing. Stat Model 11:46–69
Lindgren F, Rue H, Lindström J (2011) An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach. J R Stat Soc Ser B (Stat Methodol) 73(4):423–498
Mathéron G (1963) Principles of geostatistics. Econ Geol 58:1246–1266
Mathéron G (1981) Splines and griging: their formal equivalence. Down to earth statistics: solutions looking for geological problems. Syracuse Univ Geol Contrib 8:77–95
Nychka D (2000) Spatial process estimates as smoothers. In: Schimek M (ed) Smoothing and regression. Springer, Heidelberg
Nychka D, Bandyopadhyay S, Hammerling D, Lindgren F, Sain S (2015) A multiresolution Gaussian process model for the analysis of large spatial datasets. J Comput Gr Stat 24(2):579–599
Pebesma EJ (2004) Multivariable geostatistics in S: the gstat package. Comput Geosci 30:683–691
R Core Team (2014) R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna
Ribeiro Jr PJ, Diggle PJ (2001) geoR: a package for geostatistical analysis. R-NEWS 1(2):14–18. ISSN 1609-3631
Ruppert D, Wand MP, Carroll RJ (2003) Semiparametric regression. University Press, Cambridge
Salkauskas K (1982) Some relationships between surface splines and kriging. In: Schempp W, Zeller K (eds) Multivariate approximation theory II. Birkhauser, Basel, pp 313–325
Sang H, Huang JZ (2012) A full scale approximation of covariance functions for large spatial data sets. J R Stat Soc Ser B (Stat Methodol) 74(1):111–132
Schelin L, Sjöstedt-de Luna S (2010) Kriging prediction intervals based on semiparametric bootstrap. Math Geosci 42(8):985–1000
Searle SR (1982) Matrix algebra useful for statistics. Wiley, New York
Self SG, Liang KY (1987) Asymptotic properties of maximum likelihood estimators and likelihood ratio tests under nonstandard conditions. J Am Stat Assoc 82:605–610
Stein ML (1999) Interpolation of spatial data: some theory for kriging. Springer, New York
Stein ML (2014) Limitations on low rank approximations for covariance matrices of spatial data. Spat Stat 8:1–19
Stein ML, Handcock MS (1989) Some asymptotic properties of kriging when the covariance function is misspecified. Math Geol 21(2):171–190
Stewart GW (1973) Introduction to matrix computations. Academic Press, London
Wahba G, Wendelberger J (1980) Some new mathematical methods for variational objective analysis using splines and cross-validation. Mon Weather Rev 108:1122–1145
Watson GS (1984) Smoothing and interpolation by kriging and with splines. Math Geol 16:601–615
Wood SN (2006) Low-rank scale-invariant tensor product ssmooth for generalized additive mixed models. Biometrics 62:1025–1036
Zhang H, Zimmerman DL (2007) Hybrid estimation of semivariogram parameters. J Int Assoc Math Geol 39:247–260
Zimmerman DL, Zimmerman MB (1991) A comparison of spatial semivariogram estimators and corresponding ordinary kriging predictiors. Technometrics 33:77–91
Acknowledgements
Support from a doctoral Grant of Hasselt University is acknowledged (BOF11D04FAEC to YV). Support from the National Institutes of Health is acknowledged [award number R01CA172805 to CF]. Support from the University of Antwerp scientific chair in Evidence-Based Vaccinology, financed in 2009–2014 by a gift from Pfizer, is acknowledged [to NH]. Support from the IAP Research Network P7/06 of the Belgian State (Belgian Science Policy) is gratefully acknowledged. For the simulation study we used the infrastructure of the VSC - Flemish Supercomputer Center, funded by the Hercules Foundation and the Flemish Government - department EWI.
Author information
Authors and Affiliations
Corresponding author
Additional information
Handling Editor: Pierre Dutilleul.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
About this article
Cite this article
Vandendijck, Y., Faes, C. & Hens, N. Estimating the spatial covariance structure using the geoadditive model. Environ Ecol Stat 24, 341–361 (2017). https://doi.org/10.1007/s10651-017-0373-3
Received:
Revised:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10651-017-0373-3