[SciPy-Dev] Implementing fit method for discrete distributions