This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Download data from the TCGA for gene expression analyses in R

I'm trying to download data from the TCGA for gene expression analyses in R, but I'm in doubt if I should use FPKM, FPKM-UQ or counts? When the dataset is in counts, I suppose it's raw data, isn't it? So what's the best unit to compare multiple datasets? I'm planning to use limma or Dseq2 for GE analysis and found that with Dseq2 I need to use count(non-normalised???) data... is that correct? so what's the best package and working strategy?

rna-seq
• 4,104 views • link
written by Lucy • Lucy

Please I need help again. I need to know which column in my table corresponds to the FPKM-UQ values. Thank you for your help!

seqnames        start             end                       width         strand       ensembl_gene_id       external_gene_name      original_ensembl_gene_id
chrX         100627108         10063999             12884            -           ENSG00000000003            TSPAN6                   ENSG00000000003.13
• 0 views • link
updated by GenoMax GenoMax • written by Lucy • Lucy

No value in that output relates to FPKM-UQ

• 0 views • link
written by Kevin Blighe Kevin Blighe

I need a table with only FPKM-UQ and Genes values. How to identify the FPKM-UQ values ​​in an S4 matrix? thank you

• 0 views • link
written by Lucy • Lucy

Sorry, I have no information about which data you have or what you are aiming to do.

• 0 views • link
written by Kevin Blighe Kevin Blighe

Cancer data. I intend to compare the values ​​of FPKM-UQ expressed in normal tissue and primary tumor.

• 0 views • link
written by Lucy • Lucy

I see, but what data have you retrieved right now? If you are relatively new to programming, I may suggest using TCGAbiolinks in R / Bioconductor. If you have no programming experience, then perhaps use cBioPortal

• 0 views • link
written by Kevin Blighe Kevin Blighe

I used TCGAbiolinks. Yes, I am a beginner in programming. I used this code

query <- GDCquery(project = "TCGA-BRCA",
                  data.category = "Transcriptome Profiling",
                  data.type = "Gene Expression Quantification", 
                  workflow.type = "HTSeq - FPKM-UQ")
• 0 views • link
updated by Kevin Blighe Kevin Blighe • written by Lucy • Lucy

1 answer

If you plan to use limma / voom or DESeq2, then the best would be to obtain the raw counts and then follow the guidance for these programs.

You cannot, in any easy fashion, take FPKM expression units and re-process them using either of these programs.

Kevin

• 0 views • link
written by Kevin Blighe Kevin Blighe

Thank you very much

• 0 views • link
written by Lucy • Lucy

Log in to answer this question.

More posts like this