Tuning

find_hessian(point[, vars, model, negate_output])	Return Hessian of logp at the point passed.
find_MAP([start, vars, method, return_raw, ...])	Find the local maximum a posteriori point given a model.