r - 为三个输入参数的每个组合导出公式结果

Question

我一直在努力解决 R 中的一个问题，并尝试过 google 和搜索功能，但没有找到解决方案。

我想执行以下操作:

得出一组久期、金额和利率组合的利率支付收入。该集合定义如下:

1. 金额 (a) = 1,000 - 25,000，间隔为 100
2. 持续时间 (d) = 12 - 60(月)，间隔为 6
3. 利率 (i) 0.00 - 0.015，间隔为 0.001

公式如下:w = ad(((i/12)*(1+(i/12))^(d))/((1+(i/12))^(d)-1)) - 一个 ，其中 w 是利率支付。

我的解决方案是使用三个“for”循环计算 w 并使用“rbind”命令保存数据。遗憾的是这对我不起作用:

result = NULL

amount <- seq(1000, 25000, by=100)

duration <- seq(12, 60, by=6)

interest <- seq(0, 0.015, by=0.0001)

for (d in duration){
    for (a in amount){
      for (i in interest){
        w = a*d*(((i/12)*(1+(i/12))^(d))/((1+(i/12))^(d)-1)) - a
      }
    }
}

result = rbind(d, data.frame(d, a, i, w))

我做错了什么？我走在正确的轨道上吗？输出应该是一个包含所有可能组合和“w”值的列表。或者，输出可以是一个列表，其中包含特定 w 值的所有可能组合。

提前非常感谢您!

最好，亚历克斯

Answer 1

以下是调整当前方法的方法:

amount <- seq(1000, 25000, by=100)
duration <- seq(12, 60, by=6)
interest <- seq(0, 0.015, by=0.0001)
result <- expand.grid(amount, duration, interest)
names(result) <- c("amount", "duration", "interest")

for(i in 1:nrow(result)){
  result$w[i] <- result$amount[i]*result$duration[i]*(((result$interest[i]/12)*(1+(result$interest[i]/12))^(result$duration[i]))/((1+(result$interest[i]/12))^(result$duration[i])-1)) - result$amount[i]
}

但是，data.table 的速度大约快一百万倍:

library("data.table")
amount <- seq(1000, 25000, by=100)
duration <- seq(12, 60, by=6)
interest <- seq(0, 0.015, by=0.0001)
result <- expand.grid(amount, duration, interest)

setDT(result)
setnames(result, c("Var1", "Var2", "Var3"), c("amount", "duration", "interest"))
result[, w := amount*duration*(((interest/12)*(1+(interest/12))^(duration))/((1+(interest/12))^(duration)-1)) - amount]