[SciPy-User] How to fit a curve/function?

[Johannes Radinger]

Hello,

I've got following function describing any kind of animal dispersal kernel:

def pdf(x,s1,s2):
    return (p/(math.sqrt(2*math.pi*s1**2))*numpy.exp(-((x-0)**(2)/(2*s1**(2)))))+((1-p)/(s2*math.sqrt(2*math.pi))*numpy.exp(-((x-0)**(2)/(2*s2**(2)))))

On the other hand I've got data from literature with which I want to fit the function so that I get s1, s2 and x. 
Ususally the data in the literature are as follows:

Example 1: 50% of the animals are between -270m and +270m and 90%  are between -500m and + 500m

Example 2: 84% is between - 5000 m and +5000m, and 73% are between -1000m and +1000m

So far as I understand an integration of the function is needed to solve for s1 and s2 as all the literature data give percentage (area under the curve) Can that be used to fit the curve or can that create ranges for s1 and s2.

/Johannes

-- 
NEU: FreePhone - kostenlos mobil telefonieren!			
Jetzt informieren: http://www.gmx.net/de/go/freephone