我想对每个模拟求和(并求平均值),例如 如果模拟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。