Factor analysis: how to extract the low rank matrix and the diagonal matrix

Question

The factor analysis model is $x = Rz + \\epsilon$ where $z \\sim N(0, I)$, $\\epsilon \\sim N(0, \\Psi)$, and $\\Psi$ is a diagonal matrix. Let $\\Sigma$ be the covariance m