We study the estimation variance performance of the matrix pair (MP) method for estimating the impulse responses of multiple FIR channels driven by an unknown input sequence. A first-order perturbation analysis of the large-data-size performance of the MP method is presented and an explicit expression for the estimation variance is derived. Both the theoretical and simulation results are used to investigate the statistical performance of the MP method and a number of new insights are revealed