如何解决如何从包括tt项与时间有关的协变量的Cox PH模型绘制估计生存率?
我建立了一个Cox PH模型,其中包括几个协变量,包括基线生活质量(QoL)得分。结果就是生存。由于PH假设不适用于基线QoL与生存之间的关联,因此我为基线QoL添加了时间转换功能(即tt函数)。 现在,我想预测几个基线QoL分数的存活率,而我将平均值用于其余的协变量。 我的问题是survfit函数无法处理带有tt项的coxph模型。
我在stackoverflow(Plotting estimated HR from coxph object with time-dependent coefficient and splines)上找到了一个相关的线程,但是不幸的是,它不包含我的问题的答案。在该线程中,为包括样条线(不是我的样条线)的模型提供了一种解决方案,其目的是绘制估计的HR(我的目的是绘制生存期)。
非常感谢您的帮助!
相关代码:
model <- coxph(Surv(time,vitalstatus) ~ metastasis + age + stage + Bl_QoL + tt(Bl_QoL),data=ds,tt=function(x,t,...) x*t)
Bl.QoL.60 <- data.frame(Bl_QoL=60,metastasis=2,55,age=65,stage=0.5)
Bl.QoL.70 <- data.frame(Bl_QoL=70,stage=0.5)
Bl.QoL.80 <- data.frame(Bl_QoL=80,stage=0.5)
plot(survfit(model,newdata=Bl.QoL.60))
lines(survfit(model,newdata=Bl.QoL.70))
lines(survfit(model,newdata=Bl.QoL.80))
这是我得到的错误: survfit.coxph(模型,newdata = Bl.QoL.60)中的错误: survfit函数无法处理带有tt项的coxph模型
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。