我使用以下命令估计了一个随机系数离散时间风险模型:
(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 类定义 $ 运算符”。
有谁知道,错误可能来自哪里,如何解决或如何计算偏差残差?
任何帮助将不胜感激!