Estou com um problema e gostaria de saber se algum de vocês poderia me ajudar.
Estou gerando um Processo de Poisson não homogêneo com valores compreendidos entre 0 e 1.
A ideia é que estou pensando em um banco e o 0 representa 10:00 h e o 1 representa 16:00 h.
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.
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?
Atenciosamente,
Rafael Cunha
# lambda(t) = 12 - t , 0 <= t <= 1
# max(lambda) = 12
library(chron)
proc.5 <- function(T = 1){
t <- 0
I <- 0
S <- 0
U <- runif(1)
while(t < T){
t <- t - 1/12*log(U)
V <- runif(1)
if(V <= 12 - t){
I <- I + 1
S[I] <- t
}
}
return(S[-length(S)])
}
x <- proc.5()
times(x)
--
YNWA