R中的精度函数

时间:2015-02-07 14:31:48

标签: r time-series forecasting

我想计算包含12个预测的数值向量的准确性,但我得到了这个结果并发出警告。

accuracy(f,Test)
                    ME      RMSE       MAE       MPE      MAPE ACF1 Theil's U
Training set       NaN       NaN       NaN       NaN       NaN   NA        NA
Test set     0.9064933 0.9064933 0.9064933 0.4060658 0.4060658   NA        NA
Warning message:
In trainingaccuracy(f, test, d, D) : test elements must be within sample`

Test : ` Jan       Feb       Mar       Apr       May       Jun       Jul
2012 223.23801 203.51890 195.21139 169.68900 157.63439 135.79616  96.84895
           Aug       Sep       Oct       Nov       Dec
2012 100.42089  62.51129  49.18672  22.10280   5.00000

我通过以下代码获得f

f  <- list()

for (i in 1:12) { 
  demfit <- ets(Dem2) 
  f[[i]] <- forecast(demfit, 1)  #start filling the f list with values
  Dem2[length(Dem2)+1] <- f[[1]]$mean[1]
} 

f <- unlist(f)

此代码旨在提前一步预测12次(每次更新我的系列预测值)

感谢您的帮助

我的数据:

     Jan    Feb    Mar    Apr    May    Jun    Jul 
2008 514.43 508.07 505.13 486.00 482.20 470.37 447.30 
2009 394.93 381.19 381.94 363.38 361.50 352.96 334.86 
2010 336.10 320.44 324.10 311.16 303.78 305.80 285.19 
2011 260.85 258.78 271.43 271.43 267.46 271.77 254.83 
2012 223.23 203.51 195.21 169.68 157.63 135.79 96.84 

     Aug    Sep    Oct    Nov    Dec 
2008 438.83 415.91 408.44 409.67 395.52 
2009 340.41 326.41 331.14330 337.07 336.51 
2010 297.95 274.93 272.91 268.58 271.35 
2011 266.91 242.42 241.04 230.29 233.57 
2012 100.42 62.51 49.18 22.10 5

1 个答案:

答案 0 :(得分:0)

查看accuracy功能文档,它说明f

  

“预测”类的对象,或包含的数字向量   预测。如果是x,它也可以与Arima,ets和lm对象一起使用   省略 - 在这种情况下返回样本内准确度度量。

虽然您提供的是此类对象的列表 - 但不一样。