Réplication et blocs

Plan d'expériences avec R

Joanne Xiong

Data Scientist

Réplication

Il faut répéter une expérience pour bien évaluer la variabilité
Si nous testons l'efficacité d'un médicament sur une seule personne, comment généraliser les résultats ? (Impossible !)

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

    cyl  n
1   4     11
2   6      7
3   8     14

Blocage

Réduit la variabilité en rendant les groupes de traitement plus semblables
À l'intérieur des groupes, les écarts sont minimes. Entre groupes, ils sont plus grands
Exemple : regrouper les traitements par sexe

Boîtes à moustaches

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

texte de remplacement

Fonctions de modélisation

Modèles linéaires
```
lm(formula, data, na.action,...)
```
Modèle d'ANOVA à un facteur
```
aov(formula, data = NULL, ...)
```
Modèle d'ANOVA emboîtée
```
anova(object,...)
```

Passons à la pratique !

Plan d'expériences avec R