Edson,

Não entendi suas dúvidas. Você pode usar lm() quando todos os fatores são de efeito fixo, pode usar a aov() declarando fator de efeito aleatório com o termo Error() na formula do modelo (estimação por mínimos quadrados), e pode usar a nlme::lme() declarando fatores de efeito aleatório na opção random=, e pode usar a lme4::lmer(). Veja a página de Ridículas do LEG para ver o uso da lme() no caso de blocos incompletos www.leg.ufpr.br/ridiculas.

À disposição.
Walmes.

==========================================================================
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
VoIP: (3361 3600) 1053 1173
e-mail: walmes@ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================