Como obter resultados da anava da função AOV() ou lm()

Uma situação muito comum no R, é quando precisamos chamar resultados da análise de variância da função aov() ou lm() do R. Um exemplo claro, quando estamos fazendo um teste de comparação múltipla após a ANAVA, e estamos necessitando do quadrado médio do resíduo, como chamá-lo, a partir das funções aov() ou lm()? Veremos abaixo.

####################################

Funções

####################################

> #ANAVA:
> anava <-lm(PROD~TRAT+BLOCO,data=dados)
> anava <-anova(anava);anava
 

Analysis of Variance Table

Response: PROD
          Df Sum Sq Mean Sq F value   Pr(>F)   
TRAT       3 4.3820 1.46068  8.4706 0.001572 **
BLOCO      5 2.7590 0.55179  3.1999 0.036559 * 
Residuals 15 2.5866 0.17244                    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

> names(anava)
[1] "Df"      "Sum Sq"  "Mean Sq" "F value" "Pr(>F)" 

A partir daqui, o objeto anava está pronto como uma lista. Dessa forma podemos obter os resultadosda análise de variância da seguinte forma: 

• Grau de liberdade:

   > anava$"Df"[1:3]
  [1]  3  5 15

• Soma de quadrado: 

  > anava$"Sum Sq"[1:3]
  [1] 4.382033 2.758950 2.586617  
  

• Quadrado médio: 

  > anava$"Mean Sq"[1:3]
  [1] 1.4606778 0.5517900 0.1724411
    

Caso queira só acessar o quadrado médio do resíduo:

> anava$"Mean Sq"[3]
[1] 0.1724411

Até a próxima!!!