使用脆弱项测试考克斯比例风险假设时的错误

如何解决使用脆弱项测试考克斯比例风险假设时的错误

我正在尝试评估包含交互作用项和脆弱项的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,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)

谢谢！