dist_names = ['gamma', 'beta', 'rayleigh', 'norm', 'pareto']
for dist_name in dist_names:
    dist = getattr(scipy.stats, dist_name)
param = dist.fit(y)
Как найти распределение вероятностей и параметры для реальных данных? (Питон 3)