scipy.stats.binom

scipy.stats.binom()

A binom discrete random variable.

Discrete random variables are defined from a standard form and may require some shape parameters to complete its specification. Any optional keyword parameters can be passed to the methods of the RV object as given below:

Methods:

binom.rvs(n,pr,loc=0,size=1) :

• random variates

binom.pmf(x,n,pr,loc=0) :

• probability mass function

binom.cdf(x,n,pr,loc=0) :

• cumulative density function

binom.sf(x,n,pr,loc=0) :

• survival function (1-cdf — sometimes more accurate)

binom.ppf(q,n,pr,loc=0) :

• percent point function (inverse of cdf — percentiles)

binom.isf(q,n,pr,loc=0) :

• inverse survival function (inverse of sf)

binom.stats(n,pr,loc=0,moments=’mv’) :

• mean(‘m’,axis=0), variance(‘v’), skew(‘s’), and/or kurtosis(‘k’)

binom.entropy(n,pr,loc=0) :

• entropy of the RV

Alternatively, the object may be called (as a function) to fix :

the shape and location parameters returning a :

“frozen” discrete RV object: :

myrv = binom(n,pr,loc=0) :

• frozen RV object with the same methods but holding the given shape and location fixed.

You can construct an aribtrary discrete rv where P{X=xk} = pk :

by passing to the rv_discrete initialization method (through the values= :

keyword) a tuple of sequences (xk,pk) which describes only those values of :

X (xk) that occur with nonzero probability (pk). :

Examples

>>> import matplotlib.pyplot as plt
>>> numargs = binom.numargs
>>> [ n,pr ] = ['Replace with resonable value',]*numargs

Display frozen pmf:

>>> rv = binom(n,pr)
>>> x = np.arange(0,np.min(rv.dist.b,3)+1)
>>> h = plt.plot(x,rv.pmf(x))

Check accuracy of cdf and ppf:

>>> prb = binom.cdf(x,n,pr)
>>> h = plt.semilogy(np.abs(x-binom.ppf(prb,n,pr))+1e-20)

Random number generation:

>>> R = binom.rvs(n,pr,size=100)

Custom made discrete distribution:

>>> vals = [arange(7),(0.1,0.2,0.3,0.1,0.1,0.1,0.1)]
>>> custm = rv_discrete(name='custm',values=vals)
>>> h = plt.plot(vals[0],custm.pmf(vals[0]))

Binomial distribution

Counts the number of successes in n independent trials when the probability of success each time is pr.

binom.pmf(k,n,p) = choose(n,k)*p**k*(1-p)**(n-k) for k in {0,1,...,n}