r - 对每个模拟求和

Question

我想对每个模拟求和（和平均），例如，如果模拟 1 有 7 个事件，那么这些值应该加在一起。

代码：

sims <- 10

freq.mean <- 3
freq <- rpois(sims, lambda = freq.mean)

gamma_alpha <- 2
gamma_theta <- 5
freq.sev <- lapply(freq, function(x) rgamma(x, shape = gamma_alpha, scale = gamma_theta))

我想要每个模拟中的总和列表：例如总和simulation[[1]]将是：

78.3 =
sum(9.338442,   4.187821,   1.684569,   22.362444,  12.177631,  15.595493,  12.959588)

如果 freq 在第一次模拟中返回 7。

标签： rsumsimulation