我试图看看是否有办法对我执行的计算进行矢量化。我搜索了这个答案,无法找到我需要的东西。
我有一个增长率矢量。每个代表一个时期(在我的情况下为一年)。我想将这个向量应用到一些本金额。然后,在将第一个增长率应用于主体后,使用第一次迭代的结果并将第二个增长元素应用于新值。
这里有一些复制代码(全部在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
感谢大家的帮助。
答案 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