fit_trunc_normal function

Fit Truncated Normal Model on Covariate