我想生成一个矩阵的预测.. 例如
> x<-ts(...,start=c(2009,1),frequency=12)
如此。
A B . . . G
2009/1 2314 1456 3021
2009/2 2345 2351 2241
. . . .
. . . .
. . . .
2014/12 3210 3412 3199
获取第一列x[,1]
第二列x[,2]
..
现在我想用一个订单预测所有产品(A,B,C,....,G)。
预测方法是holtswinter。
所以我猜它会是这样的:
for(i in 1:7)
fit1[i]<-hw(x[,i],h=12,model="additive")
fit2[i]<-hw(x[,i],h=12,model="multiplicative")
答案 0 :(得分:0)
最简单的方法可能是将predict()
与HoltWinters()
一起使用,与您提出的其他question非常相似:
dat <- data.frame( a =rnorm(100), b=rnorm(100), c=rnorm(100))
myTs <- ts(dat, start=2000, frequency=12)
sapply(colnames(myTs), function(i,...) predict(HoltWinters(myTs[,i],...)), seasonal = "additive" )
a b c
-0.49806413 -0.35461294 0.02780248
或输出为list
:
lapply(myTs, function(x,...) predict(HoltWinters(x,...)), seasonal = "additive")
$a
May
2008 -0.4980641
$b
May
2008 -0.3546129
$c
May
2008 0.02780248