循环 function 通过 R 中的连续乘法 function 获得指数倍数

Question

I try to calculate values from the function:

And the loop function I write so far is:

However, even though it can run, but when I apply this function to my data, there is an error says "Error: unused argument (.)" Maybe because y[q,k] <- ((3+0.01*(q-1))^(-(1:k)))*sum(cumprod(px[1:k])) is wrong? But it looks like the equation I want. Anyone know how can I fix this problem? Thanks in advance.

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)

And the error shows

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

Answer 1

Maybe you can try the code below

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
}