Probability distributions like the negative binomial and hypergeometric show up constantly in statistics, biology, and data science. Both describe how often certain outcomes occur in specific counting situations: the negative binomial tracks how many trials are needed to achieve a set number of successes, while the hypergeometric describes drawing objects from a finite collection without replacement. For these distributions, the "mode" is the most likely outcome, meaning the value where the probability is highest. Understanding exactly how probabilities behave near that peak is useful for approximations, computations, and theoretical analysis.
The paper works out precise mathematical formulas describing how the probability values change as you move away from the mode, step by step across the discrete outcomes. These formulas are expressed as layered expansions, where each successive term captures finer and finer detail about the shape of the distribution near its peak. The coefficients in these formulas are built from well-known mathematical objects called Bernoulli polynomials, and the formulas carefully account for the fact that the mode does not always sit exactly at the mean, tracking the small offset between those two values precisely.
A bonus result concerns identifying the mode itself. For the negative binomial, the formulas exactly reproduce a classical rule for pinpointing which integer value is the most probable. For the hypergeometric distribution, the formulas reveal a small discrepancy in the analogous rule, but the paper shows this discrepancy is always too tiny to shift the answer to a different integer, confirming the classical mode formula is still correct in practice. Together with a companion paper covering the ordinary binomial distribution, this work completes a unified treatment of the three most fundamental discrete probability laws.