に欠損値で動作するように知っているが、私はあなたに私が使用してデータセットを与えることができないすべてのものを試してみました。.. 。 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
データとコードを投稿できますか?私は 'usconsumption'データセットを例として使ってこれを試してみました。 –