Relative log expression normalization without negative values
1
0
Entering edit mode
2.3 years ago

Hello. I would like to analyze a public RNA-seq data set which the author claims to be normalized with relative log expression normalization. However, in my experience, this generates negative expression values, but all values are positive in his matrix. Am I mistaken or is there a upstream processing not mentioned by the author?

rna-seq • 1.2k views
ADD COMMENT
2
Entering edit mode
2.3 years ago

RLE does not return negative values, it divides the raw counts by the size or normalization factors.

I expect you are thinking of log2 normalization, which will return negative values for counts below 1. This is why it's usually performed as log2(normalized counts + 1), e.g. the normTransform function in DESeq2.

ADD COMMENT
1
Entering edit mode

I think OP assumes that logcounts were standardizes, e.g. Z-scored so that would be relative expression, but I agree that here (based on how the counts look in the file of this GEO entry) that it is probably RLE-normalized values (that is the DESeq2 normalization method @OP).

ADD REPLY

Login before adding your answer.

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