r - 使用脆弱项测试 cox 比例风险假设时出错

Question

我正在尝试评估 Cox PH 模型的比例风险假设，该模型包括交互项和脆弱项。当我调用 cox.zph 时，我收到此错误：

imatr[kk, kk] 中的错误：下标越界

当我从模型中删除脆弱项时，我没有得到错误，并且在简化模型中，我的所有固定因子都符合比例风险假设。我是否也需要在我的完整模型上运行 cox.zph？如果是这样，我怎样才能绕过我收到的错误消息？

此代码应重新创建我收到的错误消息：

library(survival)

set.seed(80)

fake.data <- data.frame(event = c(rbinom(75, 1, 0.7), rbinom(75, 1, 0.4)),
                        time = runif(150, min=0, max=3),
                        time.period = rep(c("pre", "peri", "post")),
                        condition = rep(c("fed", "not fed")),
                        momID = floor(runif(150, min=1, max=35)))



phmodel <- coxph(Surv(time, event) ~ time.period*condition + frailty(momID) , data = fake.data )

assump.check <- cox.zph(phmodel)

谢谢！

标签： rsurvival-analysiscox-regressionsurvival