[R-br] Falha em ANAVA

[Emmanuel Arnhold]

# ajuste seu modelo como exemplo abaixo e faça a anova tipo 3

analise2<-lm(PVINICIAL~Gest*Manej,data=GANHO1, contrast=list(Manej=contr.sum, Gest=contr.sum))
Anova(analise2,type="III")

# ou seja, a análise do multcomp esta correta

# se quiser pode fazer

require(easyanova)

resposta=ea2(GANHO1[,-1])

resposta

names(resposta)

resposta[1]
resposta[11]

# Emmanuel



________________________________
 De: Fernando Antonio de souza <nandodesouza em gmail.com>
Para: r-br em listas.c3sl.ufpr.br 
Enviadas: Domingo, 9 de Dezembro de 2012 11:24
Assunto: Re: [R-br] Falha em ANAVA
 

Ah, é mesmo! ´Segue o código com correções. Quanto ao tipo de ANOVA ambos "II" ou "III" deram o mesmo resultado. e a glht deu resultados diferentes.

Caros colegas, Estou realizando uma análise de variânica do peso vazio 
inicial de ovelhas em função da idade gestacional e manejo. Após rodar a
 análise a análise de variância tipo 2 não indicou haver diferenças para
 o nível do manejo (0,05) entretanto ao utilizar a função glht 
(multcomp) a análise mostrou haver diferenças entre os pesos. Em qual 
análise confiar? segue CMR

library(multcomp)
library(car)

analise<-lm(PVINICIAL~Gest*Manej,data=GANHO1)

Anova(analise,type="II")
anova(analise)
summary(glht(analise,linfct=mcp(Manej="Tukey",interaction_average=TRUE,covariate_average=TRUE)))
summary(glht(analise,linfct=mcp(Gest="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 10:51, Walmes Zeviani <walmeszeviani em gmail.com> escreveu:

Seu CMR não corresponde a descrição do problema que você fez.
>
>À 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
>skype: walmeszeviani
>twitter: @walmeszeviani
>homepage: http://www.leg.ufpr.br/~walmes
>linux user number: 531218
>==========================================================================
>
>_______________________________________________
>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.
>

_______________________________________________
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.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121209/2952f69c/attachment.html>