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

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@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.