Hi,
I am using "seqlm" package to find DMRs. Seqlm returns statistical parameters such as p.value, Bonferroni, FDR, ...for each segment (deferentially methylated region) .In one of these analysis, most segments had meaningful p.value(<0.05) but their FDR are more than 0.5 and Bonferroni is 1 or near 1. How can i filter correct value based on statistic? Is there a problem with this run?
1 answer
I don't know seqlm package but If you have comparison between multiple segments you need to use multiple correction FDR or Bonferroni, Statistics says that 1 per every 20 segments in multiple comparison will always result with significant raw p-value (e.g. 5/100 t-tests will end up with p<0.05).
If you are checking for significance difference between two segments you can use raw p-value, when you have multiple segments - use p-value from FDR or Bonferroni.
Best,
Agata
Log in to answer this question.
-
Bismark Methylation calling
written by Keem Uarren Eiryll •
-
Bias in FDR distribution
written by Daniel •
-
Does gene raw read count mean and dispersion preserve batch effects if I used raw counts from diffe…
written by aldas02027 •
-
Comparing the whole genelist by GSEA GO BP terms to adjusted pvalue genelist by DAVID GO BP terms
written by minifoog •
-
DMRcate: extractRanges error
written by macsuejacques •
-
Forum: Choosing the appropriate test for differential analysis: Limma or Mann-Whitney (with BH corr…
written by innyus •
-
InParanoid: score and cut off for indentify orthologs
written by Izal •
-
No significant differential methylated regions in control and eutopic samples comparison
written by niutster
-
Filtering DMRs
written by niutster
-
Problem With Fdr Test : Q-Values All < 0.05 For A List Of P-Values That Are Not All Significant.
written by Francois Olivier Hébert •
For omics data most people use FDR corrected p-values, Bonferroni is much too stringent. Don't use raw p-values with omics data.
I think you should learn about multiple testing correction. It is perfectly possible that after correction, you don't get any significant test. This just means that at the level of confidence you selected, you can't reject the null hypothesis (usually of no difference between samples). However, don't make the mistake of thinking that statistical significance means biological relevance.
Nice way to put it. It is important to understand how one performs the multiple testing. The line sentence speaks volumes about correlation of statistical testing with biological significance.