1.对p值进行从小到大的排序,标记上序号1~n; 2.其中,最大的FDR(不考虑重复则为第n位)等于最大的p值; 3.对于n-1位的FDR,取下面两者的较小值: 上一步(第n位)计算得出的FDR值; p-value*n/(n-1) 4.不断迭代第三步(n-2,n-3...),直至计算到最小p值对应的FDR。 下面直接在R里实现: ###例...
它的基本原理是将整体显著性水平平均分配给每个检验,从而降低单个检验的显著性水平。这意味着在Bonferroni检验中,为了拒绝原假设,单个检验的p值必须小于整体显著性水平除以检验数量。FDR校正方法则是基于FDR值的控制来调整p值的阈值。FDR校正相对于Bonferroni来说更为温和,它试图在假阳性和假阴性之间取得平衡,将假阳性的...
【统计】对 P 值进行 FDR 校正的原理(附代码演示) 做了多年的生物信息处理,和统计显著已经打过很多交道。然而 限于生物出身背景,对于统计知识了解不是很深刻。细细想来,每次 对外说出各种检验方式的名称(通常都是以人名命名的那种),自己 却鲜有细细了解到底这些检验的原理以及具体操作过程。近日听了一 场答辩,评...
2.其中,最大的FDR(不考虑重复则为第n位)等于最大的p值; 3.对于n-1位的FDR,取下面两者的较小值: 上一步(第n位)计算得出的FDR值; p-value*n/(n-1) 4.不断迭代第三步(n-2,n-3...),直至计算到最小p值对应的FDR。 下面直接在R里实现: ###例如这10个p值进行校正 temp <- c(0.01,0.11,0....
1.对p值进行从小到大的排序,标记上序号1~n; 2.其中,最大的FDR(不考虑重复则为第n位)等于最大的p值; 3.对于n-1位的FDR,取下面两者的较小值: 上一步(第n位)计算得出的FDR值; p-value*n/(n-1) 4.不断迭代第三步(n-2,n-3...),直至计算到最小p值对应的FDR。 下面...
2.其中,最大的FDR(不考虑重复则为第n位)等于最大的p值; 3.对于n-1位的FDR,取下面两者的较小值: 上一步(第n位)计算得出的FDR值; p-value*n/(n-1) 4.不断迭代第三步(n-2,n-3...),直至计算到最小p值对应的FDR。 下面直接在R里实现: ###例如这10个p值进行校正 temp 根据上述定义,我们知道...
3.对于n-1位的FDR,取下面两者的较小值: 上一步(第n位)计算得出的FDR值; p-value*n/(n-1) 4.不断迭代第三步(n-2,n-3...),直至计算到最小p值对应的FDR。 下面直接在R里实现: ###例如这10个p值进行校正 temp <->0.01,0.11,0.21,0.31,0.41,0.51,0.61,0.71,0.81,0.91)p.adjust(temp, method...