Understanding and Mitigating Errors with MASS::glm.nb Package in R for Negative Binomial Regression