1

我有 7 天期间测量的数据。部分数据如下:

start wk    end wk      X1
2/1/2004    2/7/2004    89
2/8/2004    2/14/2004   65
2/15/2004   2/21/2004   64
2/22/2004   2/28/2004   95
2/29/2004   3/6/2004    79
3/7/2004    3/13/2004   79

我想使用 X1 的加权平均值将此每周(7 天)数据转换为每月数据。请注意,某些 7 天 X1 数据将在一个月之间重叠(X1=79,2004 年 2/29 至 3/6 期间)。

具体来说,我将通过以下方式获得 2004 年 2 月的月度数据(例如 Y1)

(7*89 + 7*65 + 7*64 + 7*95 + 1*79)/29 = 78.27

R 是否具有可以正确执行此操作的功能?(xts 库中的 to.monthly 不做我需要的)如果,不是在 R 中做到这一点的最佳方法是什么?

4

2 回答 2

3

将数据转换为每日数据,然后聚合:

Lines <- "start end X1
2/1/2004    2/7/2004    89
2/8/2004    2/14/2004   65
2/15/2004   2/21/2004   64
2/22/2004   2/28/2004   95
2/29/2004   3/6/2004    79
3/7/2004    3/13/2004   79
"

library(zoo)

# read data into data frame DF
DF <- read.table(text = Lines, header = TRUE)

# convert date columns to "Date" class
fmt <- "%m/%d/%Y"
DF <- transform(DF, start = as.Date(start, fmt), end = as.Date(end, fmt))

# convert to daily zoo series
to.day <- function(i) with(DF, zoo(X1[i], seq(start[i], end[i], "day")))
z.day <- do.call(c, lapply(1:nrow(DF), to.day))

# aggregate by month
aggregate(z.day, as.yearmon, mean)

最后一行给出:

Feb 2004 Mar 2004 
78.27586 79.00000 
于 2013-04-10T18:02:45.727 回答
0

如果你愿意从你的 DF 中去掉“结束周”,apply.monthly 会像一个魅力一样工作。

DF.xts <- xts(DF$X1, order.by=DF$start_wk)

DF.xts.monthly <- apply.monthly(DF.xts, "sum")

然后,如果您绝对需要它们,则始终可以通过添加 30 来重新创建结束日期。

于 2014-10-30T00:38:48.293 回答