Bears = read.table("c:/413/beardata.txt", header=T)
Bears$SEAS=ifelse(Bears$MONTH >= 4 & Bears$MONTH <= 9,1,0)
attach(Bears)
names(Bears)

# Does order matter??
bear1<-lm(WEIGHT~AGE+SEAS)
bear2<-lm(WEIGHT~SEAS+AGE)
summary(bear1)
summary(bear2)
anova(bear1)
anova(bear2)

bear3<-lm(WEIGHT~AGE+SEAS+CHEST+NECK)
summary(bear5)
anova(bear5)

# Making coplots to look for interactions

coplot(WEIGHT~CHEST|AGE)
coplot(WEIGHT~Chest|AGE,number=4,overlap=0)
coplot(WEIGHT~NECK|AGE,number=4,overlap=0)

detach(Bears)