lmRobMM
.
scale.lmRobMM(x, estimate="initial")
"lm.RobMM"
.
The default is to return the initial S-estimate of the scale. To obtain the
scale estimate corresponding to the final M-estimates, set
estimate
to
"final"
.
stack.df <- data.frame(Loss=stack.loss,stack.x) stack.robust <- lmRobMM(Loss~Air.Flow+Water.Temp+Acid.Conc.,data=stack.df) scale.lmRobMM(stack.robust)