# KOMMANDOER TIL FORELESNINGENE FOR UKE 3 

# EKSEMPEL: 11.1 fra boken
# (se boken for beskrivelse av dataene)

# Leser inn dataene 
d = read.table("http://www.uio.no/studier/emner/matnat/math/STK2120/v14/exmp11-01.txt",header=TRUE,sep=",")

# Enveis variansanalyse (jf. side 548 i D&B):
d$box = as.factor(d$box)
fit = aov(strength~box,data=d)
summary(fit) 

# Normalfordelingsplott av residualene:
qqnorm(fit$res);qqline(fit$res) 

# Levenes test for lik varians (jf. side 549 i D&B). 
# Kommadoen bruker R-pakken "car"
# (som maa installeres hvis du bruker egen PC, install.packages("car"))
library(car)
fit2 = leveneTest(d$strength,d$box,center=mean)
print(fit2)
 
# Simultane konfidensintervall [jf. (11.4) paa side 553 i D&B]:
tukey.fit=TukeyHSD(fit)
print(tukey.fit)
plot(tukey.fit)
abline(v=0)