source("https://bioconductor.org/biocLite.R")
biocLite("TCGAbiolinks")
library(TCGAbiolinks)
library(dplyr)
library(DT)
query.exp <- GDCquery(project = "TCGA-BRCA",
data.category = "Transcriptome Profiling",
data.type = "Gene Expression Quantification",
workflow.type = "HTSeq - FPKM-UQ")
r = query.exp$results[[1]]
dat = r[,c(4,8,9,15,26)]
'Data analysis > TCGA' 카테고리의 다른 글
[TCGA]RNA-seq download + edgeR(analysis) 완벽 분석 코드 !!! (0) | 2018.07.24 |
---|---|
[TCGA] RNA-seq + clinical data download (0) | 2018.07.24 |
TCGA barcode 해석 (0) | 2018.07.05 |
[TCGA] clinical data generation (0) | 2018.06.26 |
[TCGA] clinical data parsing code (0) | 2017.05.15 |