# A Decomposition Model for HPLC-DAD Data Set and Its Solution by Particle Swarm Optimization

Authors: Lizhi Cui; Zhihao Ling; Josiah Poon; Simon K. Poon; Junbin Gao; Paul Kwan
Journal: Applied Computational Intelligence and Soft Computing
 (2014)

---

## Abstract

This paper proposes a separation method, based on the model of Generalized Reference Curve Measurement and the algorithm of Particle Swarm Optimization (GRCM-PSO), for the High Performance Liquid Chromatography with Diode Array Detection (HPLC-DAD) data set. Firstly, initial parameters are generated to construct reference curves for the chromatogram peaks of the compounds based on its physical principle. Then, a General Reference Curve Measurement (GRCM) model is designed to transform these parameters to scalar values, which indicate the fitness for all parameters. Thirdly, rough solutions are found by searching individual target for every parameter, and reinitialization only around these rough solutions is executed. Then, the Particle Swarm Optimization (PSO) algorithm is adopted to obtain the optimal parameters by minimizing the fitness of these new parameters given by the GRCM model. Finally, spectra for the compounds are estimated based on the optimal parameters and the HPLC-DAD data set. Through simulations and experiments, following conclusions are drawn: (1) the GRCM-PSO method can separate the chromatogram peaks and spectra from the HPLC-DAD data set without knowing the number of the compounds in advance even when severe overlap and white noise exist; (2) the GRCM-PSO method is able to handle the real HPLC-DAD data set.


---

## Body

## 1. Introduction

After more than 100 years’ development, the technology of chromatography has become the collective term for a set of laboratory technique for quality control of various mixtures such as herbal medicine, grape wine, agriculture, and petroleum. With the development of the chromatographic instrument, the High Performance Liquid Chromatography with Diode Array Detector (HPLC-DAD) technology is used in many researches to generate a data set containing the chromatogram peaks and spectra for all compounds. Figure1 shows the principle of the HPLC-DAD data set. The sample is injected at the sample injection. The high pressure pump drives the solvent to carry the sample to go through the column with absorbent. Different compounds will receive different resistance when they go through the column. Given an ultraviolet detector at the bottom of the column, a chromatogram peak represented by si will be observed when one compound comes out from the column. The position and area of the peak can tell the name and the amount of the compound. If the detector is a DAD, which has more than one thousand channels to detect multiwavelength simultaneously, the spectrum for the same compound represented by ai will also be recorded as well. Di=ai×siT represents ith compound and X represents the mixture. The relationship of the variables in Figure 1 can be shown as
(1)Xw×t=∑i=1nDi=∑i=1nai×siT=a1,a2,…,an×s1Ts2T⋮snT=A×S,
where n indicates the number of the compounds.Figure 1
The principle of the HPLC-DAD data set.For the data setX in (1), there are already several methods to separate it, but with insufficiencies. The algorithm of evolving factor analysis (EFA) [1, 2] and its improvements such as evolutionary factor analysis (EVOLU) [3], fixed-size moving window evolving factor analysis (FSMWEFA) [4], heuristic evolving latent projections (HELP) [5], and orthogonal projection resolution (OPR) [6] are used for peak purity, but without full quantitative information. The method of Multivariate Curve Resolution with Alternating Least Square (MCR-ALS) [7, 8] can recover the pure species spectra and elution profiles. However, the MCR-ALS method will be unavailable when the compounds become complex (see the simulations). And the performance of the MCR-ALS method depends on two important parameters: (1) a threshold for deciding the number of the compounds; (2) the noise level of the data set for estimating initial spectra. Usually, it is not easy to decide these two parameters when noise exists (see Appendix A for explanation). The immune algorithm (IA) [9, 10] can extract the compounds from noise. But, the standard chromatogram peaks for compounds are needed from experiments in advance. The method of independent component analysis (ICA) [11] can separate the HPLC-DAD data set without knowing the number of the compounds in advance. But the cluster methods are still needed to select compounds from the obtained independent components. Our previous works proposed a model named independent components analysis constrained by reference curve (ICARC) and its solution by multiarea genetic algorithm (mGA) [12] and by multitarget Particle Swarm Optimization (mPSO) [13], respectively, which can extract the chromatogram peaks from the HPLC-DAD data set directly. However, through further analysis, we find that it is not necessary for the chromatogram peaks (source signals) to be independent from each other. So a method based on the model of Generalized Reference Curve Measurement (GRCM) and the algorithm of Particle Swarm Optimization (PSO) is proposed in this paper.The remainder of this paper is arranged as follows: Section2 introduces the principle of the GRCM-PSO method; Section 3 gives the simulations and experiments; finally, Section 4 draws the conclusions and future works.


