Parallel algorithms for downdating the least squares estimator of the regression model
Yanev, P.I. and Kontoghiorghes, Erricos J. (2008) Parallel algorithms for downdating the least squares estimator of the regression model. Parallel Computing 34 (6-8), pp. 451-468. ISSN 0167-8191.
Computationally efficient parallel algorithms for downdating the least squares estimator of the ordinary linear regression are proposed. The algorithms, which are based on the QR decomposition, are block versions of sequential Givens strategies and efficiently exploit the triangular structure of the data matrices. The first strategy utilizes only part of the orthogonal matrix which is derived from the QR decomposition of the initial data matrix. The rest of the orthogonal matrix is not updated or explicitly computed. A modification of the parallel algorithm, which explicitly computes the whole orthogonal matrix in the downdated QR decomposition, is also considered. An efficient distribution of the matrices over the processors is proposed. Furthermore, the new algorithms do not require any inter-processor communication. The theoretical complexities are derived and experimental results are presented and analyzed. The parallel strategies are scalable and highly efficient for large scale downdating least squares problems. A new parallel block-hyperbolic downdating strategy is developed. The algorithm is rich in BLAS-3 computations, involves negligible duplicated computations and requires insignificant inter-processor communication. It is found to outperform the previous downdating strategies and to be highly efficient for large scale problems. The experimental results confirm the derived theoretical complexities.
|Keyword(s) / Subject(s):||QR decomposition, least squares, downdating, parallel algorithm|
|School:||Birkbeck Schools and Departments > School of Business, Economics & Informatics > Computer Science and Information Systems|
|Date Deposited:||28 Jul 2011 14:25|
|Last Modified:||17 Apr 2013 12:21|
Additional statistics are available via IRStats2.