This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Question about how to transform RSEM expected_count of TCGA TARGET GTEX to integers?

Hi, there

I downloaded RSEM expected_count data of TCGA TARGET GTEX from UCSC Xena to perform differential expressed analysis. I noticed the unit is log2(x+1), but I failed to turn the data into integers. The input of DESeq2 should be a matrix with integers. Any suggestions? Thanks in advance!

Best,

Issac

rna-seq
• 3,979 views • link
updated by i.sudbery i.sudbery • written by dz2353 dz2353

1 answer

log2(x+1) data is not suitable for analysis with DESeq2 (nor edgeR/voom). So you would at least have to convert it back to counts by doing (2^x) -1. You could then simply round to the nearest integer. Alternatively you can download un-tarnsformed RSEM output from Firehose at https://gdac.broadinstitute.org/

• 0 views • link
written by i.sudbery i.sudbery

Yep, I did (2^x)-1 but the result is not an integer format cause the scale. I will round them. Thanks a lot!

• 0 views • link
written by dz2353 dz2353

Log in to answer this question.

More posts like this