Closed tBuLi closed 5 years ago
Mostly it was just because this question made me realize that our API also works for multivariate loglikelihood problems nowadays, but so far we only had tests and examples of single variable problems.
I agree with your second idea, I was actually thinking of extending this example by integrating out the y
variable using sympy to get the marginal distribution, and see what happens then. But more elaborate examples of how sympy can make your life easier should definitely be included.
I'll leave the marginal distribution for another time, it has some unforeseen complications because it is not a clean integral. It contains some mathematical object that cannot be transformed into numerical equivalents.
It seems like this marginal distribution should in fact be a simple Gaussian again, but for some reason sympy does not fully simplify if no further hints are given. http://www.maths.qmul.ac.uk/~ig/MTH5118/Notes11-09.pdf
After a recent stackoverflow question, I decided to include a more elaborate example of likelihood fitting to the docs, which does this feature more justice.