Double hierarchical generalized linear models (DHGLM) are a family of models that are flexible enough as to model hierarchically the mean and scale parameters. In a Bayesian framework, fitting highly parameterized hierarchical models is challenging when this problem is addressed using typical Markov chain Monte Carlo (MCMC) methods due to the potential high correlation between different parameters and effects in the model. The integrated nested Laplace approximation (INLA) could be considered instead to avoid dealing with these problems. However, DHGLM do not fit within the latent Gaussian Markov random field (GMRF) models that INLA can fit. In this paper we show how to fit DHGLM with INLA by combining INLA and importance sampling (IS) algorithms. In particular, we will illustrate how to split DHGLM into submodels that can be fitted with INLA so that the remainder of the parameters are fit using adaptive multiple IS (AMIS) with the aid of the graphical representation of the hierarchical model. This is illustrated using a simulation study on three different types of models and two real data examples.
翻译:双级双级通用线性模型(DHGLM)是一组非常灵活的模型,足以按等级对中值和比例参数进行模型。在巴伊西亚框架内,在使用典型的Markov链Monte Carlo(MCMC)方法解决这个问题时,安装高度参数等级化的等级模型具有挑战性,因为模型中不同参数和效应之间可能存在高度关联。可以考虑综合的嵌巢式Laplace近似值(INLA),以避免处理这些问题。但是,DHGLM并不适合INLA适合的潜在Gossian Markov随机字段(GMRF)模型。在本文中,我们展示了如何将DHGLM与INLA结合,将INLA和重要取样算法(IS)结合起来,如何将DHGLM分为可与INLA相适应的子模型,以便其余参数能够使用适应性多重IS(AMIS)与等级模型的图形表示法。我们用三种不同模型和两个真实数据实例进行模拟研究来说明。