Tuning

find_hessian(point[, vars, model, negate_output])	Return Hessian of logp at the point passed.
find_MAP([method, vars, use_grad, use_hess, ...])	Find the local maximum a posteriori point of a model with scipy.optimize.