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
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
