Imputing relative abundances via Gaussian random fields
- 1Uppsala University, Science for Life Laboratory
Abstract
Longitudinal studies with human participants frequently face the challenge of missing samples, which lead to incomplete trajectories of participant biomarkers and a reduction in data size. Existing approaches focus on single-variable imputation or model features one at a time, while neural networks typically require a large number of parameters and substantial training data. To address these challenges, we present an explainable model for multiple imputation of compositional samples based on a Gaussian random field (GRF). The model is specifically designed for longitudinal data: it incorporates participant-specific information and accounts for both temporal dependence across repeated observations and within-sample correlations, while using a moderate number of parameters. By coupling the GRF with sparse activation functions, the framework outputs sum-constrained compositions and explicitly models zero inflation. Under the Aitchison and Jensen--Shannon distances, our framework achieves imputation accuracy comparable to that of an observation-carrying baseline, while overcoming its static nature by modeling temporal trajectories and multivariate correlations. The resulting multiple imputations allow us to reconstruct complete datasets for downstream analysis and uncertainty quantification. We demonstrate the practical utility of our framework on gut microbiome datasets, restoring incomplete time series and achieving a balanced accuracy of 0.72 in predicting food allergy in infants (sensitivity 0.8) and of 0.64 in forecasting postpartum depression.
Related
- Imputing missing counts for waterbird species: A comparison of methods for zero-inflated biodiversity monitoring data
- Spatially heterogenous distributions of dipteran abundances identify priority areas for haemosporidian parasite surveillance
- An information-theoretic perspective on feed-forward loop abundances in transcriptional networks
- Hudson River Estuary sediments contain high fractional abundances of preserved globothalamean foraminiferal organic linings at Piermont, New York
- Joint species distributions reveal crop type and field-specific assembly rules and idiosyncrasies in carabid beetle occurrences and abundances
BibTeX
@misc{shternshis2026imputing,
title = {{Imputing relative abundances via Gaussian random fields}},
author = {A. Shternshis and P. Singh and L. W. Hugerth},
howpublished = {bioRxiv},
year = {2026},
month = oct,
doi = {10.64898/2026.10.12.758512},
url = {https://www.biorxiv.org/content/10.64898/2026.10.12.758512v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - Imputing relative abundances via Gaussian random fields
AU - A. Shternshis
AU - P. Singh
AU - L. W. Hugerth
PY - 2026
DA - 2026/10/12
PB - bioRxiv
DO - 10.64898/2026.10.12.758512
UR - https://www.biorxiv.org/content/10.64898/2026.10.12.758512v1
AB - Longitudinal studies with human participants frequently face the challenge of missing samples, which lead to incomplete trajectories of participant biomarkers and a reduction in data size. Existing approaches focus on single-variable imputation or model features one at a time, while neural networks typically require a large number of parameters and substantial training data. To address these challenges, we present an explainable model for multiple imputation of compositional samples based on a Gaussian random field (GRF). The model is specifically designed for longitudinal data: it incorporates participant-specific information and accounts for both temporal dependence across repeated observations and within-sample correlations, while using a moderate number of parameters. By coupling the GRF with sparse activation functions, the framework outputs sum-constrained compositions and explicitly models zero inflation. Under the Aitchison and Jensen--Shannon distances, our framework achieves imputation accuracy comparable to that of an observation-carrying baseline, while overcoming its static nature by modeling temporal trajectories and multivariate correlations. The resulting multiple imputations allow us to reconstruct complete datasets for downstream analysis and uncertainty quantification. We demonstrate the practical utility of our framework on gut microbiome datasets, restoring incomplete time series and achieving a balanced accuracy of 0.72 in predicting food allergy in infants (sensitivity 0.8) and of 0.64 in forecasting postpartum depression.
ER -