Fit continuous data with 89 different distributions (scipy.stats) to find which one fit the best.
Inspired and adapted from: https://stackoverflow.com/a/37616966
import scipy.stats as st
from fit_distribution import FitDistributions
data = st.norm.rvs(1, 2, size=5000)
DF = FitDistributions(plot_all = True)
DF.distributions = [st.norm, st.maxwell, st.uniform]
dist_name, params = DF.fit_distribution(data, bins = 100)
norm (0.9603151469277776, 1.997823976707984)