前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >TCGA数据下载分析(1-2):RTCGA包PCA和生存分析

TCGA数据下载分析(1-2):RTCGA包PCA和生存分析

作者头像
Y大宽
发布2018-09-10 12:00:50
1.1K0
发布2018-09-10 12:00:50
举报
文章被收录于专栏:Y大宽

1 PCA plot:pcaTCGA

Plots Two Main Components of Principal Component Analysis

用法如下: pcaTCGA(x, group.names, title = "", return.pca = FALSE, scale = TRUE, center = TRUE, var.scale = 1, obs.scale = 1, ellipse = TRUE, circle = TRUE, var.axes = FALSE, alpha = 0.8, add.lines = TRUE, ...)

代码语言:javascript
复制
expressionsTCGA(BRCA.rnaseq, OV.rnaseq, LIHC.rnaseq) %>%
  dplyr::rename(cohort = dataset) %>%  
  filter(substr(bcr_patient_barcode, 14, 15) == "01") -> BRCA.OV.LIHC.rnaseq.cancer
pcaTCGA(BRCA.OV.LIHC.rnaseq.cancer, "cohort") -> pca_plot
代码语言:javascript
复制
plot(pca_plot)

Rplot.jpeg

2 生存分析kmTCGA()

  • Kaplan-Meier 生存曲线评估乳腺癌和卵巢癌病人中TP53基因突变与生存关系
代码语言:javascript
复制
library(RTCGA.mutations)
# library(dplyr) if did not load at start
library(survminer)
mutationsTCGA(BRCA.mutations, OV.mutations) %>%
   filter(Hugo_Symbol == 'TP53') %>%
   filter(substr(bcr_patient_barcode, 14, 15) ==
   "01") %>% # cancer tissue
   mutate(bcr_patient_barcode =
   substr(bcr_patient_barcode, 1, 12)) ->
  BRCA_OV.mutations

library(RTCGA.clinical)
survivalTCGA(
  BRCA.clinical,
  OV.clinical,
  extract.cols = "admin.disease_code"
  ) %>%
   dplyr::rename(disease = admin.disease_code) ->
  BRCA_OV.clinical

BRCA_OV.clinical %>%
   left_join(
     BRCA_OV.mutations,
     by = "bcr_patient_barcode"
     ) %>%
   mutate(TP53 =
   ifelse(!is.na(Variant_Classification), "Mut","WILDorNOINFO")) ->
  BRCA_OV.clinical_mutations

BRCA_OV.clinical_mutations %>%
select(times, patient.vital_status, disease, TP53) -> BRCA_OV.2plot

kmTCGA(
  BRCA_OV.2plot,
  explanatory.names = c("TP53", "disease"),
  break.time.by = 400,
  xlim = c(0,2000),
  pval = TRUE) -> km_plot
代码语言:javascript
复制
print(km_plot)

Rplot01.jpeg

更多的看这里

本文参与 腾讯云自媒体同步曝光计划,分享自作者个人站点/博客。
原始发表:2018.08.06 ,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
目录
  • 1 PCA plot:pcaTCGA
  • 2 生存分析kmTCGA()
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档