2
$\begingroup$

Let $ Ф(x) = \frac{1}{\sqrt{2\pi}} \int_{-\infty}^{x} e^{-t^2/2} \, dt $ be the cumulative distribution function of the standard normal distribution.

Numerical calculations suggest the following inequality is true: $$ 1 - Ф(x) > \frac{2}{\pi} \log( e^{x^2/2} -1 ) e^{-x^2} $$ for $ x > 0 $. The question is how to prove it?

There are some known lower bounds on $ 1 - Ф(x) $ like $$ 1 - Ф(x) > \frac{1}{\sqrt{2\pi}} \frac{x}{x^2+1} e^{-x^2/2} $$ and $$ 1 - Ф(x) > \frac{1}{2\sqrt{2\pi}} (\sqrt{x^2+4}-x) e^{-x^2/2}. $$ In both cases the difference of left and right sides is monotonically decreasing, but in my inequality the difference isn't monotonic.

$\endgroup$
  • $\begingroup$ As Robert noted due to a different exponential your inequality is worse than known for large x. It makes sense for x near zero. But for such values it is more natural to seek for improvements in terms of powers of x, not exponentials. But on this field it seems impossible to compete with known excellent inequalities via continuous fractions and Pade approximants. $\endgroup$ – Sergei Sep 10 '14 at 6:55
  • $\begingroup$ Yes, I know that my inequality is worse than the second lower bound everywhere (though there is an improvement of this inequality that is better than the second bound on some interval). But it came up in my study, so I need to prove it. $\endgroup$ – user47855 Sep 10 '14 at 11:30
3
$\begingroup$

For sufficiently large $x$, the inequality is true, as $$1 - \Phi(x) \sim \dfrac{e^{-x^2/2}}{\sqrt{2\pi} x}$$
while your expression (let's call it $R(x)$) is $O(x^2 e^{-x^2})$. In principle one can get explicit bounds corresponding to these asymptotics and thereby determine a finite value $B$ such that $1 - \Phi(x) > R(x)$ for all $x \ge B$. You then want to look at the interval $[0,B]$. A finite number of computations using interval arithmetic should be able to handle this. Thus if $1 - \Phi(x_0) - R(x_0) \ge a$ and you have an estimate $-\Phi'(x) - R'(x) \le b$ in an interval $[a,c]$, then $1 - \Phi(x) - R(x) \ge a + b (x-x_0)$ in that interval.

$\endgroup$
  • $\begingroup$ Thank you for your answer! "A finite number of computations using interval arithmetic should be able to handle this" Could you please explain it a little more? $\endgroup$ – user47855 Sep 10 '14 at 11:32
  • $\begingroup$ In practice you could use e.g. the Maple function Optimization[Minimize] with method=branchandbound. Or you could use the Maple function gmin which I wrote: <math.ubc.ca/~israel/advisor/advisor6/advisor6.html> $\endgroup$ – Robert Israel Sep 10 '14 at 15:44
  • $\begingroup$ I've already checked this inequality in Mathematica, but now I need to prove it by hand, and it seems that the finite part is a harder one... I would very appreciate your further advice. $\endgroup$ – user47855 Sep 10 '14 at 17:56
  • $\begingroup$ I don't really see the point of doing this sort of computation by hand. If you want something more "proofy", you might try <cl.cam.ac.uk/~jrh13/hol-light> $\endgroup$ – Robert Israel Sep 10 '14 at 21:29

Your Answer

By clicking "Post Your Answer", you acknowledge that you have read our updated terms of service, privacy policy and cookie policy, and that your continued use of the website is subject to these policies.

Not the answer you're looking for? Browse other questions tagged or ask your own question.