2

我正在尝试使用 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 天的数据（因此频率很高）

谢谢您的帮助！

4

1 回答 1

7

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

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

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

于 2013-07-19T07:08:19.870 回答