Inferencia en modelos de regresión lineal múltiple

Seguimos caminando

Autor/a

Román Salmerón Gómez y Catalina García García

Ejercicios

1.- Dado el modelo de regresión lineal simple \(\mathbf{S} = \beta_{1} + \beta_{2} \cdot \mathbf{E} + \mathbf{u}\), donde \(\mathbf{S}\) es el salario medio por hora (medido en euros) y \(\mathbf{E}\) son los años de escolarización, interprete los coeficientes de las variables independientes que son significativamente distintos de cero e indique si el modelo es globalmente válido a partir de la siguiente información:

Salario medio (euros/hora) Escolarización (años)
4.45 6
5.77 7
7.31 10
7.81 12
10.67 15
13.61 17

2.- Dado el modelo de regresión lineal múltiple \(\mathbf{y} = \beta_{1} + \beta_{2} \cdot \mathbf{V} + \beta_{3} \cdot \mathbf{Z} + \mathbf{u}\) con \(n=15\) observaciones y \(k=3\) variables independientes, en los siguientes casos, responda de forma razonada a las siguientes cuestiones:

  • Indique qué coeficientes son significativamente distintos de cero usando los intervalos de confianza de los coeficientes de las variables independientes: \(\widehat{\beta}_{i} \pm t_{n-k} \cdot \sqrt{\widehat{var} \left(\widehat{\beta}_{i} \right)}\). Interprételos.
  • Indique qué coeficientes son significativamente distintos de cero usando los contrastes de significación individual de los coeficientes de las variables independientes: \(t_{exp} = \left| \frac{\widehat{\beta}_{i}}{\sqrt{\widehat{var} \left(\widehat{\beta}_{i} \right)}} \right|\).
  • ¿Es el modelo válido? \(F_{exp} = \frac{SCE/(k-1)}{SCR/(n-k)}\).

Responda de forma razonada todas las cuestiones planteadas usando, al menos, cuatro decimales.

(Caso A)

\(\left( \mathbf{y}, \mathbf{X} \right) =\)

     y cte  V  Z
1  -16   1 -7  2
2   -2   1 -6 -3
3  -10   1 -5  1
4    1   1 -4 -2
5  -16   1 -3  4
6   -3   1 -2 -1
7   -9   1 -1  3
8    1   1  0  0
9   14   1  1 -3
10   3   1  2  1
11  20   1  3 -4
12   4   1  4  2
13  16   1  5 -1
14   6   1  6  3
15  18   1  7 -2

\(\left( \mathbf{X}^{t} \mathbf{X} \right)^{-1} =\)

           cte            V            Z
cte 0.06666667 0.0000000000 0.0000000000
V   0.00000000 0.0035924233 0.0004898759
Z   0.00000000 0.0004898759 0.0114304376

\(\mathbf{X}^{t} \mathbf{y} =\)

       y
cte   27
V    571
Z   -271

(Caso B)

\(\left( \mathbf{y}, \mathbf{X} \right) =\)

     y cte  X  Z
1   -8   1 -7  0
2   -5   1 -6 -1
3  -10   1 -5  1
4   -4   1 -4 -1
5   -5   1 -3  1
6    1   1 -2 -1
7   -2   1 -1  1
8    0   1  0  0
9    5   1  1  1
10  10   1  2 -1
11  10   1  3  1
12  17   1  4 -1
13  10   1  5  1
14  19   1  6 -1
15  20   1  7  0

\(\left( \mathbf{X}^{t} \mathbf{X} \right)^{-1} =\)

           cte           X          Z
cte 0.06666667 0.000000000 0.00000000
X   0.00000000 0.003571429 0.00000000
Z   0.00000000 0.000000000 0.08333333

\(\mathbf{X}^{t} \mathbf{y} =\)

      y
cte  58
X   594
Z   -30

(Caso C)

\(\left( \mathbf{y}, \mathbf{X} \right) =\)

     y cte  X  Z
