This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Using TCGABioLinks to query multiple cancer studies

I have a list of cancer samples that I wish to download transcriptome data for using TCGABiolinks. I use the GDCquery() function for downloading data, however my list spans multiple TCGA projects. Here is my code 

query <- GDCquery(project = "", data.category = "Transcriptome Profiling",
                  data.type = "Gene Expression Quantification", 
                  workflow.type = "HTSeq - Counts",
                  barcode = acc_numbers)

I am unsure what to put in the project argument. I do not think there is an argument for the TCGA pan cancer atlas, so does anyone know what I should do here?

rna-seq r tcga
• 2,018 views • link
updated by GenoMax GenoMax • written by ivykosater • ivykosater

1 answer

https://rdrr.io/bioc/TCGAbiolinks/man/GDCquery.html

See example from the manual to perform query on more than one study:

query.met <- GDCquery(project = c("TCGA-GBM","TCGA-LGG"),
                      legacy = TRUE,
                      data.category = "DNA methylation",
                      platform = "Illumina Human Methylation 450")

If you want all to perform across all projects then you can replace list of projects with TCGAbiolinks:::getGDCprojects()$project_idor filter it as needed.

• 0 views • link
updated by GenoMax GenoMax • written by Ashutosh Pandey Ashutosh Pandey

Log in to answer this question.

More posts like this