我有时间序列数据，我怎么能做每月报告和平均为一个变量?(使用R)

Question

我的时间序列数据看起来像是6个月的8/18/2012 11:18:00 PM，我如何将它们按月子集并在一个月内对变量进行平均?(使用R)

Thank you so much

Answer 1

您可以使用xts包。首先我生成你的数据。这里我创建了一个6个月，每天一半的数据。长数据

dat <- data.frame(date = as.POSIXct('8/18/2012 11:18:00',
                                    format='%m/%d/%Y %H:%M:%S') + 
                         seq(0,by = 60*60*12,length.out=365),
                  value = rnorm(365))

然后创建一个xts对象

library(xts)
dat.xts <- xts(x= dat$value,order.by = dat$date)

最后，我使用方便的函数apply.monthly相当于lapply，得到这样的东西:

apply.monthly(dat.xts,mean)
2012-08-31 23:18:00  0.03415933
2012-09-30 23:18:00  0.02884122
2012-10-31 22:18:00 -0.27767240
2012-11-30 22:18:00 -0.15614567
2012-12-31 22:18:00 -0.02595911
2013-01-31 22:18:00 -0.23284335
2013-02-16 10:18:00  0.14537790

Answer 2

您可以format日期并使用aggregate计算平均值(感谢示例数据的@agstudy):

aggregate(value~format(date,"%Y-%m"),dat,FUN=mean)
  format(date, "%Y-%m")       value
1               2012-08 -0.31409786
2               2012-09 -0.37585310
3               2012-10 -0.04552703
4               2012-11 -0.05726177
5               2012-12  0.04822608
6               2013-01  0.03412790
7               2013-02 -0.10157931