A Bayesian approach to approximate joint diagonalization of square matrices

Mingjun Zhong, Mark Girolami

Research output: Chapter in Book/Report/Conference proceedingConference contribution


We present a Bayesian scheme for the approximate diagonalisation of several square matrices which are not necessarily symmetric. A Gibbs sampler is derived to simulate samples of the common eigenvectors and the eigenvalues for these matrices. Several synthetic examples are used to illustrate the performance of the proposed Gibbs sampler and we then provide comparisons to several other joint diagonalization algorithms, which shows that the Gibbs sampler achieves the state-of-the-art performance on the examples considered. As a byproduct, the output of the Gibbs sampler could be used to estimate the log marginal likelihood, however we employ the approximation based on the Bayesian information criterion (BIC) which in the synthetic examples considered correctly located the number of common eigenvectors. We then succesfully applied the sampler to the source separation problem as well as the common principal component analysis and the common spatial pattern analysis problems.
Original languageEnglish
Title of host publicationProceedings of the 29 th International Conference on Machine Learning, Edinburgh, Scotland, UK, 2012
Number of pages8
Publication statusPublished - 2012


Dive into the research topics of 'A Bayesian approach to approximate joint diagonalization of square matrices'. Together they form a unique fingerprint.

Cite this