# pdf("F:\\sta4210\\CH10TA01.pdf") life_ins <- read.table("http://www.stat.ufl.edu/~winner/sta4210/data/CH10TA01.txt", header=F, col.names=c("X1","X2","Y")) attach(life_ins) par(mfrow=c(1,2)) plot(X1,residuals(lm(Y~X1+X2))) # Figure 10.3a plot(residuals(lm(X1~X2)),residuals(lm(Y~X2))) # Figure 10.3b par(mfrow=c(1,1)) summary(lm(Y~X1+X2)) # Equation 10.3 summary(lm(Y~X2)) # Equation 10.4a summary(lm(X1~X2)) # Equation 10.4b # dev.off()