dropModel.R
dropModel <- function(model, drop) {
model2 <- update(model,
dropFormula(model, drop))
plot(resid(model), resid(model2),
xlab = "Original Residuals",
ylab = paste("Residuals after dropping", drop))
abline(0,1)
model2
}