If `sscipy.stats.ncx2` is usable for the test statistics the code would be much nice since no chi2 dist has to be implemented by hand