Since you are applying multiple univariate analyses, you are testing multiple hypotheses and you should correct for multiple testing. However, Bonferroni is just one way to correct, which controls the family-wise error rate, FWER. The FWER is the probability of at least one false positive. Even among FWER controlling procedures, it is very conservative. (Applying the Bonferroni correction, your adjusted p-value would be min{4p, 1}
, where pp
is your nominal p value.)
Another optionA less conservative error rate would be to correct the p-values with the false discovery rate (FDR). The false discovery rate is the expected proportion of false positives out of those mediators you call significant.
If you use R, you can get information with references on different correction methods with
?p.adjust
This documentation is also available online.
Incidentally, it so happens that I actually just posted a freely available article on more involved methods for this exact topic recently, but these approaches probably aren't worth it with only 4 hypotheses.