2016-07-12 10 views
0

ちょっとRにArimaモデルをフィッティングするのに簡単な質問があります。私は削除できない複数のNA値を持つ大きなデータセットを持っています。NA値を持つArimaモデル

このような値のないarimaでデータをモデル化する方法を考え出す必要があります。 AIC値の助けを借りて、私は最高のモデルがArima 2,0,2であると考えましたが、このモデルに合うようにしようとすると結果が得られます:Null

私はどのように働くことができますか?それと?

モデルに合うように私のコードは

fitted(final.arima) 

である私は私が私が使用するコード厥

+0

データとコードを投稿できますか?私は 'usconsumption'データセットを例として使ってこれを試してみました。 –

答えて

0

に欠損値で動作するように知っているが、私はあなたに私が使用してデータセットを与えることができないすべてのものを試してみました。.. 。 daten1 <テスト

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 
関連する問題