Quiero saber el resultado del error como RMSE, etc. de ARIMA. Tengo 45 meses de datos sobre las luces nocturnas. tengo un modelo ARIMA como este

fitARIMA <- arima(newdata, order=c(0,0,0),seasonal = list(order = c(1,0,0), period = 12))
summary(fitARIMA)
### here is the result
Call:
arima(x = newdata, order = c(0, 0, 0), seasonal = list(order = c(1, 0, 0), period = 12))
            
Coefficients:
       sar1  intercept
      0.4770   572.1038
s.e.  0.1608    38.5140
            
sigma^2 estimated as 26880:  log likelihood = -294.88,  aic = 593.76
            
Training set error measures:
              ME RMSE MAE MPE MAPE
Training set NaN  NaN NaN NaN  NaN
Warning message:
In trainingaccuracy(object, test, d, D) :
test elements must be within sample

¿Alguien sabe por qué sucedió esto y cómo resolver este problema? gracias

Aquí están los datos que utilicé

         Jan     Feb     Mar     Apr     May     Jun     Jul
2015  467.38  441.67  579.30  600.41  793.38  576.80  741.21
2016  516.02  241.41  443.20  502.98  497.31  668.08  596.89
2017  325.89  253.30  737.37  462.75  609.31  559.05  581.16
2018  428.74  584.53  508.92  655.63  867.83 1059.98  509.34
         Aug     Sep     Oct     Nov     Dec
2015  634.66  582.00  661.35  249.46  482.33
2016  686.76  598.28  598.23  391.71  492.66
2017  680.36  753.18  476.41    3.12  608.01
2018  820.85  825.13

                    
0
Ahmad Rizqullah 21 ene. 2021 a las 21:40

1 respuesta

La mejor respuesta

Reemplaza arima con Arima. La función Arima() del paquete de pronóstico guarda información adicional que es utilizada por la función summary().

Suponiendo que su objeto newdata es de la clase ts con una frecuencia establecida en 12, el siguiente código más simple debería funcionar.

fitARIMA <- Arima(newdata, order=c(0,0,0), seasonal = c(1,0,0))
summary(fitARIMA)
0
Rob Hyndman 21 ene. 2021 a las 23:30