overrepresentationを計算するときに、p(X>=x)を計算したい。
しかし、Rのphyper
ではp(X<=x)
を計算するため、工夫が必要。
x # white balls drawn
m # nunmber of white balls in bag
n # number of black balls in bag
k # number of balls drawn
p <- 1 - phyper(x-1, m, n, k) # P[X >= x] = 1 - P[X <= x-1]