# 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/v13/appelsin-juice.txt", header=T) # Plott av dataene: plot(c(0,3,7),appelsin.juice$cvit[appelsin.juice$merke==1],ylim=c(40,55),xlab="Dager",ylab="C vitamin",type='l') lines(c(0,3,7),appelsin.juice$cvit[appelsin.juice$merke==2],lty=2) lines(c(0,3,7),appelsin.juice$cvit[appelsin.juice$merke==3],lty=3) # Toveis variananalyse (jf. side 575 i D&B): fit.appelsin=aov(cvit~factor(merke)+factor(dager),data=appelsin.juice) summary(fit.appelsin) # Normalfordelingsplott (jf. side 576 i D&B): qqnorm(fit.appelsin$residuals) # Plott av residualer mot tilpassede verdier (jf. side 576 i D&B): plot(fit.appelsin$fitted.values,fit.appelsin$residuals) # ============================= # EKSEMPEL: ENERGIFORBRUK VED ULIKE FYSISKE AKTIVITETER (randomisert blokkfors?k) # Les inn dataene: energi=read.table("http://www.uio.no/studier/emner/matnat/math/STK2120/v13/energi.txt", header=T) # Variansanalyse: fit.energi=aov(kcal~factor(aktivitet)+factor(person),data=energi) summary(fit.energi)