Olá, tenho o código de uma PDF a partida da qual deveria ser possível, integrando, obter a CDF, certo pessoal? E esta CDF deveria somar 1 sobre todo o intervalo... O problema é que a CDF que obtive não soma 1. Seguem os códigos da PDF e da CDF: dpoix <- function(y, frac, rate, log=FALSE) { b <- y*log(frac) m <- log(rate) n <- (y+1)*log(rate+frac) exp(b+m-n) } ppoix <- function(y,frac,rate) { a <- function(y,frac,rate) { #a1 <- (frac+rate)^(-y-1)*frac^(y+1) a1 <- exp((-y-1)*log(frac+rate) + (y+1)*log(frac)) #a2 <- (frac+rate)^(-y-1)*frac^y*rate a2 <- exp((-y-1)*log(frac+rate) + y*log(frac) + log(rate)) rate*(-1 + a1 + a2) } b <- function(y,frac,rate) { (log(frac) - log(frac+rate))*(frac+rate) } a(y,frac,rate)/b(y,frac,rate) } Alguém poderia dar alguma dica? Desde já obrigado.