Extract gene expression values from GDCprepare object
1
0
Entering edit mode
enigmargs • 0
@enigmargs-12267
Last seen 7.2 years ago

Hello, I am trying to download gene expressions values from TCGA-GBM project from GDC using TCGAbiolinks using the example code below;

query <- GDCquery(project = "TCGA-GBM", legacy = TRUE,
        data.category = "Gene expression", data.type = "Gene expression quantification",
        platform = "Illumina HiSeq", file.type = "normalized_results", experimental.strategy = "RNA-Seq")

GDCdownload(query); GDCprepare(query, save=T)

But, I am unable to save the object from GDCprepare. Should the object from GDCprepare result in a dataframe? (It's not a dataframe when I checked it). If so how to access the gene expression values within them? Is it possible to extract only gene expression and patient ID from S4 object resulting from GDCprepare

Any help is much appreciated.

 

tcgabiolinks tcga tcgadownload gdcprepare gdc • 3.6k views
ADD COMMENT
0
Entering edit mode
@tiago-chedraoui-silva-8877
Last seen 3.6 years ago
Brazil - University of São Paulo/ Los A…

Hi,

The default GDCprepare output is a SummarizedExperiment (SE) object, but you can also set to output the data frame.

The main difference is the the SE  has gene information (start, end and others) and sample information.

There are some examples here: https://www.bioconductor.org/packages/devel/bioc/vignettes/TCGAbiolinks/inst/doc/download_prepare.html

Best regards,

Tiago

ADD COMMENT

Login before adding your answer.

Traffic: 954 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6