## 2. Mathematical Methods

It is difficult to extract theai and si in (1) only based on the data set X without any other knowledge. Fortunately, the fact that the shape of a chromatogram peak looks like a Gaussian curve [14] can help. Based on this “a priori” knowledge, the GRCM-PSO method is proposed as shown in Figure 2. Firstly, a reference curve r(θ) with parameter θ is constructed based on the general shape of the chromatogram peak, according to which the initial population θi, i=1,2,…,n are generated. Then, the GRCM model calculates the errors εi, i=1,2,…,n for the parameters. Following, a search category is used to obtain the rough solutions θir, i=1,2,…,m (m≪n). In the dashed rectangular box in Figure 2, a step called reinitialization generates t parameters randomly around one rough solution, for example, θ1ir, i=1,2,…,t in Figure 2 for the first rough solution. The GRCM model calculates the errors ε1ir for these θ1ir. Based on these errors ε1ir, the PSO algorithm is adopted to obtain the optimal parameter θ1* around θ1r. Similarly, other optimal parameters θi*, i=2,3,…,m can be found. Finally, the approximated chromatogram peaks can be constructed by the reference curve, and the spectra can be obtained by an estimator.Figure 2
Principle of the GRCM-PSO method.The structure of the parameter used in this paper is the same as that in literature [13], which is shown by (2) and (3). Equation (2) is the Gaussian curve which will be used in the simulations to demonstrate the performance of the GRCM-PSO method; (3) is a 5-parameter curve which will be used in the experiments to show the practicability of our method:
(2)r1θ=r1μ,σ=exp⁡-x-μ22*σ2,μ∈1,colX,σ∈1,μ2-μ16,(3)r2θ=r2(μ,σL,σR,hL,hR)=unitexp⁡-x-μ22*σL2+hL1-hL,x≤μ,unitexp⁡-x-μ22*σR2+hR1-hR,x>μ,μ∈1,colX,σL,σR∈1,μ2-μ13,
where col(X) is the column number of X. r2(θ) is the combination of two Gaussian curves at the peak position with σL and σR for each side’s width and hL and hR for each side’s deviation from zero. The ranges of the σ, σL, σR are limited in order to guarantee that every peak has an integral shape. hL∈[0,0.0001], hR∈[0,0.03] in the experiments due to the profile of the data set. unit{·} is the function to limit the amplitude at 1.In order to obtain initial parameters with small errors, four times of initialization with the same population ofn=2000 are implemented to generate 8000 parameters totally and only the top 2000 parameters according to their errors are chosen as the initialized parameters.


### 2.1. The Model of GRCM

The function of the GRCM model is to assess the parameters by calculating their errors, which indicate the distance between the reference curves constructed by these parameters and the chromatogram peaks existing inX. As shown in Figure 3, the GRCM model is composed of five elements: reference curve r(θ), data set X, Reference Curve Measurement (RCM) model, predicted curve (PC) y, and measurement operator (MO) ·MO.Figure 3
The GRCM model.The RCM model is designed by introducing a vectorm so that
(4)y=mT×X=mT×x1,x2,…,xt≈s⟶rθ.Equation (4) means that let y approximate s and look like r(θ). Then, we have the objective function as
(5)min⁡yT-rTθ22.Solving (5), we obtain the RCM model as
(6)yT=rTθ×1t×X~*T×X~*,
where X~* is a matrix generated from X, which will be introduced in Appendix B as well as the deducing process from (5) to (6).The MO·MO is designed as
(7)ε=yT-rTθ22rTθ22.


