It is easy to loss precision or run into overflow in the computation of `sumSquares`. Is it only for computing the sample variance? If true, we can simplify this interface to accept variance directly. We have a more stable implementation of variance computation in `DoubleRDDFunctions`.