pdf("muscle2.pdf")

muscle2 <- read.table("C:\\data\\muscle2.dat", header=F,
col.names=c("M","W","H"))

pairs(muscle2)

muscle2.reg <- lm(H ~ M + W, data=muscle2)

summary(muscle2.reg)

aov(muscle2.reg)

par(mfrow=c(2,2))

plot(muscle2.reg, which=1:4)

muscle2.rstandard <- rstandard(muscle2.reg)

muscle2.rstudent <- rstudent(muscle2.reg)

muscle2.inf <- influence.measures(muscle2.reg)

muscle2.rstandard

muscle2.rstudent

muscle2.inf


## You must have downloaded DAAG to get Variance Inflation Factors

library(DAAG) 

muscle2.vif <- vif(muscle2.reg)

muscle2.vif

dev.off()