1.需要校正p_value的原因
我粗浅的理解就是,因为样本量过大,即使是很小的一个p值(即出错的可能),也会有很多的假阳性出现,所以需要通过矫正来减小假阳性率。我就简单介绍其中一种:多重检验校正方法(BH法),在R里边,有一个自带的包stats中的p.adjust函数。
2.R中的实现
例子:
> # demo
> p <- c(0.03,0.051,0.016,0.123)
> p_adj <- round(p.adjust(p,"BH"),3)
> p_adj
[1] 0.060 0.068 0.060 0.123
计算方法:
3.参考资料
1.多重比较P值校正
2.p值校正小汇总(知乎)
里头有更详细的一些解释~~
问题:
样本量多大就需要矫正呢,根据什么来确定样本量是算大还是算小呢?争取日后搞懂呜呜呜。