补充：08.任意基因癌症表达量分组的生存分析

# 从网页工具下载感兴趣的基因，保存格式为.csv
# 一个基因在TCGA数据库中的各个癌症的生存分析
# 网页中输入感兴趣的基因，分组(高低表达量)，下载数据，
# Rstudio中读入数据,ggstatsplot包中的代码进行分析
a=read.table('08.任意基因.csv',header = T,sep = ',',fill=T)
colnames(a)
install.packages("ggstatsplot")#安装这个R包
library(ggstatsplot)
dat=a#只是把自己的数据放在这里而已
ggbetweenstats(data=dat,
               x=Group,y=Expression)

# 更改代码中的X值

install.packages("survival")#安装这个R包
install.packages("survminer")#再安装这个R包

# 自定义生存分析图，理解参数，修改参数，应用高级代码
library(ggplot2)
library(survival)
library(survminer)
table(dat$Status)
dat$Status=ifelse(dat$Status=='Dead',1,0)
sfit=survfit(Surv(Days,Status)~Group,data = dat)
sfit
summary(sfit)
ggsurvplot(sfit,conf.int=F,pval = TRUE)