问题 In scipy there is no support for fitting discrete distributions using data. I know there are a lot of subject about this. For example if i have an array like below: x = [2,3,4,5,6,7,0,1,1,0,1,8,10,9,1,1,1,0,0] I couldn' t apply for this array; from scipy.stats import nbinom param = nbinom.fit(x) But i would like to ask you up to date, is there any way to fit for these three discrete distributions and then choose the best fit for the discrete dataset? 回答1: You can use Method of Moments to fit