[R] Warning messages after auto.arima
Chua Siang Li
siang.li.chua at acceval-intl.com
Tue Sep 16 11:42:27 CEST 2008
Dear R-helpers.
Would appreciate if someone can explain the warning messages below, after
auto.arima. I couldn't find any clue in the archived help.
Also, how do I retrieve the AICs of each tried model in auto.arima? The
purposes are (1) to output to a text file, and (2) to find the 2nd best
model by finding 2nd lowest AIC instead of eyeballing thru the value at the
console screen.
For those models with drift, how can I find out what are the drifts? For
example, eyeballing thru, my 2nd best model is ARIMA(1,1,1) with drift but
it didn't state what is the drift.
Many thanks.
> modelarima <- auto.arima(Price[1:nTrain], trace=TRUE)
ARIMA(2,1,2) with drift : 4417.541
ARIMA(0,1,0) with drift : 4538.817
ARIMA(1,1,0) with drift : 4424.534
ARIMA(0,1,1) with drift : 4457.507
ARIMA(1,1,2) with drift : 4416.070
ARIMA(1,1,1) with drift : 4414.328
ARIMA(1,1,1) : 4413.551
ARIMA(0,1,1) : 4458.788
ARIMA(2,1,1) : 4415.887
ARIMA(1,1,0) : 4424.668
ARIMA(1,1,2) : 4415.277
ARIMA(0,1,0) : 4542.161
ARIMA(2,1,2) : 4416.889
Best model: ARIMA(1,1,1)
Warning message:
In if (class(fit) != "try-error") offset <- -2 * fit$loglik - length(x) * :
the condition has length > 1 and only the first element will be used
----
Chua Siang Li
Consultant - Operations Research
Acceval Pte Ltd
Tel: 6297 8740
Email: siang.li.chua at acceval-intl.com
Website: www.acceval-intl.com
This message and any attachments (the "message"...{{dropped:13}}
More information about the R-help
mailing list