0

我使用以下命令估计了一个随机系数离散时间风险模型:

(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 类定义 $ 运算符”。

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

任何帮助将不胜感激!

4

1 回答 1

1

这与lme4您使用的版本有关。您可以检查您拥有的版本:

packageVersion("lme4")

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

是否是开发版本,似乎是这样,然后使用以下方法获取最新的开发版本:

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

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

于 2013-09-16T08:47:40.747 回答