Estimating the chemical rank of three-way data arrays by a simple linear transform incorporating Monte Carlo simulation

Estimating an appropriate chemical rank of a three-way data array is very important to second-order calibration. In this paper, a simple linear transform incorporating Monte Carlo simulation approach (LTMC) to estimate the chemical rank of a three-way data array was suggested. The new method determines the chemical rank through performing a simple linear transform procedure on the original cube matrix to produce two subspaces by singular value decomposition. One of two subspaces is derived from the original three-way data array itself and the other is derived from a new three-way data array produced by the linear transformation of the original one. Projection technique incorporating the Monte Carlo approach acts as distinguishing criterion to choose the appropriate component number of the system. Simulated three-way trilinear data arrays with different noise types (homoscedastic and heteroscedastic), various noise level as well as high collinearity are used to illustrate the feasibility of the new method. The results have shown that the new method could yield accurate results with different conditions appended. The feasibility of the new method is also confirmed by two real arrays, HPLC-DAD data and excitation-emission fluorescent data. All the results are compared with the other three factor-determining methods: factor indicator function (IND), core consistency diagnostic (CORCONDIA) and two-mode subspace comparison (TMSC) approach. It shows that the newly proposed algorithm can objectively and quickly determine the chemical rank to fit the trilinear model.