Add stan model that estimates b parameter
Add estimating b model to script