19 Experimentos em Parcelas Subdivididas

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

Paulo Justiniano Ribeiro Jr