(Sums of Independent Poisson Random Variables) Let X and Y be independent Poisson random variables with respective means λ_{1} and λ_{2}. Calculate the distribution of X + Y.
Since the event {X + Y = n} may be written as the union of the disjoint events {X = k, Y = n − k}, 0 ≤ k ≤ n, we have
P\{X+Y=n\}=\sum_{k=0}^{n}P\{X=k,\ Y=n-k\}=\sum_{k=0}^{n}P\{X=k\}P\{Y=n-k\}
=\sum_{k=0}^{n}e^{-\lambda_{1}}{\frac{\lambda_{1}^{k}}{k!}}e^{-\lambda_{2}}{\frac{\lambda_{2}^{n-k}}{(n-k)!}}
=e^{-(\lambda_{1}+\lambda_{2})}\sum_{k=0}^{n}\frac{\lambda_{1}^{k}\lambda_{2}^{n-k}}{k!(n-k)!}
={\frac{e^{-(\lambda_{1}+\lambda_{2})}}{n!}}\sum_{k=0}^{n}{\frac{n!}{k!(n-k)!}}\lambda_{1}^{k}\lambda_{2}^{n-k}
={\frac{e^{-(\lambda_{1}+\lambda_{2})}}{n!}}(\lambda_{1}+\lambda_{2})^{n}
In words, X_{1}+X_{2} has a Poisson distribution with mean \lambda_{1}+\lambda_{2}.