Calculates W = Y - X*beta_est. It is used in the initialization step of the algorithm, to initialise the factorstructures.
Useful links