rv_continuous#
- class scipy.stats.rv_continuous(momtype=1, a=None, b=None, xtol=1e-14, badvalue=None, name=None, longname=None, shapes=None, seed=None)[source]#
A generic continuous random variable class meant for subclassing.
rv_continuous
is a base class to construct specific distribution classes and instances for continuous random variables. It cannot be used directly as a distribution.- Parameters:
- momtypeint, optional
The type of generic moment calculation to use: 0 for pdf, 1 (default) for ppf.
- afloat, optional
Lower bound of the support of the distribution, default is minus infinity.
- bfloat, optional
Upper bound of the support of the distribution, default is plus infinity.
- xtolfloat, optional
The tolerance for fixed point calculation for generic ppf.
- badvaluefloat, optional
The value in a result arrays that indicates a value that for which some argument restriction is violated, default is np.nan.
- namestr, optional
The name of the instance. This string is used to construct the default example for distributions.
- longnamestr, optional
This string is used as part of the first line of the docstring returned when a subclass has no docstring of its own. Note: longname exists for backwards compatibility, do not use for new subclasses.
- shapesstr, optional
The shape of the distribution. For example
"m, n"
for a distribution that takes two integers as the two shape arguments for all its methods. If not provided, shape parameters will be inferred from the signature of the private methods,_pdf
and_cdf
of the instance.- seed{None, int,
numpy.random.Generator
,numpy.random.RandomState
}, optional If seed is None (or np.random), the
numpy.random.RandomState
singleton is used. If seed is an int, a newRandomState
instance is used, seeded with seed. If seed is already aGenerator
orRandomState
instance then that instance is used.
- Attributes:
random_state
Get or set the generator object for generating random variates.
Methods
rvs
(*args, **kwds)Random variates of given type.
pdf
(x, *args, **kwds)Probability density function at x of the given RV.
logpdf
(x, *args, **kwds)Log of the probability density function at x of the given RV.
cdf
(x, *args, **kwds)Cumulative distribution function of the given RV.
logcdf
(x, *args, **kwds)Log of the cumulative distribution function at x of the given RV.
sf
(x, *args, **kwds)Survival function (1 -
cdf
) at x of the given RV.logsf
(x, *args, **kwds)Log of the survival function of the given RV.
ppf
(q, *args, **kwds)Percent point function (inverse of
cdf
) at q of the given RV.isf
(q, *args, **kwds)Inverse survival function (inverse of
sf
) at q of the given RV.moment
(order, *args, **kwds)non-central moment of distribution of specified order.
stats
(*args, **kwds)Some statistics of the given RV.
entropy
(*args, **kwds)Differential entropy of the RV.
expect
([func, args, loc, scale, lb, ub, ...])Calculate expected value of a function with respect to the distribution by numerical integration.
median
(*args, **kwds)Median of the distribution.
mean
(*args, **kwds)Mean of the distribution.
std
(*args, **kwds)Standard deviation of the distribution.
var
(*args, **kwds)Variance of the distribution.
interval
(confidence, *args, **kwds)Confidence interval with equal areas around the median.
__call__
(*args, **kwds)Freeze the distribution for the given arguments.
fit
(data, *args, **kwds)Return estimates of shape (if applicable), location, and scale parameters from data.
fit_loc_scale
(data, *args)Estimate loc and scale parameters from data using 1st and 2nd moments.
nnlf
(theta, x)Negative loglikelihood function.
support
(*args, **kwargs)Support of the distribution.
Notes
Public methods of an instance of a distribution class (e.g.,
pdf
,cdf
) check their arguments and pass valid arguments to private, computational methods (_pdf
,_cdf
). Forpdf(x)
,x
is valid if it is within the support of the distribution. Whether a shape parameter is valid is decided by an_argcheck
method (which defaults to checking that its arguments are strictly positive.)Subclassing
New random variables can be defined by subclassing the
rv_continuous
class and re-defining at least the_pdf
or the_cdf
method (normalized to location 0 and scale 1).If positive argument checking is not correct for your RV then you will also need to re-define the
_argcheck
method.For most of the scipy.stats distributions, the support interval doesn’t depend on the shape parameters.
x
being in the support interval is equivalent toself.a <= x <= self.b
. If either of the endpoints of the support do depend on the shape parameters, then i) the distribution must implement the_get_support
method; and ii) those dependent endpoints must be omitted from the distribution’s call to therv_continuous
initializer.Correct, but potentially slow defaults exist for the remaining methods but for speed and/or accuracy you can over-ride:
_logpdf, _cdf, _logcdf, _ppf, _rvs, _isf, _sf, _logsf
The default method
_rvs
relies on the inverse of the cdf,_ppf
, applied to a uniform random variate. In order to generate random variates efficiently, either the default_ppf
needs to be overwritten (e.g. if the inverse cdf can expressed in an explicit form) or a sampling method needs to be implemented in a custom_rvs
method.If possible, you should override
_isf
,_sf
or_logsf
. The main reason would be to improve numerical accuracy: for example, the survival function_sf
is computed as1 - _cdf
which can result in loss of precision if_cdf(x)
is close to one.Methods that can be overwritten by subclasses
_rvs _pdf _cdf _sf _ppf _isf _stats _munp _entropy _argcheck _get_support
There are additional (internal and private) generic methods that can be useful for cross-checking and for debugging, but might work in all cases when directly called.
A note on
shapes
: subclasses need not specify them explicitly. In this case, shapes will be automatically deduced from the signatures of the overridden methods (pdf
,cdf
etc). If, for some reason, you prefer to avoid relying on introspection, you can specifyshapes
explicitly as an argument to the instance constructor.Frozen Distributions
Normally, you must provide shape parameters (and, optionally, location and scale parameters to each call of a method of a distribution.
Alternatively, the object may be called (as a function) to fix the shape, location, and scale parameters returning a “frozen” continuous RV object:
- rv = generic(<shape(s)>, loc=0, scale=1)
rv_frozen object with the same methods but holding the given shape, location, and scale fixed
Statistics
Statistics are computed using numerical integration by default. For speed you can redefine this using
_stats
:take shape parameters and return mu, mu2, g1, g2
If you can’t compute one of these, return it as None
Can also be defined with a keyword argument
moments
, which is a string composed of “m”, “v”, “s”, and/or “k”. Only the components appearing in string should be computed and returned in the order “m”, “v”, “s”, or “k” with missing values returned as None.
Alternatively, you can override
_munp
, which takesn
and shape parameters and returns the n-th non-central moment of the distribution.Deepcopying / Pickling
If a distribution or frozen distribution is deepcopied (pickled/unpickled, etc.), any underlying random number generator is deepcopied with it. An implication is that if a distribution relies on the singleton RandomState before copying, it will rely on a copy of that random state after copying, and
np.random.seed
will no longer control the state.Examples
To create a new Gaussian distribution, we would do the following:
>>> from scipy.stats import rv_continuous >>> class gaussian_gen(rv_continuous): ... "Gaussian distribution" ... def _pdf(self, x): ... return np.exp(-x**2 / 2.) / np.sqrt(2.0 * np.pi) >>> gaussian = gaussian_gen(name='gaussian')
scipy.stats
distributions are instances, so here we subclassrv_continuous
and create an instance. With this, we now have a fully functional distribution with all relevant methods automagically generated by the framework.Note that above we defined a standard normal distribution, with zero mean and unit variance. Shifting and scaling of the distribution can be done by using
loc
andscale
parameters:gaussian.pdf(x, loc, scale)
essentially computesy = (x - loc) / scale
andgaussian._pdf(y) / scale
.