Gilson, São poucos NA, e pelo que vi nos dados não compromete a estimação dos efeitos. Se interação for significativa, dedobre. Se não for, despreze e trabalhe os efeitos principais. Você pode usar a contrast::contrast() e a multcomp::glht(). Fazer imputação dos dados para evitar os NA é um procedimento que desaconselho. (1), substituir pelo valor médio das outras repetições disponíveis vai gerar um desvio de zero, o que vai aumentar os graus liberdade e não aumentar a variância, ou seja, aumentando a precisão sem que exista. Na minha opinião, visto que há recursos disponíveis para o tratamentos dos NA, deve trabalhar apenas com os dados observados, sem imputação. Para uso da contrast() e glht() olha na ridículas do LEG. À 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 ==========================================================================