Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Convert log_eFC to log_2FC in differential gene expression chapter #211

Open
ngvananh2508 opened this issue Jun 13, 2023 · 1 comment
Open
Assignees
Labels
enhancement New feature or request

Comments

@ngvananh2508
Copy link

The formula to convert coef from log_eFC to log_2FC is not equal to coef/log_e(FC).

Let x be the fold change -> coef = log_e(x) -> e^(coef) = x -> coef.log_2(e) = log_2(x)
-> coef := coef.log_2(e)

And when using MAST, if I want to input all the cell types, do I only need to change "summaryCond <- summary(zlmCond, doLRT='groupstim.CD14_Monocytes')" to "....doLRT = list('groupstim.CD14_Monocytes', 'groupstim.B_cells', ...) and expand to all cell types in merge function to still be able to use MAST?

Thank you very much for helping me!

@ngvananh2508 ngvananh2508 added the enhancement New feature or request label Jun 13, 2023
@simonmfr
Copy link

Just FYI, here the MAST developer says "the base of the logarithm is whatever you set it to be when you log transform your counts", which means that you would wrongly do log2(FC)/ln(2) instead of ln(FC)/ln(2).
Best,
Simon

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
enhancement New feature or request
Projects
None yet
Development

No branches or pull requests

3 participants