rfsoneway(formula, data, location=mean, spread=function(x) sqrt(var(x)))y ~ x where y is the
    numeric response and x is the grouping factormedian, 0 (to
    avoid centering).sd, 1 (to
    avoid scaling).y - fitted.values).spread for their
    grouprfs, Lattice