On the fitting of multivariate autoregressions, and approximate factorization of a spectral density matrix