A statistical test for Nested Sampling algorithms

被引:171
作者
Buchner, Johannes [1 ]
机构
[1] Max Planck Inst Extraterr Phys, Giessenbachstr, D-85748 Garching, Germany
关键词
Nested sampling; MCMC; Bayesian inference; Evidence; Test; Marginal likelihood; MONTE-CARLO; INFERENCE; EFFICIENT; MODEL;
D O I
10.1007/s11222-014-9512-y
中图分类号
TP301 [理论、方法];
学科分类号
081202 ;
摘要
Nested sampling is an iterative integration procedure that shrinks the prior volume towards higher likelihoods by removing a "live" point at a time. A replacement point is drawn uniformly from the prior above an ever-increasing likelihood threshold. Thus, the problem of drawing from a space above a certain likelihood value arises naturally in nested sampling, making algorithms that solve this problem a key ingredient to the nested sampling framework. If the drawn points are distributed uniformly, the removal of a point shrinks the volume in a well-understoodway, and the integration of nested sampling is unbiased. In this work, I develop a statistical test to check whether this is the case. This "Shrinkage Test" is useful to verify nested sampling algorithms in a controlled environment. I apply the shrinkage test to a test-problem, and show that some existing algorithms fail to pass it due to over-optimisation. I then demonstrate that a simple algorithm can be constructed which is robust against this type of problem. This RADFRIENDS algorithm is, however, inefficient in comparison to MULTINEST.
引用
收藏
页码:383 / 392
页数:10
相关论文
共 15 条
[1]   Nested sampling for parameter inference in systems biology: application to an exemplar circadian model [J].
Aitken, Stuart ;
Akman, Ozgur E. .
BMC SYSTEMS BIOLOGY, 2013, 7 :72
[2]  
Beaujean F., 2013, ARXIV E PRINTS
[3]  
Betancourt Michael, 2011, AIP Conference Proceedings, V1305, P165, DOI 10.1063/1.3573613
[4]  
Cameron E., 2013, ARXIV E PRINTS
[5]  
Evans M., 2007, Bayesian Statistics, V8, P491
[6]   Multimodal nested sampling: an efficient and robust alternative to Markov Chain Monte Carlo methods for astronomical data analyses [J].
Feroz, F. ;
Hobson, M. P. .
MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY, 2008, 384 (02) :449-463
[7]   MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics [J].
Feroz, F. ;
Hobson, M. P. ;
Bridges, M. .
MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY, 2009, 398 (04) :1601-1614
[8]  
Feroz F., 2013, ARXIV E PRINTS
[9]   A nested sampling algorithm for cosmological model selection [J].
Mukherjee, P ;
Parkinson, D ;
Liddle, AR .
ASTROPHYSICAL JOURNAL, 2006, 638 (02) :L51-L54
[10]  
Robert C., 2007, BAYESIAN STAT, V8, P491