1  -15   1 -7  1
2   -9   1 -6  0
3   -8   1 -5  0
4   -4   1 -4  0
5   -3   1 -3  0
6    1   1 -2  0
7   23   1 -1 -7
8    1   1  0  0
9  -15   1  1  7
10   5   1  2  0
11   7   1  3  0
12  13   1  4  0
13  11   1  5  0
14  16   1  6  0
15  18   1  7 -1

\(\left( \mathbf{X}^{t} \mathbf{X} \right)^{-1} =\)

           cte           X    Z
cte 0.06666667 0.000000000 0.00
X   0.00000000 0.003571429 0.00
Z   0.00000000 0.000000000 0.01

\(\mathbf{X}^{t} \mathbf{y} =\)

       y
cte   41
X    544
Z   -299

3.- Un concesionario de coches desea analizar la cuantía de las ventas realizadas (en euros) a partir del número de ventas realizadas y si la venta se ha realizado en metálico. ¿Qué variables independientes influyen en la cuantía de las ventas realizadas? ¿Cómo influyen? ¿Cómo interpretaría el término indpendiente? ¿El modelo es válido?

Responda de forma razonada teniendo en cuenta que a partir de 210 observaciones se ha obtenido que \(\boldsymbol{\beta} = (66'85, 66'1, -89'97)\), \(\widehat{\sigma} = 164'27\), \(R^{2} = 0'6248\) y \[\widehat{var} \left( \widehat{\boldsymbol{\beta}} \right) = \left( \begin{array}{ccc} 23'396^{2} & \diamond &\diamond \\ \diamond & 3'814^{2} & \diamond \\ \diamond & \diamond & 22'96^{2} \\ \end{array} \right).\]

4.- Se desea lanzar una campaña de publicidad sobre cierto vehículo, por lo que se analiza qué factores influyen en su consumo. Concretamente se dispone de 150 observaciones sobre el gasto promedio (\(\textbf{G}\), en litros cada 100 km), la portencia del motor (\(\textbf{M}\), en caballos de vapor) y velocidad máxima (\(\textbf{V}\), en km/h). Conteste de forma razonada a las siguientes cuestiones:

  • Especifique el modelo econométrico del modelo económico planteado.
  • Interprete los coeficientes significativamente distintos de cero sabiendo que \(\boldsymbol{\beta} = (-123'72, -2'08, 3'03)\) y \[\widehat{var} \left( \widehat{\boldsymbol{\beta}} \right) = \left( \begin{array}{ccc} 11'69 & \diamond & \diamond \\ \diamond & 0'00048 & \diamond \\ \diamond & \diamond & 0'0022 \\ \end{array} \right).\]
  • ¿El modelo es conjuntamente válido? Tenga en cuenta que \(\overline{R}^{2} = 0'9998\) y \(\widehat{\sigma}^{2} = 0'005\).

5.- El bar de mi barrio desea analizar la relación entre las ventas diarias en metálico (\(\textbf{V}\), en euros) a partir del número de compras en metálico diarias realizadas \(\textbf{NC}\). A partir de 105 observaciones se tiene que:

\[\widehat{\textbf{V}} = 16'4291 + 54'746 \cdot \textbf{NC}, \quad \overline{R}^{2} = 0'6083, \ SCE = 2242535.\]

Responda de forma razona las siguientes cuestiones:

  • Sabiendo que las desviaciones típicas estimadas del término independiente y del número de compras son, respectivamente, 18’852 y 4’294, interprete los coeficientes significativamente distintos de cero.
  • ¿Se puede afirmar, al 90% de confianza, que cada compra en metálico suponen unos ingresos de 55 euros?
  • ¿Es el modelo conjuntamente válido?

6.- Para 200 empresas se tiene información sobre su producción (\(\mathbf{Y}\)), trabajo (\(\mathbf{L}\)) y capital (\(\mathbf{K}\)), de forma que se obtienen los siguientes resultados: \[\widehat{log \textbf{Y}} = 0'796 + 0'602 \cdot log \textbf{L} + 0'348 \cdot log \textbf{K}, \quad SCR = 0'051, \ R^{2} = 0'9509.\] Además, se tiene que: \[\left( \mathbf{X}^{t} \mathbf{X} \right)^{-1} = \left( \begin{array}{ccc} 14'31 & \diamond & \diamond\\ \diamond & 0'605 & \diamond \\ \diamond & \diamond & 0'36 \\ \end{array} \right).\]