### 2.2. Search Category and Reinitialization

After the initialization, every parameterθi will search within a small hypersphere to find one parameter with the smallest error as its target Ti. It is possible for θi to find θj as its target, that is, Ti=θj, and for θj to find θk as its target, that is, Tj=θk. In order to accelerate the searching speed, we directly set Ti=Tj=θk. Finally, only limited parameters have been chosen as targets for others, which are the rough solutions θir, i=1,2,…,m (m≪n). It is assumed that all the real solutions are around the rough solutions because the intensively and randomly distribution of the initializing parameters. So a step named reinitialization only around rough solutions will reduce the search area significantly. The areas for reinitialization are hyperspheres, whose radii are half of the smallest distance between the centre rough solution and other solutions in order to cover all the possible spaces. An example of such a hypersphere is illustrated in Figure 4. There are five rough solutions θir, (i=1,2,…,5) in the two-dimensional space. The distance between θ3r and other rough solutions is d13<d35<d34<d23. So the hypersphere for θ3r is shown by the circle in Figure 4, where R3=0.5×d13. The population in every hypersphere is set to 10.Figure 4
The hypersphere for the reinitialization.


### 2.3. Algorithm of PSO

PSO is swarm intelligence that emulates social interaction and individual cognition of bird flocks foraging [15, 16]. Equation (9) gives the algorithm of PSO:
(8)vi+1=ω×Δvi+c1×r1×pi-θi+c2×r2×pg-θi,θi+=θi+vi+1,
where θi and vi represent the position and velocity of the ith particle, respectively; ω is the inertia weight; c1 and c2 are acceleration constants; r1 and r2 are two random numbers in [0,1]; pi is the personal best position for θi; and pg is the global best position. Please see relative literatures for the values of the parameters in (8).In this paper, all the parameters are divided in several different groups within certain hyperspheres. And every group updates these particles according to (8), respectively, until the value of the best particle in every group does not change for 500 steps, or the maximum step is reached.


### 2.4. Other Processes

During the process fromθi, i=1,2,…,n to θi*, i=2,3,…,m in Figure 2, the random initialization of θi may cause inaccuracy in the results. So this process is executed multiple times to eliminate the influence of the random initialization. Through observation, ten times was chosen. Ten executions will generate 10 candidate solutions. There can be difference among the value even the number of the optimal parameters among these candidates. Firstly, select one candidate with the maximum number of parameters as a reference. Then, select one parameter from every candidate to be grouped with one parameter in the reference according to the Euclidean distance and count the number of parameters in every group. Only the groups with the number of parameters lager than 6 are selected as valid groups. Finally, choose one parameter with lowest error from every valid group to form the final result.Finally, the estimator is designed as the following equation to calculate the spectra for all the compounds:(9)A^=X×pinvS^,
where S^=[s1,s2,…,sm]T≈[y1,y2,…,ym]T is the approximation of chromatogram peaks; pinv(·) is the pseudoinverse function. Equation (9) is derived from (1) directly.


## 2.1. The Model of GRCM

The function of the GRCM model is to assess the parameters by calculating their errors, which indicate the distance between the reference curves constructed by these parameters and the chromatogram peaks existing inX. As shown in Figure 3, the GRCM model is composed of five elements: reference curve r(θ), data set X, Reference Curve Measurement (RCM) model, predicted curve (PC) y, and measurement operator (MO) ·MO.Figure 3
The GRCM model.The RCM model is designed by introducing a vectorm so that
(4)y=mT×X=mT×x1,x2,…,xt≈s⟶rθ.Equation (4) means that let y approximate s and look like r(θ). Then, we have the objective function as
(5)min⁡yT-rTθ22.Solving (5), we obtain the RCM model as
(6)yT=rTθ×1t×X~*T×X~*,
where X~* is a matrix generated from X, which will be introduced in Appendix B as well as the deducing process from (5) to (6).The MO·MO is designed as
(7)ε=yT-rTθ22rTθ22.


