R中偏差残差计算的误差

R中偏差残差计算的误差,r,R,我使用以下命令估计了一个随机系数离散时间危险模型: (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),

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

(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"))

(logit.model这与您正在使用的
lme4
版本有关。您可以使用以下工具检查您的版本:

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

这是一个开发版本(似乎是这样),然后使用以下工具获取最新的开发版本:

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

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

好的,非常感谢您的回答,这似乎就是问题所在。现在我有一个后续问题。不幸的是,估计模型中的观测数量小于原始数据集。因此,我得到了错误消息:“data.frame(…,check.names=FALSE)中的错误:参数意味着不同的行数:70416621。“你知道解决这个问题的优雅方法吗?我想这与
glmer
na.action
参数有关,默认情况下是
na.ommit
。选中
?不适用。排除