Geostatistical analysis of binomial data: Generalised linear or transformed Gaussianmodelling?

Michelle C. Stanton, Peter J. Diggle

Research output: Contribution to journalArticlepeer-review

Abstract

When fitting a binomial geostatistical model to data obtained by spatially discrete sampling, techniques such as Markov chain Monte Carlo, which are computationally intensive and require careful tuning to each application, need to be employed. As a result, this approach is often infeasible when computational resources are scarce or statistical expertise is unavailable. A practical solution to this problem is to transform the binomial samples such that the conditional distribution of the transformed data can be assumed to be approximately Gaussian. Likelihood-based inference can then be performed analytically, whereas Bayesian inference requires only routine Monte Carlo simulation. In this paper, the predictive performance of the binomial model-fitting method is compared with that of the approximate transformed Gaussian method. A simulation experiment is undertaken in which data are simulated from a binomial geostatistical model, and both methods are then used to make predictive inference. Predictions are assessed using the empirical root mean square prediction error and empirical coverage probabilities of nominal prediction intervals obtained for the approximate method. As expected, the comparative performance of the approximate method deteriorates as the denominator decreases and, to a lesser extent, as the overall proportion of successes in the simulated data deviates from 0.5. Our results provide guidance on when the approximate method can safely be used. We demonstrate the two methods on village-level prevalence data pertaining to the tropical eye disease Loa loa.

Original languageEnglish (US)
Pages (from-to)158-171
Number of pages14
JournalEnvironmetrics
Volume24
Issue number3
DOIs
StatePublished - May 2013
Externally publishedYes

Keywords

  • Model-based geostatistics
  • Prevalence surveys
  • Spatial variation
  • Trans-Gaussian kriging

ASJC Scopus subject areas

  • Ecological Modeling
  • Statistics and Probability

Fingerprint

Dive into the research topics of 'Geostatistical analysis of binomial data: Generalised linear or transformed Gaussianmodelling?'. Together they form a unique fingerprint.

Cite this