I am trying to show that the function
$f(x,\vec{y})=\alpha\ln(1+\exp(x+\vec{y}\cdot \vec{z}))+(1-\alpha)\ln(1+\exp(-x-\vec{y}\cdot\vec{z}))$
is a convex function of $(x,\vec{y})$ (where $x\in\mathbb{R}$,$\vec{y}\in\mathbb{R}^n$, $\alpha\in\{0,1\}$, and $\vec{z}$ is a pre-defined vector in $\mathbb{R}^n$).
I know that it's possible to show that a function of multiple variables is convex if its Hessian matrix is positive semi-definite, but that method seems to be very calculation-intensive, so I was wondering if there is a simpler method that I have overlooked.
I know this would be trivial if the natural log function was convex (as the composition of convex functions is convex, as is the sum of convex functions).
Thanks in advance for any help received!