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

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

R简单实现

  1. p_value # a vector
  2. p_adjust <- p.adjust(p_value,method = "BH")