if np.all(np.abs(diff) <= diffthreshold) or phat.pvalue > 0.05: