Characterizing the preferences of a decision maker in a multicriteria decision is a complex task that becomes even harder if the information available is limited. This paper addresses a particular case of project portfolio selection; in this case, the measures of project impacts are not assumed, and the available information is only projects’ ranking and costs. Usually, resource allocation follows the ranking priorities until they are depleted. This action leads to a feasible solution, but not necessarily to a good portfolio. In this paper, a good portfolio is found by solving a multiobjective problem. To effectively address such dimensionality, the decision maker’s preferences in the form of a fuzzy relational system are incorporated in an ant-colony algorithm. The Region of Interest is approached by solving a surrogate triobjective problem. The results show that the reduction of the dimensionality supports the decision maker in choosing the best portfolio.

1. Introduction

Resource allocation in institutions should address the proper distribution of a given budget among a set of available projects [1]. The type of projects can vary according to the area, and they might be involved in a wide range of tasks or activities, such as improving the skills of a professional sports team, the selection of R&D project portfolios in enterprises, the founding of projects by a government program, and supporting environmental regulations [2].

The construction of the best portfolio that accomplishes a certain balance among the selected projects and that is subject to a budget has been approached in the scientific literature (e.g., [38]). This problem can be defined as follows: where is the space of feasible portfolios and represents the functions that characterize the impact of a portfolio over the considered criteria.

Problem (1) has evolved in different particular cases. For example, recently, the Portfolio Selection Problem on a Set of Ordered Projects, or PSPSOP, has been reported to focus on the construction of a portfolio from a set of projects subject to a limited budget (cf. Bastiani et al. [9]). This case also involves a distinctive feature, which is that the only information available about the projects is their rank; that is, they are ordered according to the decision maker’s (DM) preferences. The importance of the study of this problem arises from the fact that in most situations a DM prefers simple decision methods, and the decision process in such methods has conditions that involve limited resources such as time and information.

The task in PSPSOP is closely related to R&D projects funded by enterprises; it can be seen in the discussion presented in Cooper et al. [10]. Particularly, in government organizations, R&D projects pursue the satisfaction of the citizenship, and, for this purpose, they try to identify programs that focus on the priorities of social sectors. The participatory budgeting aids in the definition of such priorities, and it combines the efforts of both citizens and government to establish them (cf. Fernandez and Olmedo [11]). With a proper definition of priorities, all that remains is the construction of a portfolio in which the costs of the projects adjust to the approved budget. This last task remains a challenge in the general case, and it has been studied from different approaches in the literature.

The construction of portfolios from projects that include information on ranking and costs is a problem that has been addressed in the scientific literature. The Ranking method is commonly used, and it involves the ordered selection of projects, previously ranked by priorities established by a DM, according to an available budget. This method has been criticized because the final portfolio that it creates only guarantees that the most important projects will be supported without considering a balance between the priorities and costs, which usually increases the number of supported projects and their impact in society (cf. [12, 13]). Additionally, it could be possible that a DM becomes reluctant to follow such a construction if he considers that the ranking information is not reliable [14].

Bastiani et al. [9] propose a method to solve PSPSOP based on the cardinality and the discrepancies present in a portfolio, where cardinality refers to the total number of projects involved in it, and the term discrepancy is a concept that reflects the negative effect that is applied over the DM’s thinking because one of the projects, when it is compared against others, seems to have merits that belong to the portfolio but it is not in it. The strategy balances the priorities and the number of projects in the final portfolio through a model that minimizes discrepancies concerning the ranking and costs and maximizes its cardinality, all of which are defined in ten objectives. Given that this model handles a large number of objective functions, this paper proposes a method to reduce the original many-objective problem to a surrogate one with three-objective functions. The reduced problem can be easily solved by an ant-colony algorithm that incorporates a knowledge mechanism (like the one proposed by Cruz et al. [15]).

This paper is organized into six sections. Section 2 presents a criticism of previous related approaches and describes the many-objective optimization model proposed by Bastiani et al. [9] for PSPSOP. Section 3 details the decision-support mechanisms used for the PSPSOP; this mechanism is based on DM preferences. Section 4 describes the metaheuristic optimization method. Finally, the experiments and results and conclusions are shown in Sections 5 and 6, respectively.

2. A Brief Outline of Previous Approaches

