我在R Studio 3.3.2上使用R Studio上的预测包。我正在用黄土方法(stl)预测时间序列。问题是当我在R版本3.4.1上运行相同的代码时,结果与我在3.3.2版本上得到的结果非常不同。以下是我面临问题的代码部分。
wcfit <- stl(wcts, t.window=4, s.window="periodic", robust=TRUE)
fcast <- forecast(wcfit,h=8, method="naive")
fcast值在两种情况下都不同。请帮忙
以下是两种情况下相同时间序列的样本输出:
R3.3.2
Point Forecast Lo 80 Hi 80 Lo 95 Hi 95
2017 Q1 3528.904 927.43592 6130.372 -449.6980 7507.506
2017 Q2 2619.054 17.58643 5220.522 -1359.5475 6597.656
2017 Q3 3393.551 792.08334 5995.019 -585.0506 7372.153
2017 Q4 647.000 -1954.46786 3248.468 -3331.6018 4625.602
2018 Q1 3528.904 -150.12735 7207.935 -2097.6888 9155.496
2018 Q2 2619.054 -1059.97684 6298.085 -3007.5383 8245.647
2018 Q3 3393.551 -285.47993 7072.582 -2233.0414 9020.144
2018 Q4 647.000 -3032.03113 4326.031 -4979.5926 6273.593
R3.4.1
Point Forecast Lo 80 Hi 80 Lo 95 Hi 95
2017 Q1 3528.904 927.4359 6130.372 -449.698 7507.506
2017 Q2 2619.054 -1059.9768 6298.085 -3007.538 8245.647
2017 Q3 3393.551 -1112.3233 7899.426 -3497.589 10284.692
2017 Q4 647.000 -4555.9357 5849.936 -7310.204 8604.204
2018 Q1 3528.904 -2288.1552 9345.963 -5367.520 12425.328
2018 Q2 2619.054 -3753.2145 8991.323 -7126.490 12364.598
2018 Q3 3393.551 -3489.2858 10276.388 -7132.840 13919.942
2018 Q4 647.000 -6711.0623 8005.062 -10606.185 11900.185
我也意识到版本R3.2.1给出与3.4.1相同的结果
谢谢!