时间戳

首页 > 解决方案 > R - 使用 effect() 函数总结混合模型 (lmer) 估计的“数据”参数类型错误

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

问题描述

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

我的代码：

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


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

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

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

标签： rmodellme4effectmixed

解决方案

有一个包“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

推荐阅读