According to the reviewed scientific literature, there are different approaches that solve PSPSOP with a lack of available information. The works of [1618] are similar to Scoring and Ranking methods [10] or additive functions [12, 19] in that they prioritize projects according to a certain utility function to measure their importance. Alternatively, the use of proxy variables [9, 11, 14] has offered versatile and satisfactory results that extend the information derived from a ranking of projects. However, none of those approaches offer a strategy that incorporates the DM’s preferences, the key element in the present work that guides the construction of better solutions. Table 1 analyzes the state of the art related to the research proposed in these works. The first column presents the work, and the second column describes it.

From the information provided in Table 1, all the strategies followed to solve PSPSOP are based on ranking, costs, discrepancies, and cardinality. Our interest is focused on the definition and management of such elements given by the work of Bastiani et al. [9]. The remainder of this section briefly discusses their definition, but it is left up to the reader whether to get a deeper understanding of the concepts by reviewing the work in [9].

The model proposed by Bastiani et al. [9] defines three objectives per category plus one for the power indicator . Hence, the number of proxy variables is , where is the number of categories in the instance of PSPSOP that is solved. The model uses a reference portfolio to reflect the DM’s disappointment in combination with the statement “ is better than ” to define three types of discrepancies over a built portfolio C, which are initially based on the idea that a discrepancy occurs when in such a portfolio but . These discrepancies are as follows: (a) weak discrepancy that occurs when the budget of a project , which is not in C, is much higher than the average budget in category k; (b) strong discrepancy that occurs when the budget of a project , which is not in C, is considerably higher than the average budget in category k; and (c) unacceptable discrepancy that occurs whenever a project , which is not in C, has a budget that is not significantly higher than the average budget in category k, or the budget required by another project that has a lower rank than . The first two discrepancies, in combination with the cardinalities per category and the power of the set of objectives to the model, are proposed in [9]. The unacceptable discrepancy is used to constrain the feasible region of PSPSOP.

The model of Bastiani et al. [9], referred to as Bastiani’s model in this work, succeeds in providing more information to a DM related to the construction of a particular portfolio. The success is derived from the specialization per category of the information derived from the ranking. However, due to this situation, the number of objectives involved in the model is not fixed, and it varies with the number of considered categories. This situation allows growth in the dimensionality, which in turn increases the difficulty in approximating the Region of Interest, that is, the zone in the Pareto frontier that fits a DM’s preferences. For this reason, the question of whether it is possible to keep the same level of information in the construction of a portfolio through an optimization approach but reduce the dimensionality remained open. Our work in this paper is dedicated to answering this question positively. The proposed approach is presented in the following section.

3. Proposed Surrogate Model for PSPSOP

According to Cruz et al. [15], metaheuristic approaches are a viable strategy to solve multiobjective optimization problems. These approaches rely on their ability to approximate the Pareto frontier to provide a set of solutions (or populations) that might satisfy complex constraints and objectives functions, a feature rarely shared with approaches based on mathematical programming, which normally yield a single solution and are limited to a narrower range of variations in the definition of an optimization problem.

During the optimization, it is possible to find some pitfalls derived from a problem with many objectives. One of the pitfalls arises if it is difficult to generate a proper population that lies on the Pareto frontier. Another pitfall corresponds to the increments in the number of Dominance Resistant Solutions, that is, solutions that are hard to dominate even though they are not Pareto optimal and are hard to take out of the final set of solutions. Another concern is related to the appropriate selection of one solution from the set given by metaheuristics. This situation leads to a process in which a human must perform a cognitive effort to choose the desired alternative. However, according to Miller [20], the human mind’s capacity to handle information is diminished when it increases, for example, when the number of objectives that he/she must use to make a decision increases. All these drawbacks have been observed and/or discussed in several research studies, as shown in [2124].

The high dimensionality of a problem can be addressed through the incorporation of preference information. According to the work in [15], a DM that is willing to provide his/her preferences aids in the localization of the Pareto frontier zone known as the Region of Interest (RoI); this zone can be understood as a relaxation of the feasible region defined by the multiobjective optimization problem that only covers solutions that are of interest for a DM. Several research papers, such as the one by Fernandez et al. [25], present strategies that address multiobjective problems using surrogate models that approximate the RoI. The remaining part of this section presents the adaptation of such model to solve the PSPSOP.

