Published online by Cambridge University Press: 09 February 2006
The distributions of products and ratios of random variables are of interest in many areas of the sciences. In this note, the exact distributions of the product XY and the ratio X/Y are derived when X and Y are gamma and Weibull random variables distributed independently of each other.The authors thank the referee and the editor for carefully reading the paper and for their help in improving the paper.
For given random variables X and Y, the distributions of the product XY and the ratio X/Y are of interest in many areas of the sciences.
In traditional portfolio selection models certain cases involve the product of random variables. The best examples of this are in the case of investment in a number of different overseas markets. In portfolio diversification models (see, e.g., Grubel, 1968) not only are prices of shares in local markets uncertain, but also the exchange rates are uncertain, and so the value of the portfolio in domestic currency is related to a product of random variables. Similarly in models of diversified production by multinationals (see, e.g., Rugman, 1979) there are local production uncertainty and exchange rate uncertainty, and so profits in home currency are again related to a product of random variables. An entirely different example is drawn from the econometric literature. In making a forecast from an estimated equation Feldstein (1971) pointed out that both the parameter and the value of the exogenous variable in the forecast period could be considered random variables. Hence the forecast was proportional to a product of random variables.
An important example of ratios of random variables is the stress-strength model in the context of reliability. It describes the life of a component that has a random strength Y and is subjected to random stress X. The component fails at the instant that the stress applied to it exceeds the strength, and the component will function satisfactorily whenever Y > X. Thus, Pr(X < Y) is a measure of component reliability. It has many applications, especially in engineering concepts such as structures, deterioration of rocket motors, static fatigue of ceramic components, fatigue failure of aircraft structures, and the aging of concrete pressure vessels.
The distributions of XY and X/Y have been studied by several authors, especially when X and Y are independent random variables and come from the same family. With respect to products of random variables, see Sakamoto (1943) for uniform family, Harter (1951) and Wallgren (1980) for Student's t family, Springer and Thompson (1970) for normal family, Stuart (1962) and Podolski (1972) for gamma family, Steece (1976), Bhargava and Khatri (1981), and Tang and Gupta (1984) for beta family, AbuSalih (1983) for power function family, and Malik and Trudel (1986) for exponential family (for a comprehensive review of known results, see also Rathie and Rohrer, 1987). With respect to ratios of random variables, see Marsaglia (1965) and Korhonen and Narula (1989) for normal family, Press (1969) for Student's t family, Basu and Lochner (1971) for Weibull family, Shcolnick (1985) for stable family, Hawkins and Han (1986) for noncentral chi-squared family, Provost (1989) for gamma family, and Pham-Gia (2000) for beta family.
There is relatively little work of the preceding kind when X and Y belong to different families (the only work known to the authors is that by Nadarajah and Ali, 2004, and Nadarajah and Gupta, 2005). In the applications mentioned before, it is quite possible that X and Y could arise from different but similar distributions. In this note, we study the exact distributions of XY and X/Y when X and Y are independent random variables having the gamma and Weibull distributions with probability density functions (p.d.f.s)
and
respectively, for x > 0, y > 0, α > 0, β > 0, λ > 0, and μ > 0. The calculations involve the incomplete gamma function defined by
and the Gauss hypergeometric function defined by
where (e)k = e(e + 1)…(e + k − 1) denotes the ascending factorial (also known as Pochhammer's symbol). Further properties of the preceding special functions can be found in Prudnikov, Brychkov, and Marichev (1986) and Gradshteyn and Ryzhik (2000).
Theorem 1 expresses the cumulative distribution function (c.d.f.) of XY as an infinite sum of elementary functions.
THEOREM 1. Suppose X and Y are distributed according to (1) and (2), respectively. The c.d.f. of Z = XY can be expressed as
for z > 0.
Proof. The c.d.f. corresponding to (1) is 1 − Γ(α,μx)/Γ(α). Thus, one can write the c.d.f. of XY as
where we have set w = 1/y. Application of equation (2.10.1.12) in Prudnikov et al. (1986, vol. 2) shows that the integral I can be calculated as
The result of the theorem follows by substituting (5) and (4). █
Figure 1 illustrates possible shapes of the p.d.f. of XY for selected values of α and β. The four curves in each plot correspond to selected values of α. The effect of the parameters is evident.
Theorem 2 expresses the c.d.f. of X/Y as an infinite sum of elementary functions.
THEOREM 2. Suppose X and Y are distributed according to (1) and (2), respectively. The c.d.f. of Z = X/Y can be expressed as
for β > 1,
for β < 1, and
for β = 1.
Proof. Because 1 − Γ(α,μx)/Γ(α) is the c.d.f. corresponding to (1), one can write the c.d.f. of X/Y as
Application of equation (2.10.1.5) in Prudnikov et al. (1986, vol. 2) shows that the integral I can be calculated as
Note that the hypergeometric term in (10) simplifies as
The result of the theorem follows by substituting (10) into (9). █
Figure 2 illustrates possible shapes of the p.d.f. of X/Y for selected values of α and β. The four curves in each plot correspond to selected values of α. The effect of the parameters is evident.