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
==========================================================================