I wish to analyze microarray data for identification of DEGs through limma analysis at the threshold P<0.05 or a fold change >1.5. I written R script as
When I run these script, get some variation in the results, In case of p.value the number of up-regulated and down-regulated genes are found more as compare to lfc. Is my script correct?
Which threshold criteria is better for fetching out of the significant genes? Any suggestion will be highly appreciated.
Thanks in Advance.