#r #for-loop #dplyr #forecasting
#r #for-цикл #dplyr #прогнозирование
Вопрос:
У меня есть один фрейм данных, который состоит из трех столбцов (AJT, NET и SAT). Мое намерение состоит в том, чтобы сделать прогноз для каждой из этих трех временных серий с пакетом прогнозов. По этой причине я преобразую фрейм данных в объект ts и делаю прогнозирование с помощью функции snaive, поэтому я написал эти строки кода:
#CODE
library(forecast)
# Data set
DATA_SET<-data.frame(
AJT=seq(1, 48, by = 2),
NET=seq(1, 24, by = 1),
SAT=seq(1, 94, by = 4)
)
# Making TS object
TS_SALES<-ts(DATA_SET,start=c(2016,1),frequency = 12)
# Making forecasting with Forecast package
SNAIVE_AJT<-snaive(TS_SALES[, 'AJT'],h=5)
SNAIVE_NET<-snaive(TS_SALES[, 'NET'],h=5)
SNAIVE_SAT<-snaive(TS_SALES[, 'SAT'],h=5)
# Union forecast in list
SNAIVE_UNION<-mapply(SNAIVE_AJT, SNAIVE_NET,SNAIVE_SAT, FUN=list, SIMPLIFY=FALSE)
Все выходные данные из функции snaive я помещаю в SNAIVE_UNION, которые содержат все результаты прогнозирования. Итак, здесь наиболее важным компонентом является «среднее значение», которое содержит результаты прогнозирования по месяцам.
SNAIVE_UNION[["mean"]]
#
# # [[1]]
# # Jan Feb Mar Apr May
# # 2018 25 27 29 31 33
# #
# # [[2]]
# # Jan Feb Mar Apr May
# # 2018 13 14 15 16 17
# #
# # [[3]]
# # Jan Feb Mar Apr May
# # 2018 49 53 57 61 65
Итак, здесь мое намерение состоит в том, чтобы поместить результаты формы SNAIVE_UNION[[«mean»]] в таблицу данных, подобную таблице ниже, с некоторым функциональным циклом, для или другой функции
Jan Feb Mar Apr May
------------------------------
AJT 25 27 29 31 33
NET 13 14 15 16 17
SAT 49 53 57 61 65
Я спрашиваю об этом, потому что этот временной ряд является лишь небольшой частью целого ряда, и я хотел бы автоматизировать этот код.
Ответ №1:
Вы можете использовать:
pred <- do.call(rbind, SNAIVE_UNION[["mean"]])
# [,1] [,2] [,3] [,4] [,5]
# [1,] 25 27 29 31 33
# [2,] 13 14 15 16 17
# [3,] 49 53 57 61 65