## 2.2. Search Category and Reinitialization

After the initialization, every parameterθi will search within a small hypersphere to find one parameter with the smallest error as its target Ti. It is possible for θi to find θj as its target, that is, Ti=θj, and for θj to find θk as its target, that is, Tj=θk. In order to accelerate the searching speed, we directly set Ti=Tj=θk. Finally, only limited parameters have been chosen as targets for others, which are the rough solutions θir, i=1,2,…,m (m≪n). It is assumed that all the real solutions are around the rough solutions because the intensively and randomly distribution of the initializing parameters. So a step named reinitialization only around rough solutions will reduce the search area significantly. The areas for reinitialization are hyperspheres, whose radii are half of the smallest distance between the centre rough solution and other solutions in order to cover all the possible spaces. An example of such a hypersphere is illustrated in Figure 4. There are five rough solutions θir, (i=1,2,…,5) in the two-dimensional space. The distance between θ3r and other rough solutions is d13<d35<d34<d23. So the hypersphere for θ3r is shown by the circle in Figure 4, where R3=0.5×d13. The population in every hypersphere is set to 10.Figure 4
The hypersphere for the reinitialization.


## 2.3. Algorithm of PSO

PSO is swarm intelligence that emulates social interaction and individual cognition of bird flocks foraging [15, 16]. Equation (9) gives the algorithm of PSO:
(8)vi+1=ω×Δvi+c1×r1×pi-θi+c2×r2×pg-θi,θi+=θi+vi+1,
where θi and vi represent the position and velocity of the ith particle, respectively; ω is the inertia weight; c1 and c2 are acceleration constants; r1 and r2 are two random numbers in [0,1]; pi is the personal best position for θi; and pg is the global best position. Please see relative literatures for the values of the parameters in (8).In this paper, all the parameters are divided in several different groups within certain hyperspheres. And every group updates these particles according to (8), respectively, until the value of the best particle in every group does not change for 500 steps, or the maximum step is reached.


## 2.4. Other Processes

During the process fromθi, i=1,2,…,n to θi*, i=2,3,…,m in Figure 2, the random initialization of θi may cause inaccuracy in the results. So this process is executed multiple times to eliminate the influence of the random initialization. Through observation, ten times was chosen. Ten executions will generate 10 candidate solutions. There can be difference among the value even the number of the optimal parameters among these candidates. Firstly, select one candidate with the maximum number of parameters as a reference. Then, select one parameter from every candidate to be grouped with one parameter in the reference according to the Euclidean distance and count the number of parameters in every group. Only the groups with the number of parameters lager than 6 are selected as valid groups. Finally, choose one parameter with lowest error from every valid group to form the final result.Finally, the estimator is designed as the following equation to calculate the spectra for all the compounds:(9)A^=X×pinvS^,
where S^=[s1,s2,…,sm]T≈[y1,y2,…,ym]T is the approximation of chromatogram peaks; pinv(·) is the pseudoinverse function. Equation (9) is derived from (1) directly.


## 3. Simulation, Experiment, and Discussion

In this section, a group of simulations are given to demonstrate the performance of the GRCM-PSO method. Then experiments on a HPLC-DAD data set are implemented to show the practicability of the GRCM-PSO method. Two criteria are used to evaluate the results: (1) whether all the chromatogram peaks can be found; (2) whether the errors between the true/simulated spectra and estimated spectra are small enough.


### 3.1. Simulation and Discussion

