R 将中介包与lme4模型一起使用时出错:尚未实现中介模型-前面的答案对我不适用
嗨,我是R和这个论坛的新手。 我的问题和你的一样,但我的问题没有得到解决。当我尝试使用包lme4和中介运行多级中介模型时,我得到了这样一条消息:“中介模型尚未实现” 我的数据:R 将中介包与lme4模型一起使用时出错:尚未实现中介模型-前面的答案对我不适用,r,lme4,R,Lme4,嗨,我是R和这个论坛的新手。 我的问题和你的一样,但我的问题没有得到解决。当我尝试使用包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 7
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\u发病率纠正如何调节T0\u教育cat
和FI\u SRH
之间的关系。我正在为PID
(项目ID)指定一个随机截距,因为数据聚集在不同的研究/项目中。我正在控制混杂因素T0\u ASEXE+T0\u AGE+T0\u 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)
med.fit您是否正在加载软件包lmerTest
?如果是这样,class(med.fit[[1]])
返回merModLmerTest
,而不是lmerMod
——但中介包需要后者。卸载lmerTest
可以解决问题。+1。为我工作。尽管有些人可能想知道,在再次调用mediate
之前,我必须重新估计med.fit
和out.fit
(调用detach\u package(lmerTest)
)。可能会节省一些时间。我使用了detach\u package()
函数。
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)