Question

我使用以下命令估计了随机系数离散时间危险模型：

(logit.model <-
   glmer(event ~ 
       a1
     + a2
     + a3
     + a4
     + high * I(a5 - 2)
     + midup * I(a5 - 2)    
     + obsnum1 + obsnum2 + obsnum3
     + (1 + obsnum1 + obsnum2 + obsnum3 | country_cluster),
     family=binomial("logit"), data=data.final, 
     verbose=TRUE, control=list(maxIter=400)))

现在我想使用以下命令计算偏差残差：

    data.final.r <- cbind(data.final,
                          dev.res = residuals(logit.model,
                                              type="deviance"))

不幸的是，我收到以下错误命令：

“没有为此S4类定义$运算符”。

有谁知道，错误可能来自哪里，如何解决或者如何计算偏差残差？

非常感谢任何帮助！

Answer 1

这与您正在使用的lme4版本有关。您可以查看您的版本：

packageVersion("lme4")

这是您从CRAN获得的稳定版本，然后使用update.packages()或install.package(lme4)获取最新版本

它是一个开发版本，似乎就是这种情况，然后使用以下方法获取最新的开发版本：

install.packages("lme4",repos="http://lme4.r-forge.r-project.org/repos")

警告：如果您使用的是稳定版/ CRAN版，则更改为开发版可能会破坏代码中的内容。

R中偏差残差计算出错

1 个答案: