<div dir="ltr">Tura,<div><br></div><div>No seu exemplo o resultado da regressão está gritando que o intercepto é diferente de zero com p-valor <span style="font-size:12.8000001907349px"><2e-16.</span></div><div><span style="font-size:12.8000001907349px"><br></span></div><div><span style="font-size:12.8000001907349px">Se o intercepto estatisticamente indistinguível de zero o p-valor deveria dar próximo de um...</span></div></div><div class="gmail_extra"><br><div class="gmail_quote">2015-04-28 6:04 GMT-03:00 Bernardo Rangel Tura <span dir="ltr"><<a href="mailto:tura@centroin.com.br" target="_blank">tura@centroin.com.br</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Bom dia a todos<br>
<br>
Estou esbarrando com um problema que não estou conseguindo resolver sozinho a pesar de achar que é muito simples. Segue a simulação do problema<br>
<br>
set.seed(1)<br>
var1 <- sample(LETTERS[1:3],1e3,replace=T)<br>
<br>
set.seed(2)<br>
var2 <- sample(LETTERS[1:3],1e3,replace=T)<br>
<br>
S <-  0.12 + 0.15*(var1=='B') +0.32*(var1=='C') +0.12*(var2=='B') -0.25*(var2=='C') -rnorm(1e3,0,0.1)<br>
<br>
summary(lm(S~var1+var2))<br>
<br>
O resultado desta regressão é<br>
Coefficients:<br>
             Estimate Std. Error t value Pr(>|t|)<br>
(Intercept)  0.114883   0.007063   16.27   <2e-16 ***<br>
var1B        0.162707   0.007621   21.35   <2e-16 ***<br>
var1C        0.320609   0.007684   41.72   <2e-16 ***<br>
var2B        0.115608   0.007639   15.13   <2e-16 ***<br>
var2C       -0.257842   0.007542  -34.19   <2e-16 ***<br>
<br>
Vamos supor que eu queira determinar que o intercepto deste mesma regressão fosse zero, pelo que aprendi poderia usar<br>
summary(lm(S~0+var1+var2)),<br>
porém o resultado é:<br>
<br>
Coefficients:<br>
       Estimate Std. Error t value Pr(>|t|)<br>
var1A  0.114883   0.007063   16.27   <2e-16 ***<br>
var1B  0.277590   0.006786   40.91   <2e-16 ***<br>
var1C  0.435491   0.006856   63.52   <2e-16 ***<br>
var2B  0.115608   0.007639   15.13   <2e-16 ***<br>
var2C -0.257842   0.007542  -34.19   <2e-16 ***<br>
<br>
Se observarem bem o que aconteceu é que no lugar do intercepto surgiu o nivel A do factor var1 e não o intercepto=0.<br>
<br>
agradeço antecipadamente a ajuda de vocês!<br>
<br>
abraços<br>
Tura<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br>
</blockquote></div><br></div>