Fast Variational Block-Sparse Bayesian Learning

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Möderl, Jakob, Leitinger, Erik, Fleury, Bernard H., Pernkopf, Franz, Witrisal, Klaus
Format: Preprint
Published: 2023
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866908627094208512
author Möderl, Jakob
Leitinger, Erik
Fleury, Bernard H.
Pernkopf, Franz
Witrisal, Klaus
author_facet Möderl, Jakob
Leitinger, Erik
Fleury, Bernard H.
Pernkopf, Franz
Witrisal, Klaus
contents We propose a variational Bayesian (VB) implementation of block-sparse Bayesian learning (BSBL) to compute proxy probability density functions (PDFs) that approximate the posterior PDFs of the weights and associated hyperparameters in a block-sparse linear model, resulting in an iterative algorithm coined variational BSBL (VA-BSBL). The priors of the hyperparameters are selected to belong to the family of generalized inverse Gaussian distributions. This family contains as special cases commonly used hyperpriors. Inspired by previous work on classical sparse Bayesian learning (SBL), we investigate the update stage in which the proxy PDFs of a single block of weights and of its associated hyperparameter are successively updated, while keeping the proxy PDFs of the other parameters fixed. This stage defines a nonlinear first-order recurrence relation for the mean of the proxy PDF of the hyperparameter. By iterating this relation "ad infinitum" we obtain a criterion that determines whether the so-generated sequence of hyperparameter means converges or diverges. Incorporating this criterion into the VA-BSBL algorithm yields a fast implementation, coined fast-BSBL (F-BSBL), which achieves a two-order-of-magnitude runtime improvement. We further identify the range of the parameters of the generalized inverse Gaussian distribution which result in an inherent pruning procedure that switches off "weak" components in the model, which is necessary to obtain sparse results. Lastly, we show that expectation-maximization (EM)-based and VB-based implementations of BSBL are identical methods. Thus, we extend a well-known result from classical SBL to BSBL. Consequently, F-BSBL and BSBL using coordinate ascent to maximize the marginal likelihood coincide. These results provide a unified framework for interpreting existing BSBL methods.
format Preprint
id arxiv_https___arxiv_org_abs_2306_00442
institution arXiv
publishDate 2023
record_format arxiv
spellingShingle Fast Variational Block-Sparse Bayesian Learning
Möderl, Jakob
Leitinger, Erik
Fleury, Bernard H.
Pernkopf, Franz
Witrisal, Klaus
Signal Processing
We propose a variational Bayesian (VB) implementation of block-sparse Bayesian learning (BSBL) to compute proxy probability density functions (PDFs) that approximate the posterior PDFs of the weights and associated hyperparameters in a block-sparse linear model, resulting in an iterative algorithm coined variational BSBL (VA-BSBL). The priors of the hyperparameters are selected to belong to the family of generalized inverse Gaussian distributions. This family contains as special cases commonly used hyperpriors. Inspired by previous work on classical sparse Bayesian learning (SBL), we investigate the update stage in which the proxy PDFs of a single block of weights and of its associated hyperparameter are successively updated, while keeping the proxy PDFs of the other parameters fixed. This stage defines a nonlinear first-order recurrence relation for the mean of the proxy PDF of the hyperparameter. By iterating this relation "ad infinitum" we obtain a criterion that determines whether the so-generated sequence of hyperparameter means converges or diverges. Incorporating this criterion into the VA-BSBL algorithm yields a fast implementation, coined fast-BSBL (F-BSBL), which achieves a two-order-of-magnitude runtime improvement. We further identify the range of the parameters of the generalized inverse Gaussian distribution which result in an inherent pruning procedure that switches off "weak" components in the model, which is necessary to obtain sparse results. Lastly, we show that expectation-maximization (EM)-based and VB-based implementations of BSBL are identical methods. Thus, we extend a well-known result from classical SBL to BSBL. Consequently, F-BSBL and BSBL using coordinate ascent to maximize the marginal likelihood coincide. These results provide a unified framework for interpreting existing BSBL methods.
title Fast Variational Block-Sparse Bayesian Learning
topic Signal Processing
url https://arxiv.org/abs/2306.00442