10.5 Ejercicios
1. Estudio de simulación: especificación incorrecta de la forma funcional
Como se indicó en el Capítulo 10.2, la especificación incorrecta de la función de regresión viola el supuesto 1 del Concepto clave 6.3, por lo que el estimador MCO estará sesgado e inconsistente. Se ha ilustrado el sesgo de para el ejemplo de la función de regresión poblacional cuadrática
y el modelo lineal usando 100 observaciones generadas aleatoriamente. Estrictamente hablando, este hallazgo podría ser solo una coincidencia porque se considera solo una estimación obtenida usando un solo conjunto de datos.
En este ejercicio, se debe generar evidencia de simulación para el sesgo de en el modelo si la función de regresión de la población es
Instrucciones:
Asegúrese de utilizar las definiciones sugeridas en el código esqueleto en script.R para completar las siguientes tareas:
Genere 1000 estimaciones de MCO de en el modelo anterior utilizando un bucle for() donde , usando muestras de tamaño . Guardar las estimaciones en beta_hats.
Comparar la media muestral de las estimaciones con el parámetro verdadero utilizando el operador ==.
Sugerencia:
- Puede generar números aleatorios a partir de una distribución uniforme utilizando runif().
2. Estudio de simulación: Sesgo de errores en variables
Considere nuevamente la aplicación del modelo clásico de error de medición presentado en el Capítulo 10.2:
El regresor único se mide con error de modo que en su lugar se observa . Por tanto, se estima en
en lugar de
con el error medio cero no correlacionado con y . Entonces es estimado de manera inconsistente por MCO:
Deje
Recuerde de (10.2) que en este caso. Además, suponga que con .
Como se mencionó en el ejercicio 1, el capítulo 10.2 analiza las consecuencias del error de medición para el estimador de MCO de en este entorno basado en una única muestra y, por lo tanto, una sola estimación. Estrictamente hablando, la conclusión obtenida podría ser incorrecta porque el sesgo observado puede deberse a una variación aleatoria. Una simulación de Monto Carlo es más apropiada aquí.
Instrucciones:
Muestre que se estima con un sesgo utilizando un estudio de simulación. Asegúrese de utilizar las definiciones sugeridas en el código esqueleto en script.R para completar las siguientes tareas:
Generar 1000 estimaciones de en el modelo de regresión simple Usar rmvnorm() para generar muestras de 100 observaciones aleatorias de la distribución normal bivariada indicada anteriormente.
Guardar las estimaciones en beta_hats.
Calcular la media muestral de las estimaciones.