The simulation data set is shown in graphs (a) and (b) of Figure5, which contains seven compounds with severe overlap. The seven chromatogram peaks are constructed by (2) with the parameters of [45,6], [53,15], [59,9], [100,30], [141,5], [149,15], and [157,8]. And the seven spectra are constructed randomly as long as they are uncorrelated with each other. The simulation data set is added with different level of whiten noise. The results are listed in Table 1. From the results, we can see the following.(1)
The GRCM-PSO method can separate the simulation data set without knowing the compounds’ number in advance even when severe overlap and white noise exist. This is a big advantage over previous method which needs to know the compounds’ number in advance. The values of theε and the error between the calculated spectra and the simulated spectra are small. The time cost by this method is much less than that by ICARCmPSO [12], which was 13.9 seconds. However, the MCR-ALS method cannot give correct results. The results given by MCR-ALS method are illustrated in graphs (c) and (d) of Figure 5, in which no noise is added to the simulated data set.(2)
The average time cost by the ten implementations is almost the same. This means that the degree of the white noise has no significant influence on the time cost.(3)
The values of theε and the error between the calculated spectra and the simulated spectra become larger with the increase of the noise level. What should be noted is that when noise becomes severe, the “Errors” for small peaks are influenced more significantly than that for big peaks.Table 1
The calculated results for the simulations with different level of white noise. The column of  “Rate” indicates the rate for this parameter emerging among the total ten candidates. The column of “Error” represents the error between simulated spectra and calculated spectra.
Number
SNR
Parameters
Steps/times
Rate
ε
Error
μ
σ
1
No noise
45
6
1048.1/3.17 s
1
3.34
e
-
27
2.63
e
-
25
53
15
1
7.22
e
-
27
2.74
e
-
26
59
9
1
5.26
e
-
26
1.06
e
-
24
100
30
1
4.12
e
-
28
1.63
e
-
27
141
5
1
2.55
e
-
25
2.48
e
-
24
149
15
1
5.85
e
-
28
1.40
e
-
26
157
8
1
1.95
e
-
27
3.14
e
-
24
2
100
45
6
1059.3/3.23 s
1
5.29
e
-
11
2.57
e
-
09
53
15
1
9.97
e
-
11
2.64
e
-
10
59
9
1
6.86
e
-
10
5.59
e
-
09
99.99
30
1
2.00
e
-
11
1.83
e
-
11
141
4.99
1
2.96
e
-
9
2.53
e
-
09
149
15
1
9.02
e
-
13
1.42
e
-
09
157
8
1
4.97
e
-
11
5.73
e
-
08
3
50
44.99
5.99
1063.1/3.20 s
1
5.32
e
-
6
0.01
52.96
15.06
1
9.52
e
-
6
2.4
e
-
5
58.98
9.01
1
5.75
e
-
5
0.09
99.99
29.99
1
1.99
e
-
6
1.50
e
-
6
140.98
5.05
1
2.59
e
-
4
0.0003
149
14.99
1
1.22
e
-
7
0.0001
157
7.99
1
5.31
e
-
6
0.005
4
40
45
5.99
1055.7/3.26 s
1
5.38
e
-
5
0.11
52.79
15.17
1
1.07
e
-
4
5.3
e
-
4
58.95
8.98
0.9
6.65
e
-
4
1.23
100.04
30.46
1
1.99
e
-
5
2.30
e
-
5
140.99
5.03
1
3.1
e
-
3
0.008
149.01
14.99
1
8.55
e
-
7
0.008
156.98
8.01
1
5.90
e
-
5
0.09(a) The profile and simulated chromatogram peaks; (b) simulated spectra; (c) the calculated chromatogram peaks by MCR-ALS method; (d) the calculated spectra by MCR-ALS method.
(a)
(b)
(c)
(d)


### 3.2. Experiment and Discussion

