## Abstract

Conventional geostatistical methodology solves the problem of predicting the realized value of a linear functional of a Gaussian spatial stochastic process S(x) based on observations Y_{i} = S(x_{i}) + Z_{i} at sampling locations x_{i}, where the Z_{i} are mutually independent, zero-mean Gaussian random variables. We describe two spatial applications for which Gaussian distributional assumptions are clearly inappropriate. The first concerns the assessment of residual contamination from nuclear weapons testing on a South Pacific island, in which the sampling method generates spatially indexed Poisson counts conditional on an unobserved spatially varying intensity of radioactivity; we conclude that a conventional geostatistical analysis oversmooths the data and underestimates the spatial extremes of the intensity. The second application provides a description of spatial variation in the risk of campylobacter infections relative to other enteric infections in part of north Lancashire and south Cumbria. For this application, we treat the data as binomial counts at unit postcode locations, conditionally on an unobserved relative risk surface which we estimate. The theoretical framework for our extension of geostatistical methods is that, conditionally on the unobserved process S(x), observations at sample locations x_{i} form a generalized linear model with the corresponding values of S(x_{i}) appearing as an offset term in the linear predictor. We use a Bayesian inferential framework, implemented via the Markov chain Monte Carlo method, to solve the prediction problem for non-linear functionals of S(x), making a proper allowance for the uncertainty in the estimation of any model parameters.

Original language | English (US) |
---|---|

Pages (from-to) | 299-325 |

Number of pages | 27 |

Journal | Journal of the Royal Statistical Society. Series C: Applied Statistics |

Volume | 47 |

Issue number | 3 |

State | Published - 1998 |

Externally published | Yes |

## Keywords

- Generalized linear mixed model
- Geostatistics
- Kriging
- Markov chain Monte Carlo method
- Spatial prediction

## ASJC Scopus subject areas

- Mathematics(all)
- Statistics and Probability