Let $X_1, \cdots, X_n$ be $iid$ normal random variables with unknown mean $\mu$ and known variance $\sigma^2$. How to find $E[\Phi(\bar X)]$, where $\bar X:=\frac{\sum_{i=1}^nX_i}{n}$, please? I guess the answer should be $\Phi(\mu)$. Here is how I started. Note that $Y:= \bar X$ is also normal with mean $\mu$ and variance $\frac{\sigma^2}{n}$.
$$E[\Phi(\bar X)] = \int_{-\infty}^\infty \Phi(y) f_Y(y) dy.$$
Let $z=\frac{y-\mu}{\sigma/\sqrt{n}}$ and the above integration becomes
$$\int_{-\infty}^\infty \Phi(\mu+\frac{\sigma z}{\sqrt{n}}) \phi(z) dz$$
I could not proceed any further from here. Does anyone know what to do, please? Thank you!
I don't know how you got to the displayed integral in your question, but here is how to evaluate it once you get there. The result does not match what you claim is the answer.
Let $Y = \Phi(X)$ where $X \sim N(\mu, \sigma^2)$. Then, $$\begin{align} E[Y] = E[\Phi(X)] &= \int_{-\infty}^\infty \Phi(x)f_X(x)\,\mathrm dx \tag{1} \end{align}$$ But, if $Z\sim N(0,1)$ is independent of $X$, then $$P\{Z \leq x\mid X = x\} = P\{Z \leq x\} = \Phi(x)$$ and so we can recognize the integral in $(1)$ as $$\begin{align} E[\Phi(X)] &= \int_{-\infty}^\infty \Phi(x)f_X(x)\,\mathrm dx\\ &= \int_{-\infty}^\infty P\{Z \leq x\mid X = x\}f_X(x)\,\mathrm dx\\ &= P\{Z \leq X\}&{\scriptstyle{\text{by the continuous version of the law of total probability}}}\\ &= P\{Z-X \leq 0\}\\ &= \Phi\left(\frac{0-(-\mu)}{\sqrt{1 + \sigma^2}}\right) &{\scriptstyle{\text{this follows because}~ Z-X \sim N(-\mu,1+\sigma^2)}}\\ &= \Phi\left(\frac{\mu}{\sqrt{1 + \sigma^2}}\right) \end{align}$$