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