Dear Goran, I need to find an EM algorithm for finding the MLE estimates for a negative binomial distribution. I saw your web post on the R help network and wonder if you might be willing to share your code or algorithm. It would be greatly appreciated. Thank you, Dave LeBlond Sr. Statistician, Non-Clinical Statistics email: david.leblond@abbott.com voice mail: 847-935-6031 Dept R436, Bldg AP9A-1 Abbott Laboratories 100 Abbott Park Road Abbott Park, IL 60064-3500 [[alternative HTML version deleted]]