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. 

あなたのモデルには2つの要素があります。

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