R中的事件时间调解

Question

我正在研究R中的时间到事件调解分析，并且一直难以让程序运行。

我试图通过问卷调查来检查治疗前后的变化，作为干预后医疗事件中群体差异的调解者。有一个问题是有几个辍学，所以我试图在下面的代码中估算这些值。我一直收到关于中间声明的以下信息：

Error in formula.default(object, env = baseenv()) : invalid formula

代码如下。任何意见都非常感谢！

med.fit <- with(imp1, lm(posttx~ pretx+ Group)) 
out.fit<- with(imp1, survreg(Surv(eventtm,event) ~ Group + posttx, dist="lognormal"))
med.out <- mediate(med.fit, out.fit, boot=TRUE, treat = "Group", mediator = "posttx", sims = 100)
summary(med.out)