使用TCGAbiolinks查询与下载TCGA数据

#install package
if (!requireNamespace("BiocManager", quietly = TRUE))
  install.packages("BiocManager")
if (!requireNamespace("BiocManager", quietly = TRUE))
  install.packages("TCGAbiolinks")
library(TCGAbiolinks)
library(SummarizedExperiment)
library(dplyr)
library(DT)

#TCGAbiolinks 下载数据,具体结合实际来写 

query <- GDCquery(project = "TCGA-LIHC",
                  legacy = FALSE,
                  experimental.strategy = "RNA-Seq",
                  data.category = "Transcriptome Profiling",   #转录组数据
                  data.type = "Gene Expression Quantification",  #代表输入的是gene表达数据
                  workflow.type = "HTSeq - Counts")
#查询数据
datatable(getResults(query))
#下载数据
GDCdownload(query)

[1]R语言|TCGAbiolinks 包系列(1)——检索并获取数据https://www.jianshu.com/p/e405073f1642
[2]TCGAbiolinks下载TCGA数据(更新版本)
https://www.jianshu.com/p/8334f80fc248

你可能感兴趣的:(R与生信)