<font face="trebuchet ms,sans-serif">Olympio,<br><br>Quando estimação é REML ou ML não existe esse personagem soma de quadrados. Os testes de hipótese para os termos de efeito fixo não são feitos por uma estatística F de uma análise de variância. Usamos a estatística F de um teste de Wald. Quando você pede anova() de um objeto nlme, você não recebe um quadro de análise de variância, mas um quadro de teste de Wald. Isso porque anova() é uma função genérica, ela olha para classe do objeto e então retorna uma saída apropriada. Dessa forma, lm e aov retornam anovas, glm retornam quadro de deviance, lme retorna teste de Wald, e assim vai. Então, não existe anova (SQ, QM, F...) para modelos aleatórios. Mesmo assim, acredito que seja possível fazer testes particionados semelhante ao que fazemos com summary(..., split=). Uma possível solução seria o ajuste de modelos aninhados/seqüenciais e a comparação do ajuste pelo anova(modelo1, modelo2). Outra possíbilide é usar a função aod::wald.test() para fazer esses testes de "particionamento de hipótese". Consulte um estatístico ou um livro de modelos lineares e veja a documentação do pacote. Um CMR é bem vindo.<br>
<br>À disposição.<br>Walmes.<br><br clear="all"></font><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Walmes Marques Zeviani</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Departamento de Estatística - Universidade Federal do Paraná</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">fone: (+55) 41 3361 3573</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">VoIP: (3361 3600) 1053 1173</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">e-mail: <a href="mailto:walmes@ufpr.br" target="_blank">walmes@ufpr.br</a></span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">twitter: @walmeszeviani</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">homepage: <a href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">linux user number: 531218</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br>