I thought that scipy.stats.binom.log* family of functions were able to return a greater range of values because somehow, they kept all computations in log-space.