expr_gene$Overall_Survival_Status <- as.numeric(expr_gene$Overall_Survival_Status) diff <- survdiff(Surv(days_to_last_followup,Overall_Survival_Status)~low_high,data = expr_gene)
library(survival) library(rms) library(tidyverse) data <- read.csv("Liver_clinical_arrange.csv") %>% mutate(age = ageNum, riskScore = riskScoreNum) View(data) dd<-datadist(data) options(datadist=dd) f <- cph(Surv(futime_year,fustat) ~ age + stage+T+riskScore,data=data, x=T, y=T, surv=T) surv <- Survival(f) nom <- nomogram(f,fun = list(function(x)surv(1,x),function(x)surv(3,x),function(x)surv(5,x)), lp=F, funlabel = c("1-year Survival Probability","3-year Survival Probability","5-year Survival Probability")) png(filename = "nomogram.png",width = 9,height = 8,res = 300,units = "in") plot(nom) dev.off()
图片alt
http://www.360doc.com/content/18/0714/09/46405145_770258087.shtml用于术前预测结直肠癌淋巴结转移的放射组学列线图的开发和验证
使用列线图对结肠癌复发的个体化预测
m6A相关lncRNA生信挖掘套路解析