Boa noite!<div><br><div>Estou com um problema e gostaria de saber se algum de vocês poderia me ajudar.</div><div><br></div><div>Estou gerando um Processo de Poisson não homogêneo com valores compreendidos entre 0 e 1.</div>
<div>A ideia é que estou pensando em um banco e o 0 representa 10:00 h e o 1 representa 16:00 h.</div><div>Gero os valores de 0 até 1 [que representam os instantes de ocorrências do processo de Poisson] e aplico a função "times" do pacote "chron" para que transforme esses valores em horas. Só que as horas geradas estão compreendidas ao longo de um dia, de 00:00 h até 23:59.</div>
<div>Alguém sabe alguma maneira de colocar algo como "from = 10:00:00 to = 16:00:00" no pacote chron? Ou alguma outra maneira que eu possa resolver esse problema?</div><div><br></div><div>Atenciosamente,</div><div>
<br></div><div>Rafael Cunha</div><div><br></div><div><div># lambda(t) = 12 - t , 0 <= t <= 1</div><div># max(lambda) = 12</div><div><br></div><div>library(chron)</div><div><br></div><div>proc.5 <- function(T = 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>I <- 0</div><div><span class="Apple-tab-span" style="white-space:pre">    </span>S <- 0</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(t < T){</div><div><span class="Apple-tab-span" style="white-space:pre">             </span>t <- t - 1/12*log(U)</div>
<div><span class="Apple-tab-span" style="white-space:pre">              </span>V <- runif(1)</div><div><span class="Apple-tab-span" style="white-space:pre">             </span>if(V <= 12 - 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>S[I] <- t</div><div><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(S[-length(S)])</div><div>} </div><div><br></div><div>x <- proc.5()</div><div>times(x)</div><br>-- <br>YNWA<br><br>
</div></div>