library(car)
modp1 <- lm(prestige ~ income + education + women, data=Prestige)
summary(modp1)

qq.plot(modp1, simulate=TRUE, line="none")

rstud <- rstudent(modp1)
plot(density(rstud))

boxplot(rstud)

influencePlot(modp1)

av.plots(modp1, ask=FALSE)

cr.plots(modp1, ask=FALSE)

modp2 <- lm(prestige ~ log10(income) + education + poly(women, 2), data=Prestige)
summary(modp2)
Anova(modp2)

influencePlot(modp2)

cr.plots(modp2, ask=FALSE)