natural log used in Seurat for differential expression
1
0
Entering edit mode
5.8 years ago
nix • 0

Hi all, I have recently taken up some single cell analysis, and would like to use Seurat for differential expression.

I notice they use log fold change (natural log) to transform fold changes. I would have thought that log2 is more interpretable (and more often-used) than natural log in the case of fold change. Is there an advantage to using the natural log over log2 ?

Thanks in advance for any advice!

RNA-Seq seurat differential expression • 5.9k views
ADD COMMENT
1
Entering edit mode

I agree that, in bioinformatics, log2 is so often used and easily interpretable that should be the default option of any tool.

ADD REPLY
3
Entering edit mode
5.8 years ago

I guess, it's just a matter of choice of Seurat developers. And in fields other than bioinformatics, natural log is the one most widely used - one of the major reasons is that the exponential function (e^x), which is the reciprocal "cousin" of natural log, comes naturally in many calculations - like differentials and integrals and by extension in the sums of infinite series. No surprise that e is also defined as the limiting sum of the infinite series (1 + 1/n)^n as n approaches infinity

Changing the base of log amounts to multiplying by a constant factor

log2(N) = log2(e) * ln(N)

So it doesn't affect the performance of downstream analyses like clustering, PCA etc as most of these analyses are invariant to scaling transformations (i.e they do not change if all the values are multiplied by a constant number).

PS: The best way to convince is to try yourself. Try to run the analyses by multiplying with a constant factor and see if the downstream results change. You may add a different scaling factor than 10000 while normalizing

NormalizeData(object, assay.type = "RNA", normalization.method = "LogNormalize", scale.factor = 10000, display.progress = TRUE)

ADD COMMENT

Login before adding your answer.

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