階層予測（htsパッケージ）R不均一なグループとカスタム予測

Question

私はhtsパッケージのために私を読んでいます。私は、カスタム予測セクションを試してみたいと思っています。カスタム予測セクションは、各シリーズを特定の予測関数でループする必要があります。提供された例は、偶数のグループを示しています。グループの数が不均一なカスタム予測をどのように実行してから、それらをgtsオブジェクトに再結合できるかを知りたいと思います。 hts()にnodes引数で指定された階層構造と一致していませんg

bts <- ts(5 + matrix(sort(rnorm(500)), ncol=5, nrow=100)) 
y <- hts(bts, nodes=list(2, c(3, 2))) 

allts_y <- aggts(y) 
allf <- matrix(allts_y, nrow=100, ncol=ncol(allts_y)) 

for(i in 1:ncol(allts_y)) 
    allf[,i] <- forecast(auto.arima(allts_y[,i]), h=10)$mean 
allf <- ts(allf, start=1) 

# below code is where I run into a hang-up 
g <- matrix(c(rep(2, 5), rep(3, 5), rep(1:5, 2)), nrow = 2, byrow = T) 
y.f <- combinef(allf, groups = g) 

Answer 1

あなたのグループ構造。最後の行を置き換えてください

y.f <- combinef(allf, nodes=y$nodes)