R中的递归for循环

Question

我很难对这个递归程序进行编码，并且想知道是否有人可以帮助我。

我想对这个递归方程进行编码：

for k=1,2,...

beta(k)=k-sum_(i=0)^(k-1)Kchoosei*beta(i)*exp(-i(k-i))

我在R中是用手动的方式做的，但是我想把它放在一个函数中。

beta0<-0

beta1<-1-choose(1,0)*beta0*exp(-0*lambdaL*(1-0))
beta2<-2-choose(2,0)*beta0*exp(-0*lambdaL*(2-0))-choose(2,1)*beta1*exp(-1*lambdaL*(2-1))
beta3<-3-choose(3,0)*beta0*exp(-0*lambdaL*(3-0))-choose(3,1)*beta1*exp(-1*lambdaL*(3-1))-choose(3,2)*beta2*exp(-2*lambdaL*(3-2))
beta4<-4-choose(4,0)*beta0*exp(-0*lambdaL*(4-0))-choose(4,1)*beta1*exp(-1*lambdaL*(4-1))-choose(4,2)*beta2*exp(-2*lambdaL*(4-2))-choose(4,3)*beta3*exp(-3*lambdaL*(4-3))

Answer 1

您可以为和定义第二个循环。请注意，这里的索引以1而不是0开始，这将导致“索引转移”。

beta = numeric()
beta[1] <- 0

for (k in 1:10){
  beta[k+1] <- k
  for (i in 0:(k-1))
    beta[k+1] <- beta[k+1] - choose(k, i)*beta[i+1]*exp(-i*(k-i))
}

beta
# [1] 0.000000 1.000000 1.264241 2.080705 3.247551 4.528104 5.748673
# [8] 6.876234 7.941197 8.972749 9.987645

Answer 2

我认为您需要两个函数，因为您需要每个以前的beta值作为输入，但是您只需要一个输出。下面是我的建议(为了适应您的lambaL，因为您的帖子有些含糊不清，它适用于LambdaL == 1)：

beta_vec <- function(k){
  if(k == 0) 0 else {
    beta_vec_old <- beta_vec(k-1)
    c(beta_vec_old,sum(sapply(0:(k-1),function(i){1-choose(k,i)*beta_vec_old[i+1]*exp(-i*(k-i))})))
}}

beta <- function(k){
  tail(beta_vec(k),1)  
}

# > beta_vec(5)
# [1] 0.000000 1.000000 1.264241 2.080705 3.247551 4.528104
# > beta(5)
# [1] 4.528104

(编辑后输入代码)