R季节性分解

Question

我模拟了一个具有周期性和线性成分的时间序列，并尝试使用R stl函数对其进行分析

n = 1000
x = ts(0.1*rnorm(n) + sin(6*pi*(1:n)/n) + (1:n)/n,frequency=n)
plot(x)
stl(x,"per")

但收到消息

Error in stl(x, "per") : 
  series is not periodic or has less than two periods

如何在模拟时间序列上使用 stl？

Answer 1

必须有2个以上的周期，所以频率必须小于n/2

n = 1000
x = ts(0.1*rnorm(n) + sin(6*pi*(1:n)/n) + (1:n)/n,
       frequency=n/2.1)

plot(x)
stl(x,"per")

R seasonal decomposition