################################################### ### chunk number 4: LinearModelDiagnosticPlots ################################################### oldpar <- par(no.readonly = TRUE) par(mfcol = c(2, 2)) plot(fit, ask = FALSE) par(oldpar)