> > +    # Load data to compare to
> > +    ar = np.load(os.path.join(cur_dir, 'res_gstat_test_MULM_OLS.npz'))

> It comes from the first implementation Benoit performed.

That doesn't mean I know I can trust it.

> What do you think of a test performed on Normal and unrelated data, with a few
> targets, many (>= 10,000) permutations, and verifying that h0 is symetric,
> centered at 0, or even that it is a Gaussian?

Sounds like a good thing to me.