This doesn't quite make sense to me. What I think the following code is doing is finding me the value of the beta distribution evaluated when p = 0.5 with alpha and beta equal to 10 each. Why is this value greater than 1?
I've implemented the function myself, seperate to this, using the relationship between the beta and gamma functions and I end up with roughly the same value.
>>> from scipy.stats import beta
>>> beta.pdf (0.5, 10, 10)
3.5239410400390625
Any advice would be appreciated.
There is obviously something lacking in my understanding - is this not actually the PDF?
Cheers.
see http://en.wikipedia.org/wiki/Beta_distribution
You are mistaking the PDF (Probability density function) and CDF (Cumulative distribution function)
CDF will be less than one at any point