几个通用的cmets:
- 您可以使用基本 R 的
arima 函数拟合 SARIMA 模型;真的不需要加载额外的包来安装基本的 SARIMA 模型。
-
AirPassengers 数据集是基础 R 的一部分,而不是 astsa 的一部分。
在arima 中,您分别通过函数参数order 和seasonal 指定非季节性和季节性参数。如果您在 R 终端中键入 ?arima,您可以找到该函数接受哪些参数的详细信息。
例如,要指定 SARIMA(0,1,1)(0,1,1)12 模型
arima(log(AirPassengers), order = c(0,1,1), seasonal = list(order = c(0,1,1), period = 12))
请注意,seasonal 参数是 list,其中包含两个元素 order(表示季节性 AR 项的顺序、季节性差异的顺序和季节性 MA 项的顺序)和 period。
您在astsa::sarima 中看到的错误似乎源于太多参数;检查您是否真的有 3 个非季节性参数、3 个(可选)季节性参数和一个(可选)周期。
上面的模型是
sarima(log(AirPassengers), 0, 1, 1, 0, 1, 1, 12)
它产生与arima 相同的拟合结果(实际上sarima 只是调用arima“在幕后”)。