Given a random process spectral factor

, the phase approximation algorithm, initiated by Jonckheere and Helton [1], constructs a reduced-order spectral factor

such that

is small in the Hankel-norm sense. In this note, we derive the

error bound

, where the σ\´s are the canonical correlation coefficients. A similar result holds in the multivariable case.