Add new argument ordered_trafo to trafo() and deal with ordered factors in a more transparent way (via a new function of_trafo()).