R语言如何用潜类别混合效应模型（LCMM）分析抑郁症状

相关视频

模型背景

当对重复测量的标志变量进行建模时，我们通常不会把它看成是一个有误差测量的潜过程。然而，这正是混合模型理论所做的基本假设。潜过程混合模型利用这个框架将线性混合模型理论扩展到任何类型的结果（有序、二元、连续、类别与任何分布）。

潜类别混合模型

潜类别混合模型在Proust-Lima等人中介绍（2006  https://doi.org/10.1111/j.1541-0420.2006.00573.x 和2013 https://doi.org/10.1111/bmsp.12000  ）。

使用线性混合模型根据时间对定义为潜过程感兴趣的变量进行建模：

其中：

X（t） 和Z（t） 是协变量的向量（Z（t） ；

β 是固定效应（即总体均值效应）；

ui 是随机效应（即个体效应）；它们根据具有协方差矩阵B 的零均值多元正态分布进行分布；

（wi（t）） 是高斯过程，可以添加到模型中以来放宽对象内部相关结构。

同时在观察方程中定义了感兴趣的潜过程标志变量Yij （针对对象i和场合j）的观察之间的关系：

其中

tij 是主题i 和场合j 的测量时间；

ϵij 是一个独立的零均值高斯误差；

H 是链接函数，可将潜过程转换为比例和度量。

使用了不同的参数族。 当标志变量为连续时，H-1 是递增单调函数的参数族，其中：

线性变换：这简化为线性混合模型（2个参数）

Beta累积分布族重新调整（4个参数）

当标志变量是离散类别（二元或有序的）时：  H是阈值函数，即Y的每个级别对应于要估计Λ（tij）+ ϵij区间的边界。

可识别性

与任何潜变量模型一样，必须定义潜变量的度量。在lcmm中，误差的方差为1，平均截距（在β中）为0。

示例

在本文中，lcmm通过研究年龄65岁左右男性的抑郁症状（由CES-D量表测量）的线性轨迹来说明潜过程混合模型 。包括截距和age65的相关随机效应。

考虑的模型：

，

固定效应部分是

估计不同连续链接函数的模型H

我们使用65岁左右的年龄变量进行中心化，并以十年为单位。

潜过程混合模型可以用不同的链接函数进行拟合，如下所示。这是用参数链接来完成的。

线性链接函数

定义线性链接函数时，模型将简化为标准线性混合模型。默认情况下具有线性链接函数：

lcmm(CESD ~ age65*male, random=~ age65 #链接=线性

它与hlme安装的模型完全相同。与hlme对象的唯一区别是截距和残差标准误差的参数化。

hlme(CESD ~ age65*male, random=~ age65 #链接=线性

对数似然相同，但估计参数β不在同一范围内

loglik

\[1\] -7056.652

非线性链接函数1：Beta累积分布函数Beta分布的重标累积分布函数(CDF)提供了标志变量与其基本潜伏过程之间的凹、凸或sigmoïd变换。

lcmm( random=~ age65, link='beta')

非线性链接函数2：二次I样条二次I样条族近似于连续增加的链接函数。它涉及在标志变量范围内分布的节点。默认情况下，使用位于标志变量范围内的5个等距结：

lcmm(random=~ age65, subject='ID', link='splines')

可以指定结的数量及其位置。首先输入节点的数目，接着 ，再指定位置equi，quant或manual用于分别等距节点，在标志变量分布分位数或内部结在参数intnodes手动输入。例如，7-equi-splines意味着具有7个等距节点，6-quant-splinesI样条， 意味着具有6个节点的I样条，其位于标志变量分布的分位数处。

例如，在分位数处有5个结：

lcmm(link='5-quant-splines')选择最佳模型

要选择最合适的链接函数，可以比较这些不同的模型。通常，这可以通过使用AIC 或 UACV等顺着根据拟合优度对模型进行比较来实现 。

AIC（每个模型的输出中都有UACV）：

在这种情况下，根据AIC标准，由I-splines和5个分位数结点链接函数的模型提供了最佳拟合度。可以在图中比较不同的估计链接函数：

plot(mli, which="linkfunction",xlab="潜过程")

legend(x="topleft", legend=c("线性", "beta","样条曲线 (5个等距结点)","样条曲线(5个分位数结点)"))