Let users give a function that returns a vcov_mat

[nathanielolin]

Should be an option to just pass a function that takes mod and whatever other parameters and returns a vcov matrix. Need to think about whether that should also return dof.