R - 将频率为20分钟的时间序列放入函数stl()

时间:2013-07-19 05:32:23

标签: r stl time-series

我正在尝试使用R中的函数stl()来平滑一个数据集,该数据集存储传感器的时间与温度。主要目标是从功能中去除噪声,以便找到环境温度(基线曲线)。我已经按时间序列格式加载了我的数据,当我尝试使用stl()时,它给了我这个错误:

Error in stl(timeseries[[1]]) : 
  series is not periodic or has less than two periods

这是我的数据:

> head(stations[[1]])
                 Date Unit Temp
1 0013-06-30 10:00:01    C 32.5
2 0013-06-30 10:20:01    C 32.5
3 0013-06-30 10:40:01    C 33.5
4 0013-06-30 11:00:01    C 34.5
5 0013-06-30 11:20:01    C 37.0
6 0013-06-30 11:40:01    C 35.5

我已转换为时间序列类:

timeseries[[1]] = as.ts(stations[[1]]$Temp,freq=26280)

注意:频率很高,因为数据是每20分钟拍摄一次

由于频率不一致导致stl()的错误?我觉得在将数据制作成时间序列时我可能做错了什么,并且这已经抛弃了计算系列周期的能力

我确实需要所有这些数据,因为整套只能覆盖4天的数据(因此频率很高)

感谢您的帮助!

1 个答案:

答案 0 :(得分:6)

错误消息告诉您,为了估算时间序列的季节性成分,您需要至少两个季节的数据。如果您有4天的温度数据,您可能希望季节性成分为几天。因此,您应该相应地设置时间序列。你每天有24 * 3次观察,所以这应该是频率。

timeseries[[1]] <- ts(stations[[1]]$Temp, frequency=24*3)

然后stl(timeseries[[1]], "periodic")应该可以工作,尽管我无法测试,因为它需要至少2天的数据,即2小时是不够的。