Estimate large covariance matrices in approximate factor models by thresholding principal orthogonal complements.