Question: GDC_Gene expression quantification values for matched tumor normal
0
gravatar for kirti.gondkar
4.1 years ago by
kirti.gondkar0 wrote:

Hi All,

i want to download Gene expression quantification data for tumor and matched normal pairs from GDC portal. Can someone suggest how can I do d same.

Thank you in advance

rna-seq • 1.3k views
ADD COMMENTlink modified 4.0 years ago by thomaskuilman790 • written 4.1 years ago by kirti.gondkar0
0
gravatar for thomaskuilman
4.0 years ago by
thomaskuilman790
thomaskuilman790 wrote:

In R you can use the fantastic TCGAbiolinks package. You can find a thorough explanation here at the Bioconductor webpage. I have used the following script to do just what you want to do:

library("TCGAbiolinks")
library("SummarizedExperiment")

DIRPREFIX <- "/PATH/"

## Download raw counts for indicated sets
## Example of TCGA.data.sets: c("TCGA-LUAD", "TCGA-SKCM")
for (id in TCGA.data.sets){
    query <- GDCquery(project = id,
        ## get data.category using for instance TCGAbiolinks:::getProjectSummary("TCGA-LUAD")
        data.category = "Transcriptome Profiling",
        ## see link above for data.types and workflow.types
        data.type="Gene Expression Quantification",
        workflow.type = "HTSeq - Counts")
    GDCdownload(query, method = "client", directory = "~/Downloads/")
    data <- GDCprepare(query, directory = "~/Downloads/")
    save(data, file = file.path(DIRPREFIX, paste0(id, ".Rdata")), compress = "xz")
}
ADD COMMENTlink modified 4.0 years ago • written 4.0 years ago by thomaskuilman790
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1575 users visited in the last hour