Hi everyone,
I have a microarray set of 26933 genes and 6 samples (3 controls and 3 mutants). I calculated differential expressed genes using limma's topTable() function. I got a result of this format:
> head(res.new)
logFC AveExpr t P.Value adj.P.Val B 17291881 -2.157936 9.021125 -41.49770 3.218440e-09 3.829474e-05 10.82552 17357688 -2.364349 7.160051 -41.15755 3.397181e-09 3.829474e-05 10.79753 17224251 1.820522 10.665778 37.01446 6.816679e-09 3.829474e-05 10.41426 17239845 -2.025497 10.054660 -35.87667 8.366167e-09 3.829474e-05 10.29350 17304523 -1.830077 10.117184 -35.36567 9.191590e-09 3.829474e-05 10.23680 17229782 2.716749 7.911040 35.24054 9.407763e-09 3.829474e-05 10.22267
The logFC values are log2(FC). How do I convert the values (positive and negative) to normal fold changes?
I tried this:
#take 2^(absolute(log2FoldChange))
> res.new$FC<-2^(abs(res.new$logFC))
#change the sign of FC according to log2FC
> res.new$FC<-ifelse(res.new$logFC<0,res.new$FC*(-1),res.new$FC*1)
But I am not sure if I am correct. Please let me know.
Thanks!