Seasonal ARIMA(1,0,1)x(1,0,1)[4]模型fixed参数长度错误问题咨询
问题原因解析
核心原因:均值项的存在与否导致参数数量不匹配
当设置
seasonal=c(1,0,1)时,季节差分阶数D=0,模型默认会包含均值参数(未做季节差分的数据集需要拟合基准均值)。此时模型的待估参数包括:- 非季节AR(1)系数
phi - 非季节MA(1)系数
theta - 季节AR(1)系数
Phi - 季节MA(1)系数
Theta - 均值项
mean
总共5个参数,但你在fixed里只指定了4个值,参数长度不匹配,因此触发报错。
- 非季节AR(1)系数
当把季节差分阶数改为
D=1(即seasonal=c(1,1,1))时,季节差分操作会消除数据的季节均值趋势,模型默认不再包含均值项。此时待估参数正好是你指定的4个(phi、theta、Phi、Theta),fixed参数长度完全匹配,因此不会报错。
修正方案
如果要在seasonal=c(1,0,1)的情况下使用fixed参数,需要把均值项也纳入固定值,示例代码如下:
library(forecast) set.seed(0) model <- Arima(ts(rnorm(100),freq=4), order=c(1,0,1), seasonal=c(1,0,1), fixed=c(phi=0.5, theta=-0.4, Phi=0.3, Theta=-0.2, mean=0))
内容的提问来源于stack exchange,提问作者Antoni Parellada
相关产品推荐
相关产品推荐

