Ana, vc precisara' de um estatistico local para maiores esclarecimentos. A ordem com que os itens aparecem no seu ficheiro e' irrelevante para as tarefas executadas pelo lm(). O que importa e' como os itens sao ordenados alfabeticamente (por isso o R toma 'C' como referencia). Antes de executar o lm(), vc pode trocar a referencia usando o comando relevel(). Veja o exemplo na ajuda do comando digitando help('relevel') . O exemplo dado e' exatamente usando lm().... b 2012/11/28 <alanarocha@sapo.pt>
modC<-lm(pthi~dose+grupo) modC
Call: lm(formula = pthi ~ dose + grupo)
Coefficients: (Intercept) dose grupoH 752.8614 0.1791 -159.5864 summary(modC)
Call: lm(formula = pthi ~ dose + grupo)
Residuals: Min 1Q Median 3Q Max -862.97 -323.87 -60.89 222.13 1936.94
Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 752.86139 84.56014 8.903 1.22e-15 *** dose 0.17913 0.08253 2.171 0.0315 * grupoH -159.58641 92.49382 -1.725 0.0864 . --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 540.3 on 158 degrees of freedom (1 observation deleted due to missingness) Multiple R-squared: 0.03447, Adjusted R-squared: 0.02225 F-statistic: 2.82 on 2 and 158 DF, p-value: 0.0626
eu fiz isto, percebo que o grupo C aparece por defeito do R. ams como esplico no relatorio qual o valor estimado para o grupo C e para o grupo H? E se no meu ficheiro de dados o H aparece primeiro porque é que o R considera o C primeiro? Ana
______________________________**_________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-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<http://www.leg.ufpr.br/r-br-guia>) e forneça código mínimo reproduzível.