Hat jemand Beispielcode von haben scipy.stats.distributions verwenden?
Frage
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.
Lösung
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.
Andere Tipps
Hier einige Beispiel-Code: Wahrscheinlichkeitsverteilungen in SciPy