[R-br] Duvida sobre como alterar o contraste no options de "contr.poly" para "contr.sum" e entender o que acontece.

Augusto Ribas ribas.aca em gmail.com
Quinta Fevereiro 7 16:09:18 BRST 2013


Muito obrigado pela explicação, eu que não estou sabendo montar as
comparações então, com esses exemplos eu acho que entendi melhor.

2013/2/7 Walmes Zeviani <walmeszeviani em gmail.com>:
> O problema do seu código foi só a ordem das colunas, veja
>
> set.seed(123)
> tratamento <- factor(rep(letters[1:3],each=10))
> resposta <- rnorm(30,rep(c(1,1,4),each=10))
>
> modelo1 <- lm(resposta~tratamento)
> summary(modelo1)
>
> modelo2 <- lm(resposta~tratamento,
>               contrasts=list(tratamento=contr.sum))
> summary(modelo2)
>
> contrasts(tratamento)
> C(tratamento, contr=contr.sum)
>
> modelo3 <- lm(resposta~tratamento,
>               contrasts=list(tratamento=cbind(c(1,-1,0), c(1,0,-1))))
> summary(modelo3)
>
> C(tratamento, contr=contr.sum) # ordem das colunas!!!
>
> modelo4 <- lm(resposta~tratamento,
>               contrasts=list(tratamento=cbind(c(1,0,-1), c(0,1,-1))))
> summary(modelo4)
>
> cbind(coef(modelo2), coef(modelo3), coef(modelo4))
>
> À 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.



-- 
Grato
Augusto C. A. Ribas

Site Pessoal: http://recologia.wordpress.com/
Lattes: http://lattes.cnpq.br/7355685961127056


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