用R分析时间序列

Question

我是R的新手，继Walter Zucchini与R的时间序列分析PDF之后。我有一些来自传感器的数据，特别是我可以每分钟或每5秒获得一次数据。 然后我想使用ts()命令来创建这些值的时间序列。所以语法应该是data1mints <- ts(data1min ,freq = 525600)，其中525600是常规年份的分钟数 之后，我尝试使用此命令plot(stl(log(data1min), s.window = "periodic"))进行绘图，但R说我

  

系列不是周期性的或少于两个周期

更准确地说，我有3月20日至3月28日的数据，所以我没有完整的年度数据，但我认为这是分析每分钟发生的事情的足够时期。

我错了什么？

Answer 1

错误消息告诉您错误 - 您的时间少于2个。

例如，

# this works since there are 3 periods
freq <- 100
ny <- 3 # no of years, i.e. periods
n <- ny * freq
set.seed(13)
tt <- ts(rnorm(n), freq = freq)
s <- stl(tt, "periodic")

# this issues error since there are less than 2 periods. (We have changed ny to 1.)
freq <- 100
ny <- 1 ##
n <- ny * freq
set.seed(13)
tt <- ts(rnorm(n), freq = freq)
s <- stl(tt, "periodic")