Geostatistical inference under preferential sampling

Daniel Dinsdale & Matias Salibian 2019-11-11

Geostatistical inference under preferential sampling

This repository contains R code illustrating the computations in Dinsdale, D. and Salibian-Barrera, M. (2018) Methods for preferential sampling in geostatistics, Journal of the Royal Statistical Society Series C, 68(1), 181-198.

We use TMB and INLA. For the latter we adapted the scripts publicly available here, with contributions from an anonymous referee.

Please note that the TMB code below was last tested with the R INLA package "INLA_19.09.03".

In this repository you will find:

• A script to generate synthetic data following the Poisson preferential sampling model here
• A script to obtain parameter estimates using INLA in R here
• A script to obtain parameter estimates using TMB in R here and here