[R-br] teste de Homogeneidade e Normalidade

[Walmes Zeviani]

Eduardo,

O código abaixo mostra como não é difícil analisar os seus dados pela
abordagem que discuti

#------------------------------------------------------------------------------------------
# experimento com 11 tratamentos e 1 controle, resposta binomial n=25, 4
repetições

trat <- c("controle", paste("trat", 1:11, sep=""))
da <- expand.grid(trat=trat, rept=1:4)
da$germinadas <- rbinom(nrow(da), size=25, prob=0.8)
levels(da$trat) # controle é o nível de referência

# ajuste um modelo linear generalizado
g0 <- glm(cbind(germinadas, 25-germinadas)~trat, data=da, family=binomial)
anova(g0, test="Chisq") # testa hipótese sobre os efeitos principais
summary(g0) # cada efeito é o contraste controle vs trat_i
# usar correção de bonferroni nas comparações com controle

#------------------------------------------------------------------------------------------

Se todo mundo continuar fazendo as coisas como todos já fizeram não vamos
evoluir, certo?

À 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
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/20110601/5331e054/attachment.html>