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)