<div>Entendi sua explanação. heehhe</div><div> </div><div>olhei o código bem por alto depois do almoço num tempó livre, única coisa "estranha" que notei foi esse loop que pra mim estava sem muito sentido.</div><div>
</div><div>Mas muito dessa minha opinião vai do fato de não conhecer ao fundo o problema, pois para mim, ele estaria interessado em cada uma das iterações (eu, quando faço simulações, procuro guardar cada uma delas, nem que seja pra brincar de bootstrap depois, hehehe)</div>
<div> </div><div>Quanto à copa américa, estava meio na cara que daria Uruguai, é o time com mais conjunto, Argentina e Brasil estavam em início de trabalho. No betclick eu apostei no uruguai e salvei algum, ehhehe</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 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">
Leonard,<br>
<br>
para cada 'g', ele determina uma realizacao (aleatoria) de Y_n<br>
condicionado em Y_{n-1}. Se essa realizacao por nao-nula, ele<br>
incrementa 'pp'. Como a chance de Y_m ser zero e' bem pequena, ele<br>
repete o laco 'm' vezes, contando quantas vezes (em 'm') ele observa<br>
Y_n > 0...<br>
<br>
Estou bem certo de que o codigo esta' correto... mas eu tambem estava<br>
certo que o Brasil ganharia a Copa America...<br>
<br>
b<br>
<br>
2011/8/18 Leonard Assis <<a href="mailto:assis.leonard@gmail.com">assis.leonard@gmail.com</a>>:<br>
<div><div></div><div class="h5">> Benilton, eu nem cheguei a rodar, analisei o que ele faz mesmo.<br>
><br>
> Tem que resolver esse problema daquele loop que só faz perder tempo senão o<br>
> próprio estudo de simulação dele está inválido.<br>
> lmassis <at> yahoo <dot> com <dot> br<br>
> assis.leonard <at> gmail <dot> com<br>
><br>
><br>
> 2011/8/18 Benilton Carvalho <<a href="mailto:beniltoncarvalho@gmail.com">beniltoncarvalho@gmail.com</a>><br>
>><br>
>> 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>
>> ## f: número máximo de filhos<br>
>> ## n: número de gerações<br>
>> ## m: número de iterações<br>
>> S <- foreach(p=seq(0,1,by=0.01), .combine=rbind) %dopar% {<br>
>> pp <- 0<br>
>> for (g in 1:m) {<br>
>> y <- 1<br>
>> for (k in 1:n)<br>
>> 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<br>
>> código mínimo reproduzível.<br>
><br>
><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<br>
> mínimo reproduzível.<br>
><br>
<br>
<br>
<br>
</div></div><font color="#888888">--<br>
Successful people ask better questions, and as a result, they get<br>
better answers. (Tony Robbins)<br>
</font><div><div></div><div class="h5">_______________________________________________<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>
</div></div></blockquote></div><br>