The HPLC-DAD data set of “adataset.mat” is downloaded fromhttp://www.mcrals.info/ for free. This data set is a three-compound mixture with two known pesticides and one unknown interferent [8]. The 5-parameter function r2(θ) shown by (3) is used in the experiments as the RC. The graphical results are illustrated in Figure 6 as well as the results given by the ICARCmPSO method and the MCR-ALS method. The values of the results are listed in Table 2. From the experiments, we can see the following.Table 2
The calculated results for the experiments. The column of “Rate” indicates the rate for this parameter emerging among the total ten candidates. The column of “Error” represents the error between the calculated spectra and the true spectra.
Method
Parameters
Steps/times
Rate
ε
Error
μ
σ
(
σ
L
)
σ
R
h
L
h
R
GRCM-PSO
23.93
7.11
12.70
0
0.03
1361/3.61 s
0.8
1.9
e
-
5
0.27
30.29
8.72
14.07
1.0
e
-
4
0.01
1
1.3
e
-
5
1.59
59.85
7.91
14.43
3
e
-
21
0.03
1
3.3
e
-
5
—
ICARCmPSO
23.07
7.37
12.28
0.9
e
-
4
0.08
1481.8/16.55 s
1
1.5
e
-
5
0.18
29.91
8.39
14.12
1.2
e
-
4
0.01
1
1.4
e
-
5
1.73
59.17
8.26
13.04
0.8
e
-
4
0.10
1
2.1
e
-
5
—
MCR-ALS
—
—
—
—
—
1000/0.39 s
—
—
0.033
—
—
—
—
—
—
—
0.153The calculated results. (a), (b), and (c) are for the GRCM-PSO method. (d), (e), and (f) are for the ICARCmPSO method. (g), (h), and (i) are for the MCR-ALS method. (a), (d), and (g) are the calculated chromatogram peaks. (b), (e), and (h) are the compares between first spectrum and corresponding calculated spectrum. (c), (f), and (i) are the compares between the second ones.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)(1) Comparison between the GRCM-PSO method and the ICARCmPSO [13]: the average time and average steps for the GRCM-PSO method are much less. The ε and the “Errors” are similar for both of these two methods.For the ICARCmPSO method, the parameter of scope for particles to search their local targets should be determined according to the specific application. For the GRCM-PSO method, this parameter is fixed to a small value. So, from the view of operability and speed, the GRCM-PSO method has advantage over the ICARCmPSO method.(2) Comparison between the GRCM-PSO method and the MCR-ALS method: both of them obtain the same number, 3, of the compounds. The speed of the MCR-ALS method is better than that of the GRCM-PSO method. The accuracy of the MCR-ALS method can be better than that of the GRCM-PSO method. But the parameters in the GRCM-PSO method are easier to be controlled.For the MCR-ALS method, the important two parameters are the threshold to select the valid singular values and the noise level for initial estimation of the spectra. If noise exists, it will be difficult to decide the first parameter as explained in AppendixA. The performance of the MCR-ALS method is also very sensitive to the second parameter as shown in graph (i) of Figure 6. A small change of the β, which is explained in Appendix A, will cause big error in the calculated spectra, while all the parameters for the GRCM-PSO method are fixed for all applications.So, from the view of operability and stability, the GRCM-PSO method has advantage over the MCR-ALS method.


## 3.1. Simulation and Discussion

