Yanev, P.I. and Kontoghiorghes, Erricos J. (2006) Efficient algorithms for estimating the general linear model. Parallel Computing 32 (2), pp. 195-204. ISSN 0167-8191.
Abstract
Computationally efficient serial and parallel algorithms for estimating the general linear model are proposed. The sequential block-recursive algorithm is an adaptation of a known Givens strategy that has as a main component the Generalized QR decomposition. The proposed algorithm is based on orthogonal transformations and exploits the triangular structure of the Cholesky QRD factor of the variance–covariance matrix. Specifically, it computes the estimator of the general linear model by solving recursively a series of smaller and smaller generalized linear least squares problems. The new algorithm is found to outperform significantly the corresponding LAPACK routine. A parallel version of the new sequential algorithm which utilizes an efficient distribution of the matrices over the processors and has low inter-processor communication is developed. The theoretical computational complexity of the parallel algorithms is derived and analyzed. Experimental results are presented which confirm the theoretical analysis. The parallel strategy is found to be scalable and highly efficient for estimating large-scale general linear estimation problems.
Metadata
Item Type: | Article |
---|---|
Keyword(s) / Subject(s): | General linear model, generalized QR decomposition, parallel algorithms |
School: | Birkbeck Faculties and Schools > Faculty of Science > School of Computing and Mathematical Sciences |
Depositing User: | Administrator |
Date Deposited: | 18 Aug 2011 08:46 |
Last Modified: | 09 Aug 2023 12:30 |
URI: | https://eprints.bbk.ac.uk/id/eprint/4028 |
Statistics
Additional statistics are available via IRStats2.