消除R的季节性趋势

时间:2015-03-22 11:20:25

标签: r

我有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")

enter image description here

我在剧情中做错了什么?

我是时间序列中的新手,我担心这个过程,我知道它应该没问题,但我不确定生成的时间序列(x.sa)是否正常(没有&# 34;季节性周期&#34;)用新的&#34;去趋势&#34;做多元回归数据

修改

是的,a是数字,而不是ts对象,这就是我在stl之前使用ts的原因。 当我绘制tm和x.sa时,我得到了:

enter image description here

enter image description here

我错过了什么?时间对象错了吗? 我将提出一个样本(因为数据是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
       )

0 个答案:

没有答案