Я работаю над прогнозированием временных рядов hts для того же фрагмента кода, предоставленного Робом Хиндманом, на уровне 1 (2 узла) и уровне 2 (5,1 узла).Однако значение MAPE, которое я получаю, очень высоко (более 30 для большинства прогнозов узлов)
Я пытался с параметром Tuning, используя модель Arima (для всех методов - сверху вниз, снизу вверх и т. Д.), ПредложенныхРоб: Иерархический временной ряд
full <- aggts(data)
fcfull <- matrix(0, ncol=ncol(full), nrow=15)
for(i in 1:ncol(full))
{
fcfull[,i] <- forecast(Arima(full[,i],
order=fmodels[[i]]$arma[c(1,6,2)],
seasonal=fmodels[[i]]$arma[c(3,7,4)]),
h=15)$mean
}
forecast.full <- combinef(fcfull, nodes=data$nodes)
Мне нужна высокая точность MAPE ок.быть около 10