R печать сводки прогноза с помощью auto.arima через список фреймов данных с помощью функции - PullRequest
Новая
       11

R печать сводки прогноза с помощью auto.arima через список фреймов данных с помощью функции

0 голосов
/ 17 февраля 2020

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

В настоящее время я построил модель, которая делает использование функции l oop через список фреймов данных для прогнозирования на 12 месяцев в будущем. Я доволен результатами, которые дает сценарий, однако я хотел бы также напечатать для каждого фрейма данных в списке сводку функции auto.arima, чтобы я мог видеть, какой метод прогнозирования был применен. Моя функция в том виде, в каком она есть, сопровождается кодом, который позволяет l oop просматривать список фреймов данных: 

 myFunctionForecast <- function(x) {

  results <- ts(x$Data, frequency = 12)

  unclass(x$Data)

  fit <- auto.arima(results, D = 1)

  fcast <- forecast(fit, h = 12)

  summary <- print(summary(fcast))

}

dataforecast_Data <- lapply(dat6split_Data, myFunctionForecast)

результаты выдаются следующим образом для каждого фрейма данных в списке ( Результаты также включают в себя остальные hi95, et c et c просто хотел сохранить ниже приборку):

  • Прогнозы:
  • Прогноз точки Lo 80
  • 4 ноября 13,88399 0,9072741
  • De c 4 15,41793 1,3156040
  • 5 января 14,31349 -1,5984508
  • 5 февраля 26,60206 9,8396039
  • 5 марта * 20.22581 2.7063289
  • 5 апреля 29.59404 11.5765060
  • 5 мая 12.96957 -5.4363124
  • 5 июня 21.36188 2.6750927
  • 5 июля 13.70204 -5.1976049
  • 5 сентября 16.13085 -3.0467465
  • 5 октября 32.03698 12.7695513

Однако то, что я также хотел бы получить для каждого печатного резюме кадр данных, как показано ниже:

Метод прогноза: ARIMA (2,0,0) (1,1,0) [12] * 104 3 *

Информация о модели: Серия: результаты ARIMA (2,0,0) (1,1,0) [12]

Коэффициенты: ar1 ar2 sar1 0,4254 0,3869 -0,5899 se 0,1591 0,1692 0,1452

сигма ^ 2 оценивается как 102,5: логарифмическая правдоподобность = -128,34 AIC = 264,68
AICc = 266,06 BIC = 270,78

Оценка ошибок: ME RMSE Training set -1.358891 8.312478 MAE MPE Training набор 5.605761 -13.99908 MAPE MASE Учебный комплект 24.16324 0.4216723 ACF1 Учебный комплект 0.02919079

В таком виде я могу заставить его напечатать сводку на консоль, но поскольку это автоматизированный процесс, мне нужен способ сохранения этих сводок, чтобы их можно было экспортировать вместе с данными обратно на наш SQL сервер. Любая помощь по этим вопросам будет принята с благодарностью:).

Приветствия.

...