15
$\begingroup$

As title. Can anyone supply a simple proof that

x \Phi(x) + \Phi'(x) \geq 0 \quad \forall x\in\mathbb{R}

where $\Phi$ is the standard normal CDF, i.e.

$\Phi(x) = \int_{-\infty}^x \frac{1}{\sqrt{2\pi}} e^{-y^2/2} {\rm d} y$


I have so far:

Defining f(x) = x \Phi(x) + \Phi'(x) we get

\begin{align} f'(x) & = \Phi(x) + x \Phi'(x) + \Phi''(x) \\ & = \Phi(x) + x\Phi'(x) - x\Phi'(x) \\ & = \Phi(x) \\ & >0 \end{align}

so it seems that if we can show

$\lim_{x\to-\infty} f(x) = 0$

then we have our proof - am I correct?

Clearly $f$ is the sum of two terms which tend to zero, so maybe I have all the machinery I require, and I just need to connect the parts in the right way! Assistance will be gratefully received.


In case anyone is interested in where this question comes from:

Bachelier's formula for an option struck at $K$ with time $T$ until maturity, with volatility $\sigma>0$ and current asset price $S$ is given by

V(S) = (S - K) \Phi\left( \frac{S-K}{\sigma S \sqrt{T}} \right) + \sigma S \sqrt{T} \Phi' \left( \frac{S-K}{\sigma S \sqrt{T}} \right)

Working in time units where $\sigma S\sqrt{T} = 1$ and letting $x=S-K$, we have

V(x) = x \Phi(x) + \Phi'(x)

and I wanted a simple proof that $V(x)>0$ $\forall x$, i.e. an option always has positive value under Bachelier's model.

  • 0
    I posted a comment with Mill's inequality (see my answer) under a wrong name ("Miller"). I removed it to prevent confusion in the future.2011-09-23

5 Answers 5

22

Writing $\phi(x) = (2\pi)^{-1/2}\exp(-x^2/2)$ for $\Phi^{\prime}(x)$ $ x\Phi(x) + \phi(x) = \int_{-\infty}^x x\phi(t)\mathrm dt + \phi(x) \geq \int_{-\infty}^x t\phi(t)\mathrm dt + \phi(x) = -\phi(t)\biggr\vert_{-\infty}^x + \phi(x) = 0 $

  • 0
    @Srivatsan: Yeah. Why didn't I read this before answering?2011-09-23
6

As was pointed out by J.M in the comment, our worry is to show the inequality for $x<0$. We will begin with the following observations:

\Phi''(x)\geq 0,\:\text{ in }(-\infty,0). Indeed \Phi''(x)=-xe^{-x^2/2}>0,\quad\forall x< 0.

Now, consider the function $g(x)=-e^{-x^2/2}\Phi(x).$ We rearrange the factors to obtain g(x)=-\sqrt{2\pi}\frac{1}{\sqrt{2\pi}}e^{-x^2/2}\Phi(x)=-\sqrt{2\pi}\Phi'(x)\Phi(x). Moreover if we pick $z, by our previous observations and the fact that both \Phi'(x) and $\Phi(x)$ are greater than $0$ we can conclude g(x)-g(z)=-\sqrt{2\pi}(\Phi'(x)\Phi(x)-\Phi'(z)\Phi(z))<0. But then $g$ is strictly decreasing on $(-\infty,0).$ Hence we can write g'(x)<0. This relation translates in -(x\Phi(x)+\Phi'(x))e^{-x^2/2}<0. Divide out by $-e^{-x^2/2}$ to obtain, for any $x<0$, x\Phi(x)+\Phi'(x)>0. Hope this is correct.

4

I will push the OP's attempt to a complete proof. If f(x) = x\Phi(x) + \Phi'(x), then, as the OP notes, f'(x) \geq 0 for all $x$; i.e., $f$ is monotonically increasing. So it suffices to show that $f(x) \to 0$ as $x \to -\infty$, which implies that $f(x) \geq 0$ for $x \in \mathbb R$, since: $ \lim_{u \to -\infty} f(u) \leq f(x). $

