point_estimates function

Compute point estimates