<div>Tiago,</div><div> </div><div>pra que o loop "for (g in 1:m)" se em momento algum do código vc faz uso deste. Como vc construiu, vc rodar com m = 1 ou m - 700 trilhoes produzirá sempre o mesmo resultado, pois vc somente armazena e retorna a ultima iteração.</div>

<div> </div><div>repare neste trecho de seu código:</div><div> </div><div>for (g in 1:m) {</div><div>y <- 1</div><div>for (k in 1:n) {</div><div>y <- sum(rbinom(y,f,p))</div><div>if (y == 0) break</div><div>}</div>
<div>
if (y > 0) pp <- pp + 1</div><div>}</div><div> </div><div>viu? está rodando m vezes por aqui e nada de usar o g</div><div> </div><div>mais adiante vc tem</div><div> </div><div><br>S[p/0.01+1,1] <- p</div><div>S[p/0.01+1,2] <- pp/m</div>

<div> </div><div> </div><div>e em seguida, retorna o S</div><div> </div><div>ou seja: um loop completamente inutil</div><div><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 Tiago Vieira <span dir="ltr"><<a href="mailto:tiagovieira7@yahoo.com.br">tiagovieira7@yahoo.com.br</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">

<div><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 12pt; background-color: rgb(255, 255, 255);"><div><font size="2" face="arial, helvetica, sans-serif">Estou tentando rodar o programa de simulação de <span>percolação (rede de bethe)</span>, mas esta <span>demandando</span> um tempo de processamento muito alto, pois há preciso de muitas iterações. O computador que estou utilizando tem processador i3 e 4<span>gb</span> de memória. Alguém pode me dar uma dica de como posso diminuir o tempo de processamento? Abaixo segue o código:</font></div>


<div> </div>
<div><span lang="EN">
<div><font size="2" face="arial, helvetica, sans-serif">programa <- <span>function</span>(f,n,m) {</font></div>
<div><font size="2" face="Arial"># f: número máximo de filhos</font></div>
<div><font size="2" face="Arial"># n: número de gerações</font></div>
<div><font size="2" face="Arial"># m: número de iterações</font></div>
<div><font size="2" face="arial, helvetica, sans-serif">S <- <span>matrix</span>(0,101,2)</font></div>
<div><font size="2" face="arial, helvetica, sans-serif">for (p <span>in</span> <span>seq</span>(0,1,<span>by</span>=0.01)) {</font></div>
<div><font size="2" face="arial, helvetica, sans-serif"><span>pp</span> <- 0</font></div>
<div><font size="2" face="arial, helvetica, sans-serif">for (g <span>in</span> 1:m) {</font></div>
<div><font size="2" face="arial, helvetica, sans-serif">y <- 1</font></div>
<div><font size="2" face="arial, helvetica, sans-serif">for (k <span>in</span> 1:n) {</font></div>
<div><font size="2" face="arial, helvetica, sans-serif">y <- <span>sum</span>(<span>rbinom</span>(y,f,p))</font></div>
<div><font size="2" face="arial, helvetica, sans-serif"><span>if</span> (y == 0) <span>break</span></font></div>
<div><font size="2" face="arial, helvetica, sans-serif">}</font></div>
<div><font size="2" face="arial, helvetica, sans-serif"><span>if</span> (y > 0) <span>pp</span> <- <span>pp</span> + 1</font></div>
<div><font size="2" face="arial, helvetica, sans-serif">}</font></div>
<div><font size="2" face="arial, helvetica, sans-serif">S[p/0.01+1,1] <- p</font></div>
<div><font size="2" face="arial, helvetica, sans-serif">S[p/0.01+1,2] <- <span>pp</span>/m</font></div></span></div>
<div><font size="2" face="arial, helvetica, sans-serif">}<var></var></font></div>
<div>
<div><font size="2" face="arial, helvetica, sans-serif"><span>return</span>(S)</font></div></div>
<div><font size="2" face="arial, helvetica, sans-serif">}</font></div>
<div>
<div><font size="2" face="arial, helvetica, sans-serif"></font> </div>
<div><font size="2" face="arial, helvetica, sans-serif"><span>prog</span> <- programa(2,100,1000)</font></div>
<div><font size="2" face="arial, helvetica, sans-serif"><span>plot</span>(<span>prog</span>, <span>xlab</span> = "p", <span>ylab</span> = "Probabilidade de <span>Percolação</span>")</font></div>
<div><font size="2" face="arial, helvetica, sans-serif"></font> </div>
<div><font size="2" face="arial, helvetica, sans-serif">Desde já agradeço pela atenção!</font></div>
<div><font size="2" face="arial, helvetica, sans-serif"></font> </div><font color="#888888">
<div><font size="2" face="arial, helvetica, sans-serif">Tiago Vieira.</font></div></font></div></div></div><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>