我使用以下代码对离散时间风险模型进行了建模:
(logit.full. <-
glmer(event ~
+ a1
+ a2
+ a3
+ obsnum1 + obsnum2 + obsnum3
+ (1 + obsnum1 + obsnum2 + obsnum3 | country_cluster),
family=binomial("logit"), data=data.final))
现在我想测试这个模型,看看在这些类型的离散时间风险模型中继承的比例假设。
对于连续时间风险模型(Cox 模型),这是使用 Schoenfield 残差检验进行测试的。在 R 中,这个测试是用 "cox.zph" 命令执行的,它是 "survival" 包的一部分。
有谁知道为我在 R 中的离散时间风险模型执行此测试的任何可能性?