R语言-FDR校正的原理

为什么要进行FDR校正？

假设检验的原理

Excel计算方法

step1：把P值从大到小排序
step2：公式：p * (总数/p的位次)
【但是，得注意下一条的第一点，建议还是R直接一步到位】

为什么校正后FDR会出现相同值？

1.如果某一个p值所对应的FDR值大于前一位p值（排序的前一位）所对应的FDR值，则放弃公式计算出来的FDR值，选用与它前一位相同的值。

2.P值都小于1。

R语言的实现代码

p.adjust(p, method = p.adjust.methods, n = length(p))

#p.adjust.methods
# c("holm", "hochberg", "hommel", "bonferroni", "BH", "BY",
#   "fdr", "none")

#"BH"和"fdr"都是FDR校正
#"bonferroni"是最严格的校正方法

参考
http://www.biotrainee.com/thread-2759-1-1.html