The research work in [25] proposes a surrogate model based on the credibility index of the statement “ is at least as good as .” The value of is computed using the method ELECTRE (cf. [26, 27]), and it is integrated into a relational system of preferences (as described by Roy [26]) to model the DM preferences. The preferences defined in [27] for pairs of alternatives and are as follows: (1) strict preference ; (2) indifference ; (3) weak preference ; (4) incomparability ; and (5) k-preference .

To define a surrogate model for PSPSOP, this work uses the net flow score defined through the objectives of Bastiani’s model, as shown in Table 2. The net flow score, denoted or NF, was introduced in Fernandez et al. [25] as a measure to enhance preference information towards a better characterization of the DM preferences in the nonstrictly outranked set, denoted by NS. Based on the elements , the following sets are defined: (a) the set of alternatives that strictly outrank ; (b) the set of portfolios that weakly outrank a given portfolio ; and (c) the set of alternatives with greater net flow score. All the previous elements are defined in Table 2, where is the set of feasible portfolios.

The combination of the net flow score with the definition of the sets and allows the formulation of the optimization problem defined in

With the previous problem, it is demonstrated that Bastiani’s model can be mapped to problem (2) independently of the original objective space dimension. The approach by Fernandez et al. [25] was successfully applied to Problem in [15]. Combined with an ant-colony metaheuristic, the approach in Bastiani et al. [9] is extended here to solve PSPSOP thru the optimization model defined in problem (2). To summarize, the contribution of this work is the solution of PSPSOP with lack of information using a surrogate objective model (which is of smaller dimension than those reported in previous approaches) and the incorporation of DM’s preferences in the search process to approximate the RoI. The details of the algorithm and each component of the ACO-SOP are presented in the next section.

4. An Ant-Colony Optimization Algorithm

This work proposes an approach based on an ant-colony optimization, denoted by ACO-SOP (Ant-Colony Optimization for Solving Portfolio Problems with Ordinal Information about Projects). The approach takes ideas from Dorigo’s ACS [28], Cruz et al. [15], and Bastiani et al. [9] to solve PSPSOP, first searching in the wider feasible region defined by Bastiani’s model and then searching in the smaller space defined by the model in problem (2).

The core strategy of ACO-SOP uses the pheromone representation, selection rule, and local search function defined in the algorithm ACO-SPRI proposed by Bastiani et al. [9], with a small variation in the construction of the nondominated fronts used in their computations. The general idea of ACO-SOP is depicted in Algorithm 1. The algorithm is characterized by five elements: (a) the initialization function; (b) the construction of a feasible set of solutions; (c) the improvement phase based on local search; (d) the construction of nondominated fronts; and (e) the updating of the pheromone matrix. These elements are briefly defined in this section (cf. Bastiani et al. [9] for further details).

Input: , , max_iter,
Output: The new set of portfolios NS
(1) Initialize Iter = 0, pheromone matrix and
(2) Construct an initial portfolio
(3) Repeat
(4) Initialize
(5) <- GenerateFeasibleSolutions()
(6) <- Perform local search to the set
(7) Calculate objective functions of Problem (2) on the set
(8) <- Generate non-dominated fronts on
(9) Updating the pheromone matrix with the set
(10) Assign
(11) iter = iter + 1
(12) until (iter = max_iter)
(13) <- Generate non-dominated fronts on
(14) return  

The algorithm ACO-SOP requires as input the set of ranked projects, the budget , the maximum number of iterations max_iter, and the maximum number of ants . Inside this algorithm, the initialization function constructs the reference portfolio (see Line ); it is done by the ordered selection of projects from by rank and in agreement with the available budget (ties are broken arbitrarily).

The phase of construction of feasible solutions (see Line ) involves the construction of one portfolio by each of the ants. Each ant starts with an empty portfolio and adds to it one project at a time based on the budget, global knowledge derived from the pheromone matrix, and local knowledge obtained from the preferences of the DM, as defined in the selection rule by Bastiani et al. [9]. At the end of this phase, a set of feasible portfolios is delivered.

Following the process of forming , the algorithm ACO-SOP performs a local search over each in , just as it is done by Bastiani et al. [9] in their local search scheme. The algorithm creates a new portfolio with every combination derived from the inclusion and/or exclusion of projects of chosen at random. Each portfolio is repaired such that constraints in rank and budget are held. Finally, the best portfolio is chosen as the local search improvement for . The union of all the new portfolios will form the new set of feasible portfolios .

