R-hjelp til exercise 5 i BSS

 

 

# Les dataene inn i en dataramme og gi navn til variablene:

gun<-read.table("http://www.math.uio.no/avdc/kurs/STK4900/data/gun.dat", col.names=c("method","phys","team","rounds"))

 

# Gj?r variablene i datarammen tilgjengelige

attach(gun)

 

# PUNKT 1)

 

# Beregner korrelasjoner:

cor(gun)

 

# Hva er korrelasjonen mellom kovariatene ("method","phys","team")?

# Kan du forklare hvorfor de er slik?

# Hvordan er korrelasjonen mellom kovariatene og "rounds"?

 

 

# PUNKTENE 2-4)

# Definerer som faktorer (kategoriske kovariater):

method=factor(method)

phys=factor(phys)

team=factor(team)

 

# Vil bruke sum-kontrast (som er vanlig for anova):

options(contrasts=c("contr.sum","contr.poly"))

 

# Tillpasser modell med alle hovedeffekter og interaksjoner og skrive ut anova tabell:

gfit<-lm(rounds~method*phys*team)

anova(gfit)

 

# Hva sier anova tabellen deg? Hvilke faktorer har signifikant betydning?

 

# Se p? parameterestimatene:

summary(gfit)

 

# Gi en fortolkning av de viktigste parameterestimatene