R/poisUNIF.R

In plpoisson: Prediction Limits for Poisson Distribution

Defines functions poisUNIF

Documented in poisUNIF

poisUNIF <-
function(xobs, n, s, t, alpha = 0.05)

{
  ns <- n * s
  nst <- ns + t
  prob <- ns / nst

  res <- list()
  res$lower <- qnbinom(alpha, xobs + 1, prob) + 1
  res$upper <- qnbinom(alpha, xobs + 1, prob, lower.tail = FALSE) - 1
  return(res)
}

Try the plpoisson package in your browser