Maurício, boa tarde! Analisando apenas a parte final do seu script, temos apenas um objeto 'repeticoes' que é atualizado/mudado a cada ciclo do 'for'. Ao final da execução vai sobrar apenas o resultado da última execução. "A saída de "repeticoes" é mostrada assim: [[1]] [1] 0.11926720 -0.06925387 -0.49590355" A saída mostra que é impressa uma lista a cada rodada e para o que você precisa, devemos acumular as sucessivas saídas... Uma forma não muito pensada, mas que deve funcionar: for(i in 1:ns) { repeticoes <- replicate(1, geral(), simplify = FALSE) if (i==1) rep=unlist(repeticoes) else rep=rbind(rep,unlist(repeticoes)) } rep; colMeans(rep) Na primeira rodada o objeto 'rep' é o próprio 'repeticoes'. Nas sucessivas, agrega usando rbind(). O unlist() é pra desmontar a lista. Ao final você deverá ter um objeto 'rep' do tipo 'matrix' e então poderá usar colMeans(rep) pra ter as médias por coluna. Não testei porque também estou no R 3.0. Mas acredito que vá funcionar... Atte. Éder Comunello