  In the context of estimating stochastically ordered distribution functions,
the pool-adjacent-violators algorithm (PAVA) can be modified such that the
computation times are reduced substantially. This is achieved by studying the
dependence of antitonic weighted least squares fits on the response vector to
be approximated.
