<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=windows-1252">
<META content="MSHTML 6.00.2900.3429" name=GENERATOR>
<STYLE></STYLE>
</HEAD>
<BODY bgColor=#ffffff>
<DIV><FONT face=Arial color=#0000ff>Muito obrigado pela dica. Mas como fica o
faça X=N-1</FONT></DIV>
<BLOCKQUOTE
style="PADDING-RIGHT: 0px; PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: #000000 2px solid; MARGIN-RIGHT: 0px">
<DIV style="FONT: 10pt arial"> </DIV>
<DIV style="FONT: 10pt arial"> </DIV>
<DIV style="FONT: 10pt arial"> </DIV>
<DIV style="FONT: 10pt arial">for não é para isso, da uma olhada no
while</DIV>
<DIV><BR></DIV>
<DIV>Exemplo basico:</DIV>
<DIV><BR></DIV>
<DIV>
<DIV>
<DIV>## Gera X até que a soma deles seja 1.5</DIV>
<DIV>x <- runif(2) </DIV>
<DIV>x</DIV>
<DIV>n<- 1</DIV>
<DIV>while (sum(x) < 1.5) { x <- runif(2) ; n <- n+1 }</DIV>
<DIV>x</DIV>
<DIV>n # numero de tentativas</DIV></DIV></DIV>
<DIV><BR></DIV>
<DIV>
<DIV><BR></DIV>
<DIV><BR>
<DIV class=gmail_quote>2011/8/23 Crysttian Arantes <SPAN dir=ltr><<A
href="mailto:crysttian@gmail.com">crysttian@gmail.com</A>></SPAN><BR>
<BLOCKQUOTE class=gmail_quote
style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">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: <A
href="tel:35%208876-9187" target=_blank value="+553588769187">35
8876-9187</A> (Oi)<BR> <A
href="tel:35%209189-7272" target=_blank value="+553591897272">35
9189-7272</A> (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"
target=_blank>sznelwar@uol.com.br</A>></SPAN> escreveu:<BR>
<BLOCKQUOTE class=gmail_quote
style="PADDING-LEFT: 1ex; MARGIN: 0pt 0pt 0pt 0.8ex; BORDER-LEFT: rgb(204,204,204) 1px solid">
<DIV>
<DIV></DIV>
<DIV class=h5><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></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.<BR></BLOCKQUOTE></DIV><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 mínimo
reproduzível.<BR></BLOCKQUOTE></DIV><BR></DIV></DIV>
<P>
<HR>
<P></P>_______________________________________________<BR>R-br mailing
list<BR>R-br@listas.c3sl.ufpr.br<BR>https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br<BR>Leia
o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo
reproduzível.</BLOCKQUOTE></BODY></HTML>