[SciPy-Dev] Subversion scipy.stats irregular problem with source code example

James Phillips zunzun@zunzun....
Fri Oct 1 13:41:48 CDT 2010


On Fri, Oct 1, 2010 at 8:24 AM,  <josef.pktd@gmail.com> wrote:
> I used loc = min(data) - 1e-6  to avoid evaluating the log likelihood
> at the lower bound, in some cases this might be infinite.
> (I picked 1e-6 arbitrarily)
>
>>>> scipy.stats.gamma.logpdf(np.linspace(0,1,3), 300.5609591140931425, loc=0, scale=0.25)
> array([          -Inf, -1205.37511721,  -999.735283  ])

Ah - understood.  Thank you.


> One more: vonmises doesn't define the bounds because it can be used as
> a circular distribution. If you want to include it as a regular
> distribution, then you should set the bound a and b.
> stats.distributions.vonmises.a = -np.pi
> stats.distributions.vonmises.b = np.pi

I shall try this.


     James


More information about the SciPy-Dev mailing list