Spec-Zone.ru › NumPy 1.16

numpy.random.RandomState.dirichlet

метод

RandomState.dirichlet(alpha, size=None)

Генерирует выборки из распределения Дирихле.

Генерирует size выборок размерности k из распределения Дирихле. Дирихле-распределение является многомерным обобщением бета-распределения. Плотность распределения Дирихле — это сопряженный априорный закон для многономиального распределения в байесовском выводе.

Параметры:
alpha : array

Параметр распределения (k-мерный для выборки размерности k).

size : int or tuple of ints, optional

Форма выходных данных. Если заданная форма, например, (m, n, k), тогда генерируется m * n * k выборок. По умолчанию None, в этом случае возвращается единственное значение.

Возвращает:
samples : ndarray,

Сгенерированные выборки, имеющие форму (размер, alpha.ndim).

Исключения:
ValueError

Если любое значение в alpha меньше или равно нулю

Примечания

X \approx \prod_{i=1}^{k}{x^{\alpha_i-1}_i}

Используется следующее свойство для вычисления: для каждой размерности генерируется случайное число y_i из стандартного генератора гамма-распределения с параметром формы alpha_i, затем X = \frac{1}{\sum_{i=1}^k{y_i}} (y_1, \ldots, y_n) имеет распределение Дирихле.

Ссылки

[1] David McKay, “Информация, вывод и алгоритмы обучения”, глава 23, http://www.inference.org.uk/mackay/itila/
[2] Википедия, “Распределение Дирихле”, https://en.wikipedia.org/wiki/Dirichlet_distribution

Примеры

Рассмотрим пример, приведенный в Википедии. Это распределение можно использовать, если необходимо разделить строки (каждая начальной длиной 1,0) на K частей с различными длинами, где каждая часть в среднем имеет заданную среднюю длину, но допускается некоторое изменение относительных размеров частей.

>>> s = np.random.dirichlet((10, 5, 3), 20).transpose()
>>> import matplotlib.pyplot as plt
>>> plt.barh(range(20), s[0])
>>> plt.barh(range(20), s[1], left=s[0], color='g')
>>> plt.barh(range(20), s[2], left=s[0]+s[1], color='r')
>>> plt.title("Lengths of Strings")
../../_images/numpy-random-RandomState-dirichlet-1.png

© 2005–2019 NumPy Developers
Licensed under the 3-clause BSD License.
https://docs.scipy.org/doc/numpy-1.16.1/reference/generated/numpy.random.RandomState.dirichlet.html

Spec-Zone.ru

Настройки Оффлайн Что нового Помощь О нас
Spec-Zone .ru
спецификации, руководства, описания, API