用coxph函数构建cox模型绘制列线图后,如何计算每个患者的列线图得分
我的代码如下
##列线图=========================================================================
library(mstate)
set.seed(1000)
df=crprep("Survival_months","CVD", #对数据train进行加权
data=train,trans=c(1,2), #trans = (终点事件, 竞争风险事件)
cens=0,id="id", #cens删失值/截尾,设置id
keep=c("Age","Marital","Race","Year_of_diagnosis","PSA.10","Gleason_score")) #数据集中需要保留的协变量(纳入的自变量)
df$Time=df$Tstop-df$Tstart
library("dplyr")
df_c <- df %>% filter(.,failcode == 1)
m.crr<- coxph(Surv(Time,status==1)~Gleason_score+PSA.10+Year_of_diagnosis+Race+Marital+Age,
data=df_c,
weight=weight.cens,
subset=failcode==1)
library(regplot)
nom_c <- regplot(m.crr,
failtime = c(36,60,96), prfail = T, droplines=T,points=T,
plots=c("density","no plot"),
title = "Nomogram")
我想要做风险分层,但是我不知道这个方法如何计算患者各自总分