У кого-нибудь есть пример кода использования scipy.stats.distributions?
Вопрос
Я изо всех сил пытаюсь понять, как использовать пакет scipy.distributions, и подумал, не мог бы кто-нибудь опубликовать для меня какой-нибудь пример кода.Кажется, он делает все, что мне нужно, я просто не могу понять, как им пользоваться.
Мне нужно сгенерировать два распределения, одно логарифмически нормальное и одно пуассоновское.Я знаю дисперсию и лямбду для каждого из них.
Ссылки на ресурсы работали бы так же хорошо.
Решение
Я предполагаю, что вы имеете в виду дистрибутивы в scipy.stats
.Чтобы создать дистрибутив, сгенерируйте случайные значения и вычислите PDF-файл:
Python 2.5.1 (r251: 54863, 4 февраля 2008, 21:48:13) [GCC 4.0.1 (Apple Inc.сборка 5465)] на darwin Введите "справка", "авторское право", "кредиты" или "лицензия" для получения дополнительной информации.
>>> 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 Сайт.
Другие советы
Вот несколько примеров кода: Распределения вероятностей в SciPy