pdf("rpd11_2r.pdf") fitness <- read.fwf(file="FITNESS_DTA.txt",width=c(3,9,10,9,8,9,9), col.names=c("age","weight","o2","ex_time","hr_rest","hr_run","hr_max")) attach(fitness) reg11_2.reg <- lm(o2 ~ ex_time + hr_rest + hr_run + hr_max) reg11_2.rstandard <- rstandard(reg11_2.reg) reg11_2.yhat <- predict(reg11_2.reg) plot(reg11_2.yhat,reg11_2.rstandard) abline(h=0) dev.off()