Premium
Using the simultaneous generalized Schur decomposition as a Candecomp/Parafac algorithm for ill‐conditioned data
Author(s) -
Stegeman Alwin
Publication year - 2009
Publication title -
journal of chemometrics
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.47
H-Index - 92
eISSN - 1099-128X
pISSN - 0886-9383
DOI - 10.1002/cem.1232
Subject(s) - degeneracy (biology) , algorithm , residual , degenerate energy levels , product (mathematics) , mathematics , decomposition , series (stratigraphy) , least squares function approximation , schur decomposition , qr decomposition , computer science , schur complement , geometry , physics , statistics , chemistry , quantum mechanics , eigenvalues and eigenvectors , paleontology , bioinformatics , organic chemistry , estimator , biology
The Candecomp/Parafac (CP)method decomposes a three‐way array into a prespecified number R of outer product arrays, by minimizing the sum‐of‐squares of the residual array. The practical use of CP is sometimes complicated by the occurrence of so‐called ‘degenerate’ sequences of solutions, in which several outer product arrays become highly correlated in all three modes and some elements of the outer product arrays become very large in magnitude. It is known that for I× J× 2 arrays, fitting a simultaneous generalized Schur decomposition (SGSD) avoids the problems of ‘degeneracy’ due to the non‐existence of an optimal CP solution. In this paper, we consider the application of the SGSD method also for other array formats, when the array has a best fitting CP decomposition with ill‐conditioned component matrices, in particular such that it resembles the pattern of a ‘degeneracy’. For these cases, we compare the performance of two SGSD algorithms and the alternating least squares (ALS) CP algorithm in a series of numerical experiments. Copyright © 2009 John Wiley & Sons, Ltd.