An Analysis of High-Dimensional Land-Water Transitions in the Great Lakes via a Bayesian Discrepancy Model for Multinomial Data
Abstract
Global water-level systems have undergone large changes in recent years. Moderate Resolution Imaging Spectroradiometer (MODIS) land cover data can be used to monitor and understand these changes. However, such data are extremely high-dimensional making them difficult to process while accounting for all sources of uncertainty. Bayesian hierarchical models are commonly implemented using Markov Chain Monte Carlo (MCMC) methods which can be computationally expensive for large datasets. Exact Posterior Regression (EPR) provides an efficient alternative for inference in latent Gaussian hierarchical models under exponential families without requiring MCMC or approximate Bayesian strategies. EPR has not been developed for multinomial data, and we extend this class of models to the multinomial setting. Additionally, we explore signal-to-noise covariance properties of discrepancy models, which we expect to be present in our data. We extend EPR to multinomial likelihoods via a stick-breaking construction to analyze water-level changes. Through simulation studies, we compare the proposed multinomial EPR approach with alternative models implemented using Gibbs sampling and the Integrated Nested Laplace Approximation (INLA). We evaluate each approach based on regression parameter estimation, category probability estimation, and computational efficiency. Overall, EPR and INLA achieve substantial reductions in computation time while EPR shows improved metrics in estimation and prediction.
Related
- Spectral Costs of Multinomial and Residual Resampling in High-Dimensional Particle Ensembles: A Free-Probability Analysis
- Bayesian calibration with model discrepancy of a two-group interfacial area transport model for upward bubbly flow in a square duct
- Logarithmic transport at topological phase transitions of one-dimensional chiral quantum systems
- Maximum Diversity and Support Transitions on Disks and Three-Dimensional Balls
- On exact solutions to three-dimensional Yukawa lattice model and BCC to FCC phase transitions
BibTeX
@misc{lee2026analysis,
title = {{An Analysis of High-Dimensional Land-Water Transitions in the Great Lakes via a Bayesian Discrepancy Model for Multinomial Data}},
author = {Seulki Lee and Jonathan R. Bradley},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.38916},
eprint = {2610.38916},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.38916},
}
RIS
TY - UNPB
TI - An Analysis of High-Dimensional Land-Water Transitions in the Great Lakes via a Bayesian Discrepancy Model for Multinomial Data
AU - Seulki Lee
AU - Jonathan R. Bradley
PY - 2026
DA - 2026/10/28
PB - arXiv
DO - 10.48550/arxiv.2610.38916
UR - https://arxiv.org/abs/2610.38916
AB - Global water-level systems have undergone large changes in recent years. Moderate Resolution Imaging Spectroradiometer (MODIS) land cover data can be used to monitor and understand these changes. However, such data are extremely high-dimensional making them difficult to process while accounting for all sources of uncertainty. Bayesian hierarchical models are commonly implemented using Markov Chain Monte Carlo (MCMC) methods which can be computationally expensive for large datasets. Exact Posterior Regression (EPR) provides an efficient alternative for inference in latent Gaussian hierarchical models under exponential families without requiring MCMC or approximate Bayesian strategies. EPR has not been developed for multinomial data, and we extend this class of models to the multinomial setting. Additionally, we explore signal-to-noise covariance properties of discrepancy models, which we expect to be present in our data. We extend EPR to multinomial likelihoods via a stick-breaking construction to analyze water-level changes. Through simulation studies, we compare the proposed multinomial EPR approach with alternative models implemented using Gibbs sampling and the Integrated Nested Laplace Approximation (INLA). We evaluate each approach based on regression parameter estimation, category probability estimation, and computational efficiency. Overall, EPR and INLA achieve substantial reductions in computation time while EPR shows improved metrics in estimation and prediction.
ER -