Maximum likelihood estimation of a Bayesian model

I dont think the symbolic expression is available (at least not from theano). If you know sympy that would be a good direction to try.