Estimation of log-Gaussian gamma processes with iterated posterior linearization and Hamiltonian Monte Carlo
Fuente:
arXiv
Saved in:
| Main Authors: | , |
|---|---|
| Format: | Preprint |
| Published: |
2026
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866912887332667392 |
|---|---|
| author | Härkönen, Teemu Särkkä, Simo |
| author_facet | Härkönen, Teemu Särkkä, Simo |
| contents | Stochastic processes are a flexible and widely used family of models for statistical modeling. While stochastic processes offer attractive properties such as inclusion of uncertainty properties, their inference is typically intractable, with the notable exception of Gaussian processes. Inference of models with non-Gaussian errors typically involves estimation of a high-dimensional latent variable. We propose two methods that use iterated posterior linearization followed by Hamiltonian Monte Carlo to sample the posterior distributions of such latent models with a particular focus on log-Gaussian gamma processes. The proposed methods are validated with two synthetic datasets generated from the log-Gaussian gamma process and a multiscale biocomposite stiffness model. In addition, we apply the methodology to an experimental Raman spectrum of argentopyrite. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2602_07454 |
| institution | arXiv |
| publishDate | 2026 |
| record_format | arxiv |
| spellingShingle | Estimation of log-Gaussian gamma processes with iterated posterior linearization and Hamiltonian Monte Carlo Härkönen, Teemu Särkkä, Simo Methodology Computation 62M20 (Primary), 62-08, 60-08 (Secondary) Stochastic processes are a flexible and widely used family of models for statistical modeling. While stochastic processes offer attractive properties such as inclusion of uncertainty properties, their inference is typically intractable, with the notable exception of Gaussian processes. Inference of models with non-Gaussian errors typically involves estimation of a high-dimensional latent variable. We propose two methods that use iterated posterior linearization followed by Hamiltonian Monte Carlo to sample the posterior distributions of such latent models with a particular focus on log-Gaussian gamma processes. The proposed methods are validated with two synthetic datasets generated from the log-Gaussian gamma process and a multiscale biocomposite stiffness model. In addition, we apply the methodology to an experimental Raman spectrum of argentopyrite. |
| title | Estimation of log-Gaussian gamma processes with iterated posterior linearization and Hamiltonian Monte Carlo |
| topic | Methodology Computation 62M20 (Primary), 62-08, 60-08 (Secondary) |
| url | https://arxiv.org/abs/2602.07454 |