TCGA数据下载困难

Question

我正在尝试下载TCGA数据，但出现此错误：

summaryMaf（maf = maf，anno = ClinicalData，chatty =详细）：在提供的临床中找不到Tumor_Sample_Barcode列数据。如果包含样品名称的列重命名为Tumor_Sample_Barcode必要。

这是我的代码：

library("TCGAbiolinks")
library("tidyverse")
library(maftools)
query <- GDCquery(   project = "TCGA-LIHC",  
                     data.category = "Clinical", 
                     file.type = "xml", 
                     legacy = FALSE)
GDCdownload(query,directory = ".")

clinical <- GDCprepare_clinic(query, clinical.info = "patient",directory = ".")
#getting the survival time of event data
survival_data <- as_tibble(clinical[,c("days_to_last_followup","days_to_death","vital_status","bcr_patient_barcode","patient_id")]) 
survival_data <- filter(survival_data,!is.na(days_to_last_followup)|!is.na(days_to_death))  #not both NA
survival_data <- filter(survival_data,!is.na(days_to_last_followup)|days_to_last_followup>0 &is.na(days_to_death)|days_to_death > 0 ) #ensuring positive values
survival_data <- survival_data[!duplicated(survival_data$patient_id),]  #ensuring no duplicates


dim(survival_data) #should be 371


maf <- GDCquery_Maf("LIHC", pipelines = "muse")
#maf <- GDCquery_Maf("LIHC", pipelines = "somaticsniper")

#clin <- GDCquery_clinic("TCGA-LIHC","clinical")
#print(clin )



laml = read.maf(
  maf,
  clinicalData = clinical,
  removeDuplicatedVariants = TRUE,
  useAll = TRUE,
  gisticAllLesionsFile = NULL,
  gisticAmpGenesFile = NULL,
  gisticDelGenesFile = NULL,
  gisticScoresFile = NULL,
  cnLevel = "all",
  cnTable = NULL,
  isTCGA = TRUE,
  vc_nonSyn = NULL,
  verbose = TRUE
)

Answer 1