title: Limma moderated t test VS stat_compare_means: why do the values differ?
body: <p>I've done differential expression on a huge list of genes and I'm selecting only the significant ones among an external list of genes that I want to study.
  When I represent my significant ones on the barplot, the p-values I obtained from the moderated t-test (BH correction) with limma, are not the same; for the barplot, I used the ggpubr's function stat_compare_means() and the p-values appearing are totally different and abberrant from the ones obtained with limma.</p>
  <p>Does anyone know if this is normal?</p>
  <p>The genes I plotted should be correct, I checked multiple times.</p>
  <p>Thank you</p>