r - 使用嵌套循环从 ARIMA 中提取 AIC

Question

我正在尝试从具有不同的 ARIMA 估计中提取 AIC 和 BIC
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 ，并将它们应用于您的 arima 模型，即

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