<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=utf-8">
<META content="MSHTML 6.00.2900.3429" name=GENERATOR>
<STYLE></STYLE>
</HEAD>
<BODY>
<DIV><FONT face=Arial size=2><FONT color=#0000ff size=4>Como faz para
executar?</FONT><BR></FONT></DIV><FONT face=Arial size=2>>O código dentro do
while não está sendo executado:<BR>> <BR>>> exp(-2)<BR>> [1]
0.1353353<BR>>> sum(runif(100))<BR>> [1] 53.06867<BR>> <BR>> logo
sum(u)<c dificilmente vai dar TRUE e n vai continuar sendo 1, o<BR>> que
faz com que x = n - 1 = 0.<BR>> <BR>> 2011/8/24 Mauro Sznelwar
<</FONT><A href="mailto:sznelwar@uol.com.br"><FONT face=Arial
size=2>sznelwar@uol.com.br</FONT></A><FONT face=Arial size=2>>:<BR>>>
Por que o X dá sempre zero, mesmo que eu substitua o lambda por
qualquer<BR>>> valor? Onde estou errando?p=1<BR>>> n=0<BR>>>
lambda=2<BR>>> c=exp(-lambda)<BR>>> u=runif(100)<BR>>>
p=p*u<BR>>> n=1<BR>>> #while (p<c){
u=runif(100);n=n+1}<BR>>> #p<BR>>> while (sum(u)<c){
u=runif(100);n=n+1}<BR>>> u<BR>>> x=n-1<BR>>>
x<BR>>><BR>>> É só trocar o nome das variaveis (na verdade nem
precisava trocar, era só<BR>>> adicionar x<-n-1 no
final)...<BR>>> ## Gera 'batata' até que a soma seja 1.5<BR>>>
batata <- runif(2)<BR>>> batata<BR>>> n<- 1<BR>>> while
(sum(batata) < 1.5) { batata <- runif(2) ; n <- n+1 }<BR>>>
batata<BR>>> x<-n-1<BR>>> 2011/8/23 Mauro Sznelwar <</FONT><A
href="mailto:sznelwar@uol.com.br"><FONT face=Arial
size=2>sznelwar@uol.com.br</FONT></A><FONT face=Arial
size=2>><BR>>>><BR>>>> Muito obrigado pela dica. Mas como
fica o faça
X=N-1<BR>>>><BR>>>><BR>>>><BR>>>><BR>>>>
for não é para isso, da uma olhada no while<BR>>>> Exemplo
basico:<BR>>>> ## Gera X até que a soma deles seja 1.5<BR>>>>
x <- runif(2)<BR>>>> x<BR>>>> n<- 1<BR>>>>
while (sum(x) < 1.5) { x <- runif(2) ; n <- n+1 }<BR>>>>
x<BR>>>> n # numero de
tentativas<BR>>>><BR>>>><BR>>>> 2011/8/23 Crysttian
Arantes <</FONT><A href="mailto:crysttian@gmail.com"><FONT face=Arial
size=2>crysttian@gmail.com</FONT></A><FONT face=Arial
size=2>><BR>>>>><BR>>>>>
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 é<BR>>>>> 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>>>>>
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<BR>>>>> outro talento.”
[Isaac Newton]<BR>>>>><BR>>>>> "Otimismo é esperar pelo
melhor. Confiança é saber lidar com o pior."<BR>>>>> [Roberto
Simonsen]<BR>>>>><BR>>>>> "Mesmo as noites totalmente
sem estrelas podem anunciar a aurora de uma<BR>>>>> grande
realização." [Martin Luther
King]<BR>>>>><BR>>>>><BR>>>>><BR>>>>><BR>>>>><BR>>>>><BR>>>>>
Em 23 de agosto de 2011 22:59, Mauro Sznelwar <</FONT><A
href="mailto:sznelwar@uol.com.br"><FONT face=Arial
size=2>sznelwar@uol.com.br</FONT></A><FONT face=Arial
size=2>><BR>>>>>
escreveu:<BR>>>>>><BR>>>>>> Tenho este
algoritmo<BR>>>>>><BR>>>>>> p=1, N=0,
c=e^-lambda<BR>>>>>><BR>>>>>>
Repita<BR>>>>>><BR>>>>>>
N=N+1<BR>>>>>><BR>>>>>> gere
Ui<BR>>>>>><BR>>>>>> atualize
p<c.e^-lambda<BR>>>>>><BR>>>>>> Faça
X=N-1<BR>>>>>><BR>>>>>><BR>>>>>>
Estou fazendo esta rotina abaixo, mas o resultado não está
consistente!<BR>>>>>> Qual o melhor
jeito?<BR>>>>>> 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<BR>>>>>><BR>>>>>>
_______________________________________________<BR>>>>>> R-br
mailing list<BR>>>>>> </FONT><A
href="mailto:R-br@listas.c3sl.ufpr.br"><FONT face=Arial
size=2>R-br@listas.c3sl.ufpr.br</FONT></A><BR><FONT face=Arial
size=2>>>>>> </FONT><A
href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br"><FONT face=Arial
size=2>https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</FONT></A><BR><FONT
face=Arial size=2>>>>>> Leia o guia de postagem (</FONT><A
href="http://www.leg.ufpr.br/r-br-guia"><FONT face=Arial
size=2>http://www.leg.ufpr.br/r-br-guia</FONT></A><FONT face=Arial size=2>) e
forneça<BR>>>>>> código mínimo
reproduzível.<BR>>>>><BR>>>>><BR>>>>>
_______________________________________________<BR>>>>> R-br mailing
list<BR>>>>> </FONT><A href="mailto:R-br@listas.c3sl.ufpr.br"><FONT
face=Arial size=2>R-br@listas.c3sl.ufpr.br</FONT></A><BR><FONT face=Arial
size=2>>>>> </FONT><A
href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br"><FONT face=Arial
size=2>https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</FONT></A><BR><FONT
face=Arial size=2>>>>> Leia o guia de postagem (</FONT><A
href="http://www.leg.ufpr.br/r-br-guia"><FONT face=Arial
size=2>http://www.leg.ufpr.br/r-br-guia</FONT></A><FONT face=Arial size=2>) e
forneça<BR>>>>> código mínimo
reproduzível.<BR>>>><BR>>>>
________________________________<BR>>>><BR>>>>
_______________________________________________<BR>>>> R-br mailing
list<BR>>>> </FONT><A href="mailto:R-br@listas.c3sl.ufpr.br"><FONT
face=Arial size=2>R-br@listas.c3sl.ufpr.br</FONT></A><BR><FONT face=Arial
size=2>>>> </FONT><A
href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br"><FONT face=Arial
size=2>https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</FONT></A><BR><FONT
face=Arial size=2>>>> Leia o guia de postagem (</FONT><A
href="http://www.leg.ufpr.br/r-br-guia"><FONT face=Arial
size=2>http://www.leg.ufpr.br/r-br-guia</FONT></A><FONT face=Arial size=2>) e
forneça<BR>>>> código mínimo
reproduzível.<BR>>>><BR>>>>
_______________________________________________<BR>>>> R-br mailing
list<BR>>>> </FONT><A href="mailto:R-br@listas.c3sl.ufpr.br"><FONT
face=Arial size=2>R-br@listas.c3sl.ufpr.br</FONT></A><BR><FONT face=Arial
size=2>>>> </FONT><A
href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br"><FONT face=Arial
size=2>https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</FONT></A><BR><FONT
face=Arial size=2>>>> Leia o guia de postagem (</FONT><A
href="http://www.leg.ufpr.br/r-br-guia"><FONT face=Arial
size=2>http://www.leg.ufpr.br/r-br-guia</FONT></A><FONT face=Arial size=2>) e
forneça<BR>>>> código mínimo reproduzível.<BR>>><BR>>>
________________________________<BR>>><BR>>>
_______________________________________________<BR>>> R-br mailing
list<BR>>> </FONT><A href="mailto:R-br@listas.c3sl.ufpr.br"><FONT
face=Arial size=2>R-br@listas.c3sl.ufpr.br</FONT></A><BR><FONT face=Arial
size=2>>> </FONT><A
href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br"><FONT face=Arial
size=2>https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</FONT></A><BR><FONT
face=Arial size=2>>> Leia o guia de postagem (</FONT><A
href="http://www.leg.ufpr.br/r-br-guia"><FONT face=Arial
size=2>http://www.leg.ufpr.br/r-br-guia</FONT></A><FONT face=Arial size=2>) e
forneça código<BR>>> mínimo reproduzível.<BR>>><BR>>>
_______________________________________________<BR>>> R-br mailing
list<BR>>> </FONT><A href="mailto:R-br@listas.c3sl.ufpr.br"><FONT
face=Arial size=2>R-br@listas.c3sl.ufpr.br</FONT></A><BR><FONT face=Arial
size=2>>> </FONT><A
href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br"><FONT face=Arial
size=2>https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</FONT></A><BR><FONT
face=Arial size=2>>> Leia o guia de postagem (</FONT><A
href="http://www.leg.ufpr.br/r-br-guia"><FONT face=Arial
size=2>http://www.leg.ufpr.br/r-br-guia</FONT></A><FONT face=Arial size=2>) e
forneça código<BR>>> mínimo reproduzível.<BR>>><BR>>
_______________________________________________<BR>> R-br mailing
list<BR>> </FONT><A href="mailto:R-br@listas.c3sl.ufpr.br"><FONT face=Arial
size=2>R-br@listas.c3sl.ufpr.br</FONT></A><BR><FONT face=Arial size=2>>
</FONT><A href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br"><FONT
face=Arial
size=2>https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</FONT></A><BR><FONT
face=Arial size=2>> Leia o guia de postagem (</FONT><A
href="http://www.leg.ufpr.br/r-br-guia"><FONT face=Arial
size=2>http://www.leg.ufpr.br/r-br-guia</FONT></A><FONT face=Arial size=2>) e
forneça código mínimo reproduzível.</FONT></BODY></HTML>