<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Se de cara, <tt>u<lambda_t(T)/lambda</tt>, ele retornará NULL<br>
<pre class="moz-signature" cols="72">[]s
Leonard de Assis
assis <dot> leonard <at> gmail <dot> com
</pre>
<br>
Em 29/11/2011 18:49, Victor Eduardo escreveu:
<blockquote
cite="mid:CAA=hdgYzMsRZkER++KXp_GD2hBKLR6wWuhS8GGQP4ZTyQ-a1JQ@mail.gmail.com"
type="cite">Pessoal, estou com dúvida no script abaixo, pois as
vezes fica me dando como resposta NULL e não estou conseguindo
identificar o erro. O meu objetivo é Implementar uma funçao que,
para um dado T , retorna os instantes das ocorrências de um
processo de Poisson nao-homogeneo com funcao de intensidade
lambda(t) = 3/(t + 1)
<div>
<br>
</div>
<div>
<div>lambda_t<-function(t){</div>
<div><span class="Apple-tab-span" style="white-space:pre"> </span>if(t<0){</div>
<div><span class="Apple-tab-span" style="white-space:pre"> </span>print("tempo
negativo")</div>
<div><span class="Apple-tab-span" style="white-space:pre"> </span>break()</div>
<div><span class="Apple-tab-span" style="white-space:pre"> </span>}</div>
<div><span class="Apple-tab-span" style="white-space:pre"> </span>lambdat<-3/(t+1)</div>
<div><span class="Apple-tab-span" style="white-space:pre"> </span>return(lambdat)</div>
<div>}</div>
<div><br>
</div>
<div>pois_naohomog<-function(){</div>
<div> lambda<-3</div>
<div><span class="Apple-tab-span" style="white-space:pre"> </span>i<-1</div>
<div><span class="Apple-tab-span" style="white-space:pre"> </span>T<-0</div>
<div><span class="Apple-tab-span" style="white-space:pre"> </span>saida<-NULL</div>
<div><span class="Apple-tab-span" style="white-space:pre"> </span>T<-T+exponencial(lambda)</div>
<div><span class="Apple-tab-span" style="white-space:pre"> </span>u<-runif(1)</div>
<div><span class="Apple-tab-span" style="white-space:pre"> </span>while(u<lambda_t(T)/lambda){</div>
<div><span class="Apple-tab-span" style="white-space:pre"> </span>saida[i]<-T</div>
<div><span class="Apple-tab-span" style="white-space:pre"> </span>i<-i+1</div>
<div><span class="Apple-tab-span" style="white-space:pre"> </span>T<-T+exponencial(lambda)<span
class="Apple-tab-span" style="white-space:pre"> </span></div>
<div><span class="Apple-tab-span" style="white-space:pre"> </span>}</div>
<div><span class="Apple-tab-span" style="white-space:pre"> </span>return(saida)<span
class="Apple-tab-span" style="white-space:pre"> </span></div>
<div>
}</div>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
</body>
</html>