Deprecate order kwarg in factorize
ENH: Use Welford's method in stats.moments.rolling_var

This PR implements a modified version of Welford's method to compute
the rolling variance. Instead of keeping track of the sum and sum of
the squares of the items in the window, it tracks the mean and the sum
of squared differences from the mean. This turns out to be (much) more
numerically stable.

The formulas to update these two variables when adding or removing an
item from the sequence are well known, see e.g.

http://en.wikipedia.org/wiki/Algorithms_for_calculating_variance

The formulas used when both adding one and removing one item I have
not seen explicitly worked out anywhere, but are not too hard to come
up with if you put pen to (a lot of) paper.
Merge pull request #6930 from TomAugspurger/factorize-kwarg

API/DEPR: Deprecate order kwarg in factorize
Merge pull request #6817 from jaimefrio/rolling_var-welford

ENH: Use Welford's method in stats.moments.rolling_var