我正在library("mediation")
中使用R中的mediation包。我有一个数据集,正在尝试使用以下代码进行中介分析(变量名已更改):
mod.y <- glmer(congr ~ age + gender + age*gender + wage + ( 1 | FileID), family = binomial(link="logit"), data = data, na.action="na.omit", control=glmerControl(optimizer="bobyqa"))
mod.m <- glmer(wage ~ age + gender + age*gender + ( 1 | FileID), family = binomial(link="logit"), data = data, na.action="na.omit")
fit <- mediate(mod.m, mod.y, sims = 50, treat="age", mediator="wage")
test.TMint(fit)
中介分析可以很好地完成,我可以使用summary(fit)
获取ACME和ADE值,但是当我尝试运行test.TMint
时,出现以下错误:
Error in test.TMint.default(fit) :
currently no test.TMint method exists for the input object.
documentation指出我需要提供一个mediate
对象。过去有没有人经历过?