Regresión Lineal Múltiple y i = b 0 + b 1 x 1i + b 2 x 2i + . . . b k x ki + u i A. Estimación

1. Regresión Lineal Múltipleyi = b0 + b1x1i + b2x2i + . . . bkxki + uiA. Estimación Javier Aparicio División de Estudios Políticos, CIDE javier.aparicio@cide.edu Primavera 2011 http://www.cide.edu/investigadores/aparicio/metodos.html

2. Similitudes con regresión simple • b0 es el intercepto • b1 a bk son k parámetros de pendiente • u es el término de error o residual • El supuesto de media condicional cero se mantiene: E(u|x1,x2, …,xk) = 0 • Igual que antes, minimizamos la suma de residuales cuadrados, de modo que tenemos k+1 condiciones de primer orden (o k+1 parámetros a estimar)

3. Interpretación de la regresión múltiple

4. b interpretada como una “derivada parcial”

5. …“derivada parcial” • La ecuación anterior implica que “regresar y en x1 y x2” tiene el mismo estimador para x1 que regresar y en los residuales de una regresión de x1 en x2 • Es decir, al relacionar x1 con y, solamente capturamos la información de xi1 que no está relacionada con xi2. • Estimamos el efecto de x1 en y después de controlar o aislar el efecto de x2

6. Estimación simple vs. múltiple

7. Suma de cuadrados: Terminología SST es la suma de “desviaciones al cuadrado” de las observaciones de la muestra: es proporcional, más no igual, a VAR(y).

8. Bondad de ajuste: R2 • ¿Cómo saber qué tan bueno es el ajuste entre la regresión y los datos de la muestra? • Podemos calcular la proporción de la Suma de cuadrados totales (SST) que es “explicada” por el modelo. • Esto es la llamada R-cuadrada de una regresión:R2 = SSE/SST = 1 – SSR/SST

9. Bondad de ajuste: R2

10. R-cuadrada: discusión • R2 nunca decrecerá conforme incluyamos más variables explicativas a la regresión, y por lo general aumentará (así sea marginalmente). • ¿Por qué? Incluir variables adicionales aumenta la SSE aunque no sean significativas. • Dado que R2típicamente aumenta con el número de variables independientes, no es por sí sola un buen criterio para comparar modelos.

11. b no sesgadas: supuestos Gauss-Markov Modelo poblacional es lineal en sus parámetros: y = b0 + b1x1 + b2x2 +…+ bkxk+ u Muestra aleatoria de tamaño n, {(xi1, xi2,…, xik, yi): i=1, 2, …, n}, representativa de la población, de modo que el modelo muestral es:yi = b0 + b1xi1 + b2xi2 +…+ bkxik+ ui E(u|x1, x2,… xk) = 0, lo cual implica que todas las variables explicativas son exógenas (no endogeneidad). Ninguna variable x es constante ni tiene una correlación lineal exacta con otra (no multicolinealidad).

12. Demasiadas vs. pocas variables • ¿Si incluimos variables que “no pertenecen al modelo poblacional” en nuestra especificación o modelo? • No tiene impacto en el resto de las b estimadas: MCO permanece sin sesgo. • ¿Si excluimos variables que “sí pertenecen al modelo”? • En general, los estimadores MCO tendrán un sesgo de variable omitida.

13. Sesgo de variable omitida ...ie, la estimación del “modelo incorrecto”. Comparémoslo con la bdel “modelo correcto”

14. Sesgo de variable omitida(continuación)

15. Sesgo de variable omitida(continuación)

16. Sesgo de variable omitida(continuación)

17. Sesgo positivo o negativo en b1

18. Sesgo de variable omitida: resumen • Dos casos donde el sesgo es igual a cero: • b2 = 0, es decir, x2 no pertenecía al modelo poblacional • x1 y x2 no están correlacionados en la muestra • Si la correlación entre (x2, x1) y entre (x2, y) es del mismo signo, el sesgo es positivo. • Si omites una variable x2 que se mueve en el mismo sentido que x1, y ésta afecta positivamente a y, b1 capturará parte de dicho impacto (sobre- estimada). • Si la correlación entre (x2, x1) y entre (x2, y) es de signo opuesto, el sesgo es negativo.

