4

我有一个时间序列数据集，我想获得 12 个月的滚动累积回报。下面是我的代码和数据的样子：

df <- data.frame(
  A = c(-0.0195,    0.0079, 0.0034, 0.0394, -0.0065,    0.0034, 0.0136, 0.0683, -0.0063,    -0.0537,    -0.0216,    -0.0036,    0.0659, -0.0377,    -0.0568,    0.0039, -0.0191,    0.0028),
  B = c(-0.0211,    0.0021, 0.0014, 0.0358, 0.0009, 0.0153, 0.0071, 0.0658, 0.0033, -0.0542,    -0.0261,    0.0064, 0.0665, -0.0304,    -0.0507,    0.0089, NA, NA),
  C= c(-0.0176, 0.0144, 0.0057, 0.0442, -0.0152,    -0.0105,    0.0213, 0.0712, -0.0176,    -0.0531,    -0.0163,    -0.0154,    0.0652, NA, NA, NA, NA, NA)
)
row.names(df) <- c("2016-10-31",    "2016-09-30",   "2016-08-31",   "2016-07-31",   "2016-06-30",   "2016-05-31",   "2016-04-30",   "2016-03-31",   "2016-02-29",   "2016-01-31",   "2015-12-31",   "2015-11-30",   "2015-10-31",   "2015-09-30",   "2015-08-31",   "2015-07-31",   "2015-06-30",   "2015-05-31")

> df
                  A       B       C
2016-10-31  -0.0195 -0.0211 -0.0176
2016-09-30   0.0079  0.0021  0.0144
2016-08-31   0.0034  0.0014  0.0057
2016-07-31   0.0394  0.0358  0.0442
2016-06-30  -0.0065  0.0009 -0.0152
2016-05-31   0.0034  0.0153 -0.0105
2016-04-30   0.0136  0.0071  0.0213
2016-03-31   0.0683  0.0658  0.0712
2016-02-29  -0.0063  0.0033 -0.0176
2016-01-31  -0.0537 -0.0542 -0.0531
2015-12-31  -0.0216 -0.0261 -0.0163
2015-11-30  -0.0036  0.0064 -0.0154
2015-10-31   0.0659  0.0665  0.0652
2015-09-30  -0.0377 -0.0304      NA
2015-08-31  -0.0568 -0.0507      NA
2015-07-31   0.0039  0.0089      NA 
2015-06-30  -0.0191      NA      NA
2015-05-31   0.0028      NA      NA

prod(1+R)-1我想使用公式（12 个单独期间收益的乘积减去 1）获得每个月的 12 个月累积收益。结果应该是：

             A(1-Y)  B(1-Y)  C(1-Y)
2016-10-31   0.0198  0.0322  0.0052
2016-09-30   0.1086  0.1246  0.0898
2016-08-31   0.0585  0.0881 
2016-07-31  -0.0050  0.0316 
2016-06-30  -0.0390  0.0048 
2016-05-31  -0.0512     
2016-04-30  -0.0517

B(1-Y)只有 5 次累积收益，因为B在 之前没有数据2015-07-31。因此，它不满足 12 个月的条件（2015-07-31和之间只有 11 个月2016-05-31）。

我试过Return.cumulative(df)了，但是这个函数从一开始就给出了累积回报，这不是我想要的。任何建议将不胜感激！

4

1 回答 1

2

正如 A.Webb 提到的，rollapply给出了预期的结果。您也可以使用 for 循环。

a <- b <- c <- 0
for(i in (nrow(df)-11):1){
    a[i] <- prod(1 + df$A[i:(i+11)]) - 1
    b[i] <- prod(1 + df$B[i:(i+11)]) - 1
    c[i] <- prod(1 + df$C[i:(i+11)]) - 1
}

cum.returns <- data.frame(a,b,c)
于 2017-02-08T23:15:29.037 回答