0

我正在尝试在模型中使用动态协变量（取自用户输入）lmer。这工作正常，但应用emmeans不再工作。

例子：

fit1 = lmer(get(uservar) ~ Intervention + period + sequence + age + (1|ID),
            data = mydata, REML=T)

和

emmeans(fit1, list(pairwise ~ Intervention), adjust = "tukey")

运行没有问题。

lmer(get(uservar) ~ Intervention + period + sequence + get(covariable) + (1|ID),
     data = mydata, REML=T)

"age"从用户输入弹出框中获取协变量也可以。

但：

emmeans(fit1, list(pairwise ~ Intervention), adjust = "tukey")

不再运行。它在说：

变量有不同的长度（为“协变量”找到）

但事实并非如此。如果我检查mydata[,covariable]它的长度与其他变量相同。

有什么建议该怎么做吗？

4

2 回答 2

1

emmeans()必须能够重建预测变量的数据，以获得诸如协变量的平均值和因子水平之类的东西。如果您将用户指定的数据保存在数据框中并使用该数据框拟合模型，它将起作用。但请记住，如果您随后更改该数据框，例如使用另一个用户输入，emmeans()即使您不重新调整模型，它也会更改结果，因为它无法知道数据已更改。

于 2019-05-12T02:36:46.850 回答
0

我找到了一个非常简单的解决方案：

fixed.part <- "get(usevar) ~ Intervention + Periode + Sequenz +  (1|ID)"
random.part <- covariable

fit1 <- lmer(formula(paste(fixed.part, random.part, sep = " + ")), 
 data =testdat_long_cc, REML=TRUE)
emmeans(fit1)
于 2019-05-13T14:09:38.100 回答