lamda=5
theta=1000
size=1000
N <- rpois(size,lambda)
Y.new <- rexp(sum(N),1/theta)
Y.adj <- pmin(pmax(Y.new-500,0))
gg <- function(i) {sum(Y.adj[(1+sum(N[1:(i-1)])):(sum(N[1:i]))])}
我的“gg”功能可以很好地输入1000个样本中的单个值:
gg(44)
[1] 1611.536
gg(145)
[1] 407.7765
但是当我尝试创建一个包含gg函数的所有1000个值的矩阵时(这就是我需要的),我得到:
gg(1:1000)
Warning messages:
1: In 1:(i - 1) :
numerical expression has 1000 elements: only the first used
2: In 1:i : numerical expression has 1000 elements: only the first used
我只是想知道如何解决这个问题。谢谢!