Journal
JOURNAL OF STATISTICAL SOFTWARE
Volume 78, Issue 8, Pages 1-29Publisher
JOURNAL STATISTICAL SOFTWARE
DOI: 10.18637/jss.v078.i08
Keywords
Bayesian analysis; geostatistics; low-rank approximations; Monte Carlo maximum likelihood; prevalence data; R
Funding
- MRC [MR/M015297/1] Funding Source: UKRI
- Medical Research Council [MR/M015297/1] Funding Source: researchfish
Ask authors/readers for more resources
In this paper we introduce a new R package, PrevMap, for the analysis of spatially referenced prevalence data, including both classical maximum likelihood and Bayesian approaches to parameter estimation and plug-in or Bayesian prediction. More specifically, the new package implements fitting of geostatistical models for binomial data, based on two distinct approaches. The first approach uses a generalized linear mixed model with logistic link function, binomial error distribution and a Gaussian spatial process as a stochastic component in the linear predictor. A simpler, but approximate, alternative approach consists of fitting a linear Gaussian model to empirical-logit-transformed data. The package also includes implementations of convolution-based low-rank approximations to the Gaussian spatial process to enable computationally efficient analysis of large spatial datasets. We illustrate the use of the package through the analysis of Loa loa prevalence data from Cameroon and Nigeria. We illustrate the use of the low rank approximation using a simulated geostatistical dataset.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available