Fisher-Pitman Permutation Tests

Computes the p-values of a Fisher-Pitman test of two independent samples or matched-pairs

These two functions, one for two independent samples and the other for matched-pairs, compute the one-sided and two-sided Fisher-Pitman permutation test p-values. Each function gives the option to use Monte Carlo simulations to compute the p-values and accepts an optional argument specifying the number of simulations to run.
The functions are quite fast for a reasonable number of simulations (<1m), but in it's current form, the code cannot handle a very large number of simulations because of the large matrices it produces.
The following articles describe how these tests are implemented in STATA:
http://www.stata-journal.com/sjpdf.html?articlenum=st0134

http://ageconsearch.umn.edu/bitstream/122702/2/sjart_st0158.pdf