У кого-нибудь есть пример кода использования scipy.stats.distributions?

Question

Я изо всех сил пытаюсь понять, как использовать пакет scipy.distributions, и подумал, может ли кто-нибудь опубликовать для меня какой-нибудь пример кода. Кажется, он делает все, что мне нужно, я просто не могу понять, как его использовать.

Мне нужно сгенерировать два распределения: одно логарифмическое и одно пуассоновское. Я знаю дисперсию и лямбду для каждого.

Ссылки на ресурсы будут работать так же хорошо.

Answer 1

Полагаю, вы имеете в виду распределения в scipy.stats. Чтобы создать распределение, сгенерируйте случайные переменные и вычислите pdf:

Python 2.5.1 (r251: 54863, 4 февраля 2008 г., 21:48:13) [GCC 4.0.1 (Apple Inc., сборка 5465)] на Дарвине Введите «помощь», «авторское право», «кредиты» или «лицензия» для получения дополнительной информации.

>>> from scipy.stats import poisson, lognorm
>>> myShape = 5;myMu=10
>>> ln = lognorm(myShape)
>>> p = poisson(myMu)
>>> ln.rvs((10,)) #generate 10 RVs from ln
array([  2.09164812e+00,   3.29062874e-01,   1.22453941e-03,
         3.80101527e+02,   7.67464002e-02,   2.53530952e+01,
         1.41850880e+03,   8.36347923e+03,   8.69209870e+03,
         1.64317413e-01])
>>> p.rvs((10,)) #generate 10 RVs from p
array([ 8,  9,  7, 12,  6, 13, 11, 11, 10,  8])
>>> ln.pdf(3) #lognorm PDF at x=3
array(0.02596183475208955)

Другие методы (и остальная часть документации scipy.stats) можно найти на новом сайте Документация SciPy .

Answer 2

Вот пример кода: Распределения вероятностей в SciPy