[R] The auto.arima function in forecast package

John mi@ojpm @ending from gm@il@com
Fri Aug 10 10:59:45 CEST 2018


Hi,

   I am wondering if I am doing correctly, but my auto.arima usually (if
not always) give me (0,1,0), whatever portion of the series I take. In the
following instances, only the last one yields ARIMA(0,1,1), and all the
other cases yield ARIMA(0,1,0). I would like to do forecast based on ARIMA.
The forecast with ARIMA(0,1,0) is exactly the random walk forecast for any
forecasting horizon h, so it does not make much sense.... Am I doing
anything wrong? Is the data very unusual?

   The length of the df_arima sequence is 68. My data is seasonally
adjusted, so I use seasonal = FALSE.

> auto.arima(df_arima, d=1, ic="aic", seasonal = FALSE)
Series: df_arima
ARIMA(0,1,0)

sigma^2 estimated as 0.9323:  log likelihood=-221.43
AIC=444.85   AICc=444.88   BIC=447.93

> auto.arima(df_arima[2:30], d=1, ic="aic", seasonal = FALSE)
Series: df_arima[2:30]
ARIMA(0,1,0)

sigma^2 estimated as 1.486:  log likelihood=-45.28
AIC=92.56   AICc=92.71   BIC=93.89
> auto.arima(df_arima[30:50], d=1, ic="aic", seasonal = FALSE)
Series: df_arima[30:50]
ARIMA(0,1,0)

sigma^2 estimated as 2.065:  log likelihood=-35.63
AIC=73.26   AICc=73.48   BIC=74.25
> auto.arima(df_arima[30:65], d=1, ic="aic", seasonal = FALSE)
Series: df_arima[30:65]
ARIMA(0,1,1)

Coefficients:
         ma1
      0.3684
s.e.  0.1591

sigma^2 estimated as 1.328:  log likelihood=-54.19
AIC=112.39   AICc=112.76   BIC=115.5


> dput(df_arima)
c(0.206955966502065, 0.572310398166964, 0.730924932354315,
1.41842200551883,
3.46103800972619, 3.61249423115895, 5.82223437508589, 4.16838241173942,
3.41346768155648, 0.862064978610255, 0.998745364906517, 0.519274399474301,
-1.11283353732583, -0.735660435061958, -0.780332562883557,
0.835729609943847,
2.72188993637388, 2.28434816118952, 1.89419627788212, 0.808622045623819,
0.628929267637468, -0.605142332415043, -1.50590253313166,
-0.288872689519148,
0.779724144195981, 1.77156391706856, 1.54853781756041, 3.31419159924959,
1.93558418905826, 1.2733744986462, 1.08821585241274, 2.62853576193935,
3.73925259835317, 4.36605858718784, 4.94030870469473, 4.95025405015856,
3.62293670685028, 2.66083271862916, 2.63088962945459, 2.64079870265661,
2.24308372288402, 4.60965981893589, 5.61665338477875, 4.13854855208327,
2.74923595127379, 1.95704807193253, 1.39700634338462, -3.2051540517546,
-4.14645815871689, -3.55954963127346, -2.75935917560429, -2.35979888878128,
-1.60997264659325, 0.107359485533354, 0.452171464181572, 1.83710430300006,
2.05388942618911, 1.42071194194746, 0.62281728667859, 0.66619920044817,
1.91899026506288, 0.876017396693274, 0.767997179685187, 0.960763098531148,
1.07168884361519, 1.7179971391635, 0.916107067445848, 0.488814956401007,
0.445930327726551, 0.700746383716555, 1.85843822164806, 0.892473264860727,
0.614568090119794, 2.4071141581713, 1.95656784524512, 1.97416030945137,
1.50143287537079, 2.16462319254254, 2.07797672518777, 0.272208427739407,
0.798064671535004, 0.754477049206304, 1.15345347401843, 1.2569270188491,
1.04667488857908, 1.61549969176582, 1.61397592409376, 2.15802281566817,
1.77514416862246, 1.37441120562698, 2.02489513295532, 2.58699935920175,
2.75245197313523, 2.23920076522477, 1.5517277192193, 2.1968256454034,
2.56612088350139, 1.58214013021238, 0.224868329892303, 0, 0.75815516785358,
0.306075421735685, -0.911618439021877, -0.851059960629386,
0.920475843398605,
0.204019124781452, 0.265049070791767, 0.693380201646732, 1.76190014304951,
1.7409678665274, 2.30875574967095, 3.1616985460498, 3.03416937102321,
2.51100317138793, 1.74651451299397, 2.1955027765292, 0.885663954102278,
0.42162278205693, -0.919995110735239, -1.09905745071531, -2.68791519080633,
-3.351439, -2.10834457783508, -2.20921682103322, -0.74375338663687,
-0.402576080931594, 1.20253696682322, 1.49061195023896, 0.526155117506621,
1.48289074384766, 0.627624111314717, 0.262639414167709, 0.101025401449073,
0.667507628973207, 0.86863891914688, 0.603620359927048, 1.53734977063493,
1.63989774985629, 2.06996035293847, 1.23150306020205, 1.00602988533123,
1.46753831297777, 1.28396719372157, 1.22372099999999, 0.42010456383097,
1.68739489046037, 1.0009947634714, 0.218850966640427, 0.218698677177565,
0.457733617597422, 1.42090579656218, 0.436897503947153, 1.05592733175737,
1.67544706250284, 2.69758588203255, 2.87640149597799, 2.13003529065778,
2.42946437023286, 2.32608159085508, 1.97956200691809, 0.274402066259904
)

	[[alternative HTML version deleted]]



More information about the R-help mailing list