estimate_dm function

Fit a DDM to Observed Data