Is the variance concave?

Let $X$ be a discrete random variables with values in the set $\{x_1,\ldots, x_n\}\subset\mathbb{R}$. Denote by $p_i$ the probability that $X=x_i$. We can then regard the variance $Var(X)$ as a function of the vector $p\in \Delta^{n-1}\subseteq\mathbb{R}^n$.

My question is: will it be a concave function of $p$? With $n=2$, we get $Var(X)=p_1p_2(x_1-x_2)^2=p_1(1-p_1)(x_1-x_2)^2$ which is concave. I'm not sure how to generalize this beyond two dimensions though.

We have $\mathrm{var}_p\left(X\right)=\mathbb{E}_p\left[X^2\right]-\mathbb{E}_p\left[X\right]^2$. Furthermore, for $\lambda\in\left[0,1\right]$ we can write $$\mathbb{E}_{\lambda p+(1-\lambda)q}\left[X^2\right]=\lambda\mathbb{E}_p\left[X^2\right]+(1-\lambda)\mathbb{E}_q\left[X^2\right].$$ Convexity of $\left(\cdot\right)^2$ and Jensen's inequality then yield $$\lambda \mathbb{E}_p\left[X\right]^2+(1-\lambda) \mathbb{E}_q\left[X\right]^2\geq \left(\lambda\mathbb{E}_p\left[X\right]+(1-\lambda)\mathbb{E}_q\left[X\right]\right)^2=\mathbb{E}_{\lambda p+(1-\lambda)q}\left[X\right]^2.$$
Thus we deduce $$\mathrm{var}_{\lambda p+(1-\lambda)q}\left(X\right)=\mathbb{E}_{\lambda p+(1-\lambda)q}\left[X^2\right]-\mathbb{E}_{\lambda p+(1-\lambda)q}\left[X\right]^2\\\geq\lambda\mathrm{var}_p\left(X\right)+(1-\lambda)\mathrm{var}_q\left(X\right),$$ that means the variance is concave.