用速率R的向量计算复合利息

时间:2018-05-29 21:46:59

标签: r loops vectorization

我试图看看是否有办法对我执行的计算进行矢量化。我搜索了这个答案,无法找到我需要的东西。

我有一个增长率矢量。每个代表一个时期(在我的情况下为一年)。我想将这个向量应用到一些本金额。然后,在将第一个增长率应用于主体后,使用第一次迭代的结果并将第二个增长元素应用于新值。

这里有一些复制代码(全部在base中):

# Vector of interest or inflation rates
RateVector <- c(0.02, 0.03, 0.04, 0.05, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01) # forecasted rates
Principal <- data.frame(Principal = 1000000) # actual value of indicator in most recent period as data frame (list)

这是我尝试进行矢量化的方法:

sapply(Principal, "*", 1 + cumsum(RateVector))

问题在于sapply函数不保存新的金额,而是将汇率向量应用于相同的初始本金。这实际上是我对此代码的期望。我不知道如何在每次迭代之后从元素到元素保存新值。

这就是我使用循环解决问题的方法:

AmountVector <- Principal # initialize output vector

# Compound growth growth calculation loop
for(i in 1:length(RateVector)){
  Principal = Principal * (1 + RateVector)[i]
  AmountVector <- rbind(AmountVector,Principal)
}

# Print result
AmountVector

感谢大家的帮助。

1 个答案:

答案 0 :(得分:3)

这是一个&#34;累积产品&#34;,因此您需要?cumprod

1000000 * cumprod(1+RateVector)
# [1] 1020000 1050600 1092624 1147255 1216091 1276895 1327971 1367810 1395166
#[10] 1409118

cbind(AmountVector, newresult = 1000000 * c(1,cumprod(1+RateVector)))
#   Principal newresult
#1    1000000   1000000
#2    1020000   1020000
#3    1050600   1050600
#4    1092624   1092624
#5    1147255   1147255
#6    1216091   1216091
#7    1276895   1276895
#8    1327971   1327971
#9    1367810   1367810
#10   1395166   1395166
#11   1409118   1409118