Caro amigo sem nome, O fato de ser glm não complica as coisas. O nível de referência segue ordem alfanumérica e não de entrada na planilha. A densidade está declarada como contínua e presente apenas com efeito de primeira ordem? Você vai ter que comparar as estimativas das médias para um valor fixo de densidade. Se usar o pacote contraste, isso é bem simples. Veja abaixo um exemplo simulado bem mediocre mas que é útil para entender o procedimento tr <- gl(3,4) # níveis de tratamento ds <- runif(length(tr), 1, 1.2) # densidade do solo y <- rbinom(length(tr), prob=0.5, size=50) # germinação # tr é categórica, ds é númerica m0 <- glm(cbind(y, 50-y)~tr*ds, family=binomial) anova(m0, test="Chisq") summary(m0) # compara os níveis de tratamento entre sí com ds=1.1 require(contrast) # tr 1 contra tr 2 sendo a densidade do solo igual a 1.1 contrast(m0, list(tr="1", ds=1.1), list(tr="2", ds=1.1)) # diferença em termos de preditor linear predict(m0, data.frame(tr="1", ds=1.1))- predict(m0, data.frame(tr="2", ds=1.1)) # diferença em termos de resposta predict(m0, data.frame(tr="1", ds=1.1), type="response")- predict(m0, data.frame(tr="2", ds=1.1), type="response") # fazer as demais comparações. usar método Bonferroni para # manutenção do erro tipo 1. Recomendo que você participe da nossa lista de discussão sobre R. Suas dúvidas serão bem recebidas e as respostas bem rápidas. Neste link estão informações para cadastro, envio e etiqueta na lista http://www.leg.ufpr.br/doku.php/software:rbr À 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 ========================================================================== 2011/9/9 <frdor85@hotmail.com>
Olá Walmes,
É a primeira vez que pergunto algo aqui e gostaria muito de resolver um dúvida nas minhas análises. Faço mestrado trabalhando com emergência de plântulas herbáceas em diferentes tratamentos de serrapilheira(Adição, controle e exclusão de serrapilheira). O meu modelo inclui o teor de argila no solo e a interaçãoentre argila e os 3 tratamentos. No summary do teste não aparece os resultados referentes ao tratamento de adição de serrapilheira. Me falaram que o R compara o primeiro nível do tratamento com os outros (na minha planilha original a adição de serrapilheira é o primeiro nível dos tratamentos). Eu inverti a ordem de entrada dos tratamentos e a comparação foi a mesma!!! Como não é possível usar HDS Tukey, como posso obter todos os resultados da comparação entre os tratamentos de uma ancova, usando GLM???