2015 journal article

CONDITIONING OF LEVERAGE SCORES AND COMPUTATION BY QR DECOMPOSITION

SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS, 36(3), 1143–1163.

By: J. Holodnak*, I. Ipsen* & T. Wentworth

author keywords: principal angles; stable rank; condition number; row-scaling; componentwise perturbations
Source: Web Of Science
Added: August 6, 2018

The leverage scores of a full-column rank matrix $A$ are the squared row norms of any orthonormal basis for $\mathrm{range}\,(A)$. We show that corresponding leverage scores of two matrices $A$ and $A+\Delta A$ are close in the relative sense if they have large magnitude and if all principal angles between the column spaces of $A$ and $A+\Delta A$ are small. We also show three classes of bounds that are based on perturbation results of QR decompositions. They demonstrate that relative differences between individual leverage scores strongly depend on the particular type of perturbation $\Delta A$. The bounds imply that the relative accuracy of an individual leverage score depends on its magnitude and the two-norm condition of $A$ if $\Delta A$ is a general perturbation; the two-norm condition number of $A$ if $\Delta A$ is a perturbation with the same normwise row-scaling as $A$; (to first order) neither condition number nor leverage score magnitude if $\Delta A$ is a componentwise row-scaled perturbation. Numerical experiments confirm the qualitative and quantitative accuracy of our bounds.