Obrigado mais uma vez pela informação. Em 1 de junho de 2011 00:52, Walmes Zeviani <walmeszeviani@gmail.com>escreveu:
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@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- José Eduardo Vargas Lopes de Aráujo Engenheiro Agrônomo Contato:(38)8811-6985 duvargas@gmail.com