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