Se pide contestar de forma razonada las siguientes cuestiones:

  • Interprete los coeficientes significativamente distintos de cero.
  • ¿Es el modelo conjuntamente válido?
  • ¿Podría comprobar si los rendimientos son constantes? En caso negativo, ¿qué información le falta?

7.- Para los tres casos del segundo ejercicio contrastar, al 5% de significación, las siguientes hipótesis: a) \(\beta_{2} = \beta_{3}\), b) \(\beta_{1} + \beta_{2} = 1\) y c) \(\beta_{1} + 2 \cdot \beta_{2} = 0\), \(\beta_{2} - 0'5 \cdot \beta_{3} = 2\).

Soluciones

Ejercicio 1

$XtX
            cte Escolaridad
cte           6          67
Escolaridad  67         843

$`inversa XtX`
                   cte Escolaridad
cte          1.4815466 -0.11775044
Escolaridad -0.1177504  0.01054482

$Xty
            Salario
cte           49.62
Escolaridad  625.33

$beta
               Salario
cte         -0.1185413
Escolaridad  0.7512127

$`Desviaciones estimadas`
[1] 1.00644786 0.08490888

$`Intervalo CTE`
[1] -2.912889  2.675806

$`Intervalo E`
[1] 0.5154678 0.9869575

$`Valor experimental CTE`
[1] -0.1177819

$`Valor experimental E`
[1] 8.84728

$`Valor teórico`
[1] 2.776445

$`Valor experimental ANOVA`
[1] 78.27436

$`Valor teórico`
[1] 7.708647

Ejercicio 2

(Caso A)

$XtX
    cte   V   Z
cte  15   0   0
V     0 280 -12
Z     0 -12  88

$`inversa XtX`
           cte            V            Z
cte 0.06666667 0.0000000000 0.0000000000
V   0.00000000 0.0035924233 0.0004898759
Z   0.00000000 0.0004898759 0.0114304376

$Xty
       y
cte   27
V    571
Z   -271

$beta
            y
cte  1.800000
V    1.918517
Z   -2.817929

$`Desviaciones estimadas`
[1] 0.4550197 0.1056257

$`Intervalo CTE`
[1] 0.8085972 2.7914028

$`Intervalo E`
[1] 1.688379 2.148656

$`Valor experimental CTE`
[1] 3.955872

$`Valor experimental E`
[1] 18.16336

$`Valor teórico`
[1] 2.178813

$`Valor experimental ANOVA`
[1] 299.3151

$`Valor teórico`
[1] 3.885294

(Caso B)

$XtX
    cte   X  Z
cte  15   0  0
X     0 280  0
Z     0   0 12

$`inversa XtX`
           cte           X          Z
cte 0.06666667 0.000000000 0.00000000
X   0.00000000 0.003571429 0.00000000
Z   0.00000000 0.000000000 0.08333333

$Xty
      y
cte  58
X   594
Z   -30

$beta
            y
cte  3.866667
X    2.121429
Z   -2.500000

$`Desviaciones estimadas`
[1] 0.5302241 0.1227230

$`Intervalo CTE`
[1] 2.711408 5.021926

$`Intervalo E`
[1] 1.854038 2.388819

$`Valor experimental CTE`
[1] 7.292514

$`Valor experimental E`
[1] 17.28631

$`Valor teórico`
[1] 2.178813

$`Valor experimental ANOVA`
[1] 158.3007

$`Valor teórico`
[1] 3.885294

(Caso C)

$XtX
    cte   X   Z
cte  15   0   0
X     0 280   0
Z     0   0 100

$`inversa XtX`
           cte           X    Z
cte 0.06666667 0.000000000 0.00
X   0.00000000 0.003571429 0.00
Z   0.00000000 0.000000000 0.01

$Xty
       y
cte   41
X    544
Z   -299

$beta
            y
cte  2.733333
X    1.942857
Z   -2.990000

$`Desviaciones estimadas`
[1] 0.4216966 0.0976038

$`Intervalo CTE`
[1] 1.814535 3.652131

$`Intervalo E`
[1] 1.730197 2.155518

$`Valor experimental CTE`
[1] 6.481753

$`Valor experimental E`
[1] 19.90555

$`Valor teórico`
[1] 2.178813

$`Valor experimental ANOVA`
[1] 365.6949

$`Valor teórico`
[1] 3.885294

Ejercicio 7

(Caso A)

$`Varianza de la perturbación aleatoria`
         [,1]
[1,] 3.105644

$`inversa XtX`
           cte            V            Z
cte 0.06666667 0.0000000000 0.0000000000
V   0.00000000 0.0035924233 0.0004898759
Z   0.00000000 0.0004898759 0.0114304376

$`Estiamción coeficientes (beta)`
         [,1]
cte  1.800000
V    1.918517
Z   -2.817929

$`R (a)`
     [,1] [,2] [,3]
[1,]    0    1   -1

$`r (a)`
[1] 0

$`Fexp (a)`
[1] 514.3874

$`F teórico (a)`
[1] 4.747225

$`R (b)`
     [,1] [,2] [,3]
[1,]    1    1    0

$`r (b)`
[1] 1

$`Fexp (b)`
[1] 33.86959

$`F teórico (b)`
[1] 4.747225

$`R (c)`
     [,1] [,2] [,3]
[1,]    1    2  0.0
[2,]    0    1 -0.5

$`r (c)`
     [,1]
[1,]    0
[2,]    2

$`Fexp (c)`
[1] 85.24315

$`F teórico (c)`
[1] 3.885294

(Caso B)

$`Varianza de la perturbación aleatoria`
         [,1]
[1,] 4.217063

$`inversa XtX`
           cte           X          Z
cte 0.06666667 0.000000000 0.00000000
X   0.00000000 0.003571429 0.00000000
Z   0.00000000 0.000000000 0.08333333

$`Estiamción coeficientes (beta)`
         [,1]
cte  3.866667
X    2.121429
Z   -2.500000

$`R (a)`
     [,1] [,2] [,3]
[1,]    0    1   -1

$`r (a)`
[1] 0

$`Fexp (a)`
[1] 58.27721

$`F teórico (a)`
[1] 4.747225

$`R (b)`
     [,1] [,2] [,3]
[1,]    1    1    0

$`r (b)`
[1] 1

$`Fexp (b)`
[1] 84.00142

$`F teórico (b)`
[1] 4.747225

$`R (c)`
     [,1] [,2] [,3]
[1,]    1    2  0.0
[2,]    0    1 -0.5

$`r (c)`
     [,1]
[1,]    0
[2,]    2

$`Fexp (c)`
[1] 98.46635

$`F teórico (c)`
[1] 3.885294

(Caso C)

$`Varianza de la perturbación aleatoria`
         [,1]
[1,] 2.667421

$`inversa XtX`
           cte           X    Z
cte 0.06666667 0.000000000 0.00
X   0.00000000 0.003571429 0.00
Z   0.00000000 0.000000000 0.01

$`Estiamción coeficientes (beta)`
         [,1]
cte  2.733333
X    1.942857
Z   -2.990000

$`R (a)`
     [,1] [,2] [,3]
[1,]    0    1   -1

$`r (a)`
[1] 0

$`Fexp (a)`
[1] 672.1714

$`F teórico (a)`
[1] 4.747225

$`R (b)`
     [,1] [,2] [,3]
[1,]    1    1    0

$`r (b)`
[1] 1

$`Fexp (b)`
[1] 72.13263

$`F teórico (b)`
[1] 4.747225

$`R (c)`
     [,1] [,2] [,3]
[1,]    1    2  0.0
[2,]    0    1 -0.5

$`r (c)`
     [,1]
[1,]    0
[2,]    2

$`Fexp (c)`
[1] 126.5615

$`F teórico (c)`
[1] 3.885294