Bayesian Methods for Variance Component Models
研究了正态理论方差分量模型的精确贝叶斯分析，使用重要性抽样和拉普拉斯T近似，比较了贝叶斯估计量与Lindley-Stein收缩估计量，发现联合众数可能过度收缩，而分层贝叶斯分析表现同样好且允许更详细的有限样本推断。
Abstract An exact Bayesian analysis can be performed for normal theory variance components models, using importance sampling. But remarkably accurate approximations are available using the Laplacian T-approximation introduced by Leonard, Hsu, and Ritter. Instead of maximizing the joint posterior density, conditional upon the parameter of interest, a device described by O'Hagan is used first. The Bayesian estimators are compared to the Lindley—Stein shrinkage estimators and the Lindley—Smith joint modal estimators. It is confirmed that joint modes can overcollapse toward prior hypotheses, when compared with more sensible Bayesian procedures. This is referred to as a “collapsing phenomenon.” A numerical example from the one-way random-effects model is considered, and the risks of the different estimators are simulated under a variety of loss functions. It is concluded that although the Lindley—Stein estimator performs well, a full hierarchical Bayesian analysis performs at least equally well, while permitting more detailed finite-sample inference regarding any parameter of interest.