### Eksempler Kap. 11

## Eksempel (opptak av fettsyrer)

# Leser dataene
fett = read.table("https://www.uio.no/studier/emner/matnat/math/STK1110/data/fettsyrer.txt",header=TRUE)

# Plott av dataene:
plot(fett$produkt,fett$opptak, ylim=c(0,300),xlab="Produkt",ylab="Opptak")

fett$produkt = as.factor(fett$product)
# Enveis variansanalyse (jf. tabell 11.2 paa side 559 i D&B):
aov.fit = aov(opptak~factor(produkt),data=fett)
summary(aov.fit)

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

# Bartletts test for lik varians (som essensielt er likelihood ratio-testen)
bartlett.test(opptak~factor(produkt),data=fett)

# Levenes test for lik varians (jf. side 562 i laereboka). 
# Kommadoen bruker R-pakken "car" 
library(car)
leveneTest(opptak~factor(produkt),data=fett,center=mean)
for(i in 1:5){show(sd(fett$opptak[fett$produkt==i]))}

#Eksempel 11.1
x = read.table("https://www.uio.no/studier/emner/matnat/math/STK1110/data/exmp11-1.txt",
               header=FALSE)
names(x) = c("Strength","Box")
plot(x$Box,x$Strength,xlab="Box",ylab="Strength")

aov.fit = aov(Strength~factor(Box),data=x)
summary(aov.fit)

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

# Bartletts test for lik varians (som essensielt er likelihood ratio-testen)
bartlett.test(Strength~factor(Box),data=x)

library(car)
leveneTest(Strength~factor(Box),data=x,center=mean)

for(i in 1:4){show(sd(x$Strength[x$Box==i]))}