Só transformar a tua variavel y<- x*6/24 + 10/24 times(y) 2011/6/22 Rafael Garcia Cunha <rafareds@gmail.com>
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 R-br@listas.c3sl.ufpr.br 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.