Estimación de un modelo AR
Vas a estimar el parámetro AR(1), \(\small \phi\), de una de las series simuladas que generaste en el ejercicio anterior. Como en una serie simulada se conocen los parámetros, es una buena forma de entender las rutinas de estimación antes de aplicarlas a datos reales.
Para simulated_data_1, cuyo \(\small \phi\) verdadero es 0.9, imprimirás la estimación de \(\small \phi\). Además, también mostrarás toda la salida que se produce al ajustar una serie temporal, para que puedas hacerte una idea de qué otras pruebas y estadísticas descriptivas ofrece statsmodels.
Este ejercicio forma parte del curso
Análisis de series temporales en Python
Instrucciones del ejercicio
- Importa la clase
ARIMAdel módulostatsmodels.tsa.arima.model. - Crea una instancia de la clase
ARIMAllamadamodusando los datos simuladossimulated_data_1y el orden (p,d,q) del modelo (en este caso, para un AR(1)) esorder=(1,0,0). - Ajusta el modelo
modcon el método.fit()y guarda el resultado en un objeto de resultados llamadores. - Imprime el resumen completo de resultados con el método
.summary(). - Imprime solo la estimación de \(\small \phi\) usando el atributo
.params[1](sin paréntesis).
ejercicio interactivo práctico
Prueba este ejercicio completando este código de ejemplo.
# Import the ARIMA module from statsmodels
from statsmodels.tsa.arima.model import ARIMA
# Fit an AR(1) model to the first simulated data
mod = ARIMA(___, order=___)
res = mod.___
# Print out summary information on the fit
print(res.___)
# Print out the estimate for phi
print("When the true phi=0.9, the estimate of phi is:")
print(res.___)