# les inn data fra eksempel 11.11 i boken
d = read.table("http://www.uio.no/studier/emner/matnat/math/STK2120/v14/exmp11-11.txt",header=TRUE,sep=",")

# sett som faktorer
d$Treat = as.factor(d$Treat)
d$Brand = as.factor(d$Brand)

# plott linjer
plot(as.vector(d$Treat),d$Response)
for(j in 1:3)
  lines(1:4,d$Response[d$Brand==j])
        
# toveis anova
fit = aov(Response~Brand+Treat,data=d)
summary(fit)

# sjekk normalfordeling 
qqnorm(fit$res);qqline(fit$res)
plot(fit)

#multippel sammenligning for vaskemidlene
plot(TukeyHSD(fit,which="Treat"))