Question

N00B问题，但有什么更有效的方法呢？我只是试图计算模型中所有36个参数排列的测试统计数据。

ARIMAP0Q<-c(BIC(arima(italy.1$gap,order=c(0,0,0))),
            BIC(arima(italy.1$gap,order=c(1,0,0))),
            BIC(arima(italy.1$gap,order=c(2,0,0))),
            BIC(arima(italy.1$gap,order=c(3,0,0))),
            BIC(arima(italy.1$gap,order=c(4,0,0))),
            BIC(arima(italy.1$gap,order=c(5,0,0))),
            BIC(arima(italy.1$gap,order=c(0,0,1))),
            BIC(arima(italy.1$gap,order=c(1,0,1))),
            BIC(arima(italy.1$gap,order=c(2,0,1))),

这样做感觉非常愚蠢

Answer 1

使用expand.grid获取所有可能的订单。

然后将apply用于您心中的内容。

ARIMAP0Q <- apply(expand.grid(0:5, 0:2, 0:3), 1L,
                  function(rw) BIC(arima(italy.1$gap, order = rw)))

有效的Schwartz-Bayes矩阵R

1 个答案: