[R-br] Duvida no summary da ancova, usando GLM

[Walmes Zeviani]

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 em ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================


2011/9/9 <frdor85 em 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???
>
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110909/3791fcb0/attachment.html>