Среднее значение всех прогнозов в emmeans - PullRequest
0 голосов
/ 06 августа 2020

У меня проблемы с использованием emmeans для оценки среднего (или средневзвешенного) всех прогнозов. Например, смешанная модель:

library(emmeans)
library(lme4)
m1 <- lmer(mpg ~ 1 + wt + (1|cyl),data=mtcars)

Фиксированные эффекты «wt» успешны:

emmeans(m1,specs="wt")
   wt emmean   SE   df lower.CL upper.CL
 3.22   20.2 1.71 1.83     12.1     28.3

Однако для расчета среднего прогнозов использовались следующие ранее работавшие (~ 12 месяцев a go), но теперь не удается:

emmeans(m1,specs="1")
NOTE: Results may be misleading due to involvement in interactions
Error in `[[<-.data.frame`(`*tmp*`, ".wgt.", value = 1) : 
  replacement has 1 row, data has 0

Та же ошибка возникает для простых линейных моделей. Большое спасибо за любую помощь.

1 Ответ

0 голосов
/ 06 августа 2020

Я думал, что использую текущую версию emmeans (1.4.8), когда у меня были проблемы, описанные в вопросе. Однако на самом деле я мог использовать emmeans 1.4.6 (см. Комментарий Расс Лент ниже). Я вернулся к emmeans v1.4.3, и код заработал. Затем я обновился до текущей версии emmeans (1.4.8), и код продолжал работать. Скорее всего, причина в том, что я использовал emmeans 1.4.6, в котором была известная ошибка. Пожалуйста, см. эту запись github для получения дополнительной информации.

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...