# KOMMANDOER TIL FORELESNINGENE FOR UKE 3 

# EKSEMPEL: OPPTAK AV FETTSYRER
# (se under "Data" p? kurssiden for beskrivelse av dataene)

# Leser inn dataene 
fett=read.table("http://www.uio.no/studier/emner/matnat/math/STK2120/v13/fettsyrer.txt",header=T) 

# Plott av dataene:
plot(fett$produkt,fett$opptak,ylim=c(0,300))

# Enveis variansanalyse (jf. side 548 i D&B):
aov.fit=aov(opptak~factor(produkt),data=fett)
summary(aov.fit)

# Normalfordelingsplott av residualene:
qqnorm(aov.fit$residuals)

# Levenes test for lik varians (jf. side 549 i D&B). 
# Kommadoen bruker R-pakken "car"
# (som m? installeres hvis du bruker egen PC)
library(car)
leveneTest(opptak~factor(produkt),data=fett,center=mean)

# Simultane konfidensintervall [jf. (11.4) p? side 553 i D&B]:
tukey.fit=TukeyHSD(aov.fit,ordered=T)
print(tukey.fit)
plot(tukey.fit)