included 'pars' argument in fit()

included 'pars' argument in fit(data,pars) to allow better control of
which covariates from the formula should be included