2013-08-22 188 views
1

我在尝试使用R软件包dlm来预测月度时间序列。我正在通过小插图阅读,并有几个问题。关于UKgas数据的预测我看到了解r中的dlm软件包

dlmGas <- dlmModPoly() + dlmModSeas(4) 

4是该系列的季刊格式,所以按月系列,我应该用dlmModSeas(12)呢?

而且对于功能

diag(W(dlmGas))[2:3] <- exp(x[1:2]) 
V(dlmGas) <- exp(x[3]) 

什么是这些数字,它们是怎么选择?这对于每月系列来说又会有什么不同?

+0

我不知道这个包,但看看*预测*包,这是非常非常好的。 – Fernando

回答

2

第一个问题:我认为dlmModSeas(12)对于每月数据是正确的。如果你在R中的命令:

dlmModPoly() + dlmModSeas(12) 

,并期待在矩阵W你看,这仍然只是第二和第三斜元素等于零。

我假设diag(W(dlmGas))[2]指趋势部件(dlmModPoly()) 和diag(W(dlmGas))[3]季节性组件的系统噪声方差的系统噪声方差(dlmModSeas(12)

所以,据我了解,增加季节数量并不会增加W的对角线元素。

+0

感谢您的帮助,所以我想除了dlmModSeas(4)之外,我不需要为每年的数据更改任何内容? –

+0

是的,我会这么说。 – DatamineR

+0

如果我的回答对你有任何用处,那么如果你接受了它就会很好。谢谢。 – DatamineR