The simulation data set is shown in graphs (a) and (b) of Figure5, which contains seven compounds with severe overlap. The seven chromatogram peaks are constructed by (2) with the parameters of [45,6], [53,15], [59,9], [100,30], [141,5], [149,15], and [157,8]. And the seven spectra are constructed randomly as long as they are uncorrelated with each other. The simulation data set is added with different level of whiten noise. The results are listed in Table 1. From the results, we can see the following.(1)
The GRCM-PSO method can separate the simulation data set without knowing the compounds’ number in advance even when severe overlap and white noise exist. This is a big advantage over previous method which needs to know the compounds’ number in advance. The values of theε and the error between the calculated spectra and the simulated spectra are small. The time cost by this method is much less than that by ICARCmPSO [12], which was 13.9 seconds. However, the MCR-ALS method cannot give correct results. The results given by MCR-ALS method are illustrated in graphs (c) and (d) of Figure 5, in which no noise is added to the simulated data set.(2)
The average time cost by the ten implementations is almost the same. This means that the degree of the white noise has no significant influence on the time cost.(3)
The values of theε and the error between the calculated spectra and the simulated spectra become larger with the increase of the noise level. What should be noted is that when noise becomes severe, the “Errors” for small peaks are influenced more significantly than that for big peaks.Table 1
The calculated results for the simulations with different level of white noise. The column of  “Rate” indicates the rate for this parameter emerging among the total ten candidates. The column of “Error” represents the error between simulated spectra and calculated spectra.
Number
SNR
Parameters
Steps/times
Rate
ε
Error
μ
σ
1
No noise
45
6
1048.1/3.17 s
1
3.34
e
-
27
2.63
e
-
25
53
15
1
7.22
e
-
27
2.74
e
-
26
59
9
1
5.26
e
-
26
1.06
e
-
24
100
30
1
4.12
e
-
28
1.63
e
-
27
141
5
1
2.55
e
-
25
2.48
e
-
24
149
15
1
5.85
e
-
28
1.40
e
-
26
157
8
1
1.95
e
-
27
3.14
e
-
24
2
100
45
6
1059.3/3.23 s
1
5.29
e
-
11
2.57
e
-
09
53
15
1
9.97
e
-
11
2.64
e
-
10
59
9
1
6.86
e
-
10
5.59
e
-
09
99.99
30
1
2.00
e
-
11
1.83
e
-
11
141
4.99
1
2.96
e
-
9
2.53
e
-
09
149
15
1
9.02
e
-
13
1.42
e
-
09
157
8
1
4.97
e
-
11
5.73
e
-
08
3
50
44.99
5.99
1063.1/3.20 s
1
5.32
e
-
6
0.01
52.96
15.06
1
9.52
e
-
6
2.4
e
-
5
58.98
9.01
1
5.75
e
-
5
0.09
99.99
29.99
1
1.99
e
-
6
1.50
e
-
6
140.98
5.05
1
2.59
e
-
4
0.0003
149
14.99
1
1.22
e
-
7
0.0001
157
7.99
1
5.31
e
-
6
0.005
4
40
45
5.99
1055.7/3.26 s
1
5.38
e
-
5
0.11
52.79
15.17
1
1.07
e
-
4
5.3
e
-
4
58.95
8.98
0.9
6.65
e
-
4
1.23
100.04
30.46
1
1.99
e
-
5
2.30
e
-
5
140.99
5.03
1
3.1
e
-
3
0.008
149.01
14.99
1
8.55
e
-
7
0.008
156.98
8.01
1
5.90
e
-
5
0.09(a) The profile and simulated chromatogram peaks; (b) simulated spectra; (c) the calculated chromatogram peaks by MCR-ALS method; (d) the calculated spectra by MCR-ALS method.
(a)
(b)
(c)
(d)


## 3.2. Experiment and Discussion

