Geostatistical inference under preferential sampling

Research output: Contribution to journalArticlepeer-review


Geostatistics involves the fitting of spatially continuous models to spatially discrete data. Preferential sampling arises when the process that determines the data locations and the process being modelled are stochastically dependent. Conventional geostatistical methods assume, if only implicitly, that sampling is non-preferential. However, these methods are often used in situations where sampling is likely to be preferential. For example, in mineral exploration, samples may be concentrated in areas that are thought likely to yield high grade ore. We give a general expression for the likelihood function of preferentially sampled geostatistical data and describe how this can be evaluated approximately by using Monte Carlo methods. We present a model for preferential sampling and demonstrate through simulated examples that ignoring preferential sampling can lead to misleading inferences. We describe an application of the model to a set of biomonitoring data from Galicia, northern Spain, in which making allowance for preferential sampling materially changes the results of the analysis.
Original languageEnglish
Pages (from-to)191-232
Number of pages42
JournalRoyal Statistical Society. Journal. Series C: Applied Statistics
Publication statusPublished - Mar 2010


Dive into the research topics of 'Geostatistical inference under preferential sampling'. Together they form a unique fingerprint.

Cite this