Alexandre,

Tenho em arquivo um conjunto de scripts (ver data) de um Curso realizado na Embrapa Arroz e Feijão: http://leg.ufpr.br/~walmes/cursoR/cnpaf2/ Eles podem servir de ponto de partida.

Esse curso http://leg.ufpr.br/~walmes/cursoR/geneticaEsalq/ considerou modelos mistos a partir do script15. Pode ser útil também.

Lembrando que devido a atualização de pacotes ou do R, talvez não seja reprozível na integra.

À disposição.

Walmes.