# Partial derivative of bivariate normal cdf and pdf

How do I take the partial derivative of bivariate normal cdf and bivariate normal pdf with its arguments (i.e. $x_{1}$ ,$x_{2}$ , and $\rho$ in the following equations)?

$$y=\Phi(x_{1},x_{2},\rho)$$

$$z=\phi(x_{1},x_{2},\rho)$$

where $x_{1}$ is normally distributed with mean 0 and variance 1 and $x_{2}$ is normally distributed with mean 0 and variance 1. $\rho$ is the correlation between $x_{1}$ and $x_{2}$.

• If you were given just the formulas for these functions and not told that they were cdfs or pdfs (in fact, suppose that you had no knowledge whatsoever of probability and/or statistics), could you find the partial derivatives using standard calculus techniques such as the chain rule? Commented Oct 5, 2013 at 13:40
• Thanks. Yes, I know how to do that from high school calculus. My question is how to do when we have cdf and pdf. Commented Oct 5, 2013 at 15:20
• I know how to do that from high school calculus. My question is how to do when we have cdf and pdf. Ignore the information that these are pdfs or cdfs and proceed. The bivariate cdf will be given as a double integral with integrand the bivariate pdf and upper limits $x_1$ and $x_2$. So the partial derivative w.r.t. $x_i$ will be a single integral of the bivariate pdf, etc. Commented Oct 5, 2013 at 16:27

\begin{align} y &= \Phi(x_1,x_2,\rho) = \int_{-\infty}^{x_1}\left[\int_{-\infty}^{x_2} \phi(a,b,\rho)\,\mathrm db\right]\,\mathrm da\\ \frac{\partial y}{\partial x_1} &= \frac{\partial}{\partial x_1}\Phi(x_1,x_2,\rho) = \frac{\partial}{\partial x_1}\int_{-\infty}^{x_1} \left[\int_{-\infty}^{x_2} \phi(a,b,\rho)\,\mathrm db\right]\,\mathrm da\\ &= \int_{-\infty}^{x_2} \phi(x_1,b,\rho)\,\mathrm db \end{align} via the rule for differentiating under the integral sign. Similarly, $$\frac{\partial y}{\partial x_2} = \int_{-\infty}^{x_1} \phi(a,x_2,\rho)\,\mathrm da.$$ If you don't recall the rule for differentiating integrals, see for example, the comments following this answer on math.SE.

The derivative with respect to $\rho$ is straightforward to find but messy in its details. We have that $$\phi(x_1,x_2,\rho)=\frac{1}{2\pi\sqrt{1-\rho^2}} \exp\left[-\frac{x^2 -2\rho xy + y^2}{2(1-\rho^2)}\right]$$ whose partial derivative with respect to $\rho$ is left to the OP to find. If $g(x_1,x_2,\rho)$ denotes this partial derivative, then $$\frac{\partial}{\partial \rho}\Phi(x_1,x_2,\rho) = \frac{\partial}{\partial \rho}\int_{-\infty}^{x_1} \left[\int_{-\infty}^{x_2} \phi(a,b,\rho)\,\mathrm db\right]\,\mathrm da = \int_{-\infty}^{x_1} \int_{-\infty}^{x_2} g(a,b,\rho)\,\mathrm db\,\mathrm da$$

• An anonymous user just proposed an edit: "I'm interested in the derivative with respect to $\rho$", which I declined, because it should have been a comment, which I am posting. Commented Oct 2, 2015 at 8:23
• @StephanKolassa Thanks. I added a few lines to the answer describing how to go about finding the derivative with respect to $\rho$. Commented Oct 2, 2015 at 19:12

While the accepted answer provides great tips for somebody looking to re-derive these quantities, for some work I was doing I actually just needed to know the closed-form partial derivatives of the bivariate normal cdf. Further, the integral for the partial with respect to $$\rho$$ is non-trivial. Since my google search brought me here (and then down a rabbit hole of statistics papers from decades ago), I thought I'd share my findings for the rest of the internet community:

\begin{align*} \frac{\partial}{\partial x_1} \Phi(x_1, x_2, \rho) &= \phi(x_1)\Phi\bigg(\frac{x_2-\rho x_1}{\sqrt{1-\rho^2}}\bigg), \\ \frac{\partial}{\partial x_2} \Phi(x_1, x_2, \rho) &= \phi(x_2)\Phi\bigg(\frac{x_1-\rho x_2}{\sqrt{1-\rho^2}}\bigg),~\text{and} \\ \frac{\partial}{\partial \rho} \Phi(x_1, x_2, \rho) &= \frac{1}{2\pi\sqrt{1-\rho^2}}\exp\bigg\{\frac{-(x_1^2-2\rho x_1x_2+x_2^2)}{2(1-\rho^2)}\bigg\}. \end{align*}

Here, $$\phi(\cdot)$$ and $$\Phi(\cdot)$$ are the standard normal pdf and cdf, respectively. The literature reference for the partial with respect to $$\rho$$ is:

Z. Drezner and G.O. Wesolowsky (1990). On the computation of the bivariate normal integral. Journal of Statistical Computation and Simulation 35 (1-2): 101-107. Equation (4).