search for: residsd

...dard deviation from mean std = sd(residuals(stepmod)) mn = mean(residuals(stepmod)) sdclas = seq(-3*std+mn, 3*std+mn, by=std) #I would like to have colors from green to red, but I can't do it so I am doing only palette( rainbow(6) ) #I am dividing residuals into earlier defined classes by residsd = cut(residuals(stepmod), sdclas, labels=FALSE) #or by (THEN POINTS DON'T APPEAR ON PLOT IN NEXT STEP) residsd = cut(residuals(stepmod), sdclas) # I am creating a plot plot(expzad$X,zadrz$Y,col=residsd, pch=15 ) # I would like to add legend to plot showing values for colours from plot but...

I have from time to time seen inquiries on r-help in respect of how to obtain the estimated standard deviation from the output of fitting a linear model. And have had occasion to want to do this myself. The way I currently do it is something like summary(fit)$sigma (where fit is returned by lm()). It strikes me that it might be a good idea to have an extractor function, analogous with coef()