The construction of nondominated fronts (see Line ) is based on the objectives defined in problem (2). The output of the algorithm in this phase is the set of fronts obtained from . Each set is composed of a subset of portfolios in that are exclusively dominated by exactly portfolios, where (following the dominance criterion established in Cruz et al. [15]). Instead of computing the dominance from the objectives of Bastiani’s model, this work uses the objectives of problem (2), a key feature in this research because it makes it more manageable for algorithms to address large dimensions.

Finally, the last and most important element that defines ACO-SOP is the pheromone matrix, and the process to construct it and update it is described herein. The bidimensional matrix represents the knowledge that the ants have gained during the construction of portfolios. They represent that knowledge in the form of pairs of projects and the gain of having them together in the same portfolio (denoted as ). The portfolios constructed by the ants at each iteration are used to update τ (see Line ) based on the number of fronts constructed and in the front where each portfolio is found. This strategy and the one incorporated to prevent premature convergence are detailed in Bastiani et al. [9].

The algorithm ACO-SOP accumulates in the front constructed at each iteration (see Line ). Then, using the last set it forms the final set of fronts and returns as its solution the front . Figure 1 presents a brief example of the construction process of ACO-SOP. Step shows how the instance information is transformed with the aid of the pheromone matrix into a set of feasible solutions , where each solution is binary vector representing the projects in a portfolio. Step shows the computed values of the 10 objectives from Bastiani’s model for each solution. After that, in step 3, ACO-SOP shows an improved set of solutions from the local search. Finally, step 4 presents the transformation of into the surrogate model, that is, a set with the objective values defined in problem (2). To exemplify the transformation, let us take the values of the solution , which indicate that there is no other solution , for , that is strictly or weakly preferred over ; that is, there is no such that or or , and there is only one solution that has greater net flow score than ; that is, , taking into account the fact that is in the nonstrictly outranked set .

The following section presents the experimental design performed to evaluate the performance of ACO-SOP when solving problem (2). The following section also shows the versatility of the approach to tackle cases of large dimensions.

5. Solving Some Computer Experiments

In this section, we develop a set of experiments to verify the validity of and to validate the advantages gained by our solution method, especially addressing problems with many objectives. Based on these goals, the exposure and resolution of two different cases of study are presented in the following paragraphs.

To test the performance of ACO-SOP, the algorithm was implemented using the programming language Java. The experimental design was made using a computer with an Intel Core i7 2.8 GHz CPU, 4 GB of RAM, and Mac OS X Lion 10.7.5 (11G63b) as operating system. The instance used for the experimentation is the one defined in Bastiani et al. [9]. The tuning process involved selecting different combinations of values similar to those previously reported for these parameters. According to the results of the fine-tuning process, the best configuration of values was as follows: tot_iter = 200 iterations, = 300 ants, = 0.63,η = 0.1,ρ = 0.9, = 0.65, and = 0.75. These parameter values were used to obtain the experimental results reported in this section. The incorporation of preferences through the outranking model is basic in order to map the problem in Bastiani et al. [9] into problem (2). The parameters of the outranking model were set as provided in Table 3.

The ACO-SOP strategy requires, besides the ranking information and costs of the PSPSOP instance, the definition of a reference portfolio and the definition of the projects’ categories. The heuristic followed to construct is based on the traditional Scoring and Ranking method and takes the projects one at a time in order of their rank until the available budget is consumed. The categories considered during the experiment were (1) priority, (2) satisfactory, and (3) acceptable, and the projects were distributed uniformly among them based on their ranks.

5.1. Analysis of Results

The results from the experiments performed on the proposed approach (ACO-SOP) are analyzed from three points of view. The first analysis involves the quality of the solutions provided by ACO-SOP in relation to their meaning to a DM and his/her preferences. The second and third analyses study the performance of ACO-SOP in comparison with a similar approach of the state of the art. These experiments differ on the quality indicator used during the comparison. Let us indicate that the considered budget was 2,500 units and that the reference portfolio is formed by the first 22 best ranked projects in the instance taken from Bastiani et al. [9].

The indicators used to measure the results from the experiment were , , and . The indicator uses the three objectives associated with the preference information derived from the proposed surrogate model (see problem (2)). The indicator uses the ten objectives associated with ranking and derived from the Bastiani model (see Section 2). Finally, the indicator is the content of the portfolio, that is, the projects that are included in it.

