Markov-Chain Monte-Carlo methods and non-identifiabilities

被引:1
|
作者
Mueller, Christian [1 ]
Weysser, Fabian [2 ]
Mrziglod, Thomas [2 ]
Schuppert, Andreas [1 ]
机构
[1] Rhein Westfal TH Aachen, Joint Res Ctr Computat Biomed, Aachen, Germany
[2] Bayer AG, Appl Math, Leverkusen, Germany
来源
MONTE CARLO METHODS AND APPLICATIONS | 2018年 / 24卷 / 03期
关键词
Markov-Chain Monte-Carlo; non-identifiability;
D O I
10.1515/mcma-2018-0018
中图分类号
O21 [概率论与数理统计]; C8 [统计学];
学科分类号
020208 ; 070103 ; 0714 ;
摘要
We consider the problem of sampling from high-dimensional likelihood functions with large amounts of non-identifiabilities via Markov-ChainMonte-Carlo algorithms. Non-identifiabilities are problematic for commonly used proposal densities, leading to a low effective sample size. To address this problem, we introduce a regularization method using an artificial prior, which restricts non-identifiable parts of the likelihood function. This enables us to sample the posterior using common MCMC methods more efficiently. We demonstrate this with three MCMC methods on a likelihood based on a complex, high-dimensional blood coagulation model and a single series of measurements. By using the approximation of the artificial prior for the non-identifiable directions, we obtain a sample quality criterion. Unlike other sample quality criteria, it is valid even for short chain lengths. We use the criterion to compare the following three MCMC variants: The Random Walk Metropolis Hastings, the Adaptive Metropolis Hastings and the Metropolis adjusted Langevin algorithm.
引用
收藏
页码:203 / 214
页数:12
相关论文
共 50 条