<div>Benilton, eu nem cheguei a rodar, analisei o que ele faz mesmo.</div><div> </div><div>Tem que resolver esse problema daquele loop que só faz perder tempo senão o próprio estudo de simulação dele está inválido.<br clear="all">

lmassis <at> yahoo <dot> com <dot> br<br>assis.leonard <at> gmail <dot> com<br>
<br><br></div><div class="gmail_quote">2011/8/18 Benilton Carvalho <span dir="ltr"><<a href="mailto:beniltoncarvalho@gmail.com">beniltoncarvalho@gmail.com</a>></span><br><blockquote style="margin: 0px 0px 0px 0.8ex; padding-left: 1ex; border-left-color: rgb(204, 204, 204); border-left-width: 1px; border-left-style: solid;" class="gmail_quote">

Tiago,<br>
<br>
esse e' um problema legal... se vc derivar as distribuicoes a mao,<br>
certamente vc consegue ganhos bem maiores...<br>
<br>
FYI: usando propriedades de estatistica, vc chega na mudanca que eu<br>
fiz na geracao das suas variaveis aleatorias.<br>
<br>
Adicionalmente, eu nao tive paciencia de esperar o seu codigo original<br>
acabar de rodar... eu o cancelei depois de um pouco mais de 1h10m.<br>
<br>
Abaixo uma sugestao para voce:<br>
<br>
<br>
library(foreach)<br>
library(doMC)<br>
registerDoMC()<br>
programa2 <- function(f,n,m) {<br>
<div class="im">    ## f: número máximo de filhos<br>
    ## n: número de gerações<br>
    ## m: número de iterações<br>
</div>    S <- foreach(p=seq(0,1,by=0.01), .combine=rbind) %dopar% {<br>
<div class="im">        pp <- 0<br>
        for (g in 1:m) {<br>
            y <- 1<br>
            for (k in 1:n)<br>
</div>                y <- rbinom(1, y*f, p)<br>
            if (y > 0) pp <- pp + 1<br>
        }<br>
        c(p, pp/m)<br>
    }<br>
    return(S)<br>
}<br>
system.time(prog2 <- programa2(2,200,1000))<br>
plot(prog2, xlab = "p", ylab = "Probabilidade de Percolação")<br>
<br>
para matar sua curiosidade, no meu computador, eu obtive:<br>
<br>
> system.time(prog2 <- programa2(2,200,1000))<br>
   user  system elapsed<br>
257.072   2.562  20.835<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br>
</blockquote></div><br>