Log transformation of RPKM data for PCA analysis
1
0
Entering edit mode
9.1 years ago
AB ▴ 360

Hi everyone,

I have the RPKM values from my rna-seq dataset and I want to perform PCA on it. Taking log2 of the RPKM data gives infinity because of all the zeros which is not accepted by prcomp. Is there any other way to use the log transformation without removing the zeros?

Thanks in advance

RNA-Seq rpkm PCA log-transformation • 5.8k views
ADD COMMENT
2
Entering edit mode
9.1 years ago
Tky ★ 1.0k

You can just add a value such as 0.01 to all RPKM values for this situations. We routinely add 1 to the raw count value in our analysis.

ADD COMMENT
4
Entering edit mode

Adding 1 is more suitable as log2 of values less than 1 will give negative values.

ADD REPLY
0
Entering edit mode

Thanks for pointing out this, I thought it is the log ratio between the two conditions. Adding 1 make sense.

ADD REPLY

Login before adding your answer.

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