Vamos mostrar aqui como especificar o modelo para análise de experimentos em parcelas subdivididas. Os comandos abaixo mostram a leitura e preparação dos dados e a obtenção da análise de variância. Deixamos por conta do leitor a análise exploratória e de resíduos, desdobramento das interações e testes de comparações múltiplas.
Considere o experimento em parcelas subdivididas de dados de produção de aveia descrito na apostila do curso. Clique aqui para ver e copiar o arquivo com conjunto de dados. A obtenção da análise de variância é ilustrada nos comandos e saídas abaixo.
> ex09 <- read.table("exemplo09.txt", header=T)
> ex09
> dim(ex09)
[1] 64 4
> names(ex09)
[1] "a" "b" "bloco" "resp"
>
> ex09$a <- as.factor(ex09$a)
> ex09$b <- as.factor(ex09$b)
> ex09$bloco <- as.factor(ex09$bloco)
>
> summary(ex09)
a b bloco resp
1:16 1:16 1:16 Min. :28.30
2:16 2:16 2:16 1st Qu.:44.90
3:16 3:16 3:16 Median :52.30
4:16 4:16 4:16 Mean :52.81
3rd Qu.:62.38
Max. :75.40
>
> ex09.av <- aov(resp ~ bloco + a*b + Error(bloco/a), data=ex09)
> summary(ex09.av)
Error: bloco
Df Sum Sq Mean Sq
bloco 3 2842.87 947.62
Error: bloco:a
Df Sum Sq Mean Sq F value Pr(>F)
a 3 2848.02 949.34 13.819 0.001022 **
Residuals 9 618.29 68.70
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1
Error: Within
Df Sum Sq Mean Sq F value Pr(>F)
b 3 170.54 56.85 2.7987 0.053859 .
a:b 9 586.47 65.16 3.2082 0.005945 **
Residuals 36 731.20 20.31
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1