Premium
Computing the conditioning of the components of a linear least‐squares solution
Author(s) -
Baboulin Marc,
Dongarra Jack,
Gratton Serge,
Langou Julien
Publication year - 2009
Publication title -
numerical linear algebra with applications
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.02
H-Index - 53
eISSN - 1099-1506
pISSN - 1070-5325
DOI - 10.1002/nla.627
Subject(s) - overdetermined system , least squares function approximation , mathematics , matrix (chemical analysis) , linear least squares , rank (graph theory) , component (thermodynamics) , variance (accounting) , covariance matrix , algebra over a field , combinatorics , algorithm , linear model , pure mathematics , statistics , materials science , physics , estimator , composite material , thermodynamics , accounting , business
In this paper, we address the accuracy of the results for the overdetermined full rank linear least‐squares problem. We recall theoretical results obtained in ( SIAM J. Matrix Anal. Appl. 2007; 29 (2):413–433) on conditioning of the least‐squares solution and the components of the solution when the matrix perturbations are measured in Frobenius or spectral norms. Then we define computable estimates for these condition numbers and we interpret them in terms of statistical quantities when the regression matrix and the right‐hand side are perturbed. In particular, we show that in the classical linear statistical model, the ratio of the variance of one component of the solution by the variance of the right‐hand side is exactly the condition number of this solution component when only perturbations on the right‐hand side are considered. We explain how to compute the variance–covariance matrix and the least‐squares conditioning using the libraries LAPACK ( LAPACK Users ' Guide (3rd edn). SIAM: Philadelphia, 1999) and ScaLAPACK ( ScaLAPACK Users ' Guide . SIAM: Philadelphia, 1997) and we give the corresponding computational cost. Finally we present a small historical numerical example that was used by Laplace ( Théorie Analytique des Probabilités . Mme Ve Courcier, 1820; 497–530) for computing the mass of Jupiter and a physical application if the area of space geodesy. Copyright © 2008 John Wiley & Sons, Ltd.