[R-br] aov() e anova() no agricolae

[Alexandre Santos]

Bom Dia pessoal, 


      Estou realizando alguns testes de Tukey com ajuda do pacote agricolae para algumas variáveis que possuem distribuição normal, então fiz para uma destas:
#  
##N.ovulos
anova.o<-aov(n.ovulos~factor(tratamento)+factor(geracao)+factor(rep),data=dados6)
summary(anova.o)
#
##Tratamento 1 casal nas geracoes
with(subset(dados6,tratamento == "1casal"),HSD.test(n.ovulos,geracao, DFerror=df.residual(anova.o),MSerro=deviance(anova.o)/df.residual(anova.o), group = TRUE))
##Tratamento 10 casais nas geracoes
with(subset(dados6,tratamento == "10casal"),HSD.test(n.ovulos,geracao, DFerror=df.residual(anova.o),MSerro=deviance(anova.o)/df.residual(anova.o), group = TRUE))

    Ate ai tudo bem, porém ... tenho algumas outras variáveis explicativas dentro do mesmo experimento como razão sexual por exemplo, que é uma proporção (número de fêmeas/número de fêmeas+número de machos) que não tem distribuição normal e fiz o ajuste por GLM com distribuição de erros binomial. E agora gostaria de saber se não existe alguma forma de se utilizar HSD.test () do agricolae utilizando as informações da  anova() do meu modelo binomial de maneira análoga ao aov() que utilizo para os dados com distribuição normal ou vou ter que tratar os dados com outras distribuições de erros com o glht() do multcomp? Ou deveria fazer toda a abordagem do problema pegando as variáveis que têm normalidade e utilizar distribuição gaussiana via glm() e tratar todo o experimento através de contrastes de modelo via multcomp(), visando a uniformizar a abordagem estatística dada ao mesmo experimento?

Obrigado,




 
Alexandre dos Santos
Engenheiro Florestal, MSc.
Universidade Federal de Lavras 
Departamento de Entomologia 
Laboratório de Entomologia Florestal 
Caixa Postal 3037 
37200-000 - Lavras/MG
Fone: +55 (35) 3829-5122
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110824/e915db9a/attachment.html>