############### # Interactives Homework #7 library(car) attach(Prestige) mod1 <- lm(prestige ~ income + type + income*type) summary(mod1) plot(income, prestige, ylim=c(0,110), xlab="Income (Canadian $)", pch=as.numeric(type), col=as.numeric(type), main="Occupational Prestige by Type and Income Level") legend(18500,20, legend=c("Blue Collar", "Professional", "White Collar"), pch=c(1:3), col=c(1:3), lty=c(1:3)) abline(mod1$coef[1],mod1$coef[2], lty=1, col=1) abline((mod1$coef[1] + mod1$coef[3]),(mod1$coef[2] +mod1$coef[5]), lty=2, col=2) abline((mod1$coef[1] + mod1$coef[4]),(mod1$coef[2] +mod1$coef[6]), lty=3, col=3)