Changes from beta versions: * fix in fit function for all response values == 0 * added the option to bound the bmd * changed bmr magic number to variable bmr_scale and updated documentation and vignette to say the default is 1.349 * added … to fitcnst so it could be called in the exact way all the other models are called