高效的 Schwartz-Bayes 矩阵 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)))

高效的 Schwartz-Bayes 矩阵 R

efficient Schwartz-Bayes matrix R

statistics

r

matrix

autoregressive-models