For multiple estimation, we currently cache demeaned Y and X for identical data sets.
In case of slide deviations - i.e. a missing variable between fits, we start from scratch.
We could likely economize a lot by relying more on "warm starts", i.e. using the solution from the last iteration
as a starting point for the current demeaning step.