19. El caso más general: sesgo en todas las bi • Técnicamente, sólo podemos anticipar el signo de este sesgo cuando el resto de las variables explicativas incluidas no están correlacionadas entre sí ni con la variable omitida • Si esto no se cumple, el sesgo afecta a todas las bi estimadas, dependiendo de las covarianzas entre las variables incluidas y con la variable omitida. • Aún así, resulta útil calcular el sesgo de variable omitida asumiendo que las otras x no están correlacionadas, aún cuando este supuesto no se cumpla.

20. Varianza de los estimadores MCO • Ya vimos que la “distribución muestral” de los estimadores está centrada en torno a los “verdaderos” parámetros (insesgamiento). • ¿Qué tan dispersa será la distribución de los estimadores? • Para analizar esto, requerimos el 5º supuesto Gauss-Markov: • Var(u|x1, x2,…, xk) = s2conocido como homoscedasticidad (homoskedasticity): varianza constante.

21. Varianza de MCO (cont.) • Sea x igual al vector de variables (x1, x2,…xk) • Suponer que Var(u|x) = s2 también implica que Var(y| x) = s2 • Los 4 supuestos requeridos para insesgamiento, más el supuesto de homoscedasticidad son los llamados supuestos Gauss-Markov.

22. Varianza de MCO (cont.) Es decir, SSTj captura la varianza de xi, mientras que R2j captura la correlación entre xj y las otras x del modelo.

23. Componentes de la Varianza de MCO • Varianza del error: a mayor s2, mayor varianza de los estimadores MCO. • Varianza muestral: a mayor SSTj, menor varianza de los estimadores MCO. • A mayor tamaño de muestra, mayor SSTj y mayor precisión de los estimadores. • Correlación entre las variables explicativas: a mayor Rj2, mayor varianza de los estimadores MCO. • Si dos variables x son altamente correlacionadas, sus b serán poco precisas. • Mayor varianza de los estimadores equivale a decir menor precisióno menor eficiencia.

24. Error de especificación y eficiencia de los estimadores MCO Estimar el modelo incorrecto produce una b1 sesgada (por la variable omitida) ¡pero de menor varianza (mayor precisión)! Un modelo con variables omitidas puede ser engañosamente preciso. Este es el llamado trade-off entre sesgo y eficiencia.

25. Trade-off entre sesgo y eficiencia • La varianza del estimador es menor en el modelo “incorrecto” pero, a menos que b2 = 0, este modelo será sesgado. • Un modelo con variables omitidas puede ser engañosamente preciso y posiblemente sesgado. • Un modelo con demasiadas variables puede ser engañosamente impreciso: pierdes más grados de libertad y enfrentas mayor multicolinearidad. • Conforme el tamaño de la muestra aumenta, la varianza de cada estimador disminuye, haciendo que las diferencias en eficiencia sean relativamente menos importantes.

26. Estimación de la varianza del error • No conocemos la varianza del error, s2, porque no observamos los errores de la población, ui • Lo que observamos son los residuales (estimados) del modelo muestral: • Pero podemos usar los residuales estimados para construir un estimador de la varianza del error.

27. Varianza del error (cont) • gl= n – (k + 1), o bien gl = n – k – 1 • gl (grados de libertad) = (número de observaciones) – (número de parámetros estimados) • A mayores grados de libertad, mayor precisión de los estimadores.

28. Teorema Gauss-Markov • Dados los 5 supuestos Gauss-Markov, puede demostrarse que MCO es “MELI” (BLUE): • Mejor Estimador Lineal Insesgado • Best Linear Unbiased Estimator • De modo que, si los supuestos G-M se sostienen, usar MCO es una buena idea. • Si, además de estos 5 supuestos, u ~ N(0, s2)  MCO es el mejor estimador (lineal o no lineal) insesgado.