Add nreplicas field in data object; adjust associated code
Avoid repeated use of hstack in nbinom
Reorganised nbinom
Brought plot_fitted_model up-to-date
Modify update model for eta
Added some R code and IPython notebooks