目标:计算 ICC 并测试结果 ICC 是否足够大以解释研究中的聚类效应。
我的研究: 52 所学校的所有学生在三年内进行了纵向跟踪。学校被随机分配到治疗组或控制组。我的问题是我们需要考虑学校集群效应吗?
ICC是一种可以帮助我们判断聚类效应是否重要的度量。
如何计算ICC?
1)查看没有人使用治疗的基线年的学生分数并计算:
米小号祝你好运_ _ _ _/ (米小号祝你好运_ _ _ _+米小号乙r _)
2)另一种方法是通过使用随机效应模型框架和使用随机效应的方差(对于学校)来计算ICC。这是最常见的不同组大小的方法。在这种情况下,ICC 是
σ2b/ (σ2b+σ2e)
在哪里σ2b是随机效应的方差和σ2e是残差方差。
如何在 R 中拟合随机效应模型?
library(nlme)
schRandModel <- lme(StudScoresBaseline ~ 1, data = data, random = ~ 1 | SCHOOL)
SchRandModel.LogLiklihood <- schRandModel$logLik
上面的模型假设学校是随机截距的。现在我们拟合没有学校的模型作为随机效应(组是学生而不是学校):
simpleModel <- lme(StudScoresBaseline ~ 1, data = data, random = ~ 1 | StudentID)
simpleModel.LogLiklihood <- simpleModel$logLik
最后,使用似然比检验,我们检验似然之间的差异:
logLikDiff = SchRandModel.LogLiklihood-simpleModel.LogLiklihood
1-pchisq(2*logLikDiff,df=1)
我得到的 p 值小于 <0.0001,表明我们有足够的证据表明随机效应模型更适合。因此,学校集群效应在本研究中很重要。
感谢Macro指导我一步一步写出答案。