Hat jemand Beispielcode von haben scipy.stats.distributions verwenden?

Question

Ich kämpfe, um herauszufinden, wie die scipy.distributions Paket verwenden und fragte sich, ob jemand einige Beispiel-Code für mich schreiben könnte. Es scheint, alles zu tun, ich brauche, ich kann einfach nicht herausfinden, wie es zu benutzen.

Ich brauche zwei Verteilungen zu erzeugen, ein Log-normal und ein poisson. Ich weiß, dass die Varianz und Lambda für jeden.

Links zu Ressourcen würden funktionieren genauso gut.

Answer 1

Ich nehme an, Sie die Verteilungen in scipy.stats bedeuten. Um eine Verteilung zu erstellen, erzeugen Zufallszahlen und berechnen Sie die pdf:

Python 2.5.1 (R251: 54863, 4. Februar 2008, 21.48.13) [GCC 4.0.1 (Apple Inc. bauen 5465)] auf Darwin Geben Sie "help", "Copyright", "Credits" oder "Lizenz" für weitere Informationen.

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

Andere Methoden (und der Rest der scipy.stats Dokumentation) können an den neuen SciPy Dokumentation Website.

Answer 2

Hier einige Beispiel-Code: Wahrscheinlichkeitsverteilungen in SciPy