Research output: Contribution to Journal/Magazine › Journal article › peer-review
Research output: Contribution to Journal/Magazine › Journal article › peer-review
}
TY - JOUR
T1 - INLA or MCMC?
T2 - a tutorial and comparative evaluation for spatial prediction in log-Gaussian Cox processes
AU - Taylor, Benjamin
AU - Diggle, Peter
PY - 2014
Y1 - 2014
N2 - We investigate two options for performing Bayesian inference on spatial log-Gaussian Cox processes assuming a spatially continuous latent field: Markov chain Monte Carlo (MCMC) and the integrated nested Laplace approximation (INLA). We first describe the device of approximating a spatially continuous Gaussian field by a Gaussian Markov random field on a discrete lattice, and present a simulation study showing that, with careful choice of parameter values, small neighbourhood sizes can give excellentapproximations. We then introduce the spatial log-Gaussian Cox process and describe MCMC and INLA methods for spatial prediction within this model class. We report the results of a simulation study in which we compare the Metropolis-adjusted Langevin Algorithm (MALA) and the technique of approximating the continuous latent field by a discrete one, followed by approximate Bayesian inference via INLA over a selection of 18 simulated scenarios. The results question the notion that the latter technique is bothsignificantly faster and more robust than MCMC in this setting; 100,000 iterations of the MALA algorithm running in 20 min on a desktop PC delivered greater predictive accuracy than the default INLA strategy,which ran in 4 min and gave comparative performance to the full Laplace approximation which ran in 39 min.
AB - We investigate two options for performing Bayesian inference on spatial log-Gaussian Cox processes assuming a spatially continuous latent field: Markov chain Monte Carlo (MCMC) and the integrated nested Laplace approximation (INLA). We first describe the device of approximating a spatially continuous Gaussian field by a Gaussian Markov random field on a discrete lattice, and present a simulation study showing that, with careful choice of parameter values, small neighbourhood sizes can give excellentapproximations. We then introduce the spatial log-Gaussian Cox process and describe MCMC and INLA methods for spatial prediction within this model class. We report the results of a simulation study in which we compare the Metropolis-adjusted Langevin Algorithm (MALA) and the technique of approximating the continuous latent field by a discrete one, followed by approximate Bayesian inference via INLA over a selection of 18 simulated scenarios. The results question the notion that the latter technique is bothsignificantly faster and more robust than MCMC in this setting; 100,000 iterations of the MALA algorithm running in 20 min on a desktop PC delivered greater predictive accuracy than the default INLA strategy,which ran in 4 min and gave comparative performance to the full Laplace approximation which ran in 39 min.
KW - log-Gaussian Cox process
KW - Markov chain Monte Carlo
KW - integrated nested Laplace approximation
KW - spatial modelling
U2 - 10.1080/00949655.2013.788653
DO - 10.1080/00949655.2013.788653
M3 - Journal article
VL - 84
SP - 2266
EP - 2284
JO - Journal of Statistical Computation and Simulation
JF - Journal of Statistical Computation and Simulation
SN - 1563-5163
IS - 10
ER -