有效的施瓦茨-贝叶斯矩阵 R

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)))