Calc-Based Stats, Prove a Conditional Distribution is Gamma Distributed

Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
1 reply · 2K views
aleph_0
Messages
6
Reaction score
0

Homework Statement



An image of the assigned problem is here: http://imgur.com/aYkaM

Homework Equations



The formula for being exponential, gamma, and probably Bayes's Law. They'd take a while to type out, and presumably anyone who can help me with this already knows the formulas or will understand them by looking at my solution thus far.

The Attempt at a Solution



[itex]P(X_{1}=x_{1}, ..., X_{n}=x_{n}|W=w)P(W=w) \quad = \\\\ P(W=w|X_{1}=x_{1}, ..., X_{n}=x_{n})P(X_{1}=x_{1}, ..., X_{n}=x_{n}) \quad \Longrightarrow \\\\ P(W=w|X_{1}=x_{1}, ..., X_{n}=x_{n}) \quad = \quad \frac{P(X_{1}=x_{1}, ..., X_{n}=x_{n}|W=w)P(W=w)}{P(X_{1}=x_{1}, ..., X_{n}=x_{n})} \quad = \\\\ \frac{w^{n}e^{-w(x_{1}+...+x_{n})}\frac{\beta^{t}}{\Gamma (t)}w^{t-1}e^{-\beta w}}{P(X_{1}=x_{1}, ..., X_{n}=x_{n})} \quad = \quad \frac{\beta^{t}}{\Gamma (t) P(X_{1}=x_{1}, ..., X_{n}=x_{n})}{w^{t+n-1}e^{-w(\beta + \sum_{i=1}^{n}x_{i})}}[/itex]

Now at this point I feel like I've gone most of the way, but the only thing is the "coefficient". For this to be a true gamma distribution with the said parameters, that front factor needs to be [itex]\frac{(\beta+\sum x_{i})^{t+n}}{\Gamma(t+n)}[/itex] and I just can't see how to make that happen.
 
Physics news on Phys.org
aleph_0 said:

Homework Statement



An image of the assigned problem is here: http://imgur.com/aYkaM

Homework Equations



The formula for being exponential, gamma, and probably Bayes's Law. They'd take a while to type out, and presumably anyone who can help me with this already knows the formulas or will understand them by looking at my solution thus far.

The Attempt at a Solution



[itex]P(X_{1}=x_{1}, ..., X_{n}=x_{n}|W=w)P(W=w) \quad = \\\\ P(W=w|X_{1}=x_{1}, ..., X_{n}=x_{n})P(X_{1}=x_{1}, ..., X_{n}=x_{n}) \quad \Longrightarrow \\\\ P(W=w|X_{1}=x_{1}, ..., X_{n}=x_{n}) \quad = \quad \frac{P(X_{1}=x_{1}, ..., X_{n}=x_{n}|W=w)P(W=w)}{P(X_{1}=x_{1}, ..., X_{n}=x_{n})} \quad = \\\\ \frac{w^{n}e^{-w(x_{1}+...+x_{n})}\frac{\beta^{t}}{\Gamma (t)}w^{t-1}e^{-\beta w}}{P(X_{1}=x_{1}, ..., X_{n}=x_{n})} \quad = \quad \frac{\beta^{t}}{\Gamma (t) P(X_{1}=x_{1}, ..., X_{n}=x_{n})}{w^{t+n-1}e^{-w(\beta + \sum_{i=1}^{n}x_{i})}}[/itex]

Now at this point I feel like I've gone most of the way, but the only thing is the "coefficient". For this to be a true gamma distribution with the said parameters, that front factor needs to be [itex]\frac{(\beta+\sum x_{i})^{t+n}}{\Gamma(t+n)}[/itex] and I just can't see how to make that happen.

[tex]P(X_1=x_1,\ldots,X_n=x_n) = \int P(X_1=x_1,\ldots,X_n=x_n|w) f_W(w) \, dw.[/tex] Just do the integral.

RGV