pan-cancer泛癌单基因分析问题之合并TCGA和GTEx

#=======================================================


#=======================================================


library(GenomicDataCommons)

setwd('D:\\SCIwork\\F33\\TCGA')

rm(list=ls())


library(dplyr)

library(TCGAbiolinks)

library(dplyr)

library(DT)

library(SummarizedExperiment)

library(stringr)

#=======================================================


#=======================================================

cancer  <- TCGAbiolinks:::getGDCprojects()$project_id

cancer <- str_subset(cancer, "TCGA")

cancer <- sort(cancer)




for (i in 1:33) {
  cancer_select <- cancer[i]
  print(cancer_select)
  #下载rna-seq的counts数据
  suppressMessages({
    query <- GDCquery(
      project = cancer_select,
      data.category = "Transcriptome Profiling",
      data.type = "Gene Expression Quantification",
      workflow.type = "HTSeq - Counts")  })
  
  
  if (is.null(query)){
    print(paste0("No Counts data of solid normal tissue for ", cancer_select ))
  } else{
    
    GDCdownload(query, method = "api", 
                files.per.chunk = 300)
    expdat <- GDCprepare(query = query, save = TRUE,
                         save.filename = paste0(cancer_select,".rda"))
    count_matrix=assay(expdat)
    write.csv(count_matrix,
              file = paste( cancer_select,"Counts.csv",
                            sep = "-"))}}