Boa noite!
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
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
