2
$\begingroup$

I am trying to find out the density of the product $XY$ of two independent Gamma random variables $X \sim \mathrm{Gamma}(k_1, \theta_1)$ and $Y \sim \mathrm{Gamma}(k_2, \theta_2)$, where $k_i$'s are the shape parameters and $\theta_i$'s are the scale parameters.

I know the formula for the density of the product of two independent RV https://en.wikipedia.org/wiki/Product_distribution, but I couldn't evaluate the integral.

I know the resulting density will involve Meijer G-functions, and I found this paper: http://epubs.siam.org/doi/abs/10.1137/0118065, but the Gamma distribution density in this paper is given as $$f_i(x_i)=\frac{1}{\Gamma{\left(b_i\right)}} x_i^{b_i-1} e^{-x_i}$$ with unity scale parameter.

Is there a formula to find the density of the product of two independent Gamma random variables with arbitrary shape and scale parameters?

Any help is highly appreciated.

$\endgroup$
  • $\begingroup$ What is amma RV? $\endgroup$ – MAN-MADE Aug 17 '17 at 2:24
1
$\begingroup$

Let $X \sim \text{Gamma}(a,b)$ with pdf $f(x)$:

enter image description here

and $Y \sim \text{Gamma}(\alpha,\beta)$ be independent with pdf $g(y)$:

enter image description here

Then, the pdf of the product $Z = X Y$ can be obtained as $h(z)$:

enter image description here

where I am using the TransformProduct function from mathStatica/Mathematica to do the nitty-gritties, and BesselK[n,z] denotes the modified Bessel function of the second kind. This is much simpler than requiring MeijerG functions. I should note that I am one of the authors of the software function used.

Quick Monte Carlo check

  • against theoretical solution derived above when $a =2$, $b = 3$, $\alpha = 4$, and $\beta = 1.1$

enter image description here

Looks fine :)

$\endgroup$
  • $\begingroup$ This is perfect! Thank you so much :) Can you point to any documentation about how this function TransformProduct works? $\endgroup$ – Mohamed Ali Aug 17 '17 at 23:18

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.