Calculando a média de cada replicação
eu tenho o seguinte código
set.seed(30)
nsim <- 50 ## NUMBER OF REPLICATIONS
demand <- c(12,13,24,12,13,12,14,10,11,10)
res <- replicate(nsim, {
load <- runif(10,11,14)
diff <- load - demand ## DIFFERENCE BETWEEN DEMAND AND LOAD
return(sum(diff < 0))
})
res
[1] 6 5 7 4 4 5 4 3 6 4 5 5 5 4 2 5 3 3 3 5 3 2 4 6 5 4 4 3 5 6 4 4 3 6 5 3 5 5 4 3 3
[42] 6 4 4 4 6 6 5 4 5
Eu tenho um grande conjunto de dados e a questão é qual é a maneira mais rápida de calcular a média para cada replicação. Por exemplo, o res na primeira replicação é 6, portanto, o resultado deve ser 6/1 = 6 para o segundo 6 + 5/2 = 5,5 para o terceiro 6 + 5 + 7/3 = 6 e para a última replicação é soma (res ) /nsim=4.38