<div dir="ltr">Pode usar a fórmula de Stirling para aproximar um fatorial, n grande.<br></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr">Atenciosamente,<br><br>Prof. Clécio da Silva Ferreira<br>Depto de Estatística - UFJF<br>Tel: (32) 2102-3306<br>Página: <a href="http://www.ufjf.br/clecio_ferreira" target="_blank">www.ufjf.br/clecio_ferreira</a><br><br></div></div></div>
<br><div class="gmail_quote">Em 5 de novembro de 2014 23:27, Marcus Nunes <span dir="ltr"><<a href="mailto:marcus.nunes@gmail.com" target="_blank">marcus.nunes@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Eu faria esta conta na mão. Espero que minha notação esteja clara. Considere log(x) como o logaritmo natural de x.<div><br></div><div>Como n! = n*(n-1)*...*3*2*1, log(n!) = log(n) + log(n-1) + ... + log(2) + log(1)</div><div><br></div><div>Assim, log(factorial(365)) = sum(log(365:1))</div><div><br></div><div>Além disso, log(365^n) = n*log(365)</div><div><br></div><div>log(factorial(365-n)) = sum(log((365-n):1))</div><div><br></div><div>Portanto, o logaritmo da expressão original torna-se</div><div><br></div><div>sum(log(365:1)) - n*log(365) - sum(log((365-n):1))<br></div><div><br></div><div>e o valor procurado é</div><div><br></div><div>exp( sum(log(365:1)) - n*log(365) - sum(log((365-n):1)) ) = 0.1087682</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">2014-11-05 23:00 GMT-02:00 Andre Oliveira <span dir="ltr"><<a href="mailto:andreolsouza@yahoo.com.br" target="_blank">andreolsouza@yahoo.com.br</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div><div style="color:#000;background-color:#fff;font-family:HelveticaNeue-Light,Helvetica Neue Light,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif;font-size:10px"><div dir="ltr"><span>Pessoas,</span></div><div dir="ltr"><span>como posso otimizar este cálculo? Com 365 o R nem faz o cálculo! Alguma sugestão? </span></div><div dir="ltr"><span><br></span></div><div dir="ltr"><span>Obrigado! </span></div><div dir="ltr"><span><br></span></div><div></div><div> </div><div><div>
<div>n=40</div><div>np=factorial(365)/((365^n)*factorial(365-n))</div><div>1-np</div><div><br></div><div><br></div><div><br></div>
<br>André Oliveira Souza. 
<br>Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo.  IFES
<br>
<br></div><div> </div></div></div></div><br></div></div>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">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.<span class="HOEnZb"><font color="#888888"><br></font></span></blockquote></div><span class="HOEnZb"><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div><div dir="ltr">Marcus Nunes<br><a href="http://marcusnunes.me/" target="_blank">http://marcusnunes.me/</a><br></div></div>
</font></span></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></div>