适合具有NA值的Arima模型

时间:2016-07-12 13:10:48

标签: r na

嘿,我有一个关于用R拟合arima模型的简短问题。我有一个包含多个NA值的大型数据集,我无法删除。

因此需要找到一种方法来使用具有缺失值的arima来建模我的数据。 因此,在AIC Value的帮助下,我认为最好的模型是Arima 2,0,2,但如果我尝试适合这个模型,我会得到结果:Null

有没有人知道如何使用它?

我的代码适合模型

{{1}}

我尝试了所有我想知道的缺失值

1 个答案:

答案 0 :(得分: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