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
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.