R包“ seasonal”提供了与季节调整软件X13-ARIMA-SEATS的接口。
通过将参数x11 = ""
传递给seas
函数,使用了基于同名过滤器的算法。
根据参考manual,在没有明确指定的季节和趋势过滤器的情况下,该算法会自动选择它们:趋势过滤器为“对于每月系列,无论是9、13或23项亨德森移动平均”,而对于季节性过滤器,“ ...程序选择使用$ 3 \ times 3、3 \ times 5 $还是$ 3 \ times 9 $移动平均值”。
我找不到确定自动选择哪些过滤器及其权重的方法。我尝试了一次网络搜索,阅读了手册的明显部分,并阅读了“季节性”软件包的文档和插图。将参数传递给x11.save
参数似乎很有希望,但是我尝试过的一些显而易见的候选方法(例如 d10
,“最终季节性因素”)当然不是过滤器权重。
所以如何确定自动选择了哪些季节和趋势过滤器及其权重?
此信息包含在X13-ARIMA-SEATS运行的诊断结果中,该结果保存在扩展名为.udg
的“诊断摘要文件”中(请参阅X13-ARIMA-SEATS参考manual页。 14)。
udg
对象的seas
字符向量包含这些诊断结果-在我的运行中为351个元素。感兴趣的两个元素分别命名为sfmsr
(根据移动季节性比选择的最终季节性过滤器,或msr
,标准)和finaltrendma
(针对最终趋势移动平均亨德森过滤器)。
这些可以用通常的方式访问,但是“季节性”软件包为此目的还提供了便利功能udg
。