Hi

I know I know I know that my question may seem silly to you. And I know I'm not as professional as you are, but please note that I HAVE RESEARCHED TO FIND THE BEST ANSWER, but I don't know if my keywords weren't appropriate because the more I searched on google and specifically on biostars, the more I feel confused.

I have a list of P-values, Log2FoldChanges and Standard Deviations (SD) of Log2FoldChanges for each gene in a time-course RNA-seq study. I'm going to find the DE genes based on the info I have been provided. According to my researches, some say that we should use the adjusted p-values to find the DEG, which means I should use the command "p.adjust" in R to adjust my current p-values. But I have also seen that some people use FDR and FC. On the other hand, everyone recommends a different criteria to find DEG. One says, "find DEG with adjusted p-values<0.1". Another one recommends using adjusted p-values<0.05.

I'm getting confused that,

1- which one is the best solution to find DEG? P-values, adjusted p-values, FDR, FC? I'm not sure whether FDR and FC could be used to find DEG. I just repeat what I have read.

2- whatever you recommend me in question number 1, what criteria do you recommend? 0.05? 0.1? etc.

3- if P-value (or its adjusted value) is enough, what is log2FoldChange and its SD good for? why is it provided beside P-values by the package I'm using?

I apologize if my question is so basic and may bother you. And thanks for your help.

adj.P < 0.05 is typical; you shouldn't need to adjust the p-values if this is done by the package you are using to call differential expression, be it DESeq2 or edgeR or limma or whatever. log2 FC will give you the expected direction of change. Some people filter their DEGs based on |log2_FC| > 0.5 or something but this is usually unnecessary in my experiments.

Thank you so much for the info. The package gives me the p-values, not the adjusted ones, so I think I have to adjust them as you say adj p-values are typical. Is bonferroni adjustment appropriate?

Hi, Just for your information, the FDR

isan adjusted pvalue. Bonferroni is a bit more stringent than FDR but at this point it shouldn't really matter which method you use.Oh, I get it now. FDR is the short term for FDR-adj-pvalue? I spotted an option in "p.adjust" command in R which could use the method "FDR" to adjust P-values. Good point. Thank you so much.

You could think of the cutoff for an adjusted p-value (of 0.05 or 0.1) as a cost function. It depends on what you want to do with the data downstream, and how "expensive" it is to tolerate false positive findings.