fit.shape function

Shape parameter estimates