numpy.random.Generator.beta

метод

random.Generator.beta(a, b, размер=None)

Извлечь выборки из бета-распределения.

Бета-распределение является частным случаем распределения Дирихле и связано с гамма-распределением. Оно имеет функцию плотности вероятности

\[f(x; a,b) = \frac{1}{B(\alpha, \beta)} x^{\alpha - 1} (1 - x)^{\beta - 1},\]

где нормализация, B, является бета-функцией,

\[B(\alpha, \beta) = \int_0^1 t^{\alpha - 1} (1 - t)^{\beta - 1} dt.\]

Часто встречается в байесовском выводе и порядковой статистике.

Параметры:

afloat или array_like из float

Альфа, положительная (>0).

bfloat или array_like из float

Бета, положительная (>0).

размерint или кортеж ints, опционально

Форма вывода. Если заданная форма, например, (m, n, k), затем m * n * k образцы извлекаются. Если size равен None (по умолчанию), возвращается единственное значение, если a и b оба являются скалярами. В противном случае, np.broadcast(a, b).size выбираются образцы.

Возвращает:

выходndarray или скаляр

Выборки из параметризованного бета-распределения.

Ссылки

[1]

Википедия, «Бета-распределение», https://en.wikipedia.org/wiki/Beta_distribution

Примеры

Попробуйте в вашем браузере!

Бета-распределение имеет среднее a/(a+b). Если a == b и оба > 1, распределение симметрично со средним 0.5.

>>> rng = np.random.default_rng()
>>> a, b, size = 2.0, 2.0, 10000
>>> sample = rng.beta(a=a, b=b, size=size)
>>> np.mean(sample)
0.5047328775385895  # may vary

В противном случае распределение смещено влево или вправо в зависимости от того, a или b больше. Распределение зеркально симметрично. См., например:

>>> a, b, size = 2, 7, 10000
>>> sample_left = rng.beta(a=a, b=b, size=size)
>>> sample_right = rng.beta(a=b, b=a, size=size)
>>> m_left, m_right = np.mean(sample_left), np.mean(sample_right)
>>> print(m_left, m_right)
0.2238596793678923 0.7774613834041182  # may vary
>>> print(m_left - a/(a+b))
0.001637457145670096  # may vary
>>> print(m_right - b/(a+b))
-0.0003163943736596009  # may vary

Отобразить гистограмму двух выборок:

>>> import matplotlib.pyplot as plt
>>> plt.hist([sample_left, sample_right],
...          50, density=True, histtype='bar')
>>> plt.show()

НазадОткрыть во вкладке