如何检查混合效应模型中因子和协变量之间的相互作用?

机器算法验证 重复测量 相互作用 lme4-nlme 对比 预测器
2022-03-25 16:36:16

在主题内设计中有 2 个因子AB(5×3) 和一个协变量。X这是我指定整体模型的方式:

lme.out = lme(y~ A*B*X, random=~1|Subject, data=mydata)

我的解释是我正在查看一个图表y~x,其中斜率由于协变量而变化,并且线根据和的不同水平向上或向下移动AB截距的变化)。

我想知道的是:如果我要修复因子A(采取任何级别),然后查看线条(y~x),效果是B什么?水平线是B向上还是向下移动(截距)还是改变线的斜率(X)。

我应该进行某种对比分析吗?但我不确定因素和协变量之间的对比如何起作用。

我能想到的一种方法是获取与 A 的不同级别相对应的数据子集并创建模型,例如:lme(y~ B+X, random=~1|Subject, data=mydata[which(mydata$A = A1,]). 这样我就可以比较这些模型的截距和斜率。

谁能告诉我我所做的是否有意义?任何形式的建议将不胜感激!

3个回答

为了可视化交互术语,您可以查看sjPlot 包请参见此处的示例)。

您的函数调用将是

sjp.int(fit, type ="eff")

但是,我不确定这是否满足您的需求?

我个人认为,如果您想在控制 X 后检查 Y 与模型中因素之间的真实关系,您应该查看绘制的调整后的平均值,而不是从您最喜欢的模型计算的原始平均值。出于此目的,有 R 包,例如 lsmeans,它们非常方便且用户友好!

您可以调查模型的 VIF。VIF 代表方差膨胀因子,是一种测量共线性的方法。

https://onlinecourses.science.psu.edu/stat501/node/347

R 的 car 包中有一个 vif 函数。