Я хотел бы суммировать (и в среднем) по каждому моделированию, например,
если в симуляции 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 в первой симуляции.