source("functionsRB.R")
load("apneaData.rda")
attach(apneaData)
par(mfrow=c(3,2))
uniplot(data.frame(apneaScore , age , bmi , gender , partgasp))
par(mfrow=c(2,2))
prelim.plot(apneaScore ~ age + bmi + gender + partgasp)

lfit <- lm(apneaScore ~ age + bmi + gender + partgasp)
summary(lfit)
confint(lfit)
lfit <- lm(apneaScore ~ I(age/20) + I(bmi/10) + gender + partgasp)
summary(lfit)
confint(lfit)
lfit <- lm(apneaScore ~ gender*(age + bmi + partgasp))
summary(lfit)
confint(lfit)
anova(lm(apneaScore ~ gender + age + bmi + partgasp ),
     lm(apneaScore ~ gender*(age + bmi + partgasp)))
AIC(lm(apneaScore ~ gender + age + bmi + partgasp),
     lm(apneaScore ~ gender*(age + bmi + partgasp)))
library(effects)
lfit <- lm(apneaScore ~ gender*(age + bmi + partgasp))
plot(all.effects(lfit),1)
plot(all.effects(lfit),3)