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