This script is written for implementation of NB permutation test with Stouffer's method.
In the script, 10,000 features were generated, and they were classified into two groups by ENPP method, then p-values were estimated with NB permutation test of K=1 for the large p-value group and K=3 for the small p-value group.
Detailed algorithm is described in the link below.