如何设置ts()的frequency

一直没有弄清楚ts()里面frequency的设置，这两天查找资料发现stackoverflow的一篇文章详细说明了这个设置，翻译出来，和大家共享。

 

 

#创建一个虚拟的数据框，包含数值和时间两个变量，下面显示了不同的设置的结果。dat <- data.frame(myts = sample(10, 24, replace = T),                  Date = seq(as.Date("2008-10-11"), as.Date("2008-10-11") + 23, by = 1))print(ts(dat$myts, frequency = 7, start = c(1950, 3)), calendar = T)#      p1 p2 p3 p4 p5 p6 p7# 1950        6  9  6  9  8# 1951  6  1  9  3  5  7  4# 1952  2  9  5  9  7  8  7# 1953  2  6  6  8  1       print(ts(dat$myts, frequency = 12, start = c(1950, 3)), calendar = T)#      Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec# 1950           6   9   6   9   8   6   1   9   3   5# 1951   7   4   2   9   5   9   7   8   7   2   6   6# 1952   8   1                                         print(ts(dat$myts, frequency = 4, start = c(1950, 3)), calendar = T)#      Qtr1 Qtr2 Qtr3 Qtr4# 1950              6    9# 1951    6    9    8    6# 1952    1    9    3    5# 1953    7    4    2    9# 1954    5    9    7    8# 1955    7    2    6    6# 1956    8    1           print(ts(dat$myts, frequency = 7), calendar = T)#   p1 p2 p3 p4 p5 p6 p7# 1  6  9  6  9  8  6  1# 2  9  3  5  7  4  2  9# 3  5  9  7  8  7  2  6# 4  6  8  1      

 

 从输出结果我们可以发现：

1：ts()自动将frequency=12,4,当成月份和季度。而如果要设置成7，表示频率为每周时，显示结果并不很直接。

2：start()中的第一个参数依赖于frequency()的设置，第二个参数是frequency设置周期的第一个事件发生时间（即并不是所有序列都从周一或者周日开始）

3：如果不确定start()参数的设置，即序列的起始点，函数默认为第一个周期的第一个时间点。

 

因此，如果设置周期为7（周），想让序列更清晰表达起始范围，可以采用以下方法。以下方法使用了strfime()函数

 

startW <- as.numeric(strftime(head(dat$Date, 1), format = "%W"))startD <- as.numeric(strftime(head(dat$Date, 1) + 1, format =" %w")) print(ts(dat$myts, frequency = 7, start = c(startW, startD)), calendar = T)#   p1 p2 p3 p4 p5 p6 p7#39                    6#40  9  6  9  8  6  1  9#41  3  5  7  4  2  9  5#42  9  7  8  7  2  6  6#43  8  1   

Which means that our first observation (which occurred in 2008-10-11) was Saturday of the 39th week of 2008