Linda Altieri

A new approach to spatial entropy measures

Entropy is widely employed in many applied sciences to measure the heterogeneity of observations. Recently, many attempts have been made to build entropy measures for spatial data, in order to capture the influence of space over the variable outcomes. The main limit of these developments is that all indices are computed conditional on a single distance and do not cover the whole spatial configuration of the phenomenon under study. Moreover, most of them do not satisfy the desirable additivity property between local and global spatial measures. This work reviews some recent developments, based on univariate distributions, and compares them to a new approach which considers the properties of entropy measures linked to bivariate distributions. This perspective introduces substantial innovations. Firstly, Shannon’s entropy may be decomposed into two terms: spatial mutual information, accounting for the role of space in determining the variable outcome, and spatial global residual entropy, summarizing the remaining heterogeneity carried by the variable itself. Secondly, these terms both satisfy the additivity property, being sums of partial entropies measuring what happens at different distance classes. The proposed indices are used for measuring the spatial entropy of a marked point pattern on rainforest tree species. The new entropy measures are shown to be more informative and to answer a wider set of questions than the current proposals of the literature.

Bayesian P-splines and advanced computing in R for a changepoint analysis on spatio-temporal point processes

ABSTRACT This work presents advanced computational aspects of a new method for changepoint detection on spatio-temporal point process data. We summarize the methodology, based on building a Bayesian hierarchical model for the data and declaring prior conjectures on the number and positions of the changepoints, and show how to take decisions regarding the acceptance of potential changepoints. The focus of this work is about choosing an approach that detects the correct changepoint and delivers smooth reliable estimates in a feasible computational time; we propose Bayesian P-splines as a suitable tool for managing spatial variation, both under a computational and a model fitting performance perspective. The main computational challenges are outlined and a solution involving parallel computing in R is proposed and tested on a simulation study. An application is also presented on a data set of seismic events in Italy over the last 20 years.