[R-br] Lattice Design

[Walmes Zeviani]

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 em ufpr.br
skype: walmeszeviani
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121116/36d3dd69/attachment.html>