A new class of [formula omitted]-transformations for the spatial analysis of Compositional Data

Abstract

Georeferenced compositional data are prominent in many scientific fields and in spatial statistics. This work addresses the problem of proposing models and methods to analyze and predict, through kriging, this type of data. To this purpose, a novel class of α-transformations, named the Isometric α-transformation (α-IT), is proposed, which encompasses the traditional Isometric Log-Ratio (ILR) transformation. Similarly to other α-transformations existing in the literature, it is shown that the ILR is the limit case of the α-IT as α tends to 0 and that α=1 corresponds to a linear transformation of the data. Unlike the ILR, the proposed transformation accepts 0s in the compositions when α>0. Maximum likelihood estimation of the parameter α is established. Prediction using kriging on α-IT transformed data is validated on synthetic spatial compositional data, using prediction scores computed either in the geometry induced by the α-IT, or in the simplex. Application to land cover data shows that the relative superiority of the various approaches w.r.t. a prediction objective depends on whether the compositions contained any zero component. When all components are positive, the limit cases (ILR or linear transformations) are optimal for none of the considered metrics. An intermediate geometry, corresponding to the α-IT with maximum likelihood estimate, better describes the dataset in a geostatistical setting. When the amount of compositions with 0s is not negligible, some side-effects of the transformation gets amplified as α decreases, entailing poor kriging performances both within the α-IT geometry and for metrics in the simplex.