具有两个随机效应的多级模型中的 ICC

机器算法验证 混合模式 多层次分析 lme4-nlme 类内相关
2022-04-09 10:02:06

我的理解是,类内相关性让您了解您的二级因素可以解释因变量的总体方差有多大的方差。它应该表明数据的方差结构是否证明了多层次分析是合理的。

现在,我指定了一个空模型,作为混合效应分析的起点,其中测试项目嵌套在受试者中。因此,我将每个主题 ID 视为包含 28 个项目的二级单元。看起来像这样:

Mt01 <- lmer(APMt ~ (1|ID) + (1|content), data=mlmData, REML=FALSE)

#Output:
Random effects:
 Groups   Name        Variance Std.Dev.
 ID       (Intercept) 0.102    0.319   
 content  (Intercept) 0.300    0.548   
 Residual             0.213    0.462   
Number of obs: 2940, groups:  ID, 105; content, 28

因变量是每个项目的 log(time)。我指定了两种随机效应,一种是由于受试者,另一种是由于项目内容(所有受试者都在相同的项目上工作)。

我将如何计算ICC?我是否为每个随机因素计算一个?在这种情况下甚至有必要吗?我将如何准确地解释该值?

以下是我的一些方法:

InterceptVarSubj <- VarCorr( Mt01 )$ID
    InterceptVarItem <- VarCorr( Mt01 )$content
ResidualVar <- attr( VarCorr( Mt01 ), "sc")^2

#A
icc <- InterceptVarItem / (InterceptVarItem + ResidualVar)
icc
[1] 0.584
#B
icc <- InterceptVarItem / (InterceptVarItem + InterceptVarSubj + ResidualVar)
icc
[1] 0.488
#C
icc <- InterceptVarSubj / (InterceptVarSubj + ResidualVar)
icc
[1] 0.323
#D
icc <- InterceptVarSubj / (InterceptVarItem + InterceptVarSubj + ResidualVar)
icc
[1] 0.165
1个回答

有几种方法可以计算和解释混合模型的 ICC。一个有用的线程在这里要计算,它是某些因素的方差量除以总方差。这类似于您的 B 和 D 计算。这可以解释为来自同一分组的两个随机抽取的单元的相关性,或者解释为这些分组解释的方差量,类似于R2.

这个计算(B):

σitem2σsubj2+σitem2+σres2
可以解释为任何项目的分数的相关性,而与使用该项目的主题无关。

同样,这个计算(D):

σsubj2σsubj2+σitem2+σres2
可以解释为任何科目的分数的相关性,无论他们正在研究什么项目。

组合计算,例如:

σsubj2+σitem2σsubj2+σitem2+σres2
可以解释为同一个人使用同一项目的分数之间的相关性(在您的数据中,这是 0.654)。

转向您列出的值(分别为 0.488 和 0.165),一种天真的解释是,无论谁在使用它,您使用的项目都适度相关,而所有项目的主题分数大多不相关。

最后,关于您的两个额外计算(A 和 C),我不知道对这些值的有用解释,但是我对这个主题有一个悬而未决的问题