# KOMMANDOER TIL FORELESNINGENE FOR UKE 3 

# EKSEMPEL: OPPTAK AV FETTSYRER
# se under "Data" paa kurssiden for beskrivelse av dataene

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

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

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

# Plott av residualer mot tilpassede verdier:
plot(aov.fit$fit,aov.fit$residuals)

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

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