问题描述
我建立了一个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 (将#修改为@）