Disease Enrichment Multiple Comparison Correction

Multiple Comparison Correction

labs <- c("ADHD", "AUT", "BPD", "DEP", "SZ")
p <- c(ADHD, AUT, BPD, DEP, SZ)

bon <- (p.adjust(p, method = "bonferroni"))
names(bon) <- labs

bh <- (p.adjust(p, method = "BH"))
names(bh) <- labs

cbind(round(p,3),round(bh,3), round(bon,3))

      [,1]  [,2]  [,3]
ADHD 0.405 0.405 1.000
AUT  0.014 0.017 0.068
BPD  0.009 0.015 0.045
DEP  0.007 0.015 0.036
SZ   0.000 0.000 0.000

cbind(p, bh, bon)

                p           bh          bon
ADHD 4.046510e-01 4.046510e-01 1.000000e+00
AUT  1.352746e-02 1.690933e-02 6.763732e-02
BPD  8.945785e-03 1.490964e-02 4.472893e-02
DEP  7.298841e-03 1.490964e-02 3.649421e-02
SZ   1.586454e-05 7.932272e-05 7.932272e-05