10 votos

auto.arima no reconocer un patrón estacional

Tengo un día el clima conjunto de datos, que tiene, como era de esperar, muy fuerte efecto estacional.

enter image description here

He adaptado un modelo ARIMA para este conjunto de datos utilizando la función de auto.arima de previsión de paquete. Para mi sorpresa, la función no se aplica ninguna de temporada de las operaciones de diferenciación estacional estacional ar o ma componentes. Aquí está el modelo es estimado:

library(forecast)
data<-ts(data,frequency=365)
auto.arima(Berlin)

Series: data
ARIMA(3,0,1) with non-zero mean 

Coefficients:
         ar1      ar2     ar3      ma1  intercept
      1.7722  -0.9166  0.1412  -0.8487   283.0378
s.e.  0.0260   0.0326  0.0177   0.0214     1.7990

sigma^2 estimated as 5.56:  log likelihood=-8313.74
AIC=16639.49   AICc=16639.51   BIC=16676.7

Y también las previsiones de uso de este modelo no son realmente satisfactorio. Aquí es la trama de la previsión: enter image description here

Puede alguien darme una pista de lo que está mal aquí?

13voto

Senseful Puntos 116

R no se ajuste a un modelo ARIMA con la estacionalidad de más de 350. Ver http://robjhyndman.com/hyndsight/longseasonality/ para una discusión de este tema. La solución es el uso de términos de Fourier de la estacionalidad, y ARMA errores para el corto plazo, la dinámica.

7voto

Owen Fraser-Green Puntos 642

La solución a tu problema es como Rob señala es la combinación de los efectos deterministas (semana del año) y los efectos estocásticos (ARIMA estructura), mientras que el aislamiento inusual días y la detección de la posible presencia de uno o más nivel de turnos y/o uno o más locales de tendencias en el tiempo. AUTOBOX , el software utilizado para el análisis fue en parte desarrollado por mí automáticamente proporcionan una sólida modelado de conjuntos de datos como este.

He puesto los datos en http://www.autobox.com/weather/weather.txt.

La acf de los datos originales enter image description here que conducen a un modelo automático de selección de la forma enter image description hereenter image description hereenter image description here . El modelo estadísticas enter image description here con un residual de la trama de enter image description here La trama de las previsiones para los próximos 60 días, que se presentan aquí enter image description here El Real/Ajuste/Previsión gráfica que se muestra aquí .enter image description here

Podría ser interesante para los demás a seguir el Prof. Hyndaman del consejo y el informe de su modelo final con disgnostic comprobaciones sobre residual de diagnóstico y los parámetros de las pruebas de significación.

Personalmente, estoy incómodo con la sugerencia acerca de la primera se realiza un análisis de fourier (posiblemente/probablemente afectados por anomalías) y, a continuación, haciendo ARIMA en los residuos es inaceptable, ya que no es una solución simultánea conduce a la ecuación 1, pero en lugar de una presunta secuencia. Mi ecuación uso de la semana-de-la-mes y también se incluye un AR(1) y remedios para la inusual puntos de datos.

Todo el software tiene limitaciones y es bueno saber de ellos. De nuevo reitero ¿por qué alguien no tratar de implementar Rob sugerencias y mostrar los resultados completos.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X