Shrinkage estimation of high-dimensional factor models with structural instabilities