Show E$\big[\frac{\bar{X}(1-\bar{X})}{n}\big] = \frac{(n-1)p((1-p)}{n^2}$ I am trying to show E$\big[\frac{\bar{X}(1-\bar{X})}{n}\big] = \frac{(n-1)p((1-p)}{n^2}$, for $X_1, X_2,\ldots, X_n\sim\operatorname{b}(1,p)$.
I can get to E$\big[\frac{\bar{X}(1-\bar{X})}{n}\big]$ = $\frac{1}{n^2}(np) - (\frac{1}{n}E\big[\bar{X}^2\big])$,which I'm pretty sure is correct, but I'm having trouble simplifying E$\big[\bar{X}^2\big]$. 
I have: \begin{align} 
\frac{1}{n} E\big[\bar{X}^2\big] &= \frac{1}{n^3} E\big[\sum_{i=1}^n X_i^2\big]\\ 
&= \frac{1}{n^3} E\big[(\sum_{i=1}^n X_i)(\sum_{i=1}^{n}X_i)\big]\\ 
&= \frac{1}{n^3} \Big(\sum_{i=1}^n\big(E(X_i)E(X_i)\Big)\\
&= \frac{1}{n^3} \Big(\sum_{i=1}^n(np)^2\Big)\\
&= \frac{1}{n^3} (n^3p^2)\\
&= p^2 \end{align}   If I use $p^2$, however, then I get $(\frac{1}{n^2})np - p^2 = \frac{p(1-np)}{n}$, which is incorrect. If someone could explain where I've gone wrong I'd greatly appreciate the help. Thank you very much.
 A: If you know that $$n\bar X \sim \operatorname{Binomial}(n,p),$$ so that $$\operatorname{E}[n \bar X] = np, \quad \operatorname{Var}[n \bar X] = np(1-p),$$ then the calculation is straightforward:  $$\begin{align*} \operatorname{E}\left[\frac{\bar X (1-\bar X)}{n}\right] &= \operatorname{E}\left[\frac{n \bar X}{n^2} - \frac{(n \bar X)^2}{n^3}\right] \\
&= \frac{1}{n^2} \operatorname{E}[n \bar X] - \frac{1}{n^3}\operatorname{E}[(n \bar X)^2] \\
&= \frac{np}{n^2} - \frac{1}{n^3} \left( \operatorname{Var}[n \bar X] + \operatorname{E}[n \bar X]^2 \right) \\
&= \frac{p}{n} - \frac{np(1-p) + (np)^2}{n^3} \\
&= \frac{np - p(1-p) - n p^2}{n^2} \\
&= \frac{(n-1)p(1-p)}{n^2}.
\end{align*}$$
If we do not assume the formulas for mean and variance of a binomial distribution are known, then we observe that $$\operatorname{E}[n \bar X] = \operatorname{E}\left[\sum_{i=1}^n X_i\right] = \sum_{i=1}^n \operatorname{E}[X_i] = \sum_{i=1}^n p = np.$$  Similarly $$\begin{align*}\operatorname{E}[(n \bar X)^2] &= \operatorname{E}\left[\biggl(\sum_{i=1}^n X_i\biggr)^2\right] \\ &= \operatorname{E}\left[\sum_{i=1}^n \sum_{j=1}^n X_i X_j \right] \\ &= \sum_{i=1}^n \sum_{j=1}^n \operatorname{E}[X_i X_j].\end{align*} $$  We then note that if $i \ne j$, $X_i$ and $X_j$ are independent, thus $$\operatorname{E}[X_i X_j] \overset{\text{ind}}{=} \operatorname{E}[X_i] \operatorname{E}[X_j] = p^2.$$  When $i = j$, then $X_i = X_j = X_i^2$, and because $X_i \sim \operatorname{Bernoulli}(p)$, $X_i^2 = X_i$ (i.e. $X_i = 1$ if and only if $X_i^2 = 1$, and $X_i = 0$ otherwise) and $\operatorname{E}[X_i^2] = p$.  Thus our second moment is $$\operatorname{E}[(n \bar X)^2] = np + n(n-1)p^2 = np(1-p) + np^2,$$ since there are exactly $n$ cases where $i = j$, and $n(n-1)$ cases where $i \ne j$.
