Here's the code for these graphs
par(mfrow=c(1,3)) plot(lm.all.income.factors, which=1) qqPlot(lm.all.income.factors$residuals, id=FALSE, main="QQ Plot") plot(lm.all.income.factors$residuals, main="Residuals vs Order") Also I'm not quite sure how to fix the QQ Plot: it looks like the data is skewed. Perhaps a non-linear graph would be more appropriate?
