add function in estimator utils for extracting map value

[Dingel321]

either evaluate log prob on a dense grid and search for max, or optimize in logprob space