Используйте Croston модель R HTS пакет

Question

Я думаю, что это похоже на пост this post passing forecasting method to hts

пример, приведенный здесь для пакета MAPE. Можно ли использовать функцию croston в пакете прогнозов аналогичным образом?

Я просто попытался

all_ts <- aggts(bts) 
    allf <- matrix(NA, nrow = 3, ncol = ncol(all_ts)) 
    for(i in 1:ncol(all_ts)){ 
     allf[,i] <- croston(all_ts[,i],h = 3) 
    } 

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

Answer 1

Посмотрите на информацию об объекте, возвращаемом croston. Это не простой вектор. Этот пример является примером, который работает.

library(hts) 
nodes <- list(2, c(3, 2)) 
abc <- ts(5 + matrix(sort(rnorm(500)), ncol = 5, nrow = 100)) 
bts <- hts(abc, nodes) 
all_ts <- aggts(bts) 
allf <- matrix(NA, nrow = 3, ncol = ncol(all_ts)) 
for(i in 1:ncol(all_ts)){ 
    allf[,i] <- croston(all_ts[,i],h = 3)$mean 
} 
y.f <- combinef(allf, bts$nodes)