期刊
JOURNAL OF EDUCATIONAL AND BEHAVIORAL STATISTICS
卷 40, 期 2, 页码 136-157出版社
SAGE PUBLICATIONS INC
DOI: 10.3102/1076998615570945
关键词
Bayes modal estimation; penalized likelihood estimation; variance estimation; Heywood case; mixed-effects model; multilevel model
资金
- Institute of Education Sciences [R305D100017]
- National Science Foundation [SES-1323977, 9911374]
- Army Grant [W911NF-14-1-0020]
- Direct For Social, Behav & Economic Scie
- Divn Of Social and Economic Sciences [1323977] Funding Source: National Science Foundation
- Division Of Research On Learning
- Direct For Education and Human Resources [9911374] Funding Source: National Science Foundation
When fitting hierarchical regression models, maximum likelihood (ML) estimation has computational (and, for some users, philosophical) advantages compared to full Bayesian inference, but when the number of groups is small, estimates of the covariance matrix (sigma) of group-level varying coefficients are often degenerate. One can do better, even from a purely point estimation perspective, by using a prior distribution or penalty function. In this article, we use Bayes modal estimation to obtain positive definite covariance matrix estimates. We recommend a class of Wishart (not inverse-Wishart) priors for sigma with a default choice of hyperparameters, that is, the degrees of freedom are set equal to the number of varying coefficients plus 2, and the scale matrix is the identity matrix multiplied by a value that is large relative to the scale of the problem. This prior is equivalent to independent gamma priors for the eigenvalues of sigma with shape parameter 1.5 and rate parameter close to 0. It is also equivalent to independent gamma priors for the variances with the same hyperparameters multiplied by a function of the correlation coefficients. With this default prior, the posterior mode for sigma is always strictly positive definite. Furthermore, the resulting uncertainty for the fixed coefficients is less underestimated than under classical ML or restricted maximum likelihood estimation. We also suggest an extension of our method that can be used when stronger prior information is available for some of the variances or correlations.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据