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

[Augusto Ribas]

Por favor, perdoe minha falta de atenção, confundi tudo no script anterior.

#exemplo
set.seed(123)
tratamento<-factor(rep(letters[1:3],each=10))
resposta<-rnorm(30,rep(c(1,1,4),each=10))

#modelo usando opções default
modelo1<-lm(resposta~tratamento)
summary(modelo1)

modelo2<- lm(resposta~tratamento,contrasts=list(tratamento=contr.sum))
summary(modelo2)

modelo3<- lm(resposta~tratamento,contrasts=list(tratamento=cbind(c(1,-1,0),c(1,0,-1))))
summary(modelo3)

Mas então no caso acima, eu não deveria obter o mesmo resultado nos 3
casos? Não entendo porque esse 3 exemplos são diferentes.

Em 7 de fevereiro de 2013 10:58, Augusto Ribas <ribas.aca em gmail.com> escreveu:
> Muito obrigado pela ajuda.
>
> Uma ultima curiosidade, o R entende automaticamente no modelo3 que é
> contrastes pela soma?
>
> #exemplo
> set.seed(123)
> tratamento<-factor(rep(letters[1:3],each=10))
> resposta<-rnorm(30,rep(c(1,1,4),each=10))
>
> #modelo usando opções default
> modelo1<-lm(resposta~tratamento)
> summary(modelo1)
>
> modelo2<- lm(resposta~variavel,contrasts=list(tratamento=contr.sum))
> summary(modelo2)
>
> modelo3<- lm(resposta~variavel,contrasts=list(tratamento=cbind(c(1,-1,0),c(1,0,-1))))
> summary(modelo3)
>
> options("contrasts")
>
>
> Em 7 de fevereiro de 2013 09:58, Walmes Zeviani
> <walmeszeviani em gmail.com> escreveu:
>> Esse tipo de atribuição não vai dar certo mesmo
>>
>>
>> options("contrasts")$contrasts[2]<-"contr.sum"
>>
>> Você deve fazer assim
>>
>> options(contrasts=c(unordered="contr.sum", ordered="contr.poly"))
>>
>> Com isso você faz uma alteração global para toda sua sessão R, ou seja, toda
>> lm(), aov(), glm(), lme(), etc, todas funções de ajuste de modelo vão usar o
>> tipo de restrição estabelecido. Eu prefiro por vezes fazer a alteração
>> localmente, ou seja, só dentro da lm() que eu tô rodando, assim, imagine que
>> eu queria blocos com contraste soma mas tratamentos não, ficaria assim
>>
>> m0 <- lm(y~bloco+trat, data=dados, contrasts=list(bloco=contr.sum))
>>
>> À 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



-- 
Grato
Augusto C. A. Ribas

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