# Predicting paleoclimate from compositional data using multivariate Gaussian process inverse prediction

## TL;DR

This paper introduces a novel computational framework using Gaussian process approximation for inverse prediction of unobserved covariates from multivariate compositional count data, with applications in paleoclimate modeling.

## Contribution

The authors develop a new efficient method for inverse prediction with Gaussian processes, enabling simultaneous inference on latent processes and covariates in high-dimensional, multi-modal spaces.

## Key findings

- Model achieves competitive predictive accuracy in paleoclimate state prediction.
- Framework allows formal statistical inference on community dynamics.
- Efficient exploration of complex latent spaces in inverse problems.

## Abstract

Multivariate compositional count data arise in many applications including ecology, microbiology, genetics, and paleoclimate. A frequent question in the analysis of multivariate compositional count data is what values of a covariate(s) give rise to the observed composition. Learning the relationship between covariates and the compositional count allows for inverse prediction of unobserved covariates given compositional count observations. Gaussian processes provide a flexible framework for modeling functional responses with respect to a covariate without assuming a functional form. Many scientific disciplines use Gaussian process approximations to improve prediction and make inference on latent processes and parameters. When prediction is desired on unobserved covariates given realizations of the response variable, this is called inverse prediction. Because inverse prediction is mathematically and computationally challenging, predicting unobserved covariates often requires fitting models that are different from the hypothesized generative model. We present a novel computational framework that allows for efficient inverse prediction using a Gaussian process approximation to generative models. Our framework enables scientific learning about how the latent processes co-vary with respect to covariates while simultaneously providing predictions of missing covariates. The proposed framework is capable of efficiently exploring the high dimensional, multi-modal latent spaces that arise in the inverse problem. To demonstrate flexibility, we apply our method in a generalized linear model framework to predict latent climate states given multivariate count data. Based on cross-validation, our model has predictive skill competitive with current methods while simultaneously providing formal, statistical inference on the underlying community dynamics of the biological system previously not available.