2016-07-12 93 views
0

嘿,我有一個關於擬合一個arima模型與R的小問題。我有一個大的數據集,我不能刪除多NA值。NA值擬合的Arima模型

因此需要設計一種方法來用缺少值的arima模型化我的數據。 因此,在AIC值的幫助下,我認爲最好的模型是Arima 2,0,2,但如果我嘗試適合這個模型,我會得到結果:空

有沒有人有一個想法,我可以如何工作接着就,隨即?

我的代碼,以適應該模型是

fitted(final.arima) 

我想我知道的一切與遺漏值工作

+0

你能發佈你的數據和代碼嗎?我試着用'usconsumption'數據集作爲例子,它爲我工作。 –

回答

0

,這是我使用的代碼,但我不能給你我使用的數據集.. daten1 < -test

final.aic <- Inf 

final.order <- c(0,0,0) 

for (p in 0:2) for (q in 0:2) for(z in 0:1) { 

    current.aic <- AIC(arima(daten1, order=c(p, z, q))) 

    if (current.aic < final.aic) { 

    final.aic <- current.aic 

    final.order <- c(p,z, q) 

    final.arima <- arima(daten1, order=final.order) 
    } 
} 

final.aic 
final.order 
final.arima 

acf(resid(final.arima)) 

Box.test(resid(final.arima), lag=20, type="Ljung-Box") 


############RMSE und R² testen für ARIMA Model############### 

ffa<-fitted(final.arima) 

daten1.mean<-mean(daten1) 

R1<-sum((ffa-daten1.mean)^2) 

R2<-sum((daten1-ffa)^2) 

Rsquared<-R1/(R1+R2) 

Rsquared