plot.lmrob {robustbase} | R Documentation |
Diagnostic plots for elements of class lmrob
## S3 method for class 'lmrob': plot(x, which = 1:5, caption = c("Standardized residuals vs. Robust Distances", "Normal Q-Q vs. Residuals", "Response vs. Fitted Values", "Residuals vs. Fitted Values" , "Sqrt of abs(Residuals) vs. Fitted Values"), panel = points, sub.caption = deparse(x$call), main = "", ask = prod(par("mfcol")) < length(which) && dev.interactive(), ..., p=0.025)
x |
an object as created by lmrob |
which |
integer number between 1 and 5 to specify which plot is desired |
caption |
Caption for the different plots |
panel |
Panel |
main |
Main title |
sub.caption |
sub titles |
ask |
waits for user input before displaying each plot |
... |
optional arguments for par |
p |
threshold for distance-distance plot |
Robust diagnostic plots as in Rousseuw and van Zomeren (1990)
lmrob
, also for examples, plot.lm
.
data(starsCYG) ## Plot simple data and fitted lines plot(starsCYG) lmST <- lm(log.light ~ log.Te, data = starsCYG) RlmST <- lmrob(log.light ~ log.Te, data = starsCYG) RlmST abline(lmST, col = "red") abline(RlmST, col = "blue") op <- par(mfrow = c(2,2), mgp = c(1.5, 0.6, 0), mar= .1+c(3,3,3,1)) plot(RlmST, which = c(1:2, 4:5)) par(op)