Finding the MLE from PyMC3

mle.py

mle = {}

for var_name in ['a','h','k','sigma']:
    var_mass, var_bins = np.histogram(trace[var_name], bins=np.arange(-50000,50000, 0.1), density=True)
    var_averaged_bins = 0.5*(var_bins[1:] + var_bins[:-1])
    bin_index = np.argmax(var_mass)
    mle[var_name] = var_averaged_bins[bin_index]

print(mle);