Guía de Regresión Cuadrática
Guía de Regresión Cuadrática
Facultad de Economía
Escuela Profesional de Economía
Departamento Académico de Economía
Asignatura
Econometría 1
Guía de práctica
Modelo de regresión polinomial de grado 2:
modelo de regresión cuadrática
Arequipa–Perú
2020
Carlos Pedro Vera Ninacondor
Contenido
Bibliografía
2
Carlos Pedro Vera Ninacondor
3
Carlos Pedro Vera Ninacondor
Interpretación
En los modelos no lineales (con respecto de las 𝑋) la función de regresión se interpreta
mejor mediante su representación gráfica y mediante el cálculo del efecto esperado
sobre 𝑌 de la variación de una o más variables independientes (Stock y Watson, 2012,
p. 187).
1.2. Modelo de regresión muestral cuadrática
El modelo de regresión muestral cuadrática (Stock y Watson, 2012, p. 183) puede
expresarse como,
𝑌𝑖 = 𝛽̂0 + 𝛽̂1 𝑋𝑖 + 𝛽̂2 𝑋𝑖2 + 𝑢̂𝑖
Donde
𝑖 es el subíndice que toma valores de 1 hasta 𝑛.
𝑛 es el número de observaciones en la muestra.
𝑌 es la variab|le dependiente (o el regresando).
𝑌𝑖 es el valor de la observación i-ésima de 𝑌 en la muestra.
𝑋 es la variable independiente (o el regresor uno).
𝑋𝑖 es el valor de la observación i-ésima de 𝑋 en la muestra.
𝑋 2 es el cuadrado de la variable independiente 𝑋 (o el regresor dos).
𝑋𝑖2 es el cuadrado del valor de la observación i-ésima de 𝑋 en la muestra.
𝛽̂0 es el estimador del coeficiente 𝛽0 de la función de regresión poblacional cuadrática.
𝛽̂1 es el estimador del coeficiente 𝛽1 de la función de regresión poblacional cuadrática.
𝛽̂2 es el estimador del coeficiente 𝛽2 de la función de regresión poblacional cuadrática.
𝑢̂𝑖 es el residuo, es el estimador del término de error 𝑢𝑖 .
1.3. Estimación de los coeficientes de la función de regresión poblacional
cuadrática
Para la estimación de los coeficientes de la función de regresión poblacional cuadrática se
utiliza el método de mínimos cuadrados ordinarios.
Estimador de mínimos cuadrados ordinarios (MCO)
Iniciando con el modelo de regresión muestral cuadrática
𝑌𝑖 = 𝛽̂0 + 𝛽̂1 𝑋𝑖 + 𝛽̂2 𝑋𝑖2 + 𝑢̂𝑖
4
Carlos Pedro Vera Ninacondor
Despejando 𝑢̂𝑖
𝛽̂0 + 𝛽̂1 𝑋𝑖 + 𝛽̂2 𝑋𝑖2 + 𝑢̂𝑖 = 𝑌𝑖
𝑢̂𝑖 = 𝑌𝑖 − 𝛽̂0 − 𝛽̂1 𝑋𝑖 − 𝛽̂2 𝑋𝑖2
Elevando al cuadrado ambos lados de la ecuación
2
𝑢̂𝑖2 = (𝑌𝑖 − 𝛽̂0 − 𝛽̂1 𝑋𝑖 − 𝛽̂2 𝑋𝑖2 )
Aplicando sumatorias a ambos lados de la ecuación
𝑛 𝑛
2
∑ 𝑢̂𝑖2 = ∑(𝑌𝑖 − 𝛽̂0 − 𝛽̂1 𝑋𝑖 − 𝛽̂2 𝑋𝑖2 )
𝑖=1 𝑖=1
El principio o método de mínimos cuadrados elige 𝛽̂0 , 𝛽̂1 y 𝛽̂2 de manera que, para una
muestra o conjunto de datos determinados, ∑ 𝑢̂𝑖2 es la más pequeña posible. En otras
palabras, para una muestra dada, proporciona valores estimados únicos de 𝛽̂0 , 𝛽̂1 y 𝛽̂2 que
producen el valor más pequeño o reducido posible de ∑ 𝑢̂𝑖2 (Gujarati y Porter, 2010, p. 57).
Minimizando la suma de los residuos elevados al cuadrado
Calculando las derivadas parciales de primer orden de ∑ 𝑢̂𝑖2 respecto de β̂0, de 𝛽̂1 y de β̂2
se obtiene
𝑛 𝑛
𝜕(∑ 𝑢̂𝑖2 )
= 2 ∑(𝑌𝑖 − 𝛽̂0 − 𝛽̂1 𝑋𝑖 − 𝛽̂2 𝑋2𝑖 ) (−1) = −2 ∑(𝑌𝑖 − 𝛽̂0 − 𝛽̂1 𝑋𝑖 − 𝛽̂2 𝑋2𝑖 )
𝜕𝛽̂0
𝑖=1 𝑖=1
𝑛 𝑛
𝜕(∑ 𝑢̂𝑖2 )
= 2 ∑(𝑌𝑖 − 𝛽̂0 − 𝛽̂1 𝑋𝑖 − 𝛽̂2 𝑋2𝑖 ) (−𝑋𝑖 ) = −2 ∑(𝑌𝑖 − 𝛽̂0 − 𝛽̂1 𝑋𝑖 − 𝛽̂2 𝑋2𝑖 ) (𝑋𝑖 )
𝜕𝛽̂1
𝑖=1 𝑖=1
𝑛 𝑛
𝜕(∑ 𝑢̂𝑖2 )
= 2 ∑(𝑌𝑖 − 𝛽̂0 − 𝛽̂1 𝑋𝑖 − 𝛽̂2 𝑋2𝑖 ) (−𝑋2𝑖 ) = −2 ∑(𝑌𝑖 − 𝛽̂0 − 𝛽̂1 𝑋𝑖 − 𝛽̂2 𝑋2𝑖 ) (𝑋2𝑖 )
̂
𝜕𝛽2
𝑖=1 𝑖=1
Seguidamente igualando a cero las derivadas parciales de primer orden, y luego resolviendo
el sistema de ecuaciones se obtienen los estimadores MCO de 𝛽0 , 𝛽1 y 𝛽2
Igualando a cero la derivada parcial de primer orden de ∑ 𝑢̂𝑖2 respecto de β̂0
𝑛
Despejando ∑ 𝑌𝑖
𝑛
5
Carlos Pedro Vera Ninacondor
𝑛 𝑛 𝑛 𝑛
Despejando ∑ 𝑌𝑖 𝑋𝑖
𝑛
Despejando ∑ 𝑌𝑖 𝑋𝑖2
𝑛
6
Carlos Pedro Vera Ninacondor
∑ 𝑌𝑖 𝑛 ∑ 𝑋𝑖 ∑ 𝑋𝑖2
𝑖=1 𝑖=1 𝑖=1
𝑛 𝑛 𝑛 𝑛 𝛽̂0
∑ 𝑌𝑖 𝑋𝑖 = ∑ 𝑋𝑖 ∑ 𝑋𝑖2 ∑ 𝑋𝑖3 × [𝛽̂1 ]
𝑖=1
𝑛
𝑖=1
𝑛
𝑖=1
𝑛
𝑖=1
𝑛
𝛽̂2 3×1
∑ 𝑌𝑖 𝑋𝑖2 ∑ 𝑋𝑖2 ∑ 𝑋𝑖3 ∑ 𝑋𝑖4
[ 𝑖=1 ]3×1 [ 𝑖=1 𝑖=1 𝑖=1 ]3×3
7
Carlos Pedro Vera Ninacondor
Despejando 𝛽0 , 𝛽1 y 𝛽2
𝑛 𝑛 𝑛
𝑛 ∑ 𝑋𝑖 ∑ 𝑋𝑖2 ∑ 𝑌𝑖
𝑖=1 𝑖=1 𝑖=1
𝑛 𝑛 𝑛 𝛽̂0 𝑛
𝑛 ∑ 𝑋𝑖 ∑ 𝑋𝑖2 ∑ 𝑌𝑖
𝑖=1 𝑖=1 𝑖=1
𝛽̂0 𝑛 𝑛 𝑛 𝑛
8
Carlos Pedro Vera Ninacondor
9
Carlos Pedro Vera Ninacondor
Donde
𝑁𝑜𝑡𝑎𝑖 es la nota lograda del estudiante i-ésimo en la población.
𝐻𝑜𝑟𝑎𝑠𝑖 es el número de horas de estudio semanal fuera de clase del estudiante i-ésimo
en la población.
𝐻𝑜𝑟𝑎𝑠𝑖2 es el cuadrado del número de horas de estudio semanal fuera de clase del
estudiante i-ésimo en la población.
𝛽0 + 𝛽1 𝐻𝑜𝑟𝑎𝑠𝑖 + 𝛽2 𝐻𝑜𝑟𝑎𝑠𝑖2 es la función de regresión poblacional cuadrática.
𝛽0 es el término constante de la función de regresión poblacional, el cual es desconocido.
𝛽1 es el coeficiente de 𝐻𝑜𝑟𝑎𝑠𝑖 de la función de regresión poblacional, el cual es
desconocido.
𝛽2 es el coeficiente de 𝐻𝑜𝑟𝑎𝑠𝑖2 de la función de regresión poblacional, el cual es
desconocido.
𝑢𝑖 es el término de error, el cual es desconocido.
1.6.3. Especificación del modelo de regresión muestral cuadrática
Considerando lo presentado en el punto 1.2, para este caso, se especifica el modelo de
regresión muestral cuadrática siguiente,
𝑁𝑜𝑡𝑎𝑖 = 𝛽̂0 + 𝛽̂1 𝐻𝑜𝑟𝑎𝑠𝑖 + 𝛽̂2 𝐻𝑜𝑟𝑎𝑠𝑖2 + 𝑢̂𝑖
Donde
𝑁𝑜𝑡𝑎𝑖 es la nota lograda del estudiante i-ésimo en la muestra.
𝐻𝑜𝑟𝑎𝑠𝑖 es el número de horas de estudio semanal fuera de clase del estudiante i-ésimo
en la muestra.
𝐻𝑜𝑟𝑎𝑠𝑖2 es el cuadrado del número de horas de estudio semanal fuera de clase del
estudiante i-ésimo en la muestra.
𝛽̂0 es el estimador del término constante 𝛽0 de la función de regresión poblacional.
𝛽̂1 es el estimador del coeficiente 𝛽1 de la función de regresión poblacional, y se espera
que sea positivo.
𝛽̂2 es el estimador del coeficiente 𝛽2 de la función de regresión poblacional, y se espera
que sea negativo.
𝑢̂𝑖 es el residuo, es el estimador del término de error 𝑢𝑖 .
11
Carlos Pedro Vera Ninacondor
12
Carlos Pedro Vera Ninacondor
𝑛 ∑ 𝑋𝑖 ∑ 𝑋𝑖2 ∑ 𝑌𝑖
𝑖=1 𝑖=1 𝑖=1
𝛽̂0 𝑛 𝑛 𝑛 𝑛
Simbólicamente
𝛽̂ = (𝑋 ′ 𝑋)−1 𝑋′𝑦
Cálculo de la matriz inversa, es decir, la matriz (𝑋 ′ 𝑋)−1
Cálculo del determinante de la matriz
𝐷 = (35 × 9694 × 7021270) + (482 × 243458 × 9694) + (9694 × 482 × 243458) −
(9694 × 9694 × 9694) − (243458 × 243458 × 35) − (7021270 × 482 × 482)
𝐷 = 2 382 246 698 300 + 1 137 559 452 664 + 1 137 559 452 664 −
910 980 427 384 − 2 074 512 921 740 − 1 631 209 531 480
𝐷 = 40 662 723 024
Cálculo de los menores de la matriz
243 458 | | 482 9 694
| 9 694 243 458 | | 482 |
243 458 7 021 270 9 694 7 021 270 9 694 243 458
482 9 694 35 9 694 35 482
𝑀= | | | | | |
243 458 7 021 270 9 694 7 021 270 9 694 243 458
482 9 694 35 9 694 35 482
[ |9 694 243 458| |
482 243 458
| | |
482 9 694 ]
8 792 393 616 1 024 170 288 23 373 120
𝑀 = [1 024 170 288 151 770 814 3 848 522 ]
23 373 120 3 848 522 106 966
13
Carlos Pedro Vera Ninacondor
𝛽̂0 7,275014710
[𝛽̂1 ] = [ 0,563203622 ]
𝛽̂2 −0,010840691 3×1
3×1
14
Carlos Pedro Vera Ninacondor
15
Carlos Pedro Vera Ninacondor
Tabla 3. Predicción de la nota lograda dado un valor del número de horas de estudio
̂ 𝑖 = 7,275014710 + 0,563203622𝐻𝑜𝑟𝑎𝑠𝑖
𝑁𝑜𝑡𝑎 Predicción
𝑛𝑖 𝐻𝑜𝑟𝑎𝑠𝑖
−0,010840691𝐻𝑜𝑟𝑎𝑠𝑖2 ̂𝑖
𝑁𝑜𝑡𝑎
1 5 7,275014710 + 0,563203622(5) − 0,010840691(5)2 9,820015545
2 5 7,275014710 + 0,563203622(5) − 0,010840691(5)2 9,820015545
3 7 7,275014710 + 0,563203622(7) − 0,010840691(7)2 10,686246205
4 5 7,275014710 + 0,563203622(5) − 0,010840691(5)2 9,820015545
5 6 7,275014710 + 0,563203622(6) − 0,010840691(6)2 10,263971566
6 3 7,275014710 + 0,563203622(3) − 0,010840691(3)2 8,867059357
7 29 7,275014710 + 0,563203622(29) − 0,010840691(29)2 14,490898617
8 35 7,275014710 + 0,563203622(35) − 0,010840691(35)2 13,707295005
9 3 7,275014710 + 0,563203622(3) − 0,010840691(3)2 8,867059357
10 16 7,275014710 + 0,563203622(16) − 0,010840691(16)2 13,511055766
11 9 7,275014710 + 0,563203622(9) − 0,010840691(9)2 11,465751337
12 13 7,275014710 + 0,563203622(13) − 0,010840691(13)2 12,764585017
13 10 7,275014710 + 0,563203622(10) − 0,010840691(10)2 11,822981830
14 7 7,275014710 + 0,563203622(7) − 0,010840691(7)2 10,686246205
15 15 7,275014710 + 0,563203622(15) − 0,010840691(15)2 13,283913565
16 20 7,275014710 + 0,563203622(20) − 0,010840691(20)2 14,202810750
17 20 7,275014710 + 0,563203622(20) − 0,010840691(20)2 14,202810750
18 24 7,275014710 + 0,563203622(24) − 0,010840691(24)2 14,547663622
19 10 7,275014710 + 0,563203622(10) − 0,010840691(10)2 11,822981830
20 15 7,275014710 + 0,563203622(15) − 0,010840691(15)2 13,283913565
21 15 7,275014710 + 0,563203622(15) − 0,010840691(15)2 13,283913565
22 21 7,275014710 + 0,563203622(21) − 0,010840691(21)2 14,321546041
23 5 7,275014710 + 0,563203622(5) − 0,010840691(5)2 9,820015545
24 21 7,275014710 + 0,563203622(21) − 0,010840691(21)2 14,321546041
25 5 7,275014710 + 0,563203622(5) − 0,010840691(5)2 9,820015545
26 3 7,275014710 + 0,563203622(3) − 0,010840691(3)2 8,867059357
27 24 7,275014710 + 0,563203622(24) − 0,010840691(24)2 14,547663622
28 13 7,275014710 + 0,563203622(13) − 0,010840691(13)2 12,764585017
29 40 7,275014710 + 0,563203622(40) − 0,010840691(40)2 12,458053990
30 5 7,275014710 + 0,563203622(5) − 0,010840691(5)2 9,820015545
31 6 7,275014710 + 0,563203622(6) − 0,010840691(6)2 10,263971566
32 19 7,275014710 + 0,563203622(19) − 0,010840691(19)2 14,062394077
33 24 7,275014710 + 0,563203622(24) − 0,010840691(24)2 14,547663622
34 18 7,275014710 + 0,563203622(18) − 0,010840691(18)2 13,900296022
35 6 7,275014710 + 0,563203622(6) − 0,010840691(6)2 10,263971566
35 35
35
∑ 𝑋𝑖 1 1 ∑ 𝑌̂𝑖
𝑛 = 35 𝑖=1
𝑌̂ = × ∑ 𝑌̂𝑖 = × 421 = 12,028571429 𝑖=1
𝑛 35
𝑖=1
482 421
Fuente: Tabla 1, ítem 1.6.6.
Elaboración: Autoría propia.
16
Carlos Pedro Vera Ninacondor
17
Carlos Pedro Vera Ninacondor
Bibliografía
INTERAMERICANA EDITORES, S. A. de C. V.
PEARSON EDUCACIÓN, S. A.
18