<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>Prezados, boa tarde.<BR> <BR>Tenho um experimento em DBC (4 blocos) instalado em dois locais distintos. Usei a função<BR> <BR>mi <- lapply(split(dados, f=dados$Local), aov, formula=Y8~bloco+trat)<br>lapply(mi, summary)<BR> <BR>para fazer as análises individuais e a função<BR> <BR>m8 <- lm(terms(Y8~Local/bloco+Local*trat, keep.order=TRUE), data=dados)<br>anova(m8) <br><BR>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.<BR> <BR>Observem as saídas:<BR>## Análises individuais<BR>$NP<br> Df Sum Sq Mean Sq F value Pr(>F) <br>bloco 1 0.51 0.5062 0.729 0.39696 <br>trat 17 35.64 2.0966 3.020 0.00108 **<br>Residuals 53 36.79 0.6942 <br>---<br>Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1<BR>$SL<br> Df Sum Sq Mean Sq F value Pr(>F) <br>bloco 1 2.3 2.34 0.167 0.684317 <br>trat 17 741.3 43.61 3.120 0.000785 ***<br>Residuals 53 740.7 13.97 <br>---<br><BR>## Análise Conjunta<BR> <BR>Response: Y8<br> Df Sum Sq Mean Sq F value Pr(>F) <br>Local 1 3872.1 3872.1 527.9346 < 2.2e-16 ***<br>Local:bloco 2 2.8 1.4 0.1937 0.8241648 <br>trat 17 380.3 22.4 3.0497 0.0002392 ***<br>Local:trat 17 396.7 23.3 3.1814 0.0001381 ***<br>Residuals 106 777.5 7.3 <BR> <BR>Alguém consegue visualizar o erro? Sugerem outras funções?<BR> <BR>Desde já agradeço<BR> <BR>Att.,<BR> <BR>Ana Paula<BR> <BR> <BR> </div></body>
</html>