Skewed Normal Distribution

Hi

I have a code snippet where I use curve_fit to optimize parameters of a function to fit some data, it’s a 5 parameter logistic function. With the optimized parameters I continue with some calculations where finding the maximum probability of the probability function is one of the operations. Now I’m curious to see if there’s some ‘standard’ library where I can use a skewed normal distribution instead of the 5PL model. I’m kind of new to Python and would appreciate help regarding libraries. Is there any such library where the mode (the highest density of the distribution) of a skewed normal distribution is passed as one of the parameters?

Our skew normal distribution is not parameterized by the mode but by a shape parameter. You could solve for the value of the shape parameter that results in the desired mode. See here for an example.

Thanks a lot mdhaber! I’ll try that approach instead.