高效Schwartz-Bayes矩陣R

N00B的問題，但是這樣做的更有效的方法是什麼？我只是試圖計算模型中所有36個參數排列的測試統計量。高效Schwartz-Bayes矩陣R

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

感覺很愚蠢做這樣

來源

2016-02-25 slap-a-da-bias

您的訂單是0：5，0：3，0：2？ – MichaelChirico

使用expand.grid得到所有可能的訂單。

然後用apply來表達你內心的感受。

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

來源

2016-02-25 03:37:32 MichaelChirico

非常感謝你 –

@ slap-a-da-bias在你重複做某件事情的任何時候都有'sapply' /'lapply' /'apply'，你很少會出錯。 – MichaelChirico

也絕對細讀這個優秀的問答：http://stackoverflow.com/questions/3505701/r-grouping-functions-sapply-vs-lapply-vs-apply-vs-tapply-vs-by-vs-aggrega – MichaelChirico

高效Schwartz-Bayes矩陣R

回答

相關問題