[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