r - 使用带有 lme4 模型 : mediator model is not yet implemented - previous answer does not work for me 的中介包时出错

Question

大家好，我是 R 和这个论坛的新手。我有同样的问题 a previous poster ，但我的答案没有得到解决。当我尝试使用包 lme4 和调解运行多级调解模型时，我收到消息:“调解器模型尚未实现”。

我的数据:

data.frame':    25383 obs. of  115 variables:
$ PID                        : num  1 1 1 1 1 1 1 1 1 1 ...
$ T0_AGE                     : num  66.6 79.7 85.6 87 79.9 67.4 80 72 80.1 
68$ T0_ASEXE                 : Factor w/ 2 levels "Male","Female": 1 1...
$ T0_ALIVING_R               : num  0 0 0 0 1 0 0 0 1 0 ...
$ T0_educationcat            : Factor w/ 3 levels "high","middle",..
$ FI_morbidity_corrected     : num  0.0625 0 0.125 0.0625 0.0625 0.125 ……
$ FI_SRH                     : num  0.5 0.375 0.75 0.625 0.375…."

我正在尝试了解 FI_morbidity_corrected 如何调节 T0_educationcat 和 FI_SRH 之间的关系。我为 PID(项目 ID)指定了一个随机拦截，因为数据聚集在不同的研究/项目中。我正在控制混杂因素 T0_ASEXE + T0_AGE + T0_ALIVING_R

med.fit <- lmer(FI_morbidity_corrected~T0_educationcat + T0_ASEXE + T0_AGE + 
T0_ALIVING_R + (1|PID),data=topicsmds)
out.fit <- lmer(FI_SRH~FI_morbidity_corrected + T0_educationcat + T0_ASEXE + 
T0_AGE + T0_ALIVING_R + (1|PID),data=topicsmds)

这行得通，但随后:

med23.out <- mediate(med.fit, out.fit, treat = "T0_educationcat", mediator = 
"FI_morbidity_corrected", control.value = "high", treat.value = "middle", 
sims = 100)
summary(med23.out)

和

med24.out <- mediate(med.fit, out.fit, treat = "T0_educationcat", mediator = 
"FI_morbidity_corrected", control.value = "high", treat.value = "low", sims 
= 100)
summary(med24.out)

给我错误:

mediator model is not yet implemented

我已经按照建议使用 library(lme4) 加载了 lme4 包，但仍然出现此错误。

inherits(mediatorModel, "merMod") 返回 TRUE

getCall(mediatorModel)[[1]] 返回 lme4::lmer

谢谢!

Answer 1

您是否正在加载包 lmerTest？如果是这样，class(med.fit[[1]]) 返回 merModLmerTest 而不是 lmerMod - 但中介包需要后者。卸载 lmerTest 即可解决问题。