Quelqu'un at-il un exemple de code d'utilisation de scipy.stats.distributions?
Question
Je ne parviens pas à comprendre comment utiliser le paquet scipy.distributions et je me demandais si quelqu'un pourrait poster un exemple de code pour moi. Il semble faire tout ce dont j'ai besoin, je ne comprends pas comment l'utiliser.
Je dois générer deux distributions, une log-normale et un poisson. Je connais la variance et le lambda pour chacun.
Les liens vers les ressources fonctionneraient aussi bien.
La solution
Je suppose que vous parlez des distributions dans scipy.stats
. Pour créer une distribution, générez des variables aléatoires et calculez le pdf:
Python 2.5.1 (r251: 54863, 4 février 2008, 21:48:13) [GCC 4.0.1 (Apple Inc. build 5465)] sur darwin Tapez & Quot; help & Quot ;, & Quot. Copyright & Quot ;, & Quotents & Quot; ou " licence " pour plus d'informations.
>>> 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)
Vous trouverez d'autres méthodes (et le reste de la <=> documentation) dans le nouveau SciPy documentation .
Autres conseils
Voici un exemple de code: Distributions de probabilités dans SciPy