[R-br] Análise Conjunta

[ana paula coelho madeira]

Prezados, boa tarde.
 
Tenho um experimento em DBC (4 blocos) instalado em dois locais distintos. Usei a função
 
mi <- lapply(split(dados, f=dados$Local), aov, formula=Y8~bloco+trat)
lapply(mi, summary)
 
para fazer as análises individuais e a função
 
m8 <- lm(terms(Y8~Local/bloco+Local*trat, keep.order=TRUE), data=dados)
anova(m8) 

para a análise conjunta. O fato é que em ambas as análises o grau de liberdade para blocos e blocos dentro de locais está incorreto.
 
Observem as saídas:
## Análises individuais
$NP
            Df Sum Sq Mean Sq F value  Pr(>F)   
bloco        1   0.51  0.5062   0.729 0.39696   
trat        17  35.64  2.0966   3.020 0.00108 **
Residuals   53  36.79  0.6942                   
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
$SL
            Df Sum Sq Mean Sq F value   Pr(>F)    
bloco        1    2.3    2.34   0.167 0.684317    
trat        17  741.3   43.61   3.120 0.000785 ***
Residuals   53  740.7   13.97                     
---

## Análise Conjunta
 
Response: Y8
             Df Sum Sq Mean Sq  F value    Pr(>F)    
Local         1 3872.1  3872.1 527.9346 < 2.2e-16 ***
Local:bloco   2    2.8     1.4   0.1937 0.8241648    
trat         17  380.3    22.4   3.0497 0.0002392 ***
Local:trat   17  396.7    23.3   3.1814 0.0001381 ***
Residuals   106  777.5     7.3     
 
Alguém consegue visualizar o erro? Sugerem outras funções?
 
Desde já agradeço
 
Att.,
 
Ana Paula
 
 
 		 	   		  
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140709/9da474ed/attachment.html>