我正在按年度时间序列执行黄土的季节性和趋势性分解,在那里我从原始数据中提取长期趋势和三年一次的准周期序列。问题是它似乎不存在ts
的选项,其周期性为年。然而,我能够通过对TS进行随机频率分析来克服这一局限,最终结果令人满意:
loessyearly <- ts(yrobs, start = 1993, frequency = 2) #create a TS with random frequency of 2
stlyearly <- stl(loessyearly, s.window = 3, t.window = 9) #s.window = 3 to smooth under a quasi-periodicity of three years
显然,在绘制图形时,X轴遵循我随机选择的频率(即2年):
但是,我希望x轴描绘每次观察一年(即频率= 1)。这可能吗?
我知道一个选项会尝试使用stlplus()
进行此分解,但我想首先通过stl()
尝试,因为我发现处理{的图形选项更加困难{1}}功能。
谢谢!