Class PascalDistribution
- All Implemented Interfaces:
Serializable,IntegerDistribution
Implementation of the Pascal distribution. The Pascal distribution is a special case of the Negative Binomial distribution where the number of successes parameter is an integer.
There are various ways to express the probability mass and distribution
functions for the Pascal distribution. The present implementation represents
the distribution of the number of failures before r successes occur.
This is the convention adopted in e.g.
MathWorld,
but not in
Wikipedia.
For a random variable X whose values are distributed according to this
distribution, the probability mass function is given by
P(X = k) = C(k + r - 1, r - 1) * p^r * (1 - p)^k,
where r is the number of successes, p is the probability of
success, and X is the total number of failures. C(n, k) is
the binomial coefficient (n choose k). The mean and variance
of X are
E(X) = (1 - p) * r / p, var(X) = (1 - p) * r / p^2.
Finally, the cumulative distribution function is given by
P(X <= k) = I(p, r, k + 1),
where I is the regularized incomplete Beta function.
- Since:
- 1.2 (changed to concrete class in 3.0)
- See Also:
-
Field Summary
Fields inherited from class org.apache.commons.math3.distribution.AbstractIntegerDistribution
random, randomData -
Constructor Summary
ConstructorsConstructorDescriptionPascalDistribution(int r, double p) Create a Pascal distribution with the given number of successes and probability of success.PascalDistribution(RandomGenerator rng, int r, double p) Create a Pascal distribution with the given number of successes and probability of success. -
Method Summary
Modifier and TypeMethodDescriptiondoublecumulativeProbability(int x) For a random variableXwhose values are distributed according to this distribution, this method returnsP(X <= x).intAccess the number of successes for this distribution.doubleUse this method to get the numerical value of the mean of this distribution.doubleUse this method to get the numerical value of the variance of this distribution.doubleAccess the probability of success for this distribution.intAccess the lower bound of the support.intAccess the upper bound of the support.booleanUse this method to get information about whether the support is connected, i.e.doublelogProbability(int x) For a random variableXwhose values are distributed according to this distribution, this method returnslog(P(X = x)), wherelogis the natural logarithm.doubleprobability(int x) For a random variableXwhose values are distributed according to this distribution, this method returnsP(X = x).Methods inherited from class org.apache.commons.math3.distribution.AbstractIntegerDistribution
cumulativeProbability, inverseCumulativeProbability, reseedRandomGenerator, sample, sample, solveInverseCumulativeProbability
-
Constructor Details
-
PascalDistribution
Create a Pascal distribution with the given number of successes and probability of success.Note: this constructor will implicitly create an instance of
Well19937cas random generator to be used for sampling only (seeAbstractIntegerDistribution.sample()andAbstractIntegerDistribution.sample(int)). In case no sampling is needed for the created distribution, it is advised to passnullas random generator via the appropriate constructors to avoid the additional initialisation overhead.- Parameters:
r- Number of successes.p- Probability of success.- Throws:
NotStrictlyPositiveException- if the number of successes is not positiveOutOfRangeException- if the probability of success is not in the range[0, 1].
-
PascalDistribution
public PascalDistribution(RandomGenerator rng, int r, double p) throws NotStrictlyPositiveException, OutOfRangeException Create a Pascal distribution with the given number of successes and probability of success.- Parameters:
rng- Random number generator.r- Number of successes.p- Probability of success.- Throws:
NotStrictlyPositiveException- if the number of successes is not positiveOutOfRangeException- if the probability of success is not in the range[0, 1].- Since:
- 3.1
-
-
Method Details
-
getNumberOfSuccesses
public int getNumberOfSuccesses()Access the number of successes for this distribution.- Returns:
- the number of successes.
-
getProbabilityOfSuccess
public double getProbabilityOfSuccess()Access the probability of success for this distribution.- Returns:
- the probability of success.
-
probability
public double probability(int x) For a random variableXwhose values are distributed according to this distribution, this method returnsP(X = x). In other words, this method represents the probability mass function (PMF) for the distribution.- Parameters:
x- the point at which the PMF is evaluated- Returns:
- the value of the probability mass function at
x
-
logProbability
public double logProbability(int x) For a random variableXwhose values are distributed according to this distribution, this method returnslog(P(X = x)), wherelogis the natural logarithm. In other words, this method represents the logarithm of the probability mass function (PMF) for the distribution. Note that due to the floating point precision and under/overflow issues, this method will for some distributions be more precise and faster than computing the logarithm ofIntegerDistribution.probability(int).The default implementation simply computes the logarithm of
probability(x).- Overrides:
logProbabilityin classAbstractIntegerDistribution- Parameters:
x- the point at which the PMF is evaluated- Returns:
- the logarithm of the value of the probability mass function at
x
-
cumulativeProbability
public double cumulativeProbability(int x) For a random variableXwhose values are distributed according to this distribution, this method returnsP(X <= x). In other words, this method represents the (cumulative) distribution function (CDF) for this distribution.- Parameters:
x- the point at which the CDF is evaluated- Returns:
- the probability that a random variable with this
distribution takes a value less than or equal to
x
-
getNumericalMean
public double getNumericalMean()Use this method to get the numerical value of the mean of this distribution. For number of successesrand probability of successp, the mean isr * (1 - p) / p.- Returns:
- the mean or
Double.NaNif it is not defined
-
getNumericalVariance
public double getNumericalVariance()Use this method to get the numerical value of the variance of this distribution. For number of successesrand probability of successp, the variance isr * (1 - p) / p^2.- Returns:
- the variance (possibly
Double.POSITIVE_INFINITYorDouble.NaNif it is not defined)
-
getSupportLowerBound
public int getSupportLowerBound()Access the lower bound of the support. This method must return the same value asinverseCumulativeProbability(0). In other words, this method must return
The lower bound of the support is always 0 no matter the parameters.inf {x in Z | P(X invalid input: '<'= x) > 0}.- Returns:
- lower bound of the support (always 0)
-
getSupportUpperBound
public int getSupportUpperBound()Access the upper bound of the support. This method must return the same value asinverseCumulativeProbability(1). In other words, this method must return
The upper bound of the support is always positive infinity no matter the parameters. Positive infinity is symbolized byinf {x in R | P(X invalid input: '<'= x) = 1}.Integer.MAX_VALUE.- Returns:
- upper bound of the support (always
Integer.MAX_VALUEfor positive infinity)
-
isSupportConnected
public boolean isSupportConnected()Use this method to get information about whether the support is connected, i.e. whether all integers between the lower and upper bound of the support are included in the support. The support of this distribution is connected.- Returns:
true
-