A Spatial Autoregressive Graphical Model with Applications in Intercropping

TL;DR

This paper introduces a novel Bayesian spatial autoregressive graphical model that captures asymmetric multivariate spatial effects, demonstrated through simulations and applied to intercropping data to reveal complex spatial relationships.

Contribution

The paper develops a new model combining multivariate spatial autoregressive and Gaussian graphical models to handle asymmetric spatial effects, filling a key methodological gap.

Findings

Model effectively captures asymmetric spatial effects in simulations.

Application to intercropping data reveals complex spatial crop relationships.

The R package facilitates practical implementation of the methodology.

Abstract

Within the statistical literature, a significant gap exists in methods capable of modeling asymmetric multivariate spatial effects that elucidate the relationships underlying complex spatial phenomena. For such a phenomenon, observations at any location are expected to arise from a combination of within- and between- location effects, where the latter exhibit asymmetry. This asymmetry is represented by heterogeneous spatial effects between locations pertaining to different categories, that is, a feature inherent to each location in the data, such that based on the feature label, asymmetric spatial relations are postulated between neighbouring locations with different labels. Our novel approach synergises the principles of multivariate spatial autoregressive models and the Gaussian graphical model. This synergy enables us to effectively address the gap by accommodating asymmetric spatial relations, overcoming the usual constraints in spatial analyses. Using a Bayesian-estimation framework, the model performance is assessed in a simulation study. We apply the model on intercropping data, where spatial effects between different crops are unlikely to be symmetric, in order to illustrate the usage of the proposed methodology. An R package containing the proposed methodology can be found on https://CRAN.R-project.org/package=SAGM.