# probability generating function of a poisson random variable

• Jun 1st 2007, 06:24 AM
chogo
probability generating function of a poisson random variable
i should start by asking is anyone on the forum well versed in the coalescent theory?

i have a problem in a derivation

if the number of mutations is given by

$\displaystyle P_n(z) = E exp((z-1)\frac{\theta}{2}(nT_n + ..... + 2T_2))$

as one can see this is the probability generating function of a poisson random variable with random mean. where $\displaystyle \lambda = \frac{\theta}{2}(nT_n + ..... + 2T_2)$

now from here i can write this as a product

so

$\displaystyle P_n(z) = \prod_{j=2}^n E\left(e^{(z-1)\frac{\theta}{2}_j T_j}\right)$

Now im confused what to do. If i have the probability generating function how do i find the density function (derivative?)

the next step in the derivation should be

$\displaystyle P_n(z) = \prod_{j=2}^n \left(1 - (z-1)\cdot \frac{j \theta /2}{j(j-1)/2}\right)^{-1}$

I cant for the life of me figure out how this was found from the pgf. PLEASE ANY IDEAS ARE WELCOME, SO IF NO ONE KNOWS THE ANSWER SUGGEST ANYTHING.

many thanks

chogo
• Jun 3rd 2007, 04:16 AM
chogo

$\displaystyle P(z) = E\left(e^{(z-1)\lambda}\right)$

how would i remove the E
• Jun 3rd 2007, 06:25 AM
JakeD
Quote:

Originally Posted by chogo

$\displaystyle P(z) = E\left(e^{(z-1)\lambda}\right)$

how would i remove the E

Are you sure that is the right formula? Because there is no random variable in $\displaystyle e^{(z-1)\lambda},$ it is a constant with respect to the expectation and $\displaystyle P(z) = E\left(e^{(z-1)\lambda}\right) = e^{(z-1)\lambda}.$ When specifying an expectation, the random variable involved should be made clear, for example, by writing $\displaystyle P(z) = E_X z^X$ when the random variable is $\displaystyle X.$
• Jun 3rd 2007, 08:23 AM
chogo
thank you jakeD sorry yeah i should not have put the expectation there. My bad. Was just trying to write something simpler to see if anyone had any suggestions

The expectation is valid in my original equation, where the random variable is $\displaystyle T_j$

the thing is how did the person who did that derivation remove the expectation? and arrive at the second formula.

http://www.mathhelpforum.com/math-he...f1547c7e-1.gif -> http://www.mathhelpforum.com/math-he...4842dd2b-1.gif

any suggestions?

i thank you so much for your help
• Jun 3rd 2007, 02:47 PM
JakeD
Quote:

Originally Posted by chogo
thank you jakeD sorry yeah i should not have put the expectation there. My bad. Was just trying to write something simpler to see if anyone had any suggestions

The expectation is valid in my original equation, where the random variable is $\displaystyle T_j$

the thing is how did the person who did that derivation remove the expectation? and arrive at the second formula.

http://www.mathhelpforum.com/math-he...f1547c7e-1.gif -> http://www.mathhelpforum.com/math-he...4842dd2b-1.gif

any suggestions?

i thank you so much for your help

What is the distribution of the random variable $\displaystyle T_j ?$ You haven't said. The expectation uses that distribution.

Let $\displaystyle M_{T_j}(t)$ be the moment generating function of $\displaystyle T_j .$ Then comparing the first and second equations

$\displaystyle E_{T_j} (e^{(z-1)\theta T_j / 2j}) = M_{T_j}((z-1)\theta/2j ) = (1 - t 2j/(j-1))^{-1}$

where $\displaystyle t = (z-1)\theta / 2j.$ Then

$\displaystyle M_{T_j}(t) = (1 - t 2j/(j-1))^{-1} = \frac{(j-1)/2j}{(j-1)/2j - t}$

which is the MGF of a Gamma(1,(j-1)/2j) distribution, that is, an exponential distribution with parameter $\displaystyle \lambda = (j-1)/2j.$ So it appears that is the distribution of $\displaystyle T_j.$ Is that correct?
• Jun 3rd 2007, 03:01 PM
topsquark
My original comment won't work, but I have to wonder if "E" isn't supposed to be a " $\displaystyle \Sigma$?"

-Dan
• Jun 3rd 2007, 03:14 PM
chogo
firstly JakeD and topsquark i cant thank you enough

Top shark, yes its definitley an E

yes $\displaystyle T_j$ is assumed to be exponentially distributed

also yes this term does cancel out

http://www.mathhelpforum.com/math-he...4842dd2b-1.gif

and becomes

$\displaystyle P_n(z) = \prod^n_{j=2}\left(1-\frac{(z-1)\theta}{j-1}\right)^{-1}$

did he substitute a value for $\displaystyle T_j$?

i dont know why the original equation which is of the form of a poisson generating function becomes what is is now, which you said is a gamma(1,1).

im almost 100% certain this is not wrong, as its a very famous theory developed in mathematical biology and very well founded.

