预测glmer对象中的错误

Question

我正在尝试predict一个glmer对象，其参数为newparams，如下所示：

library(lme4)

    set.seed(123)
    df<-data.frame(id=sample(LETTERS[1:10], 50, T),
                   y=rbinom(50, 1, 0.3),
                   x1=rbinom(50, 1, 0.5),
                   x2=as.integer(rnorm(50, 40, 5)))
    df<-df[order(df$id),]

    fit<-glmer(y~x1+x2+(1|id), data=df, binomial)
    predict(fit, newparams=list(theta=getME(fit, "theta"), beta=fixef(fit)))
    Error in predict(fit, newparams = list(theta = getME(fit, "theta"), 
unused argument (newparams = list(theta = getME(fit, "theta"), beta = fixef(fit)))

愿这里有人知道是什么原因吗？

Answer 1

我认为你应该设置newdata参数。这很好用：

predict(fit, newdata=df,
  newparams=list(theta=getME(fit, "theta"), beta=fixef(fit)))

PS

使用你的代码我得到一个不同的错误：

Error: inherits(fr, "data.frame") is not TRUE