Re: [Numpy-discussion] pareto docstring

2010-05-11 Thread josef . pktd
Assuming no typos Relationship between Pareto, Pareto(II)/Lomax and Generalized Pareto,GPD >>> import sympy as sy >>> x = sy.Symbol('x') >>> k = sy.Symbol('k') >>> c = sy.Symbol('c') >>> a = sy.Symbol('a') >>> m = sy.Symbol('m') >>> mgpd = sy.Symbol('mgpd') >>> gpd0 = (1 - c*x/k)**(1/c - 1)/k #

Re: [Numpy-discussion] pareto docstring

2010-05-11 Thread josef . pktd
On Tue, May 11, 2010 at 8:11 AM, Kevin Jacobs wrote: > On Tue, May 11, 2010 at 4:14 AM, Pauli Virtanen wrote: >> >> A third option would be just to silently fix the bug. In any case the >> change should be mentioned noticeably in the release notes. >> > > I see this as two bugs: the Lomax distri

Re: [Numpy-discussion] pareto docstring

2010-05-11 Thread Kevin Jacobs
On Tue, May 11, 2010 at 4:14 AM, Pauli Virtanen wrote: > A third option would be just to silently fix the bug. In any case the > change should be mentioned noticeably in the release notes. > > I see this as two bugs: the Lomax distribution was named incorrectly and the Parato distribution was inc

Re: [Numpy-discussion] pareto docstring

2010-05-11 Thread Pauli Virtanen
Tue, 11 May 2010 00:23:52 -0700, T J wrote: [clip] > It seems reasonable that we might have to follow the deprecation route, > but I'd be happier with a "faster" fix. > > 1.5 > - Provide numpy.random.lomax. Make numpy.random.pareto raise a > DeprecationWarning and then call lomax. > > 2.0 (

Re: [Numpy-discussion] pareto docstring

2010-05-11 Thread David Goldsmith
On Tue, May 11, 2010 at 12:23 AM, T J wrote: > On Mon, May 10, 2010 at 8:37 PM, wrote: > > > > I went googling and found a new interpretation > > > > numpy.random.pareto is actually the Lomax distribution also known as > Pareto 2, > > Pareto (II) or Pareto Second Kind distribution > > > > Great

Re: [Numpy-discussion] pareto docstring

2010-05-11 Thread T J
On Mon, May 10, 2010 at 8:37 PM, wrote: > > I went googling and found a new interpretation > > numpy.random.pareto is actually the Lomax distribution also known as Pareto 2, > Pareto (II) or Pareto Second Kind distribution > Great! > > So, from this it looks like numpy.random does not have a Pa

Re: [Numpy-discussion] pareto docstring

2010-05-10 Thread josef . pktd
On Mon, May 10, 2010 at 2:14 PM, T J wrote: > On Sun, May 9, 2010 at 4:49 AM,   wrote: >> >> I think this is the same point, I was trying to make last year. >> >> Instead of renormalizing, my conclusion was the following, >> (copied from the mailinglist August last year) >> >> """ >> my conclusion

Re: [Numpy-discussion] pareto docstring

2010-05-10 Thread David Goldsmith
On Mon, May 10, 2010 at 11:14 AM, T J wrote: > On Sun, May 9, 2010 at 4:49 AM, wrote: > > > > I think this is the same point, I was trying to make last year. > > > > Instead of renormalizing, my conclusion was the following, > > (copied from the mailinglist August last year) > > > > """ > > my

Re: [Numpy-discussion] pareto docstring

2010-05-10 Thread T J
On Sun, May 9, 2010 at 4:49 AM, wrote: > > I think this is the same point, I was trying to make last year. > > Instead of renormalizing, my conclusion was the following, > (copied from the mailinglist August last year) > > """ > my conclusion: > - > What numpy.random.pareto ac

Re: [Numpy-discussion] pareto docstring

2010-05-09 Thread David Goldsmith
On Sun, May 9, 2010 at 4:49 AM, wrote: > On Sun, May 9, 2010 at 1:01 AM, T J wrote: > > The docstring for np.pareto says: > > > >This is a simplified version of the Generalized Pareto distribution > >(available in SciPy), with the scale set to one and the location set > to > >zero. M

Re: [Numpy-discussion] pareto docstring

2010-05-09 Thread josef . pktd
On Sun, May 9, 2010 at 1:01 AM, T J wrote: > The docstring for np.pareto says: > >    This is a simplified version of the Generalized Pareto distribution >    (available in SciPy), with the scale set to one and the location set to >    zero. Most authors default the location to one. > > and also:

[Numpy-discussion] pareto docstring

2010-05-08 Thread T J
The docstring for np.pareto says: This is a simplified version of the Generalized Pareto distribution (available in SciPy), with the scale set to one and the location set to zero. Most authors default the location to one. and also: The probability density for the Pareto distribut