<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Walmes, muito obrigado pela resposta.<br>
<br>
Eu sou estatístico e sabia que o REML não trabalha com o teste F
pela partição de soma de quadrados. Eu acho que me expressei mal.
Queria mostrar somente o que queria fazer e relendo a mensagem,
realmente parece que quero fazer partição de SQ... bobeira minha e
desde já, desculpe. Era para ser o exemplo.<br>
<br>
Outra coisa, eu tentei reproduzir o código do meus dados ... não
serviu? Para mim, isso era um CMR. Reproduzi a estrutura dos meus
dados e tentei mostrar o que fiz e o exemplo da forma como seriam as
partições dos testes.<br>
Quanto às dicas, vou tentar anova com modelos aninhados... não tinha
pensado nisso e é relativamente fácil. Usei o teste anova (mod1,
mod2) para ver qual melhor estrutura de correlação funcionava e tal,
mas para testar os elementos fixos, não fiz. O pacote aod não
conheço, vou estudá-lo.<br>
<br>
Peço desculpas novamente por não me fazer entender e quando você
tiver um tempo, escreve um tópico sobre dados com medidas
repetidas... ia ser uma bênção, isso!<br>
<br>
Abs e obrigado.<br>
<br>
Olympio<br>
<br>
<br>
Em 25/03/2012 17:37, Walmes Zeviani escreveu:
<blockquote
cite="mid:CAFU=Ekb2C-buAEfEu=DNON1+kfSj1emhx8Ph_bvC8Rm-sQnABg@mail.gmail.com"
type="cite"><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
moz-do-not-send="true" 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
moz-do-not-send="true" 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>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
</body>
</html>