[Numpy-discussion] how to fit a given pdf
Renato Fabbri
renato.fabbri@gmail....
Wed Aug 11 10:00:22 CDT 2010
Dear All,
help appreciated, thanks in advance.
how do you fit a pdf you have with a given pdf (say gamma).
with the file attached, you can go like:
a=open("AC-010_ED-1m37F100P0.txt","rb")
aa=a.read()
aaa=aa[1:-1].split(",")
data=[int(i) for i in aaa]
if you do pylab.plot(data); pylab.show()
The data is something like:
___|\___
It is my pdf (probability density function).
how can i find the right parameters to make that fit with a gamma?
if i was looking for a normal pdf, for example, i would just find mean
and std and ask for the pdf.
i've been playing with scipy.stats.distributions.gamma but i have not
reached anything.
we can extend the discussion further, but this is a good starting point.
any idea?
thanks!
best of luck,
rf
--
GNU/Linux User #479299
skype: fabbri.renato
-------------- next part --------------
[3, 1, 3, 2, 2, 2, 1, 1, 1, 1, 5, 4, 9, 16, 26, 51, 55, 83, 134, 228, 222, 202, 203, 180, 142, 114, 75, 55, 48, 41, 28, 26, 17, 13, 11, 14, 15, 14, 16, 9, 3, 3, 4, 3, 7, 4, 3, 3, 5, 7, 2, 2, 1, 8, 10, 14, 22, 11, 10, 12, 6, 3, 2, 1]
More information about the NumPy-Discussion
mailing list