Só transformar a tua variavel y<- x*6/24 + 10/24 times(y)
2011/6/22 Rafael Garcia Cunha <[email protected]> > 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. >
_______________________________________________ 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.
