# KOMMANDOER TIL FORELESNINGENE FOR UKE 5 

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

# EKSEMPEL: C VITAMIN (toveis variansanalyse)

# Les inn dataene:
appelsin.juice=read.table("http://www.uio.no/studier/emner/matnat/math/STK2120/v16/appelsin-juice.txt", header=T) 

# Plott av dataene:
attach(appelsin.juice)
interaction.plot(dager,merke,cvit)
detach(appelsin.juice)

# Toveis variananalyse (jf. side 588 i D&B):
fit.appelsin=aov(cvit~factor(merke)+factor(dager),data=appelsin.juice)
summary(fit.appelsin)

# Normalfordelingsplott (jf. side 588 i D&B):
qqnorm(fit.appelsin$residuals)

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

# Simultane konfidensintervall for dager [jf. side 589 i D&B]

# For dager:
tukey.dager=TukeyHSD(fit.appelsin,"factor(dager)")
print(tukey.dager)
plot(tukey.dager)



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


# EKSEMPEL: ENERGIFORBRUK VED ULIKE FYSISKE AKTIVITETER 
# (randomisert blokkforsoek)

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

# Plott av dataene:
attach(energi)
interaction.plot(aktivitet,person,kcal)
detach(energi)

# Variansanalyse:
fit.energi=aov(kcal~factor(aktivitet)+factor(person),data=energi)
summary(fit.energi)

# Normalfordelingsplott (jf. side 588 i D&B):
qqnorm(fit.energi$residuals)

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

# Simultane konfidensintervall for aktivitet [jf. side 589 i D&B]
tukey.aktivitet=TukeyHSD(fit.energi,"factor(aktivitet)")
print(tukey.aktivitet)
plot(tukey.aktivitet)

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


# EKSEMPEL 11.16 (toveis variansanalyse med interaksjon)

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

# Toveis variananalyse (jf. side 600 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 602 i D&B)
attach(exmp11.16)
interaction.plot(Density,Variety,Yield)
detach(exmp11.16)

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

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

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

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