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