We present a statistical analysis of subspace-based methods for the retrieval of sinusoids. The formalism described previously encompasses the modern nonlinear methods of MUSIC and ESPRIT and yields methods of even higher resolution. These rely on eigendecomposition of state-covariances of linear systems, as opposed to eigendecomposition of Toeplitz matrices (originating from antenna arrays or tapped delay-lines and treated). We focus on the variability of estimates when the theory is applied to sampled covariances obtained from finite observation records, and in particular, we provide an expression for the variance of the angle operator between estimated and exact signal subspaces.