[R-br] Dúvida em simulação de processos não homogêneo

[Leonard de Assis]

Se de cara, u<lambda_t(T)/lambda, ele retornará NULL

[]s
Leonard de Assis
assis <dot> leonard <at> gmail <dot> com


Em 29/11/2011 18:49, Victor Eduardo escreveu:
> Pessoal, estou com dúvida no script abaixo, pois as vezes fica me
> dando como resposta NULL e não estou conseguindo identificar o erro. O
> meu objetivo é Implementar uma funçao que, para um dado T , retorna
> os instantes das ocorrências de um processo de Poisson nao-homogeneo
> com funcao de intensidade lambda(t) = 3/(t + 1)
>
> lambda_t<-function(t){
> if(t<0){
> print("tempo negativo")
> break()
> }
> lambdat<-3/(t+1)
> return(lambdat)
> }
>
> pois_naohomog<-function(){
>         lambda<-3
> i<-1
> T<-0
> saida<-NULL
> T<-T+exponencial(lambda)
> u<-runif(1)
> while(u<lambda_t(T)/lambda){
> saida[i]<-T
> i<-i+1
> T<-T+exponencial(lambda)
> }
> return(saida)
> }
>
>
> _______________________________________________
> R-br mailing list
> R-br em 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.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111129/59c85c73/attachment.html>