Home > Research > Publications & Outputs > Bayesian inference and data augmentation scheme...

Associated organisational unit

Electronic data

  • v63i07

    Final published version, 3.62 MB, PDF document

    Available under license: CC BY


Text available via DOI:

View graph of relations

Bayesian inference and data augmentation schemes for spatial, spatiotemporal and multivariate Log-Gaussian Cox processes in R

Research output: Contribution to Journal/MagazineSpecial issuepeer-review

Article number7
<mark>Journal publication date</mark>10/02/2015
<mark>Journal</mark>Journal of Statistical Software
Number of pages48
Pages (from-to)1-48
Publication StatusPublished
<mark>Original language</mark>English


Log-Gaussian Cox processes are an important class of models for spatial and spatiotemporal point-pattern data. Delivering robust Bayesian inference for this class of models presents a substantial challenge, since Markov chain Monte Carlo (MCMC) algorithms require careful tuning in order to work well. To address this issue, we describe recent advances in MCMC methods for these models and their implementation in the R package lgcp. Our suite of R functions provides an extensible framework for inferring covariate effects as well as the parameters of the latent field.

We also present methods for Bayesian inference in two further classes of model based on the log-Gaussian Cox process. The first of these concerns the case where we wish to fit a point process model to data consisting of event-counts aggregated to a set of spatial regions: we demonstrate how this can be achieved using data-augmentation. The second concerns Bayesian inference for a class of marked-point processes specified via a multivariate log-Gaussian Cox process model. For both of these extensions, we give details of their implementation in R.