In this answer, it says that in general the sum of the variances is not equal to the variance of the sum.
I tried to work it out by myself, and I think I got a different result, namely that the variance of the sum is equal to the sum of the variances. Do I have incorrect work here for random variables in general? Am I missing some assumption I made?
\begin{eqnarray*} Var(X+Y) & = & E[(X+Y)^{2}]-(E[X+Y])^{2}\\ & = & E[X^{2}+Y^{2}+XY]-(E[X]+E[Y])^{2}\\ & = & E[X^{2}]+E[Y^{2}]+E[XY]-(E[X])^{2}-(E[Y])^{2}-E[X]E[Y]\\ & = & E[X^{2}]-(EX)^{2}+E[Y^{2}]-(EY)^{2}\\ & = & Var(X)+Var(Y) \end{eqnarray*}