Definite Integral of loglikelihood function multipled by Gaussian

asked 2021-01-27 12:51:41 +0100

user2749584820 gravatar image

I am trying to calculate

definite_integral((-ln(2*pi)-ln(sigma)-1/2*((x-mu)/sigma)^2)*f, x, -infinity, infinity)

where

f=1/(sqrt(2*pi)*sigmaprime)*exp(-1/2*((x-muprime)/sigmaprime)^2)

using sagemath.

I have done from sage.symbolic.integration.integral import definite_integral

And I have also done

assume(sigma>0)
assume(sigmaprime>0)

however I obtain this error

ValueError: Computation failed since Maxima requested additional constraints; using the 'assume' command before evaluation *may* help (example of legal syntax is 'assume(muprime>0)', see `assume?` for more details)

Is muprime positive, negative or zero?

However muprime can be positive, negative or zero. So I don't know how to proceed.

edit retag flag offensive close merge delete