Suppose I have some empirical data which I plot and I believe it be a Beta distribution, how can I make sure that I have a Beta distribution, and how can I estimate the alpha and beta parameters of that distribution? I know there's a function in Scipy that can do that scipy.beta.fit() but I'm looking for a bit more than that, can we for example do gradient descent?