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