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

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
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