# KOMMANDOER TIL FORELESNINGENE FOR UKE 6 

# ===============


# EKSEMPEL 11.16 (toveis variansanalyse med interaksjon)

# Les inn dataene:
exmp11.16=read.table("http://www.uio.no/studier/emner/matnat/math/STK2120/v13/exmp11-16.txt", header=T,sep=",")


# Toveis variananalyse (jf. side 588 i D&B):
fit.exmp11.16=aov(Yield~factor(Variety)+factor(Density)+ factor(Variety):factor(Density),data=exmp11.16)
summary(fit.exmp11.16)

# Interaksjonsplott (jf. side 589 i D&B)
attach(exmp11.16)
interaction.plot(Density,Variety,Yield)
detach(exmp11.16)

# Normalfordelingsplott (jf. side 576 i D&B):
qqnorm(fit.exmp11.16$residuals)

# Plott av residualer mot tilpassede verdier (jf. side 576 i D&B):
plot(fit.exmp11.16$fitted.values, fit.exmp11.16$residuals)

# Multippel sammenligning av sorter med 99% simultan konfidenskoeffisient (jf. side 590 i D&B):
tukey.vareity=TukeyHSD(fit.exmp11.16,"factor(Variety)",ordered=T, conf.level=0.99)
print(tukey.vareity)
plot(tukey.vareity)

# Multippel sammenligning av tettheter 99% simultan konfidenskoeffisient (jf. side 590 i D&B):
tukey.density=TukeyHSD(fit.exmp11.16,"factor(Density)",ordered=T, conf.level=0.99)
print(tukey.density)
plot(tukey.density)