[R-br] regressao linear com constante igual azero
Bernardo Rangel Tura
tura em centroin.com.br
Terça Abril 28 06:04:25 BRT 2015
Bom dia a todos
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
set.seed(1)
var1 <- sample(LETTERS[1:3],1e3,replace=T)
set.seed(2)
var2 <- sample(LETTERS[1:3],1e3,replace=T)
S <- 0.12 + 0.15*(var1=='B') +0.32*(var1=='C') +0.12*(var2=='B')
-0.25*(var2=='C') -rnorm(1e3,0,0.1)
summary(lm(S~var1+var2))
O resultado desta regressão é
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.114883 0.007063 16.27 <2e-16 ***
var1B 0.162707 0.007621 21.35 <2e-16 ***
var1C 0.320609 0.007684 41.72 <2e-16 ***
var2B 0.115608 0.007639 15.13 <2e-16 ***
var2C -0.257842 0.007542 -34.19 <2e-16 ***
Vamos supor que eu queira determinar que o intercepto deste mesma
regressão fosse zero, pelo que aprendi poderia usar
summary(lm(S~0+var1+var2)),
porém o resultado é:
Coefficients:
Estimate Std. Error t value Pr(>|t|)
var1A 0.114883 0.007063 16.27 <2e-16 ***
var1B 0.277590 0.006786 40.91 <2e-16 ***
var1C 0.435491 0.006856 63.52 <2e-16 ***
var2B 0.115608 0.007639 15.13 <2e-16 ***
var2C -0.257842 0.007542 -34.19 <2e-16 ***
Se observarem bem o que aconteceu é que no lugar do intercepto surgiu o
nivel A do factor var1 e não o intercepto=0.
agradeço antecipadamente a ajuda de vocês!
abraços
Tura
Mais detalhes sobre a lista de discussão R-br