[SciPy-User] generic expectation operator
nicky van foreest
vanforeest@gmail....
Thu Aug 23 12:11:58 CDT 2012
Hi,
I noticed that rv_frozen does not have an expect attribute. Is there a
design reason for this? It feels somewhat unnatural to me that frozen
distributions cannot be used directly to compute expectations.
bye
Nicky
On 23 August 2012 18:36, nicky van foreest <vanforeest@gmail.com> wrote:
> Hi Josef,
>
> Thanks for your answers.
>
>>>> return scipy.integrate.quad(lambda x: x*g(x), X.dist.a, X.dist.b)
>
>>>> X = stats.norm(3,sqrt(3))
>>>> print E(sqrt, X)
>
>> I don't know why this works, sqrt of a normal distributed random
>> variable has lots of nans (or complex numbers)
>
> You are right, it shouldn't work. The point is that my code above is
> this: lambda x: x*g(x), while it should have been lambda x: f(x)*g(x).
>
> I'll give expect a try.
>
> Nicky
More information about the SciPy-User
mailing list