The first analysis involved the solution of the instance with ACO-SOP, and its results are summarized in Table 4. This table provides the values of the indicators and computed from the three best nonstrictly outranked solutions that were consistently obtained in 30 runs (the average computer time required for this experiment was 130 minutes). The solutions ACO-SOP1, ACO-SOP2, and ACO-SOP3 of Table 4 are ordered according to the indicator ; in this ordering, ACO-SOP1 is the best because it reduced all the values of the indicator to zero, while ACO-SOP2 did it only for the metrics and but it improves ACO-SOP3 in the last one, . This relation of order is associated with the DM’s preferences, and evidence of that can be supported by analyzing the values of the indicator , regarding the number of projects involved in the solutions.

First, it is possible to observe a great level of dissatisfaction in the DM when the provided solutions have strong discrepancies (see value of ACO-SOP3), because it means that some important projects that are in the reference portfolio are left out of the solution when they could have been included in place of others that were less relevant. Second, the DM could have a low level of dissatisfaction if the solutions only have weak discrepancies (see value of ACO-SOP2), because it means that only a rather small number of important projects had been left out. Third, an almost null level of dissatisfaction can be noted in the DM when the solution involved a small number of weak discrepancies and a high power (see values and of ACO-SOP1), because this is an indicator that an adequate number of high-ranked projects are in the solution.

The second analysis compares the two best solutions ACO-SOP1 and ACO-SOP2 against the solutions Sol1 and Sol2, as reported in [11]. Table 5 summarizes these results, with each of its columns representing a metric used by indicator . Observe that the lack of discrepancies on Sol1 and Sol2 (see , , , , , and ) could allow the DM to make a decision using only two criteria, the power and the cardinality ; this is a convenient situation for a DM that willingly accepts a solution with only high-ranked projects, or with a large cardinality, but does not care about their balance (because high values in are normally associated with small cardinalities and high-ranked but costly projects). But what happens if the DM has a different point of view? For example, how are these solutions affected if, for the DM, the objectives are not equally important? The latter case involves considering solutions, such as ACO-SOP1 and ACO-SOP2, with an acceptable level of discrepancies and a better balance between and the cardinalities (e.g., and ); these solutions slightly increase the discrepancies (as in metric ). Hence, while approaches such as the ones presented in [11] or in Bastiani et al. [9] succeeded in managing problems with limited information, they do not include DM preferences. This analysis showed evidence that the solutions provided by the proposed model allow the inclusion of DM’s preferences through the parameters of the outranking model (see Table 3). The preference model can be used to redefine what is acceptable or not as a solution for a DM; for example, the solutions Sol1 and Sol2 do not satisfy the DM’s preferences in Table 3 because they have unacceptable discrepancies that cause great discomfort to the DM; this situation is detailed in the following analysis.

Table 6 allows the analysis of the configurations of the portfolios obtained from Sol1, Sol2, ACO-SOP1, and ACO-SOP2; it shows the projects in each solution and in bold those in Category 1. One important difference between these solutions is the situation of Project 20 (or Pj20); while it appears on the portfolio of ACO-SOP1 and ACO-SOP2, it does not in Sol1 and Sol2. The relevance of this project is that it belongs to the reference portfolio , and the fact that it does not appear in Sol1 and Sol2 constitutes an unacceptable discrepancy for the DM characterized in Table 3.

To explain this situation, let us note that Pj20 is in Category 1 and has an average cost of and the projects of Sol1 that are in the same category have an average cost of . Hence, the unacceptability level of is . The value of is smaller than the unacceptability threshold defined for the DM; that is, the projects in the portfolio are too costly to exclude from and possibly some others, a situation that occurs with Sol1 and Sol2 and that generally does not satisfy the DM’s preferences. Such solutions where the DM’s preferences are not considered can be constructed using approaches from [9, 11], but the present approach overcomes it by including preferences during the search process; as a result, the solutions contain portfolios more convenient for the DM, because they include an increment in the number of projects of up to 20%, with the guarantee that this still satisfies his/her preferences.

