There is a problem with the "Constant Leverage" plot. Try the following: y <- c(8,9,7,9,6,9,5,9) grp <- factor(rep(1:4,each=2),labels=c('A','B','C','D')) mod <- lm(y ~ grp) windows() par(mfrow=c(2,2)) plot(mod) # Plot 5 ("Constant Leverage") has X axis labels sorted # by the fitted value (as advertised), but the data are # still plotted in the original order. # This occurs in R 2.4.1 and 2.5.0 RC for windows # but not in earlier releases