Detecting Multivariate Interactions in Spatial Point Patterns with Gibbs Models and Variable Selection

point patterns
biodiversity
Gibbs models
Authors
Affiliation

Tuomas Rajala

UCL

David Murrell

UCL

Sofia Olhede

UCL

Published

April 23, 2018

Doi

Code: R-package PenGE

Abstract

We propose a method for detecting significant interactions in very large multivariate spatial point patterns. This methodology thus develops high dimensional data understanding in the point process setting. The method is based on modelling the patterns by using a flexible Gibbs point process model to characterize point-to-point interactions at different spatial scales directly. By using the Gibbs framework significant interactions can also be captured at small scales. Subsequently, the Gibbs point process is fitted by using a pseudolikelihood approximation, and we select significant interactions automatically by using the group lasso penalty with this likelihood approximation. Thus we estimate the multivariate interactions stably even in this setting. We demonstrate the feasibility of the method with a simulation study and show its power by applying it to a large and complex rainforest plant population data set of 83 species.

Figure 6: Different estimates of the interaction matrix of the adult plant BCI ’05 census. Columns 1 and 2: Non-parametric Monte Carlo method with two different MC tests and two different spatial range vectors. The species kept fixed in the test is on the x-axis, the randomized species on the y-axis. Columns 3 and 4: Group lasso with raw residual CV and inverse residual CV methods, and two range vectors. Species are arranged by point count, increasing from left/bottom to right/top.