[Numpy-discussion] Fitting a curve on a log-normal distributed data

Gökhan Sever gokhansever at gmail.com
Fri Nov 20 13:27:48 EST 2009

On Thu, Nov 19, 2009 at 9:12 PM, Ian Mallett <geometrian at gmail.com> wrote:

> Hello,
> My analysis shows that the exponential regression gives the best result
> (r^2=87%)--power regression gives worse results (r^2=77%).  Untransformed
> data gives r^2=76%.
> I don't think you want lognorm.  If I'm not mistaken, that fits the data to
> a log(normal distribution random variable).

Lognormal fitting is the motivation behind my study since aerosol in the
atmosphere typically log-normally size distributed. See for an example case

Of course this is just a simplification. There are other approaches to
represent the size-distribution besides the lognormal. So my intention is
not actually find the best fit but represent the actuality as much as I can.

> So, take the logarithm (to any base) of all the "conc" values.  Then do a
> linear regression on those values versus "sizes".
> Try (semi-pseudocode):
> slope, intercept, p, error = scipy.stats.linregress(sizes,log(conc))

linregress also returns the r_value which I am not sure if the documentation
from the web-based editor checked-in completely to the scipy trunk yet.

> Ian
> _______________________________________________
> NumPy-Discussion mailing list
> NumPy-Discussion at scipy.org
> http://mail.scipy.org/mailman/listinfo/numpy-discussion

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/numpy-discussion/attachments/20091120/fe4a781d/attachment.html>

More information about the NumPy-Discussion mailing list