Déjame hacer mi pregunta más general y clara, tal vez de dar una respuesta, gracias a @p.s, por supuesto!
Deje $f\colon K\subseteq\mathbb{R}^n\times\mathbb{R}\to\mathbb{R}$, dado por
$$
f(\mathbf{a},b)=\int_{\Omega}\! h(\mathbf{a},b,\mathbf{x})g(\mathbf{x}) \,\mathrm{d}\mathbf{x},
$$
donde $\Omega=\{\mathbf{x}\in\mathbb{R}^n\rvert\mathbf{w}\cdot\mathbf{x}+c\geq0\}$. También, vamos a $g\colon\mathbb{R}^n\to\mathbb{R}_+$ ser un no-negativo de la función.
Si $h$ es convexa con respecto a $\mathbf{a}$, $b$ sobre el dominio $\Omega$, $f$ es convexa en a $K$.
Prueba
Deje $\lambda\in[0,1]$$(\mathbf{a}_1,b_1),(\mathbf{a}_2,b_2)\in K$. A continuación,
$$
\begin{aligned}
f(\lambda\mathbf{a}_1+(1-\lambda)\mathbf{a}_2,\lambda b_1+(1-\lambda)b_2) &= \int_\Omega\!h(\lambda\mathbf{a}_1+(1-\lambda)\mathbf{a}_2,\lambda b_1+(1-\lambda)b_2,\mathbf{x})g(\mathbf{x})\,\mathrm{d}\mathbf{x} \\
&\le \int_\Omega\!\left[\lambda h(\mathbf{a}_1,b_1,\mathbf{x}) + (1-\lambda) h(\mathbf{a}_2,b_2,\mathbf{x})\right]g(\mathbf{x})\,\mathrm{d}\mathbf{x} \\
&= \lambda \int_\Omega h(\mathbf{a}_1,b_1,\mathbf{x}) g(\mathbf{x})d\mathbf{x} + (1-\lambda) \int_\Omega\!h(\mathbf{a}_2,b_2,\mathbf{x})g(\mathbf{x})\,\mathrm{d}\mathbf{x} \\
&= \lambda f(\mathbf{a}_1,b_1) + (1-\lambda) f(\mathbf{a}_2,b_2).
\end{aligned}
$$
Entonces, yo creo que esto se puede aplicar a mi problema, ya que la función de $h(\mathbf{a},b,\mathbf{x})=\mathbf{a}\cdot\mathbf{x}+b$ es convexa (en realidad es afín con respecto a $\mathbf{a}$, $b$) con respecto a $\mathbf{a}$, $b$.
Un Contraejemplo?
Deje $h(\mathbf{a},b,\mathbf{x})=\mathbf{a}\cdot\mathbf{x}+b+1$, que es afín a w.r.t. $\mathbf{a},b$ sobre el dominio $\Omega=\{\mathbf{x}\in\mathbb{R}^n \mid \mathbf{a}\cdot\mathbf{x}+b+1 \geq 0\}$. También, vamos a $g$ ser la función de densidad de probabilidad de una $n$-dimensiones Gaussiano con media de $\mu\in\mathbb{R}^n$ y matriz de covarianza $\Sigma\in\mathbb{S}_{++}^n$ (simétrica positiva definida, $n\times n$). Entonces, la función de
$$
f(\mathbf{a},b)=\int_{\Omega}\! h(\mathbf{a},b,\mathbf{x})g(\mathbf{x}) \,\mathrm{d}\mathbf{x}
$$
se argumenta que es para ser convexa de los de arriba "Teorema", pero se evalúa de la siguiente manera
$$
f(\mathbf{a},b)
=
\frac{\mathbf{a}\cdot\mu+b+1}{2}
\left[
1 - \operatorname{fer}\left(-\frac{\mathbf{a}\cdot\mu+b+1}{\sqrt{2\mathbf{a}^T\Sigma\mathbf{a}}}\right)
\right]
+
\frac{\sqrt{\mathbf{a}^T\Sigma\mathbf{a}}}{2\pi}
\exp\left(
-\frac{(\mathbf{a}\cdot\mu+b+1)^2}{2\mathbf{a}^T\Sigma\mathbf{a}}\right),
$$
que, obviamente, no es convexa en a $(\mathbf{a},b)$.
¿Dónde está el error?