给出每月ts对象,例如:
dat <- ts(c(295, 286, 300, 278, 272, 268, 308, 321, 313, 308, 291, 296,
294, 273, 300, 271, 282, 285, 318, 323, 313, 311, 291, 293, 297,
273, 294, 259, 276, 294, 316, 325, 315, 312, 292, 301), frequency = 12)
如何按月计算平均值?即我想计算1月,1年+ 1月,2年+ 1月,3年等的平均值。然后能够比较二月的......
我想到的一种方法是将其转换为12列的矩阵并使用colMeans()
,但我想有一种更好的方法可以利用time()
的{{1}}方面对象
ts()
答案 0 :(得分:9)
好的,所以我应该在进入SO之前再给Google一次搜索post is relevant。 cycle()
函数似乎对这些事情很有用:
> tapply(dat, cycle(dat), mean)
1 2 3 4 5 6 7 8 9
295.33333 277.33333 298.00000 269.33333 276.66667 282.33333 314.00000 323.00000 313.66667
10 11 12
310.33333 291.33333 296.66667
> aggregate(c(dat), list(month = cycle(dat)), mean)
month x
1 1 295.33333
2 2 277.33333
3 3 298.00000
....
我还缺少其他基本的东西吗?