如何解决使用脆弱项测试考克斯比例风险假设时的错误
我正在尝试评估包含交互作用项和脆弱项的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)
谢谢!
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。