q值（FDR）的理解和计算，并用R实现

为什么要计算 q 值或者FDR？

https://www.jianshu.com/p/d86823ecd3ac
上面文章里挺详细。在进行多重假设检验的时候，例如生信中最常见的差异基因分析的时候，每个基因进行假设检验，也要进行2万多次。我们通常会设置信区间为0.95，也就是会存在0.05的假阳性，如果筛选出p<0.001的基因1000个，会存在50个是假阳性。所以要对p值进行矫正。

R简单实现

p_value # a vector
p_adjust <- p.adjust(p_value,method = "BH")
©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

推荐阅读更多精彩内容