Um CMR na primeira mensagem teria resolvido a questão. Veja a comparação entre Anova(, type="III") e drop1() que foi minha sugestão. Verifique se "batem". da <- expand.grid(bloco=gl(4,1), A=gl(3,1), B=gl(5,1)) da$y <- rnorm(da$A) set.seed(123) da$y[sample(1:nrow(da), 5)] <- NA da <- na.omit(da) xtabs(~A+B, da) # observações perdidas m0 <- lm(y~bloco+A*B, da) anova(m0) car::Anova(m0, type="II") car::Anova(m0, type="III") drop1(m0, test="F", scope=.~.) # idem ao type="III" db <- expand.grid(bloco=gl(4,1), A=gl(3,1), B=gl(5,1)) db$y <- rnorm(db$A) db <- subset(db, A!="1" | B!="1") db <- na.omit(db) xtabs(~A+B, db) # casela pedida m1 <- lm(y~bloco+A*B, db) anova(m1) car::Anova(m1, type="II") car::Anova(m1, type="III") # erro summary(m1) # tem NA drop1(m1, test="F", scope=.~.) # sem problemas À 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 skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================