The HPLC-DAD data set of “adataset.mat” is downloaded fromhttp://www.mcrals.info/ for free. This data set is a three-compound mixture with two known pesticides and one unknown interferent [8]. The 5-parameter function r2(θ) shown by (3) is used in the experiments as the RC. The graphical results are illustrated in Figure 6 as well as the results given by the ICARCmPSO method and the MCR-ALS method. The values of the results are listed in Table 2. From the experiments, we can see the following.Table 2
The calculated results for the experiments. The column of “Rate” indicates the rate for this parameter emerging among the total ten candidates. The column of “Error” represents the error between the calculated spectra and the true spectra.
Method
Parameters
Steps/times
Rate
ε
Error
μ
σ
(
σ
L
)
σ
R
h
L
h
R
GRCM-PSO
23.93
7.11
12.70
0
0.03
1361/3.61 s
0.8
1.9
e
-
5
0.27
30.29
8.72
14.07
1.0
e
-
4
0.01
1
1.3
e
-
5
1.59
59.85
7.91
14.43
3
e
-
21
0.03
1
3.3
e
-
5
—
ICARCmPSO
23.07
7.37
12.28
0.9
e
-
4
0.08
1481.8/16.55 s
1
1.5
e
-
5
0.18
29.91
8.39
14.12
1.2
e
-
4
0.01
1
1.4
e
-
5
1.73
59.17
8.26
13.04
0.8
e
-
4
0.10
1
2.1
e
-
5
—
MCR-ALS
—
—
—
—
—
1000/0.39 s
—
—
0.033
—
—
—
—
—
—
—
0.153The calculated results. (a), (b), and (c) are for the GRCM-PSO method. (d), (e), and (f) are for the ICARCmPSO method. (g), (h), and (i) are for the MCR-ALS method. (a), (d), and (g) are the calculated chromatogram peaks. (b), (e), and (h) are the compares between first spectrum and corresponding calculated spectrum. (c), (f), and (i) are the compares between the second ones.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)(1) Comparison between the GRCM-PSO method and the ICARCmPSO [13]: the average time and average steps for the GRCM-PSO method are much less. The ε and the “Errors” are similar for both of these two methods.For the ICARCmPSO method, the parameter of scope for particles to search their local targets should be determined according to the specific application. For the GRCM-PSO method, this parameter is fixed to a small value. So, from the view of operability and speed, the GRCM-PSO method has advantage over the ICARCmPSO method.(2) Comparison between the GRCM-PSO method and the MCR-ALS method: both of them obtain the same number, 3, of the compounds. The speed of the MCR-ALS method is better than that of the GRCM-PSO method. The accuracy of the MCR-ALS method can be better than that of the GRCM-PSO method. But the parameters in the GRCM-PSO method are easier to be controlled.For the MCR-ALS method, the important two parameters are the threshold to select the valid singular values and the noise level for initial estimation of the spectra. If noise exists, it will be difficult to decide the first parameter as explained in AppendixA. The performance of the MCR-ALS method is also very sensitive to the second parameter as shown in graph (i) of Figure 6. A small change of the β, which is explained in Appendix A, will cause big error in the calculated spectra, while all the parameters for the GRCM-PSO method are fixed for all applications.So, from the view of operability and stability, the GRCM-PSO method has advantage over the MCR-ALS method.


## 4. Conclusions and Future Works

A method named GRCM-PSO was proposed in this paper to separate the chromatogram peaks and spectra for compounds from the HPLC-DAD data set. The GRCM model transformed the separation problem to a multiparameter optimization issue. The PSO algorithm was introduced to calculate the optimal parameters. Groups of simulations with different noise level were implemented. A simulated data set was constructed with severe overlap among seven compounds. The GRCM-PSO method separated the chromatogram peaks and spectra from this simulated data set without knowing the number of the compounds in advance. And the speed was fast. Groups of experiments on a real HPLC-DAD data set were implemented. And comparisons among the results by the GRCM-PSO method, the ICARCmPSO method, and the MCR-ALS method were given. The results showed that the GRCM-PSO method was an effective, efficient, and practical method to separate HPLC-DAD data set even when severe overlap and white noise existed. The speed and practicability of the GRCM-PSO method are better than that of the ICARCmPSO method. The stability and operability of the GRCM-PSO method are better than that of the MCR-ALS method.Currently, the performance of the GRCM-PSO method depends on the selection of the reference curve. So it is only suitable for the separation task with “a priori” knowledge to be known, such as the separation of HPLC-DAD data set. The accuracy of the result by the GRCM-PSO method can be improved by further research on more accurate reference curves.



---

*Source: 276741-2014-11-25.xml*