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