使用http://www.oncolnc.org/
1.打開(kāi)網(wǎng)站，輸入“favorite gene”
2.選擇cancer，P-K一下
image.png
3.輸入percentile
image.png
可見(jiàn)，高TP53表達(dá)影響B(tài)RCA患者的長(zhǎng)期生存
同樣地，我們自己畫(huà)一下吧
1.導(dǎo)出數(shù)據(jù)
image.png
2.運(yùn)行R代碼
options(stringsAsFactors = F)
a=read.table('BRCA_7157_25_25.csv',sep = ',',fill = T,header = T)
dat=a
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)
ggsave('survival_TP53_in_BRCA_TCGA.png')
image.png
沒(méi)差
參考來(lái)源：生信技能樹(shù)
友情鏈接：
課程分享
生信技能樹(shù)全球公益巡講
（https://mp.weixin.qq.com/s/E9ykuIbc-2Ja9HOY0bn_6g）
B站公益74小時(shí)生信工程師教學(xué)視頻合輯
（https://mp.weixin.qq.com/s/IyFK7l_WBAiUgqQi8O7Hxw）
招學(xué)徒：
（https://mp.weixin.qq.com/s/KgbilzXnFjbKKunuw7NVfw）
歡迎關(guān)注公眾號(hào)：青島生信菜鳥(niǎo)團(tuán)
