How to do negative binomial model to fit for time series prevalence data and forecasting?

I need to fit a negative binomial model to a data of prevalence estimates collected across several years and forecast based on this model, How to do it? what software to use? 

19 January 2015 1,871 3 View