coef.bmstdr function

Prints and returns the estimates of the coefficients