Spatial designs and properties of spatial correlation: Effects on covariance estimation

Abstract

In a spatial regression context, scientists are often interested in a physical interpretation of components of the parametric covariance function. For example, spatial covariance parameter estimates in ecological settings have been interpreted to describe spatial heterogeneity or "patchiness" in a landscape that cannot be explained by measured covariates. In this article, we investigate the influence of the strength of spatial dependence on maximum likelihood (ML) and restricted maximum likelihood (REML) estimates of covariance parameters in an exponential-with-nugget model, and we also examine these influences under different sampling designs-specifically, lattice designs and more realistic random and cluster designs-at differing intensities of sampling (n = 144 and 361). We find that neither ML nor REML estimates perform well when the range parameter and/or the nugget-to-sill ratio is large-ML tends to underestimate the autocorrelation function and REML produces highly variable estimates of the autocorrelation function. The best estimates of both the covariance parameters and the autocorrelation function come under the cluster sampling design and large sample sizes. As a motivating example, we consider a spatial model for stream sulfate concentration.

An Explicit Link between Gaussian Fields and Gaussian Markov Random Fields: The Stochastic Partial Differential Equation Approach

“…A challenge in spatial statistics is that complex models require large sample sizes to estimate model parameters reliably (Irvine et al, 2007), but many modelling procedures are too computationally complex for large sample sizes. The link between GFs and GMRFs will allow more researchers to investigate important statistical issues like model selection, precision of parameter estimates (especially spatial covariance parameters), sampling designs and more.…”

Parameter identifiability, constraint, and equifinality in data assimilation with ecosystem models

“…The sample sites were mainly located in Pennsylvania, West Virginia, Maryland, and Virginia. For more details about this example and the issues described below, see Irvine et al (2007).…”

“…When the goal of an analysis is to provide a map or some other inference across a sampling area, then additional considerations should be made when designing the study. It has been shown in a number of contexts that a cluster sampling design is appropriate for spatially correlated data (e.g., Zimmerman 2006, Irvine et al 2007, Ritter and Leecaster 2007. A cluster design includes some observations observed at close distances as well as sampling coverage over the entire sampling area.…”

Prediction of spatial soil property information from ancillary sensor data using ordinary linear regression: Model derivations, residual assumptions and model validation tests

“…In practice, the β parameter vector and Σ variance structure must be jointly estimated from the sample data, typically using maximum likelihood (ML) or restricted maximum likelihood (REML) estimation techniques (Littell et al, 1996). In such situations it is generally necessary to collect a fairly large amount of sample data in order to reasonably estimate the parameters associated with the covariance structure when even the simplest isotropic covariance functions are employed (Irvine et al, 2007).…”