我正在尝试使用R中的phyper计算超几何概率,并注意到一种奇怪的行为。我正在研究基因集重叠概率,在一种情况下,没有“成功抽奖”,所以:
x=0
m=430
n=19500
k=2
因为我一般都在寻找过度浓缩,所以我使用:
phyper(0,600,19000,2,lower.tail=FALSE)
得到0.043,这看起来很重要。 然而,选择两个“n”基因应具有
的概率19500/19930*19499/19929=0.957
所以,不应该将phyper结果大于0.957?