Discriminatory power

[aalfons]

Check the function we use for computing the eta-squared (method = "eta2" in function discriminatory_crit() and possibly elsewhere) does what we think it does and what we describe in the documentation.