On optimality properties of the power prior

被引:76
作者
Ibrahim, JG [1 ]
Chen, MH
Sinha, D
机构
[1] Univ N Carolina, Dept Biostat, Chapel Hill, NC 27599 USA
[2] Univ Connecticut, Dept Stat, Storrs, CT 06269 USA
[3] Med Univ S Carolina, Dept Biometry & Epidemiol, Charleston, SC 29425 USA
基金
美国国家卫生研究院; 美国国家科学基金会;
关键词
generalized linear model; Gibbs sampling; historical data; Kullback-Leibler divergence; power prior; prior elicitation; proportional-hazards model; random-effects model; sensitivity analyses;
D O I
10.1198/016214503388619229
中图分类号
O21 [概率论与数理统计]; C8 [统计学];
学科分类号
020208 ; 070103 ; 0714 ;
摘要
The power prior is a useful general class of priors that can be used for arbitrary classes of regression models, including generalized linear models, generalized linear mixed models semiparametric survival models with censored data, frailty models, multivariate models, and nonlinear models. The power prior specification for the regression coefficients focuses on observable quanities in that the elicitation is based on historical data, D-0, and a scalar quantity, a(0), quantifying the heterogeneity between the current data, D, and the historical data D-0. The power prior distribution is then constructed by raising the likelihood function of the historical data to the power a(0), where 0 less than or equal to a(0) < 1. The scalar a(0) is a precision parameter that can be viewed as a measure of compatibility between the historical and current data. In this article we give a formal justification of the power prior and show that it is an optimal class of informative priors in the sense that it minimizes a convex sum of Kullback-Leibler (KL) divergences between two specific posterior densities, in which one density is based on no incorporation of historical data and the other density is based on pooling the historical and current data. This result provides a strong motivation for using the power prior as an informative prior in Bayesian inference. In addition, we derive a formal relationship between this convex sum of KL divergences and the information-processing rules proposed by others. Specifically, we show that the power prior is a 100% efficient information-processing rule in the sense defined earlier. Several examples involving simulations as well as real datasets are examined to demonstrate the proposed methodology.
引用
收藏
页码:204 / 213
页数:10
相关论文
共 10 条
[1]   Power prior distributions for generalized linear models [J].
Chen, MH ;
Ibrahim, JG ;
Shao, QM .
JOURNAL OF STATISTICAL PLANNING AND INFERENCE, 2000, 84 (1-2) :121-137
[2]  
Ibrahim JG, 2000, STAT SCI, V15, P46
[3]   Using historical controls to adjust for covariates in trend tests for binary data [J].
Ibrahim, JG ;
Ryan, LM ;
Chen, MH .
JOURNAL OF THE AMERICAN STATISTICAL ASSOCIATION, 1998, 93 (444) :1282-1293
[4]  
IBRAHIM JG, 2001, BAYESIAN SURVIVAL AN
[5]   BAYESIAN AND NON-BAYESIAN METHODS FOR COMBINING MODELS AND FORECASTS WITH APPLICATIONS TO FORECASTING INTERNATIONAL GROWTH-RATES [J].
MIN, CK ;
ZELLNER, A .
JOURNAL OF ECONOMETRICS, 1993, 56 (1-2) :89-118
[6]   Information indices: unification and applications [J].
Soofi, ES ;
Retzer, JJ .
JOURNAL OF ECONOMETRICS, 2002, 107 (1-2) :17-40
[7]   BAYESIAN APPROACHES TO RANDOMIZED TRIALS [J].
SPIEGELHALTER, DJ ;
FREEDMAN, LS ;
PARMAR, MKB .
JOURNAL OF THE ROYAL STATISTICAL SOCIETY SERIES A-STATISTICS IN SOCIETY, 1994, 157 :357-387
[8]   Information processing and Bayesian analysis [J].
Zellner, A .
JOURNAL OF ECONOMETRICS, 2002, 107 (1-2) :41-50
[9]   OPTIMAL INFORMATION-PROCESSING AND BAYES THEOREM [J].
ZELLNER, A .
AMERICAN STATISTICIAN, 1988, 42 (04) :278-280
[10]  
Zellner A., 1997, BAYESIAN ANAL ECONOM