# How calculate the probability density function of $Z = X_1/X_2$

Let $X_1$ and $X_2$ be two continuous r.v., my question is: what is the p.d.f of $Z=X_1/X_2$?

-
Funny how no one has pointed out yet that this question is not even well-posed without at least one additional assumption. – cardinal Nov 1 '11 at 23:18
@cardinal OK, I edited my answer to say that it is for the case when $X$ and $Y$ are jointly continuous, which I suspect is what the OP meant when he said that $X$ and $Y$ are continuous random variables. I also added a simpler argument for the case when the joint density has circular symmetry and thus $X$ and $Y$ are not independent except in the special case when their marginal densities are $N(0,\sigma^2)$. What else am I missing? – Dilip Sarwate Nov 2 '11 at 12:16
Funny how all the answers to this post, including the accepted one, quietly exhibit/compute a density for the distribution of $Z$ although (as @cardinal briefly signaled it) nothing in the post ensures the existence of a density which in fact might very well not exist. – Did Nov 6 '11 at 20:25

Let's first assume that $X_1$ and $X_2$ are independent continuous random variables.

The cumulative distribution function for $Z$ is $$\begin{eqnarray} F_Z(z) &=& \mathbb{P}( X_1/X_2 \le z) = \mathbb{P}( X_1 \le X_2 z ; X_2 >0) + \mathbb{P}( X_1 \ge X_2 z ; X_2< 0) \\ &=& \mathbb{E}( F_{X_1}( z X_2) ; X_2 > 0) + \mathbb{E}( 1 - F_{X_1}( z X_2) ; X_2 < 0) \end{eqnarray}$$

The cumulative function $F_Z(z)$ is continuous, so the probability density for $Z$ is found by differentiation $f_Z(z) = F_Z^\prime(z)$: $$f_Z(z) = \mathbb{E}( f_{X_1}( z X_2) X_2 ; X_2 > 0) + \mathbb{E}( -f_{X_1}( z X_2) X_2 ; X_2 < 0) = \mathbb{E}\left( \vert X_2\vert \cdot f_{X_1}(z X_2)\right)$$

Example: Let's work out an example. Let $X_1$ and $X_2$ be independent random variables each from standard normal distribution. Then $$\begin{eqnarray} f_Z(z) &=& \mathbb{E}\left( \vert X_2\vert \frac{1}{\sqrt{2\pi}} \mathrm{e}^{-\frac{1}{2} z^2 X_2^2} \right) \\ &=& \int_{-\infty}^\infty \vert x\vert \frac{1}{\sqrt{2\pi}} \mathrm{e}^{-\frac{1}{2} z^2 x^2} \cdot \frac{1}{\sqrt{2\pi}} \mathrm{e}^{-\frac{1}{2} x^2} \mathrm{d} x \\ &\stackrel{\text{use parity}}{=}& \frac{1}{2 \pi} 2 \int_0^\infty x \exp\left( -\frac{x^2}{2} \left(1 + z^2 \right) \right) \mathrm{d} x \\ &\stackrel{t = x^2}{=}& \frac{1}{\pi} \int_0^\infty \exp(-t \left(1+z^2 \right)) \mathrm{d} t \\ &=& = \frac{1}{\pi} \frac{1}{1+z^2} \end{eqnarray}$$

And we obtained the expected density of the Cauchy distribution.

Dependent variables

The above derivation carries through almost the same for dependent variables, whose measure is absolutely continuous, that is when $\mathrm{d}F_{X_1,X_2}(x_1,x_2) = f_{X_1,X_2}(x_1,x_2) \mathrm{d}x_1 \mathrm{d}x_2$. Dilip worked it out to be $$f_Z(z) = \int_{-\infty}^\infty \vert x_2\vert f_{X_1,X_2}(z x_2, x_2) \mathrm{d}x_2 = \mathbb{E}\left( \vert X_2\vert \cdot f_{X_1|X_2}(z X_2, X_2)\right)$$

Absolute continuity of the joint measure is sufficient, but not necessary, condition to ensure that $Z$ is a continuous random variable, i.e. that $Z$ also has absolutely continuous measure.

It is instructive to work out few examples to see what happens when the joint measure is not absolutely continuous. Borrowing an example from this post, consider $(X_1,X_2) = (\sin(U), \cos(U))$, where $U$ follows uniform distribution on $\left(-\frac{\pi}{2}, \frac{\pi}{2}\right)$ interval. Then, although the joint density is not absolutely continuous, $Z = X_1/X_2 = \tan(U)$ is continuous random variable and follows the standard Cauchy distribution.

However, things can go awry. Consider this example from Didier Piau, and let $(X_1,X_2) = (Y,Y)$, where $Y$ follows, say, standard normal distribution. In this case $Z = X_1/X_2 = 1$, and $F_Z(z) = \mathbf{1}_{z \ge 1}$, so the measure is not absolutely continuous.

