R中的连续乘法函数的指数倍的循环函数

Question

我尝试从该函数计算值：

到目前为止，我编写的循环函数是：

但是，即使它可以运行，但是当我将此函数应用于数据时，出现错误消息“错误：未使用的参数（。）”。也许是因为 y[q,k] <- ((3+0.01*(q-1))^(-(1:k)))*sum(cumprod(px[1:k]))是错误的吗？但它看起来像我想要的等式。有人知道我该如何解决这个问题？预先感谢。

Y <- function(x){
  k <- 50-x
  Cset <- seq(0,0.2,0.01)
  y <- matrix(length(Cset)*k,ncol = k)
  for(C in seq_along(Cset)){
    for(s in 1:k){
      y[C,s] <- sum((3+C)^(-s)*cumprod(data$Bx[1:s]))
    }
  }
  y
}

Y(x=30)

错误显示

Error in `[<-`(`*tmp*`, C, s, value = sum((3 + C)^(-s) * cumprod(data$px[1:s]))) : 
  subscript out of bounds

Answer 1

也许您可以尝试下面的代码

Y <- function(x){
  k <- 50-x
  Cset <- seq(0,0.2,0.01)
  y <- matrix(length(Cset)*k,ncol = k)
  for(C in seq_along(Cset)){
    for(s in 1:k){
      y[C,s] <- (3+C)^(-s)*sum(cumprod(Bx[1:s]))
    }
  }
  y
}