SkellamDistribution

SkellamDistribution[μ1,μ2]
represents a Skellam distribution with shape parameters μ1 and μ2.

DetailsDetails

Background & Context
Background & Context

  • SkellamDistribution[μ1,μ2] represents a discrete statistical distribution defined for integer values , which is determined by the positive real parameters μ1 and μ2 and is defined as the distribution of the difference XY1-Y2, where Y1PoissonDistribution[μ1] and Y2PoissonDistribution[μ2] are independent variates with means μ1 and μ2, respectively. The Skellam distribution has a probability density function (PDF) that is discrete and unimodal and whose overall shape (its height, its spread, and the horizontal location of its maximum) is determined by the values of μ1 and μ2.
  • The Skellam distribution was first derived by J. G. Skellam in the mid-1940s as the theoretical distribution modeling the difference of two Poisson-distributed variates from different populations. Since then, the distribution has been used in image analysis to study image differences in the presence of ambient noise as well as differences in the number of accidents between two locations and/or two periods of time. The Skellam distribution has also been used to model phenomena including sports scoring, stock prices, queueing models, and gene expression.
  • RandomVariate can be used to give one or more machine- or arbitrary-precision (the latter via the WorkingPrecision option) pseudorandom variates from a Skellam distribution. Distributed[x,SkellamDistribution[μ1,μ2]], written more concisely as xSkellamDistribution[μ1,μ2], can be used to assert that a random variable x is distributed according to a Skellam distribution. Such an assertion can then be used in functions such as Probability, NProbability, Expectation, and NExpectation.
  • The probability density and cumulative distribution functions may be given using PDF[SkellamDistribution[μ1,μ2],x] and CDF[SkellamDistribution[μ1,μ2],x], though one should note that there is no closed-form expression for its PDF. The mean, median, variance, raw moments, and central moments may be computed using Mean, Median, Variance, Moment, and CentralMoment, respectively. These quantities can be visualized using DiscretePlot.
  • DistributionFitTest can be used to test if a given dataset is consistent with a Skellam distribution, EstimatedDistribution to estimate a Skellam parametric distribution from given data, and FindDistributionParameters to fit data to a Skellam distribution. ProbabilityPlot can be used to generate a plot of the CDF of given data against the CDF of a symbolic Skellam distribution and QuantilePlot to generate a plot of the quantiles of given data against the quantiles of a symbolic Skellam distribution.
  • TransformedDistribution can be used to represent a transformed Skellam distribution, CensoredDistribution to represent the distribution of values censored between upper and lower values, and TruncatedDistribution to represent the distribution of values truncated between upper and lower values. CopulaDistribution can be used to build higher-dimensional distributions that contain a Skellam distribution, and ProductDistribution can be used to compute a joint distribution with independent component distributions involving Skellam distributions.
  • SkellamDistribution is related to a number of other statistical distributions. It is a transformation (TransformedDistribution) of PoissonDistribution in the sense that SkellamDistribution[μ1,μ2] is equivalent to the distribution of x-y where xPoissonDistribution[μ1] and yPoissonDistribution[μ2]. SkellamDistribution is also closely related to PoissonConsulDistribution, CompoundPoissonDistribution, and PolyaAeppliDistribution.

ExamplesExamplesopen allclose all

Basic Examples  (3)Basic Examples  (3)

Probability mass function:

In[1]:=
Click for copyable input
Out[1]=
In[2]:=
Click for copyable input
Out[2]=

Cumulative distribution function:

In[1]:=
Click for copyable input
Out[1]=

Mean and variance:

In[1]:=
Click for copyable input
Out[1]=
In[2]:=
Click for copyable input
Out[2]=
Introduced in 2010
(8.0)
Translate this page: