We propose to couple the R2 performance measure and Particle Swarm Optimization in order to handle multi/manyobjective problems. Our proposal shows that through a welldesigned interaction process we could maintain the metaheuristic almost inalterable and through the R2 performance measure we did not use neither an external archive nor Pareto dominance to guide the search. The proposed approach is validated using several test problems and performance measures commonly adopted in the specialized literature. Results indicate that the proposed algorithm produces results that are competitive with respect to those obtained by four wellknown MOEAs. Additionally, we validate our proposal in manyobjective optimization problems. In these problems, our approach showed its main strength, since it could outperform another wellknown indicatorbased MOEA.
1. Introduction
Evolutionary Algorithms (EAs) encompass a set of bioinspired techniques that make a multidimensional search and have been found to be effective in locating solutions close to the global optimum even in highly rugged search spaces. EAs are suitable alternatives for solving problems with two or more objectives (the socalled multiobjective optimization problems or MOPs for short), since they are able to simultaneously explore different regions of the search space and obtain several points from the tradeoff surface in a single run. Since the mid1980s, the field of evolutionary multiobjective optimization (EMO) has grown and a wide variety of multiobjective EAs (MOEAs) to solve real applications have been proposed so far [1–3].
Moreover, the use of Pareto dominance (PD) to solve multiobjective optimization problems (MOPs) has been successfully used for several years. However, when the number of objectives increases, the proportion of nondominated solutions also increases but in an exponential way [4–6]. Therefore, very quickly, it becomes impossible to distinguish individuals for selection purposes. Such a behavior dilutes the selection pressure, since the choice of solutions is performed practically at a random way.
For this reason, the evolutionary multiobjective optimization (EMO) community has developed several approaches to overcome this shortcoming in the fitness assignment process [7, 8]. Several of these approaches drive the search using a quality assessment indicator. This idea has become more popular in the last few years, mainly because of the growing interest in tackling multiobjective problems having 4 or more objectives (commonly called “manyobjective optimization problems” or MaMOP for short), for which indicatorbased MOEAs seem to be particularly suitable [9]. The idea of using an indicatorbased selection is to identify the solutions that contribute the most to the improvement of the performance indicator adopted in the selection mechanism [10, 11].
The IndicatorBased Evolutionary Algorithm (IBEA) proposed by Zitzler and Künzli [12] is the most general version of an algorithm of this sort. Instead of using the problem at hand as a fitness function, these methods minimize or maximize (whichever the case) a performance indicator. Originally, IBEA was proposed to be used with two different performance measures: hypervolume [13] and the ϵindicator. Among other results, Zitzler and Künzli found that no additional diversity preservation mechanism was required when ranking the solutions with an indicator. Similarly, Beume et al. [14] proposed S Metric Selection Evolutionary Multiobjective Algorithm (SMSEMOA). However, in this case, Beume et al. replaced the Crowding Distance with the hypervolume indicator. Ishibuchi et al. [15] presented a novel approach that iteratively optimizes separately each objective and searches for the solution that contributes more to the hypervolume indicator. This approach was designed to search for a small number of nondominated solutions along the entire Pareto front. Igel et al. [16] proposed the Multiobjective Covariance Matrix Adaptation Evolution Strategy (MOCMAES). This algorithm uses a set of monoobjective optimizers as its population. Each optimizer generates new solutions that can be accepted back into the population according to their ranking with respect to PD or their contribution to the hypervolume.
Although the hypervolume’s nice theoretical property has positioned it as the most popular choice for implementing indicatorbased MOEAs [14], it is wellknown that its computational cost considerably increases as we raise the number of objectives. To overcome this drawback, some researchers [9] have opted for approximating the hypervolume. However, this sort of scheme can decrease the accuracy of the selection mechanism. Rodríguez Villalobos and Coello Coello [17] recently proposed Δpdifferential evolution in which the authors adopted Δp performance indicator [18] as an alternative to the hypervolume. The fitness assignment for each solution is performed through the contribution to Δp. Δp indicator requires a reference set to be calculated, and the authors used the nadir point and the ideal vector to create such a reference set. The authors reported that this approach could obtain competitive results with respect to other MOEAs (including SMSEMOA), having as its main advantage its very low computational cost, even when dealing with manyobjective problems.
Moreover, the PSO has been used to solve a lot of problems [19–21]. However, in the literature the works with respect to the use of an indicator to guide the search of a PSO are very limited. In [22], Padhye uses the contribution to the hypervolume to select the gbest and the pbest of the PSO. However, in this work the hypervolume is not used to select the particles that advance to the next iteration. This algorithm is used for the topology optimization of a compliant mechanism. In [23], the authors proposed hybridization of MOPSO with a local search operator. The proposed MOPSO uses an indicator to truncate an external archive of solutions. In this work, two approaches were proposed one that uses the ϵindicator and another that adopts the hypervolume performance measure. Both approaches reached similar results. However, these approaches were compared only with problems with 2 and 3 objectives.
Other suitable performance indicators are the R2 indicator [24]. In recent works, it has been reported that desirable properties (i.e., it is weakly monotonic, it produces welldistributed solutions, and it can be computed in a fast manner) make the R2 indicator a viable candidate to be incorporated into an indicatorbased MOEA [10, 25–29]. In such works, the R2 behavior has been compared with respect to that of the hypervolume and concludes that they both have a similar behavior, but R2 has a considerably lower computational cost.
In [27], the authors proposed an approach to fast ranking the population (of a genetic algorithm and a differential evolution) using the R2 indicator. Although their approach is able to work with manyobjective problems, it uses an external archive. Therefore, the metaheuristic at hand has to be highly modified in order to work with the approach. In this paper, we propose a R2based multiobjective approach that maintains the nature of PSO, while empowering it to handle manyobjective problems.
In this work, we propose to use the R2 indicator to guide the search of MOPSO. The new approach is then compared with respect to some stateoftheart MOEAs taken from the specialized literature. Furthermore, we present a scalability study to analyze the behavior of our proposed approach as the number of objectives of the problem increases. The remainder of this work is organized as follows. Details of the R2 indicator are given in Section 2. Section 3 presents our proposed approach. A comparative study with respect to other algorithms is presented in Section 4. Finally, Conclusions and Future Work are given in Section 5.
The family of R indicators [24] is based on utility functions which map a vector y→∈Rk to a scalar value u∈R in order to measure the quality of two approximations of the Pareto front.
Definition 1.
For a set U of utility functions, probability distribution p on U, and a reference set R, the R2 indicator of a solution set A is defined as(1)R2R,A,U,p=∫u∈Umaxr∈Rurpudu∫u∈Umaxa∈Auapudu.
Definition 2.
For a discrete and finite set U and uniform distribution p over U, the R2 indicator can be defined as [30](2)R2R,A,U=1U∑u∈Umaxr∈Rurmaxa∈Aua.
Since the first summand (maxr∈R{u(r)}) is constant if we assume a constant R, the first summand can be deleted in order to have unary indicator as a result (called R for simplicity) [25].
Definition 3.
For a constant reference set, the R2 indicator can be defined as unary indicator as follows:(3)R2A,U=1U∑u∈Uua.
We selected the Tchebycheff function to be used as the utility function of our approach. This function works well when optimizing different types of Pareto fronts. Also, this aggregation function is not smooth for continuous multiobjective problems. However, our algorithm does not need to compute the derivative of the aggregation function. The Tchebycheff function can be defined as u(z)=uλ(z→)=maxj∈{1,…,k}λjzj∗zj, where λ=(λ1,…,λk)∈Λ is a weight vector and z∗ is utopian point (an objective vector that is not dominated by any feasible search point).
Definition 4.
The R2 indicator of a solution set A for a given set of weight vectors Λ and utopian point z∗ is defined as(4)R2A,Λ,z∗=1Λ∑λ∈Λmina∈Amaxj∈1,…,kλjzj∗aj.
Definition 5.
Finally, we say that the contribution of one solution a∈A to the R2 indicator can be defined as(5)CR2a,A,Λ,z∗=R2A,Λ,z∗R2A∖a,Λ,z∗.
3. Proposed Approach3.1. The PSO Algorithm
PSO has been successfully used for both continuous nonlinear and discrete binary single objective optimization [31]. The pseudocode of PSO is shown in Algorithm 1.
<bold>Algorithm 1: </bold>PSO algorithm.
(1) gbest→←x0→
(2) fori=0 to nparticlesdo
(3) pbesti→←xi→←initialize_randomly()
(4) fitnessi←fxi→
(5) iffitnessi<f(gbest→)then
(6) gbest→←xi→
(7) end if
(8) end for
(9) repeat
(10) fori=0 to nparticlesdo
(11) ford=0 to ndimensionsdo
(12) vid←W∗vid+C1∗U(0,1)∗(pbestidxid)
+ C2∗U(0,1)∗gbestxid
(13) xid←xid+vid
(14) end for
(15) fitnessi←f(xi→)
(16) iffitnessi<f(pbesti→)then
(17) pbesti→←xi→
(18) end if
(19) iffitnessi<f(gbesti→)then
(20) gbesti→←xi→
(21) end if
(22) end for
(23) until Termination criterion
3.2. PSO Based on the <inlineformula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M74"><mml:mi>R</mml:mi><mml:mn>2</mml:mn></mml:math></inlineformula> Indicator (<inlineformula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M75"><mml:mi>R</mml:mi><mml:mn>2</mml:mn></mml:math></inlineformula>MOPSO)
PSO has been particularly suitable for multiobjective optimization mainly because of the high speed of convergence that the algorithm presents. Based on such behavior, one would expect a multiobjective PSO (MOPSO) based on indicators to be very efficient computationally speaking. However, it is necessary to perform two main modifications to the original PSO:
To modify the algorithm to handle multiple objectives and produce a set of nondominated solutions in a single run.
To modify the algorithm to obtain good distribution of solutions.
A natural modification to a PSO algorithm aimed to handle multiple objectives which is to replace the comparison operator in order to determine whether a solution a is better than a solution b. Most of the existing approaches use the Pareto ranking scheme to extend the PSO to handle multiobjective optimization problems. However, with a Pareto ranking scheme a set of nondominated solutions will be produced (by definition, all nondominated solutions are equally good). Having several nondominated solutions implies the inclusion into the algorithm of both additional criteria to decide whether a new nondominated solution is pbest or gbest and a strategy to select the guide particles (pbest and gbest). Therefore, in our proposal the contribution to the R2 indicator (see (5)) is used to replace the comparison operator. Additionally, an external file A of nondominated solutions is used with all the evaluated solutions found in the optimization process. If the number of nondominated solutions is greater than the limit of the file, we select the solutions that represent the best R2 contributions. The utopian point (z∗) is formed with respect to the best obtained values for each objective (of this external file). This utopian point is updated when the external file changes.
The MOPSO implemented in this paper is called the R2MOPSO. At the start of the optimization cycle, all the N particle positions are initialized randomly and their velocities are set to zero. At the onset pbest for each particle is assigned as the particle itself. A different gbest is obtained for each particle. As gbest, pbest was selected with the best contribution to the R2 indicator. Next, the velocity and position of the particle are updated. Afterwards, the new position of the particle is evaluated, and the contribution to R2 is calculated for the union of all the pbest positions with the new position of the particle, so that if the new position has a better R2 contribution than its pbest then the pbest of the particle will be replaced with the new position. Moreover, if the new position has a similar R2 contribution than its pbest then either one is chosen randomly to be the new pbest.
An important feature of R2MOPSO is that the Pareto dominance is completely removed from the evolutionary process. The Pareto dominance is only used in the external file. However, this external file is not used in the optimization process of the PSO. Thus, we can say that R2MOPSO maintains the essence of the original PSO. We show the pseudocode of our R2MOPSO in Algorithm 2.
In order to compute the R2 indicator, it is necessary to have a set of uniformly distributed weight vectors in order to obtain good distribution of solutions. For biobjective problems and for threeobjective problems, we compute each weight vector as follows:(6)ϵ,1H,2H,…,HH, where H controls the number of weight vectors and ϵ is a value close to zero (ϵ=104 is suggested), in order to prevent cancellation in the calculations. The total number of vectors is represented by Cm1H+m1. In this work, for biobjective problems, we decided to use H=99, and for threeobjective problems we adopted H=23.
However, for manyobjective problems the weight vectors were randomly initialized in such a way that the sum of each weight vector is equal to one (the random weights’ vector was generated as in MOMHLib++ [32]).
4. Performance Assessment
The proposed approach was evaluated using 15 test functions. Five functions were taken from the ZitzlerDebThiele (ZDT) test suite [33], six test problems were taken from [34], and the remaining were taken from the DebThieleLaumannsZitzler (DTLZ) test suite [35]. The main features of these test problems are shown in Table 1.
Test problems adopted.
Test problem
# of variables
# of objectives
ZDT1–3, UF1–3
30
2
ZDT4, 6
10
2
DTLZ1
7
3
DTLZ2–4
12
3
UF8–10
30
3
In order to assess the performance of the proposed approach, we decided to include the hypervolume performance measure.
(1) Hypervolume (HV). The HV computes the area covered for all the solutions in the approximated Pareto front Q with the help of a reference point W. Equation (7) shows the mathematical definition of HV:(7)HV=volume⋃i=1Qvi,where, for each solution i∈Q, a hypercube vi is constructed using the reference point W. Therefore, HV is the union of the volume of all the hypercubes.
4.1. Experiment 1 (Two and Three Objectives)
In order to make a comparative study, we chose the four following approaches: (1) NSGAII (this is, by far, the most popular Paretobased MOEA), (2) MOEA/D (a more recent MOEA, based on decomposition, which has been found to be more effective than NSGAII in a number of problems), (3) SMSEMOA (this is perhaps the most popular indicatorbased MOEA in use today), and (4) MOMBIII (an indicatorbased MOEA based on the R2 indicator). All these MOEAs adopted a population size of 100 (except for MOEA/D in problems with 3 objectives where the population size was of 105), and the remainder parameters for each algorithm were the ones suggested by their authors. Table 2 summarizes the parameters adopted in our comparative study.
Adopted parameters for each MOEA.
NSGAII
MOEA/D
SMSEMOA
MOMBIII
R2MOPSO
pc=1.0
pc=1.0
pc=1.0
pc=1.0
c1=1.4
pm=1n
pm=1n
pm=1n
pm=1n
c2=1.4
nc=15
nc=15
nc=15
nc=30
w=0.4
nm=15
nm=20
nm=20
nm=20
size_file=100
T=20
Since we would like to investigate about behavior of the algorithms, we decided to measure the online convergence. Therefore, 30 independent executions were performed, and we measured the HV every 100 evaluations during 20,000 function evaluations in biobjective problems and 30,000 function evaluations in problems with three objectives.
Figure 1 shows the results of the MOEAs in all the adopted problems. In this figures, xaxis shows the number of evaluations and yaxis shows the average performance according to the HV. A better performance in most of the problems with two objectives (ZDT1, ZDT2, ZDT3, and ZDT6 and UF1 and UF2) of R2MOPSO with respect to the other algorithms can be observed. A rapid convergence in this type of problem is shown. However, in ZDT4 the R2MOPSO produced a slower convergence than the other approaches. Moreover, in threeobjective problems the behavior of the R2MOPSO is similar than the rest of the algorithms.
Graphical results of the online convergence for all the test problems.
The application of the hypervolume performance measure to the results obtained by the four approaches is shown in Table 3. From these results, it is easy to see that all the approaches behaved similarly. The NSGAII slightly outperformed others for the ZDT test problems. However, this approach did not behave well when optimizing the threeobjective test problems. It can be observed that the R2MOPSO are competitive in most of the problems with respect to the results obtained by the other algorithms.
Comparison of the results obtained by NSGAII, MOEA/D, SMSEMOA, MOMBIII and R2MOPSO with respect to the hypervolume with 20,000 evaluations for biobjective problems and 30,000 for three objective problems.
Problem
NSGAII
MOEA/D
SMSEMOA
MOMBIII
R2MOPSO
ZDT1
120.644 ± 0.00
120.375 ± 0.40
120.643 ± 0.00
120.532 ± 0.03
120.455 ± 0.13
ZDT2
120.283 ± 0.01
119.019 ± 2.19
120.276 ± 0.03
119.001 ± 0.10
120.108 ± 0.16
ZDT3
128.752 ± 0.01
128.198 ± 1.21
128.749 ± 0.01
128.218 ± 0.18
128.453 ± 0.20
ZDT4
120.271 ± 0.55
118.496 ± 1.36
119.063 ± 1.49
119.127 ± 0.01
118.842 ± 1.02
ZDT6
116.090 ± 0.06
117.383 ± 0.03
115.938 ± 0.07
117.387 ± 0.12
116.281 ± 0.19
UF1
14.812 ± 0.01
14.789 ± 0.02
14.715 ± 0.01
14.943 ± 0.01
15.018 ± 0.02
UF2
5.359 ± 0.02
5.328 ± 0.02
5.343 ± 0.00
5.323 ± 0.02
5.437 ± 0.00
UF3
16.295 ± 0.01
16.248 ± 0.01
16.371 ± 0.02
16.705 ± 0.01
16.4009 ± 0.01
DTLZ1
0.177 ± 0.02
0.184 ± 0.00
0.188 ± 0.00
0.164 ± 0.01
0.168 ± 0.03
DTLZ2
0.690 ± 0.01
0.709 ± 0.00
0.740 ± 0.00
0.718 ± 0.01
0.721 ± 0.00
DTLZ3
0.961 ± 3.94
26.033 ± 1.14
2.032 ± 4.09
4.352 ± 3.21
7.229 ± 9.78
DTLZ4
0.646 ± 0.18
0.454 ± 0.22
0.637 ± 0.14
0.710 ± 0.12
0.721 ± 0.01
UF8
2.020 ± 0.01
2.021 ± 0.02
2.022 ± 0.01
2.016 ± 0.02
2.021 ± 0.00
UF9
1.435 ± 0.02
1.437 ± 0.01
1.435 ± 0.00
1.398 ± 0.01
1.438 ± 0.01
UF10
1.096 ± 0.01
1.099 ± 0.01
1.099 ± 0.02
1.1139 ± 0.03
1.100 ± 0.01
Additionally, a visual comparison was performed in order to help to understand the obtained results. In this experiment the median of the executions with respect to the hypervolume indicator was plotted. In Figure 2 the produced Pareto fronts are shown. From this figure, it can be noticed that the R2MOPSO is competitive with respect to the stateoftheart algorithms in most of the problems.
Graphical results of the visual comparison (Pareto fronts) for all the test problems.
4.2. Experiment 2 (Statistical Analysis)
Additionally, a statistical analysis was performed in order to verify the statistical differences of the proposed approach with respect to the SMSEMOA and MOMBIII. The comparison methodology is as follows: first, the ShapiroWilk test was performed to verify the normality of the data distributions. If both samples are normally distributed, the homogeneity of their variances is verified with Bartlett’s test. If both variances are homogeneous then an ANOVA test is performed to verify the statistical significance; otherwise Welch’s test is used. For nonnormally distributed data, the KruskalWallis test was performed to verify the statistical significance. In all the test a significance level of 0.05 was used. The results are presented in Table 4. The values marked with “+” suggest a statistically significant difference in performance in favor of our proposal, while the values marked with “−” suggest a difference in favor of the stateoftheart algorithms. Additionally, the values in blank suggest a similar performance. The obtained results confirm the above discussion. The proposed approach showed a similar statistical performance with respect to the MOMBIII. However, in some problems (UF1, UF2, DTLZ1, DTLZ3, and DTLZ4), a statistically significant difference in favor of the proposed approach with respect to the SMSEMOA can be observed.
HV: statistical analysis.
ZDT1
ZDT2
ZDT3
ZDT4
ZDT6
UF1
UF2
UF3
DTLZ1
DTLZ2
DTLZ3
DTLZ4
UF8
UF9
UF10
R2MOPSO/SMSEMOA
+
+
+
+
+
R2MOPSO/MOMBIII
+
+
+
+

4.3. Experiment 3 (ManyObjective Problems)
Since one of the aims of using an indicatorbased MOEA is its capability to perform well in the presence of many objective functions, we decided to test the behavior of the R2MOPSO in such problems. For this experiment, we focused our efforts on solving the DTLZ{1–4} test problems in order to investigate the behavior of our proposed approach as the number of objectives of the problem increases (we increased from 4 to 10 objectives). Our results are compared with respect to those obtained by SMSEMOA and MOMBIII. Each MOEA was executed 30 times and their results were evaluated using the hypervolume indicator. The reference points used by the adopted problems were of [1.5] for DTLZ1 and [2.0] for the remaining problems. The parameters were similar to those adopted in the previous experiment. Since the computational time required by the original SMSEMOA algorithm (using exact hypervolume) becomes prohibitive very quickly, the number of objectives is raised (this has also been illustrated in other works [17]). For this reason, we decided to compare our approach with a modified version of SMSEMOA. This version uses the approximation to the contribution to the hypervolume (proposed in [9]) in order to decrease the execution time of the algorithm. The number of samples used in this latter algorithm is 100,000.
Additionally, the same statistical analysis performed in experiment two was realized in this experiment. The values marked with “(+)” suggest a statistically significant difference in performance in favor of our proposal, while the values marked with “(−)” suggest a difference in favor of the stateoftheart algorithms. Moreover, the values in with “()” suggest a similar performance.
The results for the compared approaches are shown in Tables 5, 6, 7, and 8. In some cases, the hypervolume has a value of zero. This value indicates that the algorithm did not not achieve any nondominated solution with respect to the reference point. For DTLZ2 and DTLZ4 shown in Tables 6 and 8, respectively, the results of the compared approaches were similar (these two problems were easiest to solve in our benchmark). However, for DTLZ2 from 7 objectives, the R2MOPSO is slightly better than SMSEMOA, and in DTLZ4 the R2MOPSO is always better than SMSEMOA. Moreover, in the statistical analysis the results indicate that the R2MOPSO is better with respect to the MOMBIII and the SMSEMOA. On the other hand, for DTLZ1 and DTLZ3 (which are shown in Tables 5 and 7, resp.), our approach clearly outperforms the SMSEMOA and the MOMBIII. For DTLZ1 the performance of SMSEMOA decreased as we increased the number of objectives. The main problem of this algorithm is that the exact hypervolume calculation was replaced with an approximation, and when the number of objectives is raised it is necessary to increase the number of samples for the approximation as well. However, if the number of samples is increased, the computational cost would also increase. Moreover, in the DTLZ1 the numeric results are very similar for MOMBIII and R2MOPSO. However, as in DTLZ4 in the statistical analysis the R2MOPSO is better than the MOMBIII. Finally, for DTLZ3 our R2MOPSO clearly outperformed the compared approaches. In this challenging problem, the SMSEMOA was unable to converge.
Comparison of the hypervolume’s average obtained by MOMBIII and R2MOPSO for the DTLZ1 test problem with different number of objectives (M).
M
SMSEMOA
MOMBIII
R2MOPSO
4
15.740759 ( )
15.96795 ( )
15.916804
5
23.495387 (+)
31.678247 ( )
31.647112
6
15.939020 (+)
63.881022 ( )
63.871206
7
37.503944 (+)
126.443144 ( )
127.956645
8
58.822058 (+)
253.535178 (+)
255.534288
9
71.723936 (+)
509.121631 (+)
510.093877
10
388.922373 (+)
1020.781209 (+)
1022.891997
Comparison of the hypervolume’s average obtained by MOMBIII and R2MOPSO for the DTLZ2 test problem with different number of objectives (M).
M
SMSEMOA
MOMBIII
R2MOPSO
4
15.521335 ( )
15.490132 ( )
15.485958
5
31.598023 ( )
31.599134 ( )
31.581132
6
63.648437 ( )
63.639010 ( )
63.638744
7
127.672459 ( )
127.680517 ( )
127.680623
8
255.683900 ( )
255.693671 ( )
255.707244
9
511.666530 ( )
511.665803 ( )
511.660240
10
1023.617417 ( )
1023.521349 ( )
1023.471587
Comparison of the hypervolume’s average obtained by MOMBIII and R2MOPSO for the DTLZ3 test problem with different number of objectives (M).
M
SMSEMOA
MOMBIII
R2MOPSO
4
0.0 (+)
0.310916 (+)
0.529142
5
0.0 (+)
1.231819 (+)
1.389182
6
0.0 (+)
0.201825 ( )
0.201977
7
0.0 (+)
2.351874 (+)
2.803984
8
0.0 (+)
5.925318 (+)
6.112229
9
0.0 (+)
38.102465 (+)
39.989576
10
0.0 (+)
15.234512 (+)
16.484128
Comparison of the hypervolume’s average obtained by MOMBIII and R2MOPSO for the DTLZ4 test problem with different number of objectives (M).
M
SMSEMOA
MOMBIII
R2MOPSO
4
15.110649 (+)
15.221809 (+)
15.480663
5
31.489328 (+)
31.531203 (+)
31.580678
6
63.465142 (+)
63.530765 (+)
63.599144
7
127.641483 (+)
127.642368 (+)
127.678039
8
255.711370 (+)
255.721893 ( )
255.734554
9
511.741231 (+)
511.761209 (+)
511.781234
10
1023.679123 (+)
1023.680125 (+)
1023.691262
5. Conclusions and Future Work
In this work, the incorporation of an indicator to guide the search of a PSO was proposed. Our proposed approach was validated using two and threeobjective function problems. The proposed R2MOPSO obtained competitive results when compared to NSGAII, MOEA/D, and SMSEMOA using several test problems. Therefore, we can say that our approach is successful to work with multiobjective problems. However, since our main target was manyobjective problems, we decided to study our approach with respect to this sort of problem. For this sake, we adopted scalable test functions and we compared our results with respect to a variation of a wellknown hypervolumebased approach (SMSEMOA) that approximates the hypervolume contributions, in order to have a more efficient performance. Our results indicate that our approach outperforms the SMSEMOA with respect to the hypervolume. Finally, an important feature of this proposal is that it does not adopt Pareto dominance to guide the search; it is only adopted to report the found solutions.
As part of our future work, we would like to further investigate the use of the R2 indicator to restrict the size of external archives in other types of MOEAs.
Competing Interests
The authors declare that they have no competing interests.
CuevasE.DíazM.A method for estimating view transformations from image correspondences based on the harmony search algorithmGuoZ.HuangH.DengC.YueX.WuZ.An enhanced differential evolution with elite chaotic local searchSurajTiwariP.GhoshS.SinhaR. K.Classification of two class motor imagery tasks using hybrid GAPSO based Kmeans clusteringHughesE. J.Evolutionary manyobjective optimisation: many once or one many?1Proceedings of the IEEE Congress on Evolutionary ComputationSeptember 2005Edinburgh, UK22222710.1109/CEC.2005.1554688KhareV.YaoX.DebK.Performance scaling of multiobjective evolutionary algorithmsProceedings of the 2nd International Conference on Evolutionary MultiCriterion Optimization (EMO '03)2003Faro, PortugalSpringer376390PurshouseR. C.FlemingP. J.Evolutionary manyobjective optimisation: an exploratory analysis3Proceedings of the Congress on Evolutionary Computation (CEC '03)December 20032066207310.1109/cec.2003.12999272s2.084901459041HughesE. J.KnowlesJ.CorneD.KalyanmoyD.ChairD.Fitness assignment methods for manyobjective problemsSatoH.Analysis of inverted PBI and comparison with other scalarizing functions in decomposition based MOEAsBaderJ.ZitzlerE.HypE: an algorithm for fast hypervolumebased manyobjective optimizationGómezR. H.CoelloC. A. C.Improved metaheuristic based on the R2 indicator for manyobjective optimizationProceedings of the Annual Conference on Genetic and Evolutionary Computation (GECCO '15)2015ACM679686MenchacaMendezA.CoelloC. A. C.GasparCunhaA.AntunesC. H.CoelloC. C.GDMOEA: a new multiobjective evolutionary algorithm based on the generational distance indicatorZitzlerE.KünzliS.YaoX.BurkeE. K.LozanoJ. A.Indicatorbased selection in multiobjective searchZitzlerE.BeumeN.NaujoksB.EmmerichM.SMSEMOA: multiobjective selection based on dominated hypervolumeIshibuchiH.TsukamotoN.NojimaY.Iterative approach to indicatorbased multiobjective optimizationProceedings of the IEEE Congress on Evolutionary Computation (CEC '07)September 2007SingaporeIEEE3967397410.1109/cec.2007.44249882s2.079955357239IgelC.HansenN.RothS.Covariance matrix adaptation for multiobjective optimizationRodríguez VillalobosC. A.Coello CoelloC. A.A new multiobjective evolutionary algorithm based on a performance assessment indicatorProceedings of the 14th International Conference on Genetic and Evolutionary Computation (GECCO '12)July 2012Philadelphia, Pa, USAACM50551210.1145/2330163.23302352s2.084864659168SchützeO.EsquivelX.LaraA.CoelloC. A. C.Using the averaged hausdorff distance as a performance measure in evolutionary multiobjective optimizationAlamS.DobbieG.RehmanS. U.Analysis of particle swarm optimization based hierarchical data clustering approachesSinghR. P.MukherjeeV.GhoshalS. P.Particle swarm optimization with an aging leader and challengers algorithm for optimal power flow problem with FACTS devicesYeZ.WangM.HuZ.LiuW.An adaptive image enhancement technique by combining cuckoo search and particle swarm optimization algorithmPadhyeN.Topology optimization of compliant mechanism using multiobjective particle swarm optimizationProceedings of the 10th Annual Genetic and Evolutionary Computation Conference (GECCO '08)July 2008Atlanta, Ga, USAACM183118342s2.057349170287JiaS.ZhuJ.DuB.YueH.Indicatorbased particle swarm optimization with local searchProceedings of the 7th International Conference on Natural Computation (ICNC '11)July 2011Shanghai, China1180118410.1109/icnc.2011.60221682s2.080053405652HansenM. P.JaszkiewiczA.Evaluating the quality of approximations to the nondominated set1998IMM, Department of Mathematical Modelling, Technical University of DenmarkBrockhoffD.WagnerT.TrautmannH.On the properties of the R2 indicatorProceedings of the 14th International Conference on Genetic and Evolutionary Computation (GECCO '12)July 2012Philadelphia, Pa, USAACM46547210.1145/2330163.23302302s2.084864669908DíazManríquezA.ToscanoPulidoG.CoelloC. A. C.LandaBecerraR.A ranking method based on the R2 indicator for manyobjective optimizationProceedings of the IEEE Congress on Evolutionary Computation (CEC '13)June 2013Cancun, MexicoIEEE1523153010.1109/cec.2013.65577432s2.084881592325DíazManríquezA.ToscanoPulidoG.CoelloC. A. C.LandaBecerraR.A ranking method based on the R2 indicator for manyobjective optimizationProceedings of the IEEE Congress on Evolutionary Computation (CEC '13)June 2013Cancun, Mexico1523153010.1109/cec.2013.65577432s2.084881592325PhanD. H.SuzukiJ.R2IBEA: R2 indicator based evolutionary algorithm for multiobjective optimizationProceedings of the IEEE Congress on Evolutionary Computation (CEC '13)June 20131836184510.1109/cec.2013.65577832s2.084881601149WessingS.NaujoksB.Sequential parameter optimization for multiobjective problemsProceedings of the IEEE Congress on Evolutionary Computation (CEC '10)July 2010Barcelona, ItalyIEEE1810.1109/CEC.2010.5586529ZitzlerE.KnowlesJ.ThieleL.Quality assessment of pareto set approximationsKennedyJ.EberhartR. C.JaszkiewiczA.MOMH multipleobjective metaheuristicshttp://home.gna.org/momh/index.htmlZitzlerE.DebK.ThieleL.Comparison of multiobjective evolutionary algorithms: empirical resultsZhangQ.ZhouA.ZhaoS.SuganthanP. N.LiuW.TiwariS.Multiobjective optimization test instances for the cec 2009 special session and competitionDebK.ThieleL.LaumannsM.ZitzlerE.AbrahamA.JainL.GoldbergR.Scalable test problems for evolutionary multiobjective optimization