Извлечение AIC из ARIMA с использованием вложенного цикла

Question

Я пытаюсь извлечь AIC и BIC из оценки ARIMA с различными комбинациями p & q (p = 0,1,2,3 и q = 0,1.2,3). Я попытался использовать следующий код, безуспешно.

code: storage1 <- numeric(16) for (p in 0:3){

>     for (q in 0:3){
>  
>     storage1[p]  <- arima(x,order=c(p,0,q), method="ML")} } storage1$aic

Answer 1

Один из способов - использовать expand.grid(), чтобы получить все комбинации p и q и применить их к вашей модели арима, т.е.

apply(expand.grid(p, q), 1, function(i) arima(d1$cnt, order = c(i[1], 0, i[2]), method = "ML")$aic)

#[1] 47222.43 38589.14 36935.33 36118.44 42569.73 35183.35 35141.13 35143.02 39448.38 35142.96 35142.58 35142.60