To summarize, this section has presented the results from the experiment to analyze the quality of the solutions provided by ACO-SOP and their robustness. The robustness was demonstrated by providing evidence that ACO-SOP provides solutions of the same quality as the other approaches, whenever they satisfy the DM’s preferences; evidence of this is provided in Table 4, where solutions of similar quality are presented, solutions that are in the RoI for a DM. In general, whenever a solution fits the interest of the DM expressed in the parameter values of the preference model, this solution can be found in the search space, showing that the present approach maintains the same level of information. Finally, the lack of preferences on approaches such as [11] or [9] can lead to the construction of solutions that are not in the RoI for a DM, that is, solutions the DM is not interested in.

6. Some Conclusions

This paper proposes a knowledge-based decision-support system for the Portfolio Selection Problem on a Set of Ordered Projects. This problem commonly has a lack of available information, which is why the incorporation of knowledge mechanisms in the solution strategies can improve the quality of the solutions.

This work is a crucial refinement from a recent proposal that modeled the DM’s attitude through a many-objective optimization problem. The high dimensionality of this problem is a great concern for metaheuristic approaches in generating an acceptable approximation to the Pareto frontier, and it is also a great concern for the DM in making the final decision due to human cognitive limitations. Here, we have presented a method that incorporates the DM’s preferences through a fuzzy outranking-based relational model. This allows mapping the original many-objective problem into a surrogate three-objective problem. With this transformation, the process of identifying the Region of Interest (the preference privileged zone of the Pareto frontier) becomes easy. Another advantage of the proposal is to be very robust with respect to an increasing number of objective functions. Robustness is important to make a finer representation of the DM’s preferences, including project synergic effects. So this contribution should be significant in this particular field of project portfolio optimization.

The strongly nonlinear surrogate problem is solved by using an ant-colony algorithm. This algorithm includes a knowledge and learning mechanism in the form of a pheromone matrix and integrates a decision-support mechanism based on the DM’s preferences. The experimental results show its capacity to obtain consistent solutions in the Region of Interest. These solutions surpass in quality others obtained by alternative proposals taking into account the amount and impact of the supported projects.

Let us observe that the surrogate model not only reduces the dimensionality of the problem but also reduces the number of alternatives that result from the multiobjective optimization process, because it provides only nonoutranked solutions that are on the privileged zone of the known Pareto frontier. Comparisons between pair of alternatives are not so difficult as the number of objectives seems to suggest. Note that the objectives and are strongly more important than the remaining ones. So, the DM can apply a sort of lexicographic priority that allows a radical reduction of the set of solutions. Even though it is possible that ACO-SOP produces a number of portfolios that could not be easily handled by the DM, this situation can be addressed by the support of an expert analyst who, using a multicriteria decision method, provides a sufficiently small set of alternatives that can be handled by the DM.

Another point that is worthy of discussion is that the projects might involve objectives that are hard to judge, as is the case of the risk, because they are complex in both quantitative and qualitative elements. However, it is important to note that such cases are handled in the present approach through the ranking. Following this idea, let us note that this work addresses the special case of PSPSOP that has as the only available information the ranking of the projects and their costs and that it is the rank that integrates all the previous existing information about the criteria that characterized the projects. In other words, the instances of PSPSOP already capture information about criteria such as risk in the rank, and the methodology followed to do such integration is out of the scope of the present research. However, if it were the case that an alternative to solve the problem of integrating objectives into ranking should be proposed, it could be done through the family of ELECTRE methods [26], which can address qualitative and quantitative information in a unique aggregation model of preferences.

Finally, some situations that can be present in the instances of the special case of PSPSOP studied here are the existence of interdependence among criteria and/or projects. Fuzzy multicriteria methods derived from outranking approaches can handle interdependence among criteria and reflect it on the rank ordering of candidate projects (cf. [29]). So, the rank used by our proposal takes into account the criterion interdependence. On the other hand, the project interdependence remains a challenge for the proposed model; it is defined as part of the future work of this research to explore the possibility of integrating synergy negative and positive effects among projects through the incorporation of artificial projects that are related to their interdependence.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


This work has been partially supported by the following CONACyT projects: (a) Fronteras de la Ciencias Project 1340; (b) Consolidation National Lab Project 280712; (c) Projects 236154 and 269890; (d) Project 280081, Red Temática para el Apoyo a la Decisión y Optimización Inteligente de Sistemas Complejos y de Gran Escala (OPTISAD), Universidad Autónoma de Nuevo León; and (e) Project 3058 from the program Cátedras CONACyT.