Bayesian Monte Carlo testing with one-dimensional measures of evidence

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS

卷:351

Article

Silva, Ivair R.^1,2,3  Marques, Reinaldo^4,5 

[1] Harvard Med Sch, Dept Populat Med, Boston, MA 02115 USA

[2] Harvard Pilgrim Hlth Care Inst, Boston, MA USA

[3] Univ Fed Ouro Preto, Dept Stat, Ouro Preto, MG, Brazil

[4] Univ Oslo, Dept Stat, Oslo, Norway

[5] Univ Fed Alfenas, Actuarial Div, Alfenas, MG, Brazil

关键词: Bayes factor;    Expected loss;    Exact test;    Measure of evidence;    Monte Carlo;   

DOI  :  10.1016/j.cam.2018.11.016

来源: Elsevier

【 摘 要 】

Bayesian hypothesis testing procedures are constructed by means of test statistics which are functions of the posterior distribution. Usually, the whole sample vector is selected to form the sufficient empirical part of the posterior distribution. But, in certain problems, one may prefer to use well-established one-dimensional sufficient statistics in place of the sample vector. This paper introduces a Bayesian Monte Carlo procedure specially designed for such cases. It is shown that the performance of this new approach is arbitrarily close to the exact Bayesian test. In addition, for arbitrary desired precisions, we develop a theoretical rule of thumb for choosing the minimum number m(0) of Monte Carlo simulations. Surprisingly, m(0) does not depend on the shape of loss/cost functions when those are used to compound the test statistic. The method is illustrated for testing mean vectors in high dimension and for detecting spatial clusters of diseases in aggregated maps. (C) 2018 Elsevier B.V. All rights reserved.

【 授权许可】

Free