Abstract
When modeling spatially distributed normal responses Y i in terms of vectors x i of explanatory variables, one may fit a linear model assuming independence, and then use the empirical variogram of the residuals to determine an appropriate parametric form for the autocorrelation function. Suppose, however, that the responses are not normally distributed—for example, Poisson or Bernoulli. One may model spatial dependence using a hierarchical generalized linear model in which, conditional on a latent Gaussian field Z = {Z i }, the Y i have independent distributions from the exponential family, with an appropriate link function connecting their conditional means with the linear predictors x t i β + Z i . The question then is how to determine an appropriate model for the autocorrelation function of Z. The empirical variogram of the Y i is no longer appropriate, since (unless the link function is the identity) it is on the wrong scale. We propose here an alternative, the latent scale covariogram, whose graph reflects the autocorrelation structure of the underlying normal field. We illustrate its use on several real datasets, together with a simulated dataset, and obtain results quite different from those obtained using the variogram. Supplementary materials for this article are available online.
Original language | English |
---|---|
Pages (from-to) | 127-141 |
Number of pages | 15 |
Journal | Journal of Computational and Graphical Statistics |
Volume | 28 |
Issue number | 1 |
DOIs | |
State | Published - 2 Jan 2019 |
Bibliographical note
Publisher Copyright:© 2019, © 2019 American Statistical Association, Institute of Mathematical Statistics, and Interface Foundation of North America.
Keywords
- Hierarchical generalized linear model
- Latent Gaussian field
- Spatial correlation
- Variogram