Блог

Я совсем новичок в R, и недавно мне было поручено использовать ggplot для визуализации результатов модели lmer. Для этого я сначала пытаюсь суммировать и преобразовать оценки смешанной модели в фрейм данных.

мой код:

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


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

Но вот где я застрял. Я продолжаю получать это сообщение об ошибке
Ошибка в терминах.формула (formula, data = данные) : аргумент ‘data’ имеет неправильный тип

После прочтения я понял, что проблема может заключаться в классе моей «модели», но я не уверен, как исправить эту проблему. Любая помощь была бы оценена!

Существует пакет «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