Question about how to transform RSEM expected_count of TCGA TARGET GTEX to integers?
1
1
Entering edit mode
3.7 years ago
dz2353 ▴ 120

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 • 2.4k views
ADD COMMENT
3
Entering edit mode
3.7 years ago

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/

ADD COMMENT
0
Entering edit mode

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

ADD REPLY

Login before adding your answer.

Traffic: 1292 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