Four methods for estimating a dynamic factor model, the direct autoregressive factor score (DAFS) model, are evaluated and compared. The first method estimates the DAFS model using a Kalman filter algorithm based on its state space model representation. The second one employs the maximum likelihood estimation method based on the construction of a block-Toeplitz covariance matrix in the structural equation modeling framework. The third method is built in the Bayesian framework and implemented using Gibbs sampling. The fourth is the least squares method, which also employs the block-Toeplitz matrix. All 4 methods are implemented in currently available software. The simulation study shows that all 4 methods reach appropriate parameter estimates with comparable precision. Differences among the 4 estimation methods and related software are discussed. (Contains 1 figure, 5 tables and 7 footnotes.)