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: [email protected]
skype: walmeszeviani
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
