次のコマンドを使用して、ランダム係数の離散時間ハザード モデルを推定しました。
(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 クラスには $ 演算子が定義されていません」。
エラーがどこから来たのか、どのように解決できるのか、または逸脱残差を代わりに計算する方法を知っている人はいますか?
どんな助けでも大歓迎です!