script to calculate poisson p val