编辑代码

>data <- c(0.001, 0.003, 0.006, 0.01, 0.02, 0.15, 0.25, 0.46, 0.68, 0.79 )
>databo <- p.adjust(data, method="bonferroni")
>databo
[1] 0.01 0.03 0.06 0.10 0.20 1.00 1.00 1.00 1.00 1.00
> databh <- p.adjust(data, method = "BH")
>databh
[1] 0.0100000 0.0150000 0.0200000 0.0250000 0.0400000 0.2500000 0.3571429 0.5750000 0.7555556 0.7900000