Function to compute the Monte Carlo P-value for the observed Hollander-Proschan T statistic.
Useful links