使用先前计算结果的迭代计算 - R

Question

我不久前发布了这个问题，提出了很多很好的建议，但其中 none 的最终结果与我使用 Excel 获得的结果相同。

假设我们有以下参数：

Pd.sil.initial = 0.011
Pd.sul.i = 30
Pd.sul.f = ?
R.inc = 100
D.Pd = 536000

我需要使用以下公式计算 Pd.sul.f：

Pd.sul.f = (Pd.sul.i + (R.inc * Pd.sil.i)) / (1 + (R.inc / D.Pd))

这需要迭代完成，将上一次计算的Pd.sul.f结果用作下一次计算的Pd.sul.i。我提供了我的 Excel 工作表的屏幕截图，突出显示了这一点，以及前几个计算的答案。

之前向我建议的 R 代码，但不幸的是没有提供与我的 Excel 工作表相同的结果：

Pd.sul.f <- function(n) {

 return(Pd.sul.i + n * (R.inc * Pd.sil.i) / (1 + R.inc/D.Pd))

}

sapply(1:100000, Pd.sul.f)

如果能帮助我弄清楚如何实现这一点，我们将不胜感激。我是 R 的新手，所以非常感谢更简单的答案。

预先感谢您的帮助。

Answer 1

这应该有效：

Pd.sil.initial = 0.011
Pd.sul.i = 30
R.inc = 100
D.Pd = 536000

n=10

a <- data.frame(matrix(,nrow=n, ncol=1))

for (i in 1:n) {
   Pd.sul.f <- ((Pd.sul.i + (R.inc * Pd.sil.initial)) / (1 + R.inc / D.Pd))
   a[i,] <- Pd.sul.f
   Pd.sul.i <- Pd.sul.f
}

其中 n 是迭代次数。

编辑以回答 Matt 第二个问题：

Pd.sil.initial = 0.011
Pd.sul.i = 30
R.inc = 100
D.Pd = 536000

n=10

a <- data.frame(matrix(,nrow=n, ncol=2))
a[1,1] <- 30
a[1,2] <- 100

for (i in 2:n) {
   Pd.sul.f <- ((Pd.sul.i + (R.inc * Pd.sil.initial)) / (1 + R.inc / D.Pd))
   a[i,1] <- Pd.sul.f
   a[i,2] <- i * 100
   Pd.sul.i <- Pd.sul.f
}

或者简单地说：

for (i in 2:n) {
  a[i,1] <- ((a[[1]][[i-1]] + (R.inc * Pd.sil.initial)) / (1 + R.inc / D.Pd))
  a[i,2] <- a[[2]][[i-1]] + 100
}

使用先前计算结果的迭代计算 - R

Iterative calculations using result from previous calculation - R

iteration

r

calculation