Take the 2-minute tour ×
Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. It's 100% free, no registration required.

$X \sim \mathcal{P}( \lambda) $ and $Y \sim \mathcal{P}( \mu)$ meaning that $X$ and $Y$ are Poisson distributions. What is the probability distribution law of $X + Y$. I know it is $X+Y \sim \mathcal{P}( \lambda + \mu)$ but I don't understand how to derive it.

share|improve this question

4 Answers 4

up vote 14 down vote accepted

This only holds if $X$ and $Y$ are independent, so we suppose this from now on. We have for $k \ge 0$: \begin{align*} P(X+ Y =k) &= \sum_{i = 0}^k P(X+ Y = k, X = i)\\ &= \sum_{i=0}^k P(Y = k-i , X =i)\\ &= \sum_{i=0}^k P(Y = k-i)P(X=i)\\ &= \sum_{i=0}^k e^{-\mu}\frac{\mu^{k-i}}{(k-i)!}e^{-\lambda}\frac{\lambda^i}{i!}\\ &= e^{-(\mu + \lambda)}\frac 1{k!}\sum_{i=0}^k \frac{k!}{i!(k-i)!}\mu^{k-i}\lambda^i\\ &= e^{-(\mu + \lambda)}\frac 1{k!}\sum_{i=0}^k \binom ki\mu^{k-i}\lambda^i\\ &= \frac{(\mu + \lambda)^k}{k!} \cdot e^{-(\mu + \lambda)} \end{align*} Hence, $X+ Y \sim \mathcal P(\mu + \lambda)$.

share|improve this answer
Thank you! but what happens if they are not independent? –  user31280 Oct 25 '12 at 20:20
In general we can't say anything then. It depends on how they depend on another. –  martini Oct 25 '12 at 20:22
Thank you! it's very simple and I feel like a complete idiot. –  user31280 Oct 25 '12 at 20:40
Nice derivation: specifically the transformation of (a) the i/k factorials and (b) the mu/lambda polynomials into the binomial form of the polynomial power expression. –  javadba Aug 30 at 20:59

Another approach is to use characteristic functions. If $X\sim \mathrm{po}(\lambda)$, then the characteristic function of $X$ is (if this is unknown, just calculate it) $$ \varphi_X(t)=E[e^{itX}]=e^{\lambda(e^{it}-1)},\quad t\in\mathbb{R}. $$ Now suppose that $X$ and $Y$ are independent Poisson distributed random variables with parameters $\lambda$ and $\mu$ respectively. Then due to the independence we have that $$ \varphi_{X+Y}(t)=\varphi_X(t)\varphi_Y(t)=e^{\lambda(e^{it}-1)}e^{\mu(e^{it}-1)}=e^{(\mu+\lambda)(e^{it}-1)},\quad t\in\mathbb{R}. $$ As the characteristic function completely determines the distribution, we conclude that $X+Y\sim\mathrm{po}(\lambda+\mu)$.

share|improve this answer

You can use Probability Generating Function(P.G.F). As poisson distribution is a discrete probability distribution, P.G.F. fits better in this case.For independent X and Y random variable which follows distribution Po($\lambda$) and Po($\mu$). P.G.F of X is \begin{equation*} \begin{split} P_X[t] = E[t^X]&= \sum_{x=0}^{\infty}t^xe^{-\lambda}\frac{\lambda^x}{x!}\\ &=\sum_{x=0}^{\infty}e^{-\lambda}\frac{(\lambda t)^x}{x!}\\ &=e^{-\lambda}e^{\lambda t}\\ &=e^{-\lambda (1-t)}\\ \end{split} \end{equation*} P.G.F of Y is \begin{equation*} \begin{split} P_Y[t] = E[t^Y]&= \sum_{y=0}^{\infty}t^ye^{-\mu}\frac{\mu^y}{y!}\\ &=\sum_{y=0}^{\infty}e^{-\mu}\frac{(\mu t)^y}{y!}\\ &=e^{-\mu}e^{\mu t}\\ &=e^{-\mu (1-t)}\\ \end{split} \end{equation*}

Now think about P.G.F of U = X+Y. As X and Y are independent, \begin{equation*} \begin{split} P_U(t)=P_{X+Y}(t)=P_X(t)P_Y(t)=E[t^{X+Y}]=E[t^X t^Y]&= E[t^X]E[t^Y]\\ &= e^{-\lambda (1-t)}e^{-\mu (1-t)}\\ &= e^{-(\lambda+\mu) (1-t)}\\ \end{split} \end{equation*}

Now this is the P.G.F of $Po(\lambda + \mu)$ distribution. Therefore,we can say U=X+Y follows Po($\lambda+\mu$)

share|improve this answer
Seems a typo in third line (out of 4) in both PGF of X and Y. Should be exp(-u)exp(mut) instead of exp(-u)exp(-1*mut) –  javadba Aug 30 at 21:19
Separate comment/question: please explain why in second line of derivation for PGF of Y refers to exp(x) and x! (instead of referring to exp(y) and y!). I believe these are cut/paste errors - but please confirm. –  javadba Aug 30 at 21:26
Yeah.. You are absolutely correct. Those were typos. I have edited them now. If you find more, let me know. –  Ananda Sep 2 at 5:35

hint: $\sum_{k=0}^{n} P(X = k)P(Y = n-k)$

share|improve this answer
why this hint, why the sum? This is what I don't understand –  user31280 Oct 25 '12 at 20:22
adding two random variables is simply convolution of those random variables. That's why. –  jay-sun Oct 25 '12 at 20:24
gotcha! Thanks! –  user31280 Oct 25 '12 at 20:31
adding two random variables is simply convolution of those random variables... Sorry but no. –  Did Feb 13 '13 at 6:28
There is no usual sense for convolution of random variables. Either convolution of distributions or addition of random variables. –  Did Feb 13 '13 at 6:51

Your Answer


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