Fernando, Acabei de conferir meus resultados com os do Ramalho et al, página 171. Eles não são exatamente iguais apesar de eu estar seguro sobre o modelo. Isso pode ser devido a tabela do livro estar com valores arredondados e os autores terem usado número com mais casas, algum erro de digitação na tabela ou magia negra do $@$ que foi o aplicativo considerado pelos autores para obter os resultados. Vai saber o que a caixa preta está fazendo. Segue a rotina que rodei. da <- read.table("http://www.leg.ufpr.br/~walmes/data/ramalho_milho.txt", header=TRUE, sep="\t", colClasses=rep(c("factor","numeric"),c(3,1))) str(da) m0 <- lm(terms(producao~repeticao/bloco+cultivar, keep.order=TRUE), data=da) anova(m0) require(doBy) popMeans(m0, effect="cultivar") À 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 ==========================================================================