-
Please: There is no such thing as a cumulative density function. The word "cumulative" contradicts the word "density". There are probability density functions and there are cumulative probability distribution functions. There are no "cumulative density functions". – Michael Hardy Nov 1 '11 at 20:10
Yes, I mistyped. It should be cumulative distribution function. Thanks for point this out, Michael. I will correct this momentarily. – Sasha Nov 1 '11 at 20:12
Are you assuming that $X$ and $Y$ are independent? For (positive) random variables, you get $$f_Z(z) = \int_0^\infty x_2f_{X_1}(zx_2)f_{X_2}(x_2)\mathrm dx_2$$ while my calculation gives $$f_Z(z) = \int_0^\infty x_2f_{X_1, X_2}(zx_2,x_2)\mathrm dx_2$$ which reduces to yours for independent random variables. – Dilip Sarwate Nov 2 '11 at 13:06
The random variable $Z=X_1/X_2$ is well defined on a set of full measure as soon as $X_2\ne0$ almost surely, and this holds if $X_2$ is continuous, hence no, this is not the problem. Rather, as I wrote explicitely in the first comment I referred you to and, before that, as @cardinal had pointed out succinctly without ANY reaction from the OP nor from you (a fact which I find BAFFLING, compare to Dilip's attitude), nothing ensures the existence of the density $f_{X_1,X_2}$ your solution is based on, nor of $f_Z$. Here is a simple example: $X_1$ standard gaussian and $X_2=X_1$. – Did Nov 13 '11 at 8:31
The trouble begins at the first displayed equation since in general $P(X\leqslant \frac12Y)$ is not $E(F_X(\frac12Y))$. If $X=Y$ is uniform on $(0,1)$, the former is $0$ but $F_X(y)=y$ for every $y$ in $(0,1)$ hence the latter is $E(\frac12Y)=\frac14$. // If I were you I would rewrite the post entirely, with a first part assuming independence (including the Cauchy example, which is a good one) and a second part with the counterexamples .../... – Did Nov 13 '11 at 14:46

As an alternative to Sasha's answer, with $Z = Y/X$, $F_Z(z)= P\{Y/X \leq z\}$ is the total probability mass in the region of the $x$-$y$ plane where $y \leq zx$. Although the OP did not state so, I assume that he meant that $X$ and $Y$ are jointly continuous, in which case this probability can be obtained by integrating the joint density function $f_{X,Y}(x,y)$ over this region. Sketching the region in question (for ease in setting up the integrals) we have \begin{align*} F_Z(z) &= \int_{x=0}^{\infty}\int_{y=-\infty}^{zx} f_{X,Y}(x,y) \mathrm dy\ \mathrm dx + \int_{x=-\infty}^{0}\int_{y=zx}^{\infty} f_{X,Y}(x,y) \mathrm dy\ \mathrm dx,\\ f_Z(z) = \frac{\mathrm d}{\mathrm dz}F_Z(z) &= \int_{0}^{\infty} x\cdot f_{X,Y}(x,zx) \mathrm dx - \int_{-\infty}^{0} x\cdot f_{X,Y}(x,zx) \mathrm dx, \end{align*} via the formula for "differentiating under the integral sign".

When the joint density of $X$ and $Y$ has circular symmetry about the origin, then $X$ and $Y$ are in general, dependent random variables except in the special case when the marginal densities of $X$ and $Y$ are zero-mean normal densities with identical variance. For circularly symmetric joint densities, the integrals above can be evaluated readily by switching to polar coordinates, and even more simply by noting that the volume under the surface in the region of integration is a linear function of the angle $\theta$ between the line $y = zx$ and the $x$ axis, with the volume being $0$ at angle $\theta = -\pi/2$ and $1$ at $\theta = \pi/2$. Since $\theta = \arctan(z)$, we get $$F_Z(z) = \frac{1}{2} + \frac{1}{\pi}\arctan(z)$$ and thus $Z = Y/X$ has a Cauchy density if $X$ and $Y$ have a joint density that is circularly symmetric about the origin.

The more detailed calculation is as follows. If $f_{X,Y}(x,y) = g(r)$, then since $$\int_{x=-\infty}^{\infty}\int_{y=-\infty}^{\infty} f_{X,Y}(x,y) \mathrm dy\ \mathrm dx = \int_0^{\infty}\int_0^{2\pi} g(r)\cdot r\cdot \mathrm d\theta\ \mathrm dr = 1$$ so that $\int_0^{\infty} g(r)\ \mathrm dr = 1/2\pi$, we have \begin{align*} F_Z(z) &= \int_{x=0}^{\infty}\int_{y=-\infty}^{zx} f_{X,Y}(x,y) \mathrm dy\ \mathrm dx + \int_{x=-\infty}^{0}\int_{y=zx}^{\infty} f_{X,Y}(x,y) \mathrm dy\ \mathrm dx\\ &= 2\int_{\theta=-\pi/2}^{\arctan(z)}\int_{r=0}^{\infty} r \cdot g(r) \mathrm dr \ \mathrm d\theta\\ &= 2(\arctan(z) + \pi/2)/(2\pi) = \frac{1}{2} + \frac{1}{\pi}\arctan(z) \end{align*} and the Cauchy density is obtained upon differentiating.

-
I am not quite clear about the step using "differentiating under the integral sign". – Fan Zhang Nov 2 '11 at 16:46
If $$F(\alpha)=\int_{a(\alpha)}^{b(\alpha)}f(x;\alpha)\mathrm dx$$ where $a(\cdot)$ and $b(\cdot)$ are differentiable at $\alpha$ and $\frac{\partial f(x;\alpha)}{\partial \alpha}$ exists, then $$\frac{\mathrm dF(\alpha)}{\mathrm d\alpha}=\int_{a(\alpha)}^{b(\alpha)}\frac{\partial f(x;\alpha)}{\partial\alpha}\ \mathrm dx+f(b(\alpha);\alpha)\frac{\mathrm db(\alpha)}{\mathrm d\alpha}-f(a(\alpha);\alpha)\frac{\mathrm db(\alpha)}{\mathrm d\alpha}.$$ Apply this result. In outer integral, limits don't depend on $z$, giving integral of $\frac{d}{dz}$(inner integral) $= \int\pm xf_{X,Y}(x,zx)dx$ – Dilip Sarwate Nov 2 '11 at 17:07