Olá,<br><br>O problema está na funcao runif(0).<br><br>O zero indica que o numero de valores uniformes que ele tem que gerar é zero, coloca 1.<br><br>runif(1):  gera 1 valor<br>runif(2):  gera 2 valores <br><br>Entendeu:<br>
<br>Até<br><br><br><br clear="all">Crysttian Arantes Paixão<br>Universidade Federal de Lavras - UFLA<br>Cel: 35 8876-9187 (Oi)<br>       35 9189-7272 (Tim)<br>Linux User: 456734 - Debian<br><br>“Se fiz descobertas valiosas, foi mais por ter paciência do que qualquer outro talento.” [Isaac Newton]<br>
<br><font face="Arial, Helvetica, sans-serif" size="2">"Otimismo 
                      é esperar pelo melhor. Confiança é 
                      saber lidar com o pior."  [Roberto Simonsen]<br><br></font><font face="Arial, Helvetica, sans-serif" size="2">"Mesmo 
                      as noites totalmente sem estrelas podem anunciar a aurora 
                      de uma grande realização." [Martin Luther 
                      King]<br><br></font><br><font face="Arial, Helvetica, sans-serif" size="2"><br></font><br>
<br><br><div class="gmail_quote">Em 23 de agosto de 2011 22:59, Mauro Sznelwar <span dir="ltr"><<a href="mailto:sznelwar@uol.com.br">sznelwar@uol.com.br</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<u></u>





<div bgcolor="#ffffff">
<div><font face="Arial">
<p class="MsoNormal" style="margin: 0cm 0cm 0pt; text-align: justify;"><span style="font-size: 12pt;"><font size="2">Tenho este algoritmo<u></u><u></u><u></u></font></span></p>
<p class="MsoNormal" style="margin: 0cm 0cm 0pt; text-align: justify;"><span style="font-size: 12pt;"><font size="2">p=1, N=0, c=e^-lambda<u></u><u></u></font></span></p>
<p class="MsoNormal" style="margin: 0cm 0cm 0pt; text-align: justify;"><span style="font-size: 12pt;"><font size="2">Repita <u></u><u></u></font></span></p>
<p class="MsoNormal" style="margin: 0cm 0cm 0pt; text-align: justify;"><span style="font-size: 12pt;"><font size="2"><span> </span>N=N+1<u></u><u></u></font></span></p>
<p class="MsoNormal" style="margin: 0cm 0cm 0pt; text-align: justify;"><span style="font-size: 12pt;"><font size="2">gere Ui<u></u><u></u></font></span></p>
<p class="MsoNormal" style="margin: 0cm 0cm 0pt; text-align: justify;"><span style="font-size: 12pt;"><font size="2">atualize p<c.e^-lambda<u></u><u></u></font></span></p>
<p class="MsoNormal" style="margin: 0cm 0cm 0pt; text-align: justify;"><span style="font-size: 12pt;"><font size="2">Faça X=N-1<u></u><u></u></font></span></p></font></div>
<div><font face="Arial"></font> </div>
<div><font face="Arial">Estou fazendo esta rotina abaixo, mas o resultado não está 
consistente! Qual o melhor jeito?</font></div>
<div><font face="Arial" size="2"><font size="3">p=1<br>n=0<br>lambda=2<br>c=exp(-lambda)<br> for(n in 1:100){ 
<br> u=runif(n)<br> p=p*u<br>p<c <br>}<br>x=n-1</font><br> </font></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>