<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=windows-1252">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<div class="moz-text-flowed" style="font-family: -moz-fixed;
font-size: 14px;" lang="x-western">Prezados UseRs,
<br>
<br>
Usando o livro do Faraway (practical regression and anova using r)
<br>
adaptei um trecho de código (abaixo) e consigo traçar um gráfico
<br>
de modelo de regressão (reta) e com os intervalos de confiança
(IC)
<br>
da regressão.
<br>
Com uso da função <b class="moz-txt-star"><span
class="moz-txt-tag">*</span>predicit<span class="moz-txt-tag">*</span></b>
eu obtenho os IC...
<br>
<br>
Minha perguntinha é:
<br>
<br>
Tem alguma formulazinha bonitinha e elegante para se obter esses
*s.e.*
<br>
da regressão tal como são bonitinhas e elegantes as formulazinhas
<br>
de estimativas de coeficientes dos parâmetros e erros segundo
OLS?? <span class="moz-smiley-s5" title=":-D"></span> <span
class="moz-smiley-s5" title=":-D"></span>
<br>
<br>
tipo:
<br>
1) b = ( X^t * X )<sup class="moz-txt-sup"><span
style="display:inline-block;width:0;height:0;overflow:hidden">^</span>1</sup>
* X^t * y
<br>
<br>
2) Var( b ) = ( X^t * X )<sup class="moz-txt-sup"><span
style="display:inline-block;width:0;height:0;overflow:hidden">^</span>1</sup>
* Var_residual
<br>
<br>
3) Algo do gênero de 1) e 2) ???
<br>
<br>
Agradeço pelo help com antecedência!
<br>
Cleber
<br>
<br>
ps.: vi uma formula mas é muita feia!
<br>
<br>
############
<br>
x <- rnorm( 100 )
<br>
y <- x + rnorm( 100, s=5 )
<br>
reg <- lm( y ~ x )
<br>
grid <- seq( -3, 3, 0.1 )
<br>
p <- predict( reg, data.frame( x=grid ), se=T )
<br>
IC <- cbind( p$fit - 1.96*p$se, p$fit, p$fit + 1.96*p$se )
<br>
matplot( grid, IC, lty=1, col=c(2,1,2), type="l", lwd=3 )
<br>
<br>
<br>
</div>
<br /><br />
<hr style='border:none; color:#909090; background-color:#B0B0B0; height: 1px; width: 99%;' />
<table style='border-collapse:collapse;border:none;'>
<tr>
<td style='border:none;padding:0px 15px 0px 8px'>
<a href="https://www.avast.com/antivirus">
<img border=0 src="http://static.avast.com/emails/avast-mail-stamp.png" alt="Avast logo" />
</a>
</td>
<td>
<p style='color:#3d4d5a; font-family:"Calibri","Verdana","Arial","Helvetica"; font-size:12pt;'>
Este email foi escaneado pelo Avast antivírus.
<br><a href="https://www.avast.com/antivirus">www.avast.com</a>
</p>
</td>
</tr>
</table>
<br />
</body>
</html>