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

如何解决使用先前计算的结果进行迭代计算 - R

我不久前发布了这个问题，提出了很多很好的建议，但没有一个最终得到与我使用 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 比较陌生，因此非常感谢您提供更简单的答案。

预先感谢您的帮助。

解决方法

这应该有效：

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 是迭代次数。

编辑回答马特第二个问题：

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

n=10

a <- data.frame(matrix(,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
}