Gamma distribution from sample mean of Exponential distribution

Click For Summary
SUMMARY

The discussion focuses on demonstrating that the sample mean \(\overline{X}\) from an exponential distribution with mean \(\theta\) follows a Gamma distribution, specifically \(\overline{X} \sim \text{Gamma}(n, \frac{n}{\theta})\). Participants utilize the moment-generating functions (MGFs) of the exponential and Gamma distributions, noting that the random variable \(S_n\) (the sum of the sample) has an n-Erlang distribution. Key equations include the MGF of the exponential distribution \(\frac{\lambda}{\lambda - t}\) and the MGF of the Gamma distribution \((\frac{\beta}{\beta - t})^{\alpha}\).

PREREQUISITES
  • Understanding of exponential distribution and its properties
  • Familiarity with Gamma distribution and its applications
  • Knowledge of moment-generating functions (MGFs)
  • Basic calculus, particularly differentiation and integration techniques
NEXT STEPS
  • Study the derivation of the Gamma distribution from the Erlang distribution
  • Learn about moment-generating functions and their applications in probability theory
  • Explore change-of-variables techniques in probability density functions
  • Investigate the properties and applications of the Gamma function
USEFUL FOR

Statisticians, data scientists, and students studying probability theory, particularly those interested in the properties of the Gamma distribution and its relationship with the exponential distribution.

tmbrwlf730
Messages
41
Reaction score
0

Homework Statement


Let X1, X2,...,Xn be a random sample from the exponential distribution with mean θ and \overline{X} = \sum^{n}_{i = 1}X_i

Show that \overline{X} ~ Gamma(n, \frac{n}{θ})

Homework Equations



θ = \frac{1}{λ}

MGF Exponential Distribution = \frac{λ}{λ - t}

MGF Gamma Distribution = (\frac{β}{β - t})α



The Attempt at a Solution


I've tried using the generating function of the exponential distribution but I end up with

\frac{(\frac{λ}{λ-t})^{n}}{n}

I don't know what to do with the n in the denominator to get λ = \frac{n}{θ}
 
Physics news on Phys.org
tmbrwlf730 said:

Homework Statement


Let X1, X2,...,Xn be a random sample from the exponential distribution with mean θ and \overline{X} = \sum^{n}_{i = 1}X_i

Show that \overline{X} ~ Gamma(n, \frac{n}{θ})

Homework Equations



θ = \frac{1}{λ}

MGF Exponential Distribution = \frac{λ}{λ - t}

MGF Gamma Distribution = (\frac{β}{β - t})α



The Attempt at a Solution


I've tried using the generating function of the exponential distribution but I end up with

\frac{(\frac{λ}{λ-t})^{n}}{n}

I don't know what to do with the n in the denominator to get λ = \frac{n}{θ}

We have ##\bar{X} = S_n/n,## where ##S_n = \sum_{i=1}^n X_i##. The random variable ##S_n## has an n-Erlang distribution with mean ##n \theta##. (Erlang is a special case of Gamma; its density is widely available in textbooks and on line.) To get the distribution pdf ##f(x)## of ##\bar{X}##, use
f(x) =\frac{d}{dx} P (\bar{X} \leq x), \text{ and }<br /> P(\bar{X} \leq x) = P(S_n \leq nx).
 
Ray Vickson said:
We have ##\bar{X} = S_n/n,## where ##S_n = \sum_{i=1}^n X_i##. The random variable ##S_n## has an n-Erlang distribution with mean ##n \theta##. (Erlang is a special case of Gamma; its density is widely available in textbooks and on line.) To get the distribution pdf ##f(x)## of ##\bar{X}##, use
f(x) =\frac{d}{dx} P (\bar{X} \leq x), \text{ and }<br /> P(\bar{X} \leq x) = P(S_n \leq nx).

I'm sorry but I'm still not getting it.

I understand about the erlang distribution and after plugging in what should be the parameters for it from

(\frac{λ}{λ - t})n

I get ∫^{nx}_{0}\frac{s^{n-1} e^{-\frac{s}{θ}}}{θ^n \Gamma (n)}

where \Gamma(n) = ∫^{∞}_{0}sn-1 e-s ds

I'm guessing that the sn-1 cancel out but I feel like I'm missing a change of variables by a Jacobian.

Also in your response you mentioned we took \frac{d}{dx} P (\bar{X} \leq x). Why do we take a derivative?

Thank you.
 
Last edited:
Ray Vickson said:
We have ##\bar{X} = S_n/n,## where ##S_n = \sum_{i=1}^n X_i##. The random variable ##S_n## has an n-Erlang distribution with mean ##n \theta##. (Erlang is a special case of Gamma; its density is widely available in textbooks and on line.) To get the distribution pdf ##f(x)## of ##\bar{X}##, use
f(x) =\frac{d}{dx} P (\bar{X} \leq x), \text{ and }<br /> P(\bar{X} \leq x) = P(S_n \leq nx).

So think I got it.

First to make things easier I'm just going to call \frac{1}{θ} = λ

So..
∫^{nx}_{0} \frac{λ^{n} s^{n-1} e^{-λs}}{\Gamma (n)} ds

where the sn-1 cancel out right?

∫^{nx}_{0} λn e-λs + s ds after working with the e-s in the denominator of the integrand.

∫^{nx}_{0} \frac{λ^{n} s^{n-1} e^{-λs}}{\Gamma (n)} ds

\frac{-λ^{n}}{λ-1} e-s(λ-1)|^{nx}_{0}

After taking the derivative
\frac{(-λ^{n})(-n)(λ-1)e^{-nx(λ-1)}}{λ-1}

cancel out like terms and bring e^{-nx} back down to the denominator, and since \frac{(nx)^{n-1}}{(nx)^{n-1}} = 1 we can put that back in so

\frac{(nx)^{n-1}(λ^n)(n)e^{-nxλ}}{(nx)^{n-1}e^{-nx}}

How do I get the integral sign back into the denominator? Also I still don't know why we integrated then took the derivative. Clarification on that would help. Thank you.
 
tmbrwlf730 said:
So think I got it.

First to make things easier I'm just going to call \frac{1}{θ} = λ

So..
∫^{nx}_{0} \frac{λ^{n} s^{n-1} e^{-λs}}{\Gamma (n)} ds

where the sn-1 cancel out right?

∫^{nx}_{0} λn e-λs + s ds after working with the e-s in the denominator of the integrand.

∫^{nx}_{0} \frac{λ^{n} s^{n-1} e^{-λs}}{\Gamma (n)} ds

\frac{-λ^{n}}{λ-1} e-s(λ-1)|^{nx}_{0}

After taking the derivative
\frac{(-λ^{n})(-n)(λ-1)e^{-nx(λ-1)}}{λ-1}

cancel out like terms and bring e^{-nx} back down to the denominator, and since \frac{(nx)^{n-1}}{(nx)^{n-1}} = 1 we can put that back in so

\frac{(nx)^{n-1}(λ^n)(n)e^{-nxλ}}{(nx)^{n-1}e^{-nx}}

How do I get the integral sign back into the denominator? Also I still don't know why we integrated then took the derivative. Clarification on that would help. Thank you.

If ##g(w)## is the density function of some random variable ##W##, how do we find the density function of ##W/n##? One way is to do it is to differentiate the cdf of ##W/n##. However, if you prefer to use formulas for change-of-variables in probability you can do that instead. You should to it first for a general pdf ##g(w)##, then specialize this to the case of the Erlang density.

I cannot figure out what you are doing with all your 'cancellations' and whatnot.
 
Ray Vickson said:
If ##g(w)## is the density function of some random variable ##W##, how do we find the density function of ##W/n##? One way is to do it is to differentiate the cdf of ##W/n##. However, if you prefer to use formulas for change-of-variables in probability you can do that instead. You should to it first for a general pdf ##g(w)##, then specialize this to the case of the Erlang density.

I cannot figure out what you are doing with all your 'cancellations' and whatnot.



What I was thinking that something had to cancel so that I can integrate the function. I haven't worked with the gamma function so I'm not sure about some things.

So the \Gamma(n) = ∫ sn-1 e-s ds.

So we get ∫\frac{λ^{n}s^{n-1}e^{-λs}}{∫s^{n-1}e^{-s}ds}ds

Would the sn-1 in the numerator cancel out with the one in the denominator?
Could you group e-λs and e-stogether to make e-s(λ+1) ?
How do I work with the gamma function, an integral, being a denominator of another integral?
 
Last edited:

Similar threads

  • · Replies 6 ·
Replies
6
Views
2K
Replies
2
Views
2K
Replies
4
Views
1K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 15 ·
Replies
15
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 8 ·
Replies
8
Views
2K
Replies
3
Views
2K
  • · Replies 10 ·
Replies
10
Views
4K
  • · Replies 1 ·
Replies
1
Views
5K