2017-01-29 8 views
0

R-newbie質問。リストでARMA coeftest出力を得るためのR - ループ

すべてのループARMAモデルでcoeftest関数を計算し、関連する出力をリストに表示する方法を理解することにいくつかの問題があります。

次のコードを最も簡単に調整しますか?

library(lmtest) 

bchain_2012_logreturns=diff(log(prices_2012)) 
bchain_2013_logreturns=diff(log(prices_2013)) 

bchain_logreturns_Arima_coef=list() 
k=1 
for(i in 2012:2013){ 
for(p in 0:1){ 
for(q in 0:1){   
bchain_logreturns_Arima=Arima(get(paste("bchain_",i,"_logreturns",sep="")),order=c(p,0,q)) 
bchain_logreturns_Arima_coef[[k]]=round(transpose(coeftest(bchain_logreturns_Arima)),digit=3) 
} 
} 
k=k+1 
} 

答えて

0

アキュムレータkの仕組みがわかりませんでした。今それはいいです。

library(lmtest) 

bchain_2012_logreturns=diff(log(prices_2012)) 
bchain_2013_logreturns=diff(log(prices_2013)) 

bchain_logreturns_Arima_coef=list() 
k=1 
for(i in 2012:2013){ 
for(p in 0:1){ 
for(q in 0:1){   
bchain_logreturns_Arima=Arima(get(paste("bchain_",i,"_logreturns",sep="")),order=c(p,0,q)) 
bchain_logreturns_Arima_coef[[k]]=round(transpose(coeftest(bchain_logreturns_Arima)),digit=3) 
k=k+1 
} 
} 
} 
関連する問題