Replikering och blockering

Experimentell design i R

Joanne Xiong

Data Scientist

Replikering

Ett experiment måste upprepas för att variabiliteten ska kunna bedömas korrekt
Om vi bara testar ett läkemedels effekt på en person, hur kan vi generalisera resultaten? (Det kan vi inte!)

library(dplyr)
mtcars %>%
    count(cyl)

    cyl  n
1   4     11
2   6      7
3   8     14

Blockering

Hjälper till att kontrollera variabilitet genom att göra behandlingsgrupperna mer likartade
Inom grupper är skillnaderna små – mellan grupper är de större
Ett exempel är att blockera behandlingsgrupper efter kön

Lådagram

# Boxplot of MPG by Car Cylinders
ggplot(mtcars, aes(x=as.factor(cyl), 
                   y=mpg)) + 
  geom_boxplot(fill="slateblue", 
               alpha=0.2) + 
  xlab("cyl")

alt text

Funktioner för modellering

Linjära modeller
```
lm(formula, data, na.action,...)
```
Envägs ANOVA-modell
```
aov(formula, data = NULL, ...)
```
Nästlad ANOVA-modell
```
anova(object,...)
```

Nu kör vi en övning!

Experimentell design i R