¿Alguien tiene un código de ejemplo para usar scipy.stats.distributions?

Question

Estoy luchando por descubrir cómo usar el paquete scipy.distributions y me pregunté si alguien podría publicar algún código de ejemplo para mí. Parece que hace todo lo que necesito, simplemente no puedo entender cómo usarlo.

Necesito generar dos distribuciones, una logarítmica normal y una poisson. Sé la varianza y lambda para cada uno.

Los enlaces a recursos funcionarían igual de bien.

Answer 1

Supongo que te refieres a las distribuciones en scipy.stats. Para crear una distribución, genere variables aleatorias y calcule el pdf:

Python 2.5.1 (r251: 54863, 4 de febrero de 2008, 21:48:13) [GCC 4.0.1 (Apple Inc. compilación 5465)] en darwin Escriba & Quot; help & Quot ;, & Quot; copyright & Quot ;, & Quot; créditos & Quot; o " licencia " para más información.

>>> 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)

Se pueden encontrar otros métodos (y el resto de la documentación de <=>) en el nuevo documentación de SciPy sitio.

Answer 2

Aquí hay un código de muestra: Distribuciones de probabilidad en SciPy