How to calculate $\int_{-\infty}^{\infty} \Phi(y+c)^{k-1}d(\Phi(y))$

Question

How to calculate $\int_{-\infty}^{\infty} \Phi(y+c)^{k-1}d(\Phi(y))$?When I check the list of integrals of Gaussian functions, I only find k-1=2.

Answer 1

Let $X_0, X_1, \cdots, X_{k-1}$ denote $k$ independent unit-variance normal random variables, and suppose that $E[X_0]=0$, and $E[X_i]=-c$, $1 \leq i \leq k-1$. Let $C$ denote the event that $X_0 > \max_i X_i$. Then, $$P(C\mid X_0 = y) = P(\max_i X_i \leq y) = \prod_{i=1}^{k-1}P(X_i \leq y) = [\Phi(y+c)]^{k-1}.$$ Consequently, $$P(C) = \int_{-\infty}^\infty [\Phi(y+c)]^{k-1}\,\mathrm d[\Phi(y)] = \int_{-\infty}^\infty [\Phi(y+c)]^{k-1}\phi(y)\,\mathrm dy$$ where $\phi(y)$ is the standard normal density function. As whuber says in a comment, there is no closed-form expression for the value of this integral but it is nicely behaved and amenable to numerical evaluation. Evaluating integrals of this type is of great interest to communications engineers -- more truthfully, accurate evaluation of $1-P(C)$ is of great interest to communications engineers-- because $1-P(C)$ corresponds to the probability of error in orthogonal signaling. See, for example, this answer of mine for some details.