0

我有一些每日时间序列数据,我需要提取相对于一周平均值的“工作日百分比”。例如,如果第一周有mean = 100并且本周的星期日值为 20,则星期日变为 0.2。

这是一些随机数据:

set.seed(0)
y = rnorm(1705)
date = seq(as.Date('2008-01-01'), by = 'day', length = length(y))
data.df = data.frame(y, date)

我需要一个名为 的新列pecent,即上面解释的值。我尝试添加一些列然后使用tapply,但失败了。感谢任何帮助!

4

1 回答 1

3

首先week使用创建一个变量format。然后使用ddplytransform

library(plyr)
data.df$week <- format(data.df$date,'%W %Y') #week begins with Monday
data.df <- ddply(data.df,~week,transform,percent=y/mean(y))
head(data.df)

           y       date    week    percent
1  1.2629543 2008-01-01 00 2008  3.1395415
2 -0.3262334 2008-01-02 00 2008 -0.8109741
3  1.3297993 2008-01-03 00 2008  3.3057095
4  1.2724293 2008-01-04 00 2008  3.1630952
5  0.4146414 2008-01-05 00 2008  1.0307451
6 -1.5399500 2008-01-06 00 2008 -3.8281172

请注意,第 00 周通常不像一年中的最后一周那样是一整周。如果重要的话,合并随后几年的最后几周和第一周。

于 2012-11-30T14:16:03.273 回答