fit_mreg function

Fit a model for the mediator given the treatment and covariates.