How to implement the Bayesian inference correctly?

My input data is basically input array that represents continuous distribution. And yes st is scipy.stats