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
}