Прогнозируемые значения для glmer AICcmodavg

Question

Я знаю, что с помощью AICcmodavg можно получить прогнозируемые значения (в исходном масштабе ~ вероятность) и их SE для фиксированных эффектов, но я пытаюсь безуспешно... Может ли кто-нибудь помочь мне с этим? заранее спасибо

library(lme4)
(gm1 <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd),
              data = cbpp, family = binomial))
fixef(gm1)

library("AICcmodavg")

predictSE(gm1, 
          newdata=as.data.frame(period=c("period1","period2","period3","period4")), 
          type="response", 
          se.fit=TRUE, 
          level=0, 
          print.matrix=F)

Answer 1

Было бы лучше прочитать levels(cbpp$period), и не as.data.frame() но data.frame()

levels(cbpp$period)
# [1] "1" "2" "3" "4"

predictSE(gm1, 
          newdata = data.frame(period=c("1", "2", "3", "4")),
          type = "response", 
          se.fit = TRUE, 
          level = 0, 
          print.matrix = F)

[Изменено]

простой метод поиска причины ошибки

fit <- ...(..., data = df)

predictSE(fit, newdata = df)
predictSE(fit, newdata = ...)

# If 1st predictSE() doesn't run, it means the model causes error.
# If 1st runs but 2nd doesn't, it means it is due to newdata.

Если ваша модель имеет два фактора;

newd <- expand.grid(name1 = levels(df$name1), name2 = levels(df$name2))
predictSE(fit, newdata = newd)
 # pred <- predictSE(fit, newdata = newd)
 # cbind(newd, pred)            # help to interpret