[SciPy-User] how to fit a given pdf

[Renato Fabbri]

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]