Em 28/04/2015 14:10, "Robert Iquiapaza" <rbali@ufmg.br> escreveu:

Se é isso o que deseja:

dumvar1 = table(1:length(var1),as.factor(var1)) 

dumvar2 = table(1:length(var2),as.factor(var2)) 

summary(lm(S~0+dumvar1[,-1]+dumvar2[,-1]))

Coefficients:

                Estimate Std. Error t value Pr(>|t|)    

dumvar1[, -1]B  0.234001   0.007011   33.38   <2e-16 ***

dumvar1[, -1]C  0.391907   0.007097   55.22   <2e-16 ***

dumvar2[, -1]B  0.182071   0.007258   25.09   <2e-16 ***

dumvar2[, -1]C -0.190209   0.007076  -26.88   <2e-16 ***

---

Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Att

Robert

#cuidado com a interpretação dos coeficientes.