[R-br] Obtendo a CDF a partir de uma PDF

Cristiano Strieder cstrieder em gmail.com
Domingo Junho 19 23:38:12 BRT 2011


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.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110619/fbe2a257/attachment.html>


Mais detalhes sobre a lista de discussão R-br