The Laplace approximation (LA) has been proposed as a method for approximating the marginal likelihood of statistical models with latent variables. However, the approximate maximum likelihood estimators (MLEs) based on the LA are often biased for binary or spatial data, and the corresponding Hessian matrix underestimates the standard errors of these approximate MLEs. A higher-order approximation has been proposed; however, it cannot be applied to complicated models such as correlated random effects models and does not provide consistent variance estimators. In this paper, we propose an enhanced LA (ELA) that provides the true MLE and its consistent variance estimator. We study its relationship to the variational Bayes method. We also introduce a new restricted maximum likelihood estimator (REMLE) for estimating dispersion parameters. The results of numerical studies show that the ELA provides a satisfactory MLE and REMLE, as well as their variance estimators for fixed parameters. The MLE and REMLE can be viewed as posterior mode and marginal posterior mode under flat priors, respectively. Some comparisons are also made with Bayesian procedures under different priors.
翻译:提议采用拉普尔近似值(LA)作为方法,以接近潜在变量统计模型的边际可能性;然而,基于LA的大致最大可能性估计值(MLE)往往偏向二进制数据或空间数据,相应的赫森矩阵低估了这些近似 MLE的标准差。提出了更高层次近似值;但无法应用于相关随机效应模型等复杂模型,也没有提供一致的差异估计器。在本文中,我们提议加强LA(ELA),提供真实的MLE及其一致的差异估计器。我们研究了其与变异海湾方法的关系。我们还采用了新的限制最大可能性估计分散参数。数字研究结果显示,ELE提供了令人满意的MLE和REMLE,以及固定参数的差异估计器。MLE和REMLE可被视为平坦前的后传模式和边际外延模式。一些比较也分别与先前不同巴伊斯程序进行。