SHORT VERSION: how would I go about solving the following inverse problem:
$F(a,a_v,b,b_v,c,c_v)=\frac{A B}{C}$ where $A,B,C$ are gamma distributed according to A's mean=$a$, A's variance=$a_v$, etc. I'd like to determine the gamma distribution shape parameters from the shape of $F$.
LONG VERSION:
I have an output probability distribution function $F$ which is produced by taking a simple equation $\frac{a b}{c}$ and varying each of those constants according to a gamma distribution. For simplicity, I can assume that these distributions are independent of each other, though I'd be interested in the case where they are not.
Numerically, it is simple to generate $F$. I randomly select a value $a_i$ from a gamma distribution with mean=$a$ and variance=$a_v$. I do the same with $b$ and $c$. Then I compute and record $F_i=\frac{a_i b_i}{c_i}$. I do this many times until I have enough $F_i$'s to build a histogram, and ultimately the PDF of $F$ which is a function of 6 parameters, the three means and three variances of the gamma distribution inputs.
The problem is that I'm trying to derive the means and variances of the input gamma distributions from the output distribution. I've searched extensively, and have not found anyone with a close-formed result for the PDF of the product / division of gamma distributions. Is it possible to derive these functions?
Alternatively, is it possible to derive moments of the $F$ as a function of the gamma distribution parameters? This would be almost as helpful.
OPTIONAL BACKGROUND One reason I'd love to work these out analytically is that, for example, it's clear to me that there are multiple, or even an infinite number of child gamma distributions that can produce one $F$. Fortunately, I have more information about my system to help determine the parameters. Unfortunately, numerical solvers are very unhappy when the "solution" is a plane or line rather than being discrete minima. If I could set constraints on the system analytically, it would be very helpful.