dgcnt.Rd 1.93 KB
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/gcnt.R
\name{dgcnt}
\alias{dgcnt}
\title{Probabilidades do Modelo Gamma Count}
\usage{
dgcnt(y, lambda, alpha)
}
\arguments{
\item{y}{Valor da variável de contagem.}

\item{lambda}{Valor do parâmetro \eqn{\lambda} que é a média da
distribuição do tempo entre eventos.}

\item{alpha}{Valor do parâmetro \eqn{\alpha} que é o parâmetro de
dispersão.}
}
\value{
Retorna uma probabilidade, ou seja \eqn{\Pr(Y = y | \lambda,
    \alpha) = p(y, \lambda, \alpha)}.
}
\description{
Calcula as probabilidades para uma variável aleatória
    com distribuição Gamma Count:

\deqn{p(y,\lambda,\alpha) =
  \left(\int_{0}^{1}
  \frac{(\alpha\lambda)^{y\alpha}}{\Gamma(y\alpha)}\,
  u^{y\alpha-1}
  \exp\{-\alpha\lambda u\}\, \textrm{d}u \right)
  - \left(\int_{0}^{1}
  \frac{(\alpha\lambda)^{y\alpha}}{\Gamma((y+1)\alpha)}\,
  u^{(y+1)\alpha-1}
  \exp\{-\alpha\lambda u\}\, \textrm{d}u \right),}

em que \eqn{\lambda > 0} é a média da variável aleatória tempo entre
    eventos e \eqn{\alpha > 0} é o parâmetro de dispersão.
}
\examples{

dpois(5, lambda = 5)
dgcnt(5, lambda = 5, alpha = 1)

probs <- data.frame(y = 0:30)
within(probs, {
    py0 <- dpois(y, lambda = 15)
    py1 <- dgcnt(y, lambda = 15, alpha = 1)
    py2 <- dgcnt(y, lambda = 15, alpha = 0.5)
    py3 <- dgcnt(y, lambda = 15, alpha = 1.5)
    plot(py0 ~ y, type = "h",
         ylim = c(0, max(c(py0, py2, py3))),
         ylab = expression(Pr(Y == y)))
    points(y + 0.1, py1, type = "h", col = 2)
    points(y - 0.3, py2, type = "h", col = 3)
    points(y + 0.3, py3, type = "h", col = 4)
    legend("topleft", bty = "n",
           col = c(1:4), lty = 1,
           legend = expression(
               Poisson(lambda == 15),
               GC(lambda == 15, alpha == 1),
               GC(lambda == 15, alpha == 0.5),
               GC(lambda == 15, alpha == 1.5)))
})

}
\author{
Walmes Zeviani, \email{walmes@ufpr.br}.
}