In using the Bienaymé formula to find the variance of means, I do not understand why
$$\operatorname{Var}\!\left(\frac{1}{n}\sum_{i=1}^{n}{X_i}\right)=\frac{1}{n^2}\sum_{i=1}^n{\operatorname{Var}(X_i)}$$
I assume it is a matter of simple algebraic manipulation, but I do not understand it. Could someone please explain this property?
Here is a detailed derivation.
So, let's start with the "easy" part: $$ \operatorname{Var} \left(\frac{1}{n}\sum_{i=1}^n X_i\right) = \frac{1}{n^2}\operatorname{Var} \left(\sum_{i=1}^n X_i\right)\,. $$ This is just because $\operatorname{Var} ( a X) = a^2 \operatorname{Var} X$ for any constant $a$ (and r.v. $X$ such that the variance exists).
Now, the "hard" part, where the non-correlation will be used: $$\begin{align*} \operatorname{Var} \left(\sum_{i=1}^n X_i\right) &= \mathbb{E}\left[\left(\sum_{i=1}^n X_i-\sum_{i=1}^n \mathbb{E}[X_i]\right)^2\right] = \mathbb{E}\left[\left(\sum_{i=1}^n (X_i- \mathbb{E}[X_i])\right)^2\right]\\ &= \mathbb{E}\left[\sum_{i=1}^n\sum_{j=1}^n (X_i- \mathbb{E}[X_i])(X_j- \mathbb{E}[X_j])\right]\\ &= \sum_{i=1}^n\mathbb{E}\left[(X_i- \mathbb{E}[X_i])^2\right] + \sum_{i\neq j} \mathbb{E}\left[(X_i- \mathbb{E}[X_i])(X_j- \mathbb{E}[X_j])\right]\\ &= \sum_{i=1}^n\operatorname{Var}(X_i) + \sum_{i\neq j} 0 \\ &= \sum_{i=1}^n\operatorname{Var}(X_i) \end{align*}$$ as claimed. To go from the first to the second line, we expanded the square; to go from the second to the third, we divided the sum in two different sums, and then used linearity of expectation; and from the third to the fourth line, we used the definition of variance, and the assumption of pairwise uncorrelation.