13

How can one evaluate the expectation of the squared normal CDF in closed-form?

$$\mathbb{E}\left[\Phi\left(aZ+b\right)^{2}\right] = \int_{-\infty}^{\infty}\Phi\left(az+b\right)^{2}\phi(z)\,dz$$

Here, $a$, $b$ are real numbers, $Z\sim\mathcal{N}(0,1)$, and $\phi(\cdot)$ and $\Phi(\cdot)$ are the density and distribution functions of a standard normal random variable, respectively.

Andrei
  • 143
  • 7
  • Well where do you get stuck? Have you tried to evaluate it? Maybe use the fact that $\text{Var}(g(X))=E[g(X)^2]-(E[g(X)])^2$ – stoched Mar 23 '15 at 00:33
  • I tried to evaluate the integral, using integration by parts and other (simple) techniques, but that did not lead me anywhere. Also, I actually started from the variance to get here. I found a similar question (http://stats.stackexchange.com/questions/61080/how-can-i-calculate-int-infty-infty-phi-left-fracw-ab-right-phiw), but extending to the squared CDF does not seem to be trivial. – Andrei Mar 23 '15 at 00:35
  • Have you considered using polar coordinates? – StatsStudent Mar 23 '15 at 00:57
  • No I haven't, can you detail a bit? – Andrei Mar 23 '15 at 01:35
  • 1
    If $b=0$ and $a=1$, then $\Phi(Z)$ is uniformly distributed between 0 and 1. Its second moment is then $1/3$. I recall trying to calculate something like what you ask for general $a$ and $b$, but I found no closed-form solutions. – StijnDeVuyst Mar 23 '15 at 02:32
  • See here: http://stats.stackexchange.com/questions/61080/how-can-i-calculate-int-infty-infty-phi-left-fracw-ab-right-phiw?lq=1 – StatsStudent Mar 23 '15 at 02:52
  • Thank you, but I already looked there and cited it above, in a previous comment. – Andrei Mar 23 '15 at 03:21
  • Look in Wikipedia under list of Gaussian integrals. – soakley Mar 26 '15 at 23:50

1 Answers1

6

As noted in my comment above, check Wikipedia for a list of integrals of Gaussian functions. Using your notation, it gives $$\int_{-\infty}^{\infty} \Phi (az+b)^2 \phi(z) dz=\Phi \left( {{b} \over {\sqrt{1+a^2}} }\right) - 2T \left( {{b} \over {\sqrt{1+a^2}}} \ , {{1} \over {\sqrt{1+2a^2}}} \right) ,$$ where $T(h,q)$ is Owen's T function defined by $$T(h,q)=\phi(h) \int_0^q {{\phi(hx) } \over {1+x^2}} dx$$

If you plug in $a=1,b=0$ you will get $1 \over 3$ as the comments indicate you should.

soakley
  • 4,341
  • 3
  • 16
  • 27