This paper proposes maximum (quasi)likelihood estimation for high dimensional factor models with regime switching in the loadings. The model parameters are estimated jointly by EM algorithm, which in the current context only requires iteratively calculating regime probabilities and principal...