[R-br] Diferença entre duas abordagens de comparação usando GLM

[Walmes Zeviani]

Na minha opinião, o procedimento 2 carrega muito da idéia do
Skott-Knott. Ambos são válidos pois o 1 é um teste de comparação
múltipla de parâmetros e o segundo é um teste de razão de
verossimilhanças entre modelos aninhados. Fique com o que for mais
confortável mas lembre-se que em 1 você compara parâmetros e em 2 você
compara modelos. Um estudo sobre os níveis de cobertura dessas
abordagens seria um coisa interessante.

À disposição.
Walmes.

On 6/17/12, ASANTOS <alexandresantosbr em yahoo.com.br> wrote:
> Boa noite pessoal,
>
>        Gostaria de saber se do ponto de vista da interpretação
> estatística do método empregado, existe diferença entre comparar 3
> tratamento hipotéticos usando:
>
> 1)
> require(multcomp)
> y<- c(mapply(rpois, lambda=c(5,20,5), MoreArgs=list(n=100)))##Criação da
> variável resposta
> trat <- gl(3,100) ##Criação dos tratamentos
> m.1 <- glm(y~trat,family="poisson")
> anova.m.1 <- anova(m.1, test="Chi")
> anova.m.1
> summary(glht(m.1, linfct=mcp(trat="Tukey")))
> #Os tratamentos 1 e 3 são semelhantes e diferentes do tratamentos 2
>
> ou
>
> 2)
> m.1 <-glm(y~trat,family="poisson")
> m.0 <- glm(y~1,family="poisson")
> anova(m.1,m.0,test="Chi")####Existe diferença
> sort(tapply(y,trat,mean))## Ordenando as médias
> levels(trat)
> ##Modelo T1 e T3
> t1t3<-trat
> levels(t1t3)
> levels(t1t3)[1]<-"1+3"
> levels(t1t3)[3]<-"1+3"
> levels(t1t3)
> m.2<- glm(y~t1t3,family="poisson")
> anova(m.1,m.2,test="Chi")####Semelhantes
> ## Mas se faço qualquer uma das outras combinações "1+2" ou "2+3" os
> valores do teste Chi são significativos, podendo-se concluir que os
> tratamentos 1 e 3 são semelhantes e diferentes do tratamentos 2.
>
> Obrigado,
>
>
> --
> Alexandre dos Santos
> Engenheiro Florestal, Dr.
> Universidade Federal de Lavras
> Departamento de Entomologia
> Laboratório de Entomologia Florestal
> Caixa Postal 3037
> 37200-000 - Lavras/MG
> Fone: +55 (35) 9223-0304
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
> mínimo reproduzível.
>


-- 
==========================================================================
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
==========================================================================