thank you so much for you help again its really much appreciated. If you want i can write the entire derivation our for you guys, will this help?
• Jun 3rd 2007, 03:19 PM
chogo
what did you do here

http://www.mathhelpforum.com/math-he...8fe3363f-1.gif

this seems to be exactly what i need
• Jun 3rd 2007, 03:38 PM
JakeD
Quote:

Originally Posted by chogo
firstly JakeD and topsquark i cant thank you enough

Top shark, yes its definitley an E

yes $\displaystyle T_j$ is assumed to be exponentially distributed

also yes this term does cancel out

http://www.mathhelpforum.com/math-he...4842dd2b-1.gif

and becomes

$\displaystyle P_n(z) = \prod^n_{j=2}\left(1-\frac{(z-1)\theta}{j-1}\right)^{-1}$

did he substitute a value for $\displaystyle T_j$?

i dont know why the original equation which is of the form of a poisson generating function becomes what is is now, which you said is a gamma(1,1).

im almost 100% certain this is not wrong, as its a very famous theory developed in mathematical biology and very well founded.

thank you so much for you help again its really much appreciated. If you want i can write the entire derivation our for you guys, will this help?

This post appears to be looking at my post before I finished editing it. So it responds to comments of mine I resolved and edited out.
• Jun 3rd 2007, 03:50 PM
JakeD
Quote:

Originally Posted by JakeD
Let $\displaystyle M_{T_j}(t)$ be the moment generating function of $\displaystyle T_j .$ Then comparing the first and second equations

$\displaystyle E_{T_j} (e^{(z-1)\theta T_j / 2j}) = M_{T_j}((z-1)\theta/2j ) = (1 - t 2j/(j-1))^{-1}$

where $\displaystyle t = (z-1)\theta / 2j.$

Quote:

Originally Posted by chogo
what did you do here

http://www.mathhelpforum.com/math-he...8fe3363f-1.gif

this seems to be exactly what i need

The explanations before and after this equation are needed. I used the definition of moment generating function and then equated corresponding parts of your first and second equations while substituting in the variable $\displaystyle t.$
• Jun 3rd 2007, 04:14 PM
chogo
thank your help is invaluable. but one last problem. Like you before i used the definition

$\displaystyle M(t) = E(e^{t,x})$

to go from $\displaystyle E_{t_j}(e^{(z-1)\theta T_j/2j}) = M_{t_j}((z-1)\theta T_j/2)$

but when you said you compared the forms of the first and second equations to arrive at the third equation im not fully satistifed.

Is there a formal definition which says the moment generating function i have is equal to $\displaystyle (1-t2_j / (j-1))^{-1}$

where http://www.mathhelpforum.com/math-he...48dd28b5-1.gif
• Jun 3rd 2007, 04:44 PM
JakeD
Quote:

Originally Posted by JakeD
What is the distribution of the random variable $\displaystyle T_j ?$ You haven't said. The expectation uses that distribution.

Let $\displaystyle M_{T_j}(t)$ be the moment generating function of $\displaystyle T_j .$ Then comparing the first and second equations

$\displaystyle E_{T_j} (e^{(z-1)\theta T_j / 2j}) = M_{T_j}((z-1)\theta/2j ) = (1 - t 2j/(j-1))^{-1}$

where $\displaystyle t = (z-1)\theta / 2j.$ Then

$\displaystyle M_{T_j}(t) = (1 - t 2j/(j-1))^{-1} = \frac{(j-1)/2j}{(j-1)/2j - t}$

which is the MGF of a Gamma(1,(j-1)/2j) distribution, that is, an exponential distribution with parameter $\displaystyle \lambda = (j-1)/2j.$ So it appears that is the distribution of $\displaystyle T_j.$ Is that correct?

Quote:

Originally Posted by chogo
thank your help is invaluable. but one last problem. Like you before i used the definition

$\displaystyle M(t) = E(e^{t,x})$

to go from $\displaystyle E_{t_j}(e^{(z-1)\theta T_j/2j}) = M_{t_j}((z-1)\theta T_j/2)$

but when you said you compared the forms of the first and second equations to arrive at the third equation im not fully satistifed.

Is there a formal definition which says the moment generating function i have is equal to $\displaystyle (1-t2_j / (j-1))^{-1}$

where http://www.mathhelpforum.com/math-he...48dd28b5-1.gif

I was trying to deduce what the distribution of $\displaystyle T_j$ was because you didn't say what it was. So I deduced it was an exponential distribution with parameter $\displaystyle \lambda = (j-1)/2j.$ Is this correct?

Now you know what the distribution of $\displaystyle T_j$ is. So you can say what the moment generating function for $\displaystyle T_j$ is. I was working backwards to deduce the distribution; you can work forwards knowing the distribution.

I got that $\displaystyle M_T (t) = \frac{\lambda}{\lambda - t} = (1 - t/\lambda)^{-1}$ when $\displaystyle T \sim \text{Exp}(\lambda)$ from a probability text. See Exponential distribution - Wikipedia, the free encyclopedia .
• Jun 4th 2007, 02:00 AM
chogo
sheer brilliance! thanks alot for the help