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
。选中?不适用。排除。