refit.model_mean function

Refit a MEAN model