生存分析

R安装
install.packages(“surviaval”)
install.packages(“survival”)
library(KMsurv)
library(survival)

测试

library(“survival”)
library(“survminer”)
data(“lung”)
fit <- survfit(Surv(time, status) ~ sex, data = lung)
surv_diff <- survdiff(Surv(time, status) ~ sex, data = lung)
ggsurvplot(fit,ggtheme=theme_bw(),pval=TRUE,xlab=c(“Overall days”),linetype=c(“dashed”),title=c(“Lung cancer”),legend.labs = c(“Male”, “Female”),legend.title = “Sex”) + ylab(survival probility)

最后如果想对不同组的生存率进行假设检验（Log-Rank test）的话，可以用survdiff()函数，Log-Rank test是无参数检验，近似于卡方检验，零假设是组间没有差异.