[R-br] equacao de S.E. para regressao

[Cleber Borges]

Prezados UseRs,

Usando o livro do Faraway (practical regression and anova using r)
adaptei um trecho de código (abaixo) e consigo traçar um gráfico
de modelo de regressão (reta) e com os intervalos de confiança (IC)
da regressão.
Com uso da função *predicit* eu obtenho os IC...

Minha perguntinha é:

Tem alguma formulazinha bonitinha e elegante para se obter esses *s.e.*
da regressão tal como são bonitinhas e elegantes as formulazinhas
de estimativas de coeficientes dos parâmetros e erros segundo OLS??

tipo:
1)    b = (  X^t * X )^1 * X^t * y

2)  Var( b ) = (  X^t * X )^1 *  Var_residual

3)  Algo do gênero de 1) e 2)  ???

Agradeço pelo help com antecedência!
Cleber

ps.:  vi uma formula mas é muita feia!

############
x <- rnorm( 100 )
y <- x + rnorm( 100, s=5 )
reg <- lm( y ~ x )
grid <- seq( -3, 3, 0.1 )
p <- predict( reg, data.frame( x=grid ), se=T )
IC <- cbind( p$fit - 1.96*p$se, p$fit, p$fit + 1.96*p$se )
matplot( grid, IC, lty=1, col=c(2,1,2), type="l", lwd=3 )




---
Este email foi escaneado pelo Avast antivírus.
https://www.avast.com/antivirus
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20151109/e6044f82/attachment.html>