[R-br] Falha em ANAVA

Fernando Antonio de souza nandodesouza em gmail.com
Domingo Dezembro 9 13:40:40 BRST 2012


Na verdade Emanuel, estou olhando novamente aqui, e continua a mesma coisa!
A Anova indica haver diferenças entre Manejo, no entanto quando utilizo a
função glht ela diz não haver.
library(car)
library(multcomp)

analise<-lm(PVINICIAL~Gest*Manej,data=GANHO1,
contrast=list(Manej=contr.sum, Gest=contr.sum))
Anova(analise,type="III")
summary(glht(analise,
linfct=mcp(Manej="Tukey",interaction_average=TRUE,covariate_average=TRUE)))

dput(GANHO1)structure(list(Animal = c(25L, 38L, 42L, 50L, 53L, 18L, 22L,
34L, 41L, 44L, 48L, 2L, 12L, 14L, 39L, 46L, 4L, 5L, 21L, 26L,
49L, 3L, 32L, 37L, 43L, 47L, 8L, 15L, 30L, 40L, 1L, 17L, 23L,
28L, 7L, 9L, 19L, 35L), Gest = structure(c(1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L), .Label = c("0", "100", "130", "140"), class = "factor"),
    Manej = structure(c(1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L,
    2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L,
    1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L), .Label = c("1",
    "2"), class = "factor"), PVINICIAL = c(36.5, 37, 25.6, 30.3,
    30.1, 27.3, 31, 25.1, 32.5, 29.6, 35, 40, 34.3, 33.8, 33.7,
    31.4, 40, 32.3, 34.5, 29.5, 39.7, 38, 33.2, 34.5, 34.8, 31.3,
    29.2, 32, 30.5, 27.5, 37, 32, 36.5, 40.6, 35.1, 29.8, 36.2,
    28.6)), .Names = c("Animal", "Gest", "Manej", "PVINICIAL"
), class = "data.frame", row.names = c(NA, -38L))




Em 9 de dezembro de 2012 13:33, Fernando Antonio de souza <
nandodesouza em gmail.com> escreveu:

> oi Emanuel , obrigado pela ajuda. A diferença que observei é porque o glht
> (multcomp) utiliza o contraste soma (default) e o Anova o contraste
> tratamento (default)?
> att
>
> Em 9 de dezembro de 2012 12:18, Emmanuel Arnhold <
> emmanuelarnhold em yahoo.com.br> escreveu:
>
>> analise2<-lm(PVINICIAL~Gest*Manej,data=GANHO1,
>> contrast=list(Manej=contr.sum, Gest=contr.sum))
>
>
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121209/a87c44ad/attachment.html>


Mais detalhes sobre a lista de discussão R-br