from scipy.stats import rv_continuous, gamma
x = gamma.rvs(2.5, loc=0, scale=4, size=1000)
rv_continuous.fit(gamma, x, floc=0, fscale=4)
Использование python scipy для подбора гамма-распределения к данным