Функция Forecast:: precision () возвращает "Не поддерживается" и NA

Я сделал прогноз на 365 шагов вперед, используя пакет Forecast() и следующий код:

##The time series is 3650 daily observations of rainfall
x <- ts(x$obs, start=c(2007, 10), end=c(2017, 9), frequency = 365)

##create training set - first 9 years of observations
x_train <- subset(x, start = 1, end = 3285)

##test set - last year of observations
x_test <- subset(x, start = 3286, end = 3650)

##fit the model
x_train_fit <- auto.arima(x_train, seasonal=FALSE, xreg=fourier(x_train, K=1))

##forecast using the model 
x_fcast_test <- forecast(x_train_fit,h=365, xreg=fourier(x_test, K=1))
plot(x_fcast_test, col="black") 
lines(x_test,col="red")

Проблема: Когда я запускаю функцию precision () для проверки точности этого прогноза, выходные данные выглядят так:

> accuracy(x_fcast_test)
$`Models`
   Call           
1  "Not supported"
2  "Not supported"
3  "Not supported"
4  "Not supported"
5  "Not supported"
6  "Not supported"
7  "Not supported"
8  "Not supported"
9  "Not supported"
10 "Not supported"

$Fit.criteria
   Min.max.accuracy MAE MAPE MSE RMSE NRMSE.mean NRMSE.median NRMSE.mean.accuracy NRMSE.median.accuracy Efron.r.squared
1                NA  NA   NA  NA   NA         NA           NA                  NA                    NA              NA
2                NA  NA   NA  NA   NA         NA           NA                  NA                    NA              NA
3                NA  NA   NA  NA   NA         NA           NA                  NA                    NA              NA
4                NA  NA   NA  NA   NA         NA           NA                  NA                    NA              NA
5                NA  NA   NA  NA   NA         NA           NA                  NA                    NA              NA
6                NA  NA   NA  NA   NA         NA           NA                  NA                    NA              NA
7                NA  NA   NA  NA   NA         NA           NA                  NA                    NA              NA
8                NA  NA   NA  NA   NA         NA           NA                  NA                    NA              NA
9                NA  NA   NA  NA   NA         NA           NA                  NA                    NA              NA
10               NA  NA   NA  NA   NA         NA           NA                  NA                    NA              NA

Точность () указывает на то, что мне нужно вводить объект прогноза класса и, необязательно, числовой вектор или временной ряд той же длины для сравнения.

> class(x_fcast_test)
[1] "forecast"
> class(x_test)
[1] "ts"

Вот некоторые варианты, которые я попробовал, чтобы убедиться, что я достигаю этого. Ничего из этого тоже не работает.

accuracy(x_fcast_test, x_test)
accuracy(x_fcast_test$mean, x_test)
accuracy(x_fcast_test$mean, as.numeric(x_test))
accuracy(x_fcast_test, as.numeric(x_test))

Как я должен интерпретировать эти результаты? Что может вызвать этот тип ошибки?

Источник

1 ответ

Решение

Похоже, что вы загрузили пакет, который перезаписывает accuracy() функция в пакете прогноза.

Источник
Другие вопросы по тегам