flopscope.

flopscope.numpy.random.Generator.pareto

fnp.random.Generator.pareto(self, a, size=None)

Draw samples from a Pareto II (AKA Lomax) distribution with specified shape.

Adapted from NumPy docs np.random.Generator.pareto

Arearandom
Typecounted
NumPy Refnp.random.Generator.pareto
Cost
numel(output)\text{numel}(\text{output})
Flopscope Context

Pareto distribution; cost = numel(output).

Parameters

a:float or array_like of floats

Shape of the distribution. Must be positive.

size:int or tuple of ints, optional

Output shape. If the given shape is, e.g., (m, n, k), then m * n * k samples are drawn. If size is None (default), a single value is returned if a is a scalar. Otherwise, flops.array(a).size samples are drawn.

Returns

out:ndarray or scalar

Drawn samples from the Pareto II distribution.

See also

Notes

The probability density for the Pareto II distribution is

p(x)=ax+1a+1,x≥0p(x) = \frac{a}{{x+1}^{a+1}} , x \ge 0

where a>0a > 0 is the shape.

The Pareto II distribution is a shifted and scaled version of the Pareto I distribution, which can be found in scipy.stats.pareto.

References

footnote
1

Francis Hunt and Paul Johnson, On the Pareto Distribution of
Sourceforge projects.
footnote
2

Pareto, V. (1896). Course of Political Economy. Lausanne.
footnote
3

Reiss, R.D., Thomas, M.(2001), Statistical Analysis of Extreme
Values, Birkhauser Verlag, Basel, pp 23-30.
footnote
4

Wikipedia, "Pareto distribution",
https://en.wikipedia.org/wiki/Pareto_distribution

Examples

Draw samples from the distribution:

>>> a = 3.
>>> rng = flops.random.default_rng()
>>> s = rng.pareto(a, 10000)

Display the histogram of the samples, along with the probability density function:

>>> import matplotlib.pyplot as plt
>>> x = flops.linspace(0, 3, 50)
>>> pdf = a / (x+1)**(a+1)
>>> plt.hist(s, bins=x, density=True, label='histogram')
>>> plt.plot(x, pdf, linewidth=2, color='r', label='pdf')
>>> plt.xlim(x.min(), x.max())
>>> plt.legend()
>>> plt.show()