Estimativa de um modelo AR
Você estimará o parâmetro AR(1), \(\small \phi\), de uma das séries simuladas que você gerou no exercício anterior. Como os parâmetros são conhecidos para uma série simulada, essa é uma boa maneira de entender as rotinas de estimativa antes de aplicá-las a dados reais.
Para simulated_data_1
com um verdadeiro \(\small \phi\) de 0,9, você imprimirá a estimativa de \(\small \phi\). Além disso, você também imprimirá toda a saída produzida ao ajustar uma série temporal, para que possa ter uma ideia de quais outros testes e estatísticas resumidas estão disponíveis no statsmodels.
Este exercício faz parte do curso
Análise de séries temporais em Python
Instruções de exercício
- Importe a classe
ARIMA
no módulostatsmodels.tsa.arima.model
. - Crie uma instância da classe
ARIMA
chamadamod
usando os dados simuladossimulated_data_1
e a ordem (p,d,q) do modelo (nesse caso, para um AR(1)), éorder=(1,0,0)
. - Ajuste o modelo
mod
usando o método.fit()
e salve-o em um objeto de resultados chamadores
. - Imprima o resumo completo dos resultados usando o método
.summary()
. - Basta imprimir uma estimativa de \(\small \phi\) usando o atributo
.params[1]
(sem parênteses).
Exercício interativo prático
Experimente este exercício preenchendo este código de exemplo.
# 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.___)