We introduce a frequency domain version of the EM algorithm for general dynamic factor models. We consider both AR and ARMA processes, for which we develop iterative indirect inference procedures analogous to the algorithms in Hannan (1969). Although our proposed procedure allows researchers to estimate such models by maximum likelihood with many series even without good initial values, we recommend switching to a gradient method that uses the EM principle to swiftly compute frequency domain analytical scores near the optimum. We successfully employ our algorithm to construct an index that captures the common movements of US sectoral employment growth rates.