如何解决用R绘制竞争风险列线图
我试图根据Zhongheng等人的论文绘制列线图。我使用自己的数据集“在存在竞争风险的情况下进行生存分析的线型图”,该数据集与论文中使用的数据集非常相似:3307例患者,502个关注事件(骨折),264例死亡(竞争风险)。时间是骨骼检查(DXA)后的几年,状态以这种方式编码(O =删失,1 =骨折,2 =死亡)。
我的数据集的格式已转换为包含SAS的审查权重,如下所示:
我使用了这种语法(摘自本文):
library(rms)
ddist<-datadist(df)
options(datadist='ddist')
mod<-cph(Surv(Tstart,Tstop,status==1)~age3cc+htot_bmd+tot_bmd+neck_bmd+amof+BMI2C+
AlcFR+malchronFR+malchronFR*BMI2C,data=df,weight=Weightc,subset=failcode==1,surv=T)
nom.sur<-nomogram(mod,fun=function(x) 1-Surv(5,x),fun.at=c(0.05,0.1,0.15,0.2,0.25),funlabel=c("5-year incidence risk"),lp=F)
plot(nom.sur,fun.side=c(1,1,3,1))
但是我收到一条错误消息Error in Surv(5,x) : Time and status are different lengths
。
如何解决此错误?
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。