R - 使用 effect() 函数总结混合模型 (lmer) 估计的“数据”参数类型错误答案

【问题标题】：R - 'data' argument is of the wrong type using effect() function to summarise mixed models (lmer) estimateR - 使用 effect() 函数总结混合模型 (lmer) 估计的“数据”参数类型错误
【发布时间】：2020-12-17 04:44:14
【问题描述】：

我对 R 很陌生，最近我的任务是使用 ggplot 来可视化 lmer 模型的结果。为此，我首先尝试总结混合模型估计并将其转换为数据框。

我的代码：

model <- lmer (outcome ~ group*time + (1|ID), data)


model.eff <- effect("group*time", model) #which supposedly summarises the mixed model estimates.

但这就是我卡住的地方。我不断收到此错误消息 terms.formula(formula, data = data) 中的错误：'data' 参数类型错误

阅读后，我认为问题可能出在我的“模型”类上，但我不知道如何解决这个问题。任何帮助将不胜感激！

【问题讨论】：

标签： r model lme4 effect mixed

【解决方案1】：

有一个包“broom”可以更轻松地处理模型输出。它只需要您将模型输入到函数“tidy”中：

library(broom)
model <- 
 lmer (outcome ~ group*time + (1|ID), data)

model.eff <- tidy(model)

然后您将获得这种样式的输出（来自我的数据的示例，因为您没有提供任何示例数据）：

  effect   group    term            estimate std.error statistic
1 fixed    NA       (Intercept)        6.14      4.68      1.31 
2 fixed    NA       PFS_days          -0.561     0.573    -0.981
3 ran_pars sex      sd__(Intercept)    1.36     NA        NA    
4 ran_pars Residual sd__Observation    3.50     NA        NA

【讨论】：