We derive this decomposition as follows, starting with the sample variance S computed by the usual formula:
Replacing the mean vector by the sum which defines
it,
Adding an appropriate form of zero,
which is
Notice that the second sum as exactly S, so finally
or
where is the total number of distinct pairs of data positions, of which
there are
.