2

我使用以下代码对离散时间风险模型进行了建模:

(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 中的离散时间风险模型执行此测试的任何可能性?

4

0 回答 0