Uma alternativa seria realizar a análise com ambiente/bloco... e depois colocar da mesma forma ambientes/genótipos e desdobrar a interação, saindo uma entre locais e a dentro de genótipos dentro de cada local! Mas acho que essa minha parametrização só funciona com o aov...! não lme não sei! abraço, FH 2011/5/28 Fabio Mathias Corrêa <fabio.ufla@yahoo.com.br>
Emmanuel,
A variância residual para cada ambiente, vc terá somente se realizar uma análise para cada ambiente separadamente.
No seu caso, o usual a se fazer é:
Realizar uma análise para cada ambiente, depois realizar a análise conjunta dos ambientes.
Desta forma forma vc tem a variância ambiental total e a variância de cada ambiente.
Tem outra coisa, o seu modelo esta errado.
O certo seria:
Ambiente/Bloco = Bloco dentro de Ambiente. (Wilkinson & Rogers, 1973)
A notação Bloco/Ambiente, é utilizada pelo Satanic Analysis System.
Vc está dizendo ao R que realize a análise de Ambientes dentro de blocos, o que não tem sentido.
Valeu!!!
Fábio Mathias Corrêa
Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e da Terra - DCET
Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna CEP 45662-900. Ilhéus-Bahia
Tel.: 73-3680-5076 ------------------------------ *De:* Emmanuel Arnhold <emmanuelarnhold@yahoo.com.br> *Para:* R-br <r-br@listas.c3sl.ufpr.br> *Enviadas:* Sábado, 28 de Maio de 2011 9:28 *Assunto:* [R-br] Ajuda em modelo misto
Ajustei o seguinte modelo misto:
ga=interaction(genótipo,ambiente)
modelo=lmer(resposta~bloco/ambiente+ambiente+(1|genótipo)+(1|ga))
É um modelo de avaliação de genótipos em vários ambientes. Neste modelo tenho as seguintes variâncias: ga genótipo residual
Gostaria que alguém me ajuda-se a ajustar o modelo acima para obter as seguintes variâncias, supondo 3 ambientes: ga genótipo residual (ambiente 1) residual(ambiente 2) residual (ambiente 3)
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br