r语言 - 时间序列每周数据

Question

我对R很陌生，并尝试在Excel中操作我的数据，然后再将其作为csv文件移动到R中。 我想使用ts()作为每周数据。 你能做一些像第 1-12 个月、第 1-4 周、一年这样简单的事情吗？ 哗啦啦。 我在想如果我使用第 1-7 天作为第一周，依此类推，那将是统一的并且很容易达到我的目的，但我不知道如何写它。 使用这个网站和另一个教程，我想出了这个：

myts <- ts(Time2012, start = c(8/3/2013,1), end = c(9/2/2013,4), frequency = 52)

有没有简单的方法来表示日期以显示我想计算周数？

Answer 1

我会推荐一个稍微不同的工作流程，您可能会发现它具有更广泛的实用性：

> end = Sys.Date()
> start = end - 365
> class
> # create the index array comprised of date objects
> ndx = seq(start, end, by='weeks')
> class(ndx)
  [1] "Date"
> length(ndx)
  [1] 53
> # create a fake data array
> x = 1:length(ndx)
> mydata = sin(x/2)
> # import a time series library 
> require(xts)
> # create the time series
> myts = xts(mydata, order.by=ndx)
> myts[1:5]
               [,1]
  2012-09-19 3.479426
  2012-09-26 3.841471
  2012-10-03 3.997495
  2012-10-10 3.909297
  2012-10-17 3.598472
> class(myts)
  [1] "xts" "zoo"
> periodicity(myts)
  Weekly periodicity from 2012-09-19 to 2013-09-18 

或者，如果您的数据不是按周，那么您可以创建一个具有更高分辨率（例如，天）的时间序列，然后将其汇总到周：

> ndx = seq(start, end, by='days')
> x = 1:length(ndx)
> mydata = sin(x/2) + 3
> myts = xts(mydata, order.by=ndx)
> myts[1:5]  
             [,1]
2012-09-19 3.479426
2012-09-20 3.841471
2012-09-21 3.997495
2012-09-22 3.909297
2012-09-23 3.598472
> periodicity(myts)
    Daily periodicity from 2012-09-19 to 2013-09-19 
> # now roll-up this daily series to weeks
> require(xts)
> # first create the endpoints
> np = endpoints(myts, on='weeks')

> myts_weeks = period.apply(x=myts, INDEX=np, FUN=sum, na.rm=TRUE)
> myts_weeks[1:5]
               [,1]
  2012-09-23 18.82616
  2012-09-30 17.11212
  2012-10-07 24.93492
  2012-10-14 17.51811
  2012-10-21 23.58635
> periodicity(myts_weeks)
  Weekly periodicity from 2012-09-23 to 2013-09-19 

Answer 2

创建

一组年、月、周（月）的更简单方法是使用 lubridate 。

require(lubridate)
# Your starting date, plus 52 more dates at weekly intervals
xDates <- dmy("8/3/2013") + weeks(0:52)
# A data frame of the dates, the month of the year, and the week of the month
xYMW <- data.frame(date=(xDates), month=month(xDates), week=mday(xDates) %/% 7 + 1)
xYMW[1:5, ]
        date month week
1 2013-03-08     3    2
2 2013-03-15     3    3
3 2013-03-22     3    4
4 2013-03-29     3    5
5 2013-04-05     4    1