Tuning

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