# plpoisson-package: Prediction Limits for Poisson Distribution In plpoisson: Prediction Limits for Poisson Distribution

## Description

Prediction limits for Poisson distribution are useful when quantifying the uncertainty associated with predicting the occurrences of real life phenomena. The plpoisson package provides a set of functions to compute prediction limits of the inferred Poisson distribution under both, frequentist and Bayesian frameworks.

For frequentist prediction a common approach is to estimate the parameter based on the observed data firstly, then, to predict based on the estimated parameter. Different from the common approach of frequentist prediction, this approach does not require the estimation of the parameter. In a Bayesian setting, the uniform, Jeffreys and gamma distributions are used as priors when deriving the predictive posterior distribution.

## Details

 Package: plpoisson Type: Package Version: 0.2.0 Date: 2021-02-14 License: GPL-3

For a complete list of exported functions, use `library(help = "plpoisson")`.

## Author(s)

Valbona Bejleri, Luca Sartore and Balgobin Nandram

Maintainer: Luca Sartore drwolf85@gmail.com

## References

Bejleri, V., & Nandram, B. (2018). Bayesian and frequentist prediction limits for the Poisson distribution. Communications in Statistics-Theory and Methods, 47(17), 4254-4271.

Bejleri, V. (2005). Bayesian Prediction Intervals for the Poisson Model, Noninformative Priors, Ph.D. Dissertation, American University, Washington, DC.

## Examples

 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23``` ```## Loading the package library(plpoisson) ## Setting quantities of interest xobs <- rpois(1, 50) # Number of the observed occurrencies n <- 1 # Total number of the time windows of # of size 's' observed in the past s <- rgamma(1, 4, .567) # Fixed size of observed time windows t <- rgamma(1, 3, .33) # Future time window a <- 5 # Shape hyperparameter of a gamma prior b <- 1.558 # Rate hyperparameter of a gamma prior ## Frequentist prediction limits poiss(xobs, n, s, t) ## Bayesian prediction limits (with uniform prior) poisUNIF(xobs, n, s, t) ## Bayesian prediction limits (with Jeffreys prior) poisJEFF(xobs, n, s, t) ## Bayesian prediction limits (with gamma prior) poisBayes(xobs, n, s, t, a, b) ```

plpoisson documentation built on Feb. 15, 2021, 1:08 a.m.