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 <[email protected]>: > 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: [email protected] > skype: walmeszeviani > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > > > _______________________________________________ > R-br mailing list > [email protected] > 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 _______________________________________________ R-br mailing list [email protected] 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.
