我使用以下代码建立了离散时间危险模型:
(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”命令执行,该命令是“生存”包的一部分。
有人知道在R中对我的离散时间危险模型进行此测试的可行性吗?