Researchers are often interested in understanding the relationship between a set of covariates and a set of response variables. To achieve this goal, the use of regression analysis, either linear or generalized linear models, is largely applied. However, such models only allow users to model one response variable at a time. Moreover, it is not possible to directly calculate from the regression model a correlation measure between the response variables. In this article, we employed the Multivariate Generalized Linear Mixed Models framework, which allows the specification of a set of response variables and calculates the correlation between them through a random effect structure that follows a multivariate normal distribution. We used the maximum likelihood estimation framework to estimate all model parameters using Laplace approximation to integrate out the random effects. The derivatives are provided by automatic differentiation. The outer maximization was made using a general-purpose algorithm such as \texttt{PORT} and \texttt{BFGS}. We delimited this problem by studying only count response variables with the following distributions: Poisson, negative binomial (NB) and COM-Poisson. The models were implemented on software \texttt{R} with package \texttt{TMB}. Besides the full specification, models with simpler structures in the covariance matrix were considered (fixed and common variance, fixed dispersion, $\rho$ set to 0). These models were applied to a dataset from the National Health and Nutrition Examination Survey, where three underdispersed response variables were measured at 1281 subjects. The COM-Poisson model full specified overcome the other two competitors considering three goodness-of-fit indexes. Therefore, the proposed model can deal with multivariate count responses and measures the correlation between them taking into account the effects of the covariates.
翻译:研究人员通常对理解一组共变数和一组响应变量之间的关系感兴趣。 为了实现这一目标, 基本上应用回归分析( 线性或通用线性模型) 。 但是, 这些模型只允许用户一次模拟一个响应变量。 此外, 无法直接从回归模型中计算反应变量之间的相关度量。 在本篇文章中, 我们使用了多变通用线性混合模型框架, 该框架允许对一组响应变量进行规格说明, 并通过一个随机效果结构来计算它们之间的相互关系。 我们使用最大可能性估算框架来估算所有模型参数, 使用 Laplace 近似来将随机效果整合起来。 外最大化是使用一个通用算法, 如\ textt{Port} 和\ textt{BGS} 等。 我们通过只用下列分布来计算特定响应变量: Poisson, 负的双双曲线( NB) 和 COM-Poisson。 模型是在软件中应用的, 与Sqentral- developal comliversal deal deal deal deal developations 3 developmental developmental developations commotional deal deal demod the commotional demotions.