Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

I have that Y ~ Bernouill(Pi) and Pi~Beta(a,b), find the prior predictive distribution.

I know the pmf of Y is Pi^Y(1-Pi)^1-Y And the pdf of Pi is Pi^(Y+a-1)(1-Pi)^(b-Y)

I then multiplied these together and put the product in an intgral over 0,1 but now I am stuck...


share|cite|improve this question

Basic fact: For every positive $u$ and $v$, $\displaystyle\int\limits_0^1x^{u-1}(1-x)^{v-1}\,\mathrm dx=\frac{\Gamma(u)\Gamma(v)}{\Gamma(u+v)}$.

Thus the PDF of $\Pi$ is not what you wrote but a multiple of this, and, for every $y$ in $\{0,1\}$, $$ L_{a,b}^y=\int_0^1f_{Y\mid\Pi}(y\mid x)\cdot f_\Pi(x)\,\mathrm dx=\int_0^1x^y(1-x)^{1-y}\cdot\frac{\Gamma(a+b)}{\Gamma(a)\Gamma(b)}x^{a-1}(1-x)^{b-1}\,\mathrm dx, $$ that is, $$ L_{a,b}^y=\frac{\Gamma(a+b)}{\Gamma(a)\Gamma(b)}\cdot\frac{\Gamma(a+y)\Gamma(b+1-y)}{\Gamma(a+b+1)}. $$ In other words, $L_{a,b}^0=\dfrac{b}{a+b}$ and $L_{a,b}^1=\dfrac{a}{a+b}$.

share|cite|improve this answer

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.