Efficient iterative methods for hyperparameter estimation in large-scale linear inverse problems
Fuente:
arXiv
Saved in:
| Main Authors: | , , , |
|---|---|
| Format: | Preprint |
| Published: |
2023
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866911747978297344 |
|---|---|
| author | Hall-Hooper, Khalil A Saibaba, Arvind K Chung, Julianne Miller, Scot M |
| author_facet | Hall-Hooper, Khalil A Saibaba, Arvind K Chung, Julianne Miller, Scot M |
| contents | We study Bayesian methods for large-scale linear inverse problems, focusing on the challenging task of hyperparameter estimation. Typical hierarchical Bayesian formulations that follow a Markov Chain Monte Carlo approach are possible for small problems with very few hyperparameters but are not computationally feasible for problems with a very large number of unknown parameters. In this work, we describe an empirical Bayesian (EB) method to estimate hyperparameters that maximize the marginal posterior, i.e., the probability density of the hyperparameters conditioned on the data, and then we use the estimated values to compute the posterior of the inverse parameters. For problems where the computation of the square root and inverse of prior covariance matrices are not feasible, we describe an approach based on the generalized Golub-Kahan bidiagonalization to approximate the marginal posterior and seek hyperparameters that minimize the approximate marginal posterior. Numerical results from seismic and atmospheric tomography demonstrate the accuracy, robustness, and potential benefits of the proposed approach. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2311_15827 |
| institution | arXiv |
| publishDate | 2023 |
| record_format | arxiv |
| spellingShingle | Efficient iterative methods for hyperparameter estimation in large-scale linear inverse problems Hall-Hooper, Khalil A Saibaba, Arvind K Chung, Julianne Miller, Scot M Numerical Analysis 65F22, 65F55, 65M32, 62F10, 86A22 We study Bayesian methods for large-scale linear inverse problems, focusing on the challenging task of hyperparameter estimation. Typical hierarchical Bayesian formulations that follow a Markov Chain Monte Carlo approach are possible for small problems with very few hyperparameters but are not computationally feasible for problems with a very large number of unknown parameters. In this work, we describe an empirical Bayesian (EB) method to estimate hyperparameters that maximize the marginal posterior, i.e., the probability density of the hyperparameters conditioned on the data, and then we use the estimated values to compute the posterior of the inverse parameters. For problems where the computation of the square root and inverse of prior covariance matrices are not feasible, we describe an approach based on the generalized Golub-Kahan bidiagonalization to approximate the marginal posterior and seek hyperparameters that minimize the approximate marginal posterior. Numerical results from seismic and atmospheric tomography demonstrate the accuracy, robustness, and potential benefits of the proposed approach. |
| title | Efficient iterative methods for hyperparameter estimation in large-scale linear inverse problems |
| topic | Numerical Analysis 65F22, 65F55, 65M32, 62F10, 86A22 |
| url | https://arxiv.org/abs/2311.15827 |