A three-step method for obtaining asymptotically maximum likelihood AR
MA spectral estimates using only linear transformations of the data is
presented. The computational efficiency of the algorithm is comparabl
e to that of Yule-Walker algorithms, but the three-step method gives a
symptotically statistically efficient estimates. The implementation of
the algorithm is discussed in detail, and numerical examples are pres
ented to illustrate its performance.