Thanks @Kelsey_Thompson !
I am sorry to be annoying, but I still have doubts about how to do that.
That question explains how to convert the coefficient to log2fold change in case of a glm with poisson distribution. But I usually use either 'lm" or “negative binomial”. Would that be the same?
I think in case of lm, I do not need to do the following step:
fc<- exp(fit$coefficients[2]) ## Antilog coef #2
What about negative binomial?
Thanks for considering this implementation for the future.
I think it is a good idea to have it has fold change, so it is better quantifiable.
Now we have this coefficient but we do not really know how to biologically interpret it.
I thank you in advance and for your awesome work!!
Gabri