# KOMMANDOER TIL FORELESNINGENE FOR UKE 6 

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

# Lineaer regresjon

# Eksempel: Karakterer for computer science studenter

# Les inn data:

karakterer=read.table("http://www.uio.no/studier/emner/matnat/math/STK2120/v13/karakterer.txt", header=T)

# Tilpasser lineaer regresjonsmodell:
fit.kar=lm(kar~hsm+hse+satm+satv,data=karakterer)
summary(fit.kar)

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


# Fra avnitt 3 i "Variansanalyse og lineaer regresjon"

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

# Enveis variansanalyse:
fit.aov=aov(sleep~factor(treat),data=exmp11.5)
summary(fit.aov)

# Vil bruke restriksjonen at alfa-ene summerer til null:
options(contrasts=c("contr.sum","contr.poly"))

# Tilpass lineaer regresjonsmodell:
fit.lm=lm(sleep~factor(treat),data=exmp11.5)
summary(fit.lm)

# Variansanalysetabellen (fra tilpasset regresjonsmodell):
anova(fit.lm)

# Designmatrisen (dvs X-matrisen naar modellen er skrevet som
#  en multippel lineaer regresjonsmodell paa matriseform):
model.matrix(fit.lm) 


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

# Fra avnitt 5 i "Variansanalyse og lineaer regresjon"

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

# Tilpass lineaer regresjonsmodell:
fit2.lm=lm(Response~factor(Brand)+factor(Treat), 
           data=exmp11.11)
summary(fit2.lm)

# Variansanalysetabellen:
anova(fit2.lm)

# Designmatrisen:
model.matrix(fit2.lm) 

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

# Fra avnitt 7 i "Variansanalyse og lineaer regresjon"

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

# Tilpass lineaer regresjonsmodell:
fit3.lm=lm(Yield~factor(Variety)+factor(Density)+ 
        factor(Variety):factor(Density), data=exmp11.16)
summary(fit3.lm)

# Variansanalysetabellen:
anova(fit3.lm)

# Designmatrisen:
model.matrix(fit3.lm) 

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

# Fra avnitt 8 i "Variansanalyse og lineaer regresjon"

# Bruker naa restriksjonen at alfa1=0 (default):
options(contrasts=c("contr.treatment","contr.poly"))

# Ser paa nytt paa eksempelet i avsnitt 3:
fit.lm=lm(sleep~factor(treat),data=exmp11.5)
summary(fit.lm)
anova(fit.lm)


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

# Fra avnitt 9 i "Variansanalyse og lineaer regresjon"

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

# Variansanalyse der kjonn skrives opp foerst i modellen:
 fit4.lm=lm(vekst~factor(kjonn)+factor(hemming)+
         factor(kjonn):factor(hemming), data=vekst)
anova(fit4.lm)

# Variansanalyse der hemming skrives opp foerst i modellen:
fit4b.lm=lm(vekst~factor(hemming)+factor(kjonn)+
         factor(hemming):factor(kjonn),data=vekst)
anova(fit4b.lm)