predicted values for glmer AICcmodavg

Question

I know that one can get the predicted values (in original scale ~ probability) and their SE for the fixed effects with AICcmodavg, but I'm trying without any success... Can someone help me with that? Thanks in advance

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)

cuttlefish44 cuttlefish44 · Accepted Answer · 2016-09-22T12:59:18

It would be better to read levels(cbpp$period), and not as.data.frame() but 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)

[Edited]

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

predicted values for glmer AICcmodavg

1 Answers