Olá,<div><br></div><div>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...</div><div><br></div><div>O problema é que a CDF que obtive não soma 1. Seguem os códigos da PDF e da CDF:</div>

<div><br></div><div><div>dpoix <- function(y, frac, rate, log=FALSE) {</div><div>        b <- y*log(frac)</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>      m <- log(rate)</div><div><span class="Apple-tab-span" style="white-space:pre">      </span>      n <- (y+1)*log(rate+frac)</div>

<div>        exp(b+m-n)</div><div>}</div><div><br></div><div><div>ppoix <- function(y,frac,rate) {</div><div>  a <- function(y,frac,rate) {</div><div>    #a1 <- (frac+rate)^(-y-1)*frac^(y+1)</div><div>    a1 <- exp((-y-1)*log(frac+rate) + (y+1)*log(frac))</div>

<div>    #a2 <- (frac+rate)^(-y-1)*frac^y*rate</div><div>    a2 <- exp((-y-1)*log(frac+rate) + y*log(frac) + log(rate))</div><div>    rate*(-1 + a1 + a2)</div><div>  }</div><div>  b <- function(y,frac,rate) {</div>

<div>    (log(frac) - log(frac+rate))*(frac+rate)</div><div>  }</div><div>  a(y,frac,rate)/b(y,frac,rate)</div><div>}</div></div><div><br></div><div>Alguém poderia dar alguma dica?</div><div><br></div><div>Desde já obrigado.</div>


</div>