return coefficients and se as well as Bayes factors
put  covariate coefficients at end
fixed bug when using non-null q
update