在 R 版本 3.2.0 和 R 版本 3.2.2 中 d=1 的 auto.arima 的不同结果答案

【问题标题】：Different results for auto.arima with d=1 in R version 3.2.0 and R version 3.2.2在 R 版本 3.2.0 和 R 版本 3.2.2 中 d=1 的 auto.arima 的不同结果
【发布时间】：2015-08-26 18:30:59
【问题描述】：

我在一台计算机上安装了 R 版本 3.2.0，在另一台计算机上安装了 3.2.2 版本，并且它们都安装了相同版本的包“预测”(6.1)。当我将带有d=1 的auto.arima 函数应用于3.2.0 版计算机中的向量fec1（以下数据）时，它给了我一个ARIMA(1,1,1)。在版本为 3.2.2 的计算机中，它给了我一个带有漂移的 ARIMA(1,1,0)。为什么会发生这样的差异？当我省略选项 d=1 时，我在两台计算机上得到相同的结果。

fec1 <- c(9.53408021, 9.48993631, 9.44261613, 9.38585482, 9.31611858, 
          9.24743801, 9.19483663, 9.11690414, 9.05586032, 8.96781034, 
          8.85437896, 8.73956021, 8.62587246, 8.49758284, 8.41877774, 
          8.29737132, 8.02294327, 7.62604834, 7.32936532, 7.10006394, 
          6.87039464, 6.56332757, 6.27195705, 5.97089310, 5.69537808, 
          5.46074967, 5.19061116, 4.90749261, 4.61822473, 4.39635507, 
          3.96883450, 3.71230296, 3.43604354, 3.15150833, 2.82864081, 
          2.51857560, 2.27264470, 2.06785600, 1.86116599, 1.61452534, 
          1.41214847, 1.26673153, 1.15033592, 1.00646045, 0.83493055, 
          0.65399085, 0.47450499, 0.30408840, 0.18216540, 0.09006937, 
          0.00000000)

【问题讨论】：

标签： r time-series forecasting

【解决方案1】：

请参阅PR#16278，“嵌套的 arima 模型具有更高的对数似然”，如 change log for 3.2.1 中所述，会影响 d >= 1 的 arima 计算：

arima(*, xreg = .)（对于 d >= 1）计算估计方差基于有效观察的数量，如 R 版本 3.0.1 及更早版本。 (PR#16278)

【讨论】：