This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Why doesn't TCGA .rsem.isoforms.results have a length column?

Hi,

I'm analyzing TCGA data and need the rsem output files. I thought those would be the files I download from GDC, which end with .rsem.isoforms.results. However, this is their format (first 5 lines of one of them):

isoform_id  raw_count   scaled_estimate
uc011lsn.1  0.00    0
uc010unu.1  16.09   4.35780241451249e-07
uc010uoa.1  4.00    1.08327337941115e-07
uc002bgz.2  21.91   4.34097970946371e-07

First, I'm not sure why the files have the RSEM output names (.rsem.isoforms.results), but are not actually RSEM outputs, according to this post, because they don't have the length column. Is there a way I can download the correct RSEM outputs from GDC legacy data? Second, I'm not sure if I can fix this manually? From what I understood, the length column is just the length of the transcript, so can I then just find out the length of each isoform ID (e.g. from UCSC) and add a corresponding column?

These files were downloaded from GDC API using Python.

I would really appreciate the help of someone more experienced.

tcga rna-seq alternative splicing rsem r
• 1,235 views • link
updated by Biostar Biostar • written by majd.abdulghani • majd.abdulghani

0 answers

No answers yet.

Log in to answer this question.

More posts like this