scale y for better convergence moved from rfit to jaeckel. drop.test produces an error if the statistic is negative