我有10年(2000-2009)的每日数据,我想确定季节周期,删除它,并使用这些"去趋势"回归分析中的数据。为此,我使用ts和stl作为:
tm <- ts(a, start = 2000, freq = 365.25)
tm2 <- stl(tm, s.window="periodic")
x.sa <- seasadj(tm2)
当我绘制&#34;趋势&#34; ,我有一个奇怪的情节:
plot(a,col="grey")
lines(x.da,col="red")
我在剧情中做错了什么?
我是时间序列中的新手,我担心这个过程,我知道它应该没问题,但我不确定生成的时间序列(x.sa)是否正常(没有&# 34;季节性周期&#34;)用新的&#34;去趋势&#34;做多元回归数据
修改
是的,a是数字,而不是ts对象,这就是我在stl之前使用ts的原因。 当我绘制tm和x.sa时,我得到了:
我错过了什么?时间对象错了吗? 我将提出一个样本(因为数据是10年的每日数据:3653次观察),所以,这将是一个简短的样本:
mydata<- dput(a[1:40])
c(281.127075195312, 280.629150390625, 282.007202148438, 280.685577392578,
278.718994140625, 281.126098632812, 279.95458984375, 280.4013671875,
278.469146728516, 274.627807617188, 277.929779052734, 276.451690673828,
275.632934570312, 277.067443847656, 277.834289550781, 280.159820556641,
280.40087890625, 280.572998046875, 278.621704101562, 279.200439453125,
279.614135742188, 278.984436035156, 276.867004394531, 272.86865234375,
275.842346191406, 278.622436523438, 277.736022949219, 277.709503173828,
282.733032226562, 281.467529296875, 282.947448730469, 282.413330078125,
281.055633544922, 281.113983154297, 281.544097900391, 282.1435546875,
281.356384277344, 281.776458740234, 281.724639892578, 280.923461914062
)