Now, to show $f(x) \to 0$ as $x \to -\infty$, it is enough to show that $x\Phi(x) \to 0$ and \Phi'(x) \to 0 as $x \to -\infty$.

  1. \Phi'(x) = \exp\left(-\frac{1}{2} |x|^2 \right) clearly approaches $0$ as $x \to -\infty$.
  2. For $x \leq -2$, let $y = |x| \geq 2$. $ |x \Phi(x)| = y \int_{-\infty}^x e^{-t^2/2} dt = y \int_{y}^\infty e^{-t^2/2} \leq y \int_{y}^\infty e^{-t} dt = y e^{-y} , $ which approaches $0$ as $x \to -\infty$.

So we are done.


Notes on the name. The quantity $ R(x) = e^{x^2/2} \int_x^{\infty} e^{-t^2/2} dt $ is called Mill's ratio (see, for e.g., http://www.jstor.org/stable/2236360). I have also seen the inequality $ \int_{x}^\infty e^{-t^2/2} dt \leq \frac{e^{-x^2/2}}{x} \ \ \ \ \ (x > 0), $ sometimes referred to as Mill's inequality (see Theorem 6 in this lecture notes: http://www.stat.cmu.edu/~larry/=stat705/N3.pdf). Notice that by a change of variables $x \to -x$, this inequality can be seen to be equivalent to the one in the question.

  • 0
    (+1) [Here](http://math.stackexchange.com/questions/28751/proof-of-upper-tail-inequality-for-standard-normal-distribution/28754#28754) is a simple proof of the inequality you state in your addendum. L'Hopital shows that it is asymptotically tight as $x \to \infty$. This inequality (and much more) was known to Laplace.2011-09-23
2

I will concentrate at $\lim_{x \to -\infty} f(x) = 0$, where f(x) = x \Phi(x) + \Phi'(x).

Consider $ g(x) = \sqrt{2 \pi} f(-x) = \mathrm{e}^{-\frac{x^2}{2}} - x \int_{x}^\infty \mathrm{e}^{-\frac{y^2}{2}} \mathrm{d} y \qquad \text{for} \qquad x > 0. $ Then $ g(x) = \mathrm{e}^{-\frac{x^2}{2}} - x \int_{\frac{x^2}{2}}^\infty \mathrm{e}^{-t} \frac{\mathrm{d} t}{\sqrt{2t}} \, > \, \mathrm{e}^{-\frac{x^2}{2}} - \int_{\frac{x^2}{2}}^\infty \mathrm{e}^{-t} \mathrm{d} t = 0 $ where $ x \int_{\frac{x^2}{2}}^\infty \mathrm{e}^{-t} \frac{\mathrm{d} t}{\sqrt{2t}} < \frac{x}{\sqrt{2 \frac{x^2}{2}}} \int_{\frac{x^2}{2}}^\infty \mathrm{e}^{-t} \mathrm{d} t = \int_{\frac{x^2}{2}}^\infty \mathrm{e}^{-t} \mathrm{d} t = \mathrm{e}^{-\frac{x^2}{2}}$ for $x>0$.

On the other hand $g(x) < \mathrm{e}^{-\frac{x^2}{2}}$ by definition, so $\lim_{x \to \infty} g(x) = \lim_{x \to \infty} \sqrt{2 \pi} f(-x)$ vanishes being sandwiched between 0 and $\mathrm{e}^{-\frac{x^2}{2}}$ which also tends to zero for large $x$.

1

Taking the derivative of $\Phi(x)$, \Phi'(x)=\frac{1}{\sqrt{2\pi}}e^{-x^2/2}\tag{1} For $x\ge0$, everything is non-negative, so x\Phi(x)+\Phi'(x)\ge0.

For $x<0$, \begin{align} \Phi(x) &=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^xe^{-t^2/2}\;\mathrm{d}t\\ &\le\frac{1}{\sqrt{2\pi}}\int_{-\infty}^xe^{-t^2/2}\frac{t}{x}\;\mathrm{d}t\\ &=-\frac{1}{x}\frac{1}{\sqrt{2\pi}}e^{-x^2/2}\\ &=-\frac{1}{x}\Phi'(x)\tag{2} \end{align} Therefore, we also get x\Phi(x)+\Phi'(x)\ge0.

  • 0
    I removed the rest anyway.2011-09-23