hyper2binomP(3,4,5,6) # 0.38856
## The function is simply defined as
function (x, m, n, k) {
N <- m + n
p <- m/N
N.n <- N - (k - 1)/2
(m - x/2)/N.n - k * (x - k * p - 1/2)/(6 * N.n^2)
}
Run the code above in your browser using DataLab