The idea that always ##\langle f(\vec{x},\vec{p}) \rangle=f(\langle \vec{x} \rangle,\langle \vec{p}##, \rangle) is obviously WRONG.
Take an ideal gas in the classical limit. Then the phase-space distribution function is given by the Maxwell-Boltzmann distribution,
$$\frac{\mathrm{d} N}{\mathrm{d}^3 x \mathrm{d}^3 p}=\frac{1}{(2 \pi \hbar)^3} \exp\left (-\frac{\vec{p}^2}{2m k T} \right)=f(\vec{x},\vec{p}).$$
The average momentum is obviously
$$\langle \vec{p} \rangle=\frac{1}{Z} V \int_{\mathbb{R}^3} \mathrm{d}^3 p \vec{p} f(\vec{p})=0$$
but
$$\langle \vec{p}^2 \rangle=2m \langle E \rangle=3m k T \neq \langle \vec{p} \rangle^2.$$