Éder, Muito obrigado pela atenção. Pelo que observei aqui, isso soluciona o problema. Grato, Maurício
Em 26 de outubro de 2013 16:51, Éder Comunello <[email protected]>escreveu: > 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. >
_______________________________________________ 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.
