[R-br] Ajuda com pacote Chron / Horas

[Rafael Garcia Cunha]

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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110622/571827d6/attachment.html>