Stochastic Portfolio Selection Problem with Reliability Criteria

Portfolio selection focuses on allocating the capital to a set of securities such that the profit or the risks can be optimized. Due to the uncertainty of the real-world life, the return parameters always take uncertain information in the realistic environments because of the scarcity of the a priori knowledge or uncertain disturbances. This paper particularly considers a portfolio selection process in the stochastic environment, where the return parameters are characterized by sample-based correlated random variables. To decrease the decision risks, three evaluation criteria are proposed to generate the reliable portfolio selection plans, including max-min reliability criterion, percentile reliability criterion, and expected disutility criterion.The equivalent linear (mixed integer) programming models are also deduced for different evaluation strategies. A genetic algorithm with a polishing strategy is designed to search for the approximate optimal solutions of the proposedmodels. Finally, a series of numerical experiments are implemented to demonstrate the effectiveness and performance of the proposed approaches.


Introduction
Portfolio selection problems deal with how to allocate one capital to a given number of securities so that the involved return can be maximized in this process.When the return of each security is a constant, this problem can typically be formulated as a linear programming model and efficiently solved through the simplex method.Due to the uncertainty of the real-world applications, the actual return of each security usually cannot be prespecified in advance.In this case, how to effectively choose portfolio section strategy is a key problem for the investors in order to generate the least-risk plans and produce the expected return as much as possible.Along this line, Markowitz [1,2] first proposed the mean-variance models in stochastic environments, in which the variance is used to quantify the existing risks in the uncertain return.In this method, a tolerance threshold is usually given for the portfolios.A portfolio selection plan is a safe (or a low-risk) strategy if the variance of its corresponding random return is not greater than this threshold.Based on this approach, a variety of researches about the portfolio selection with either random parameters or fuzzy parameters have been proposed in the literature, such as Markowitz et al. [3], Gao et al. [4], Yi et al. [5], Xing et al. [6], H. Levy and M. Levy [7], Chiu and Wong [8], A. Palczewski and J. Palczewski [9], Castellano and Cerqueti [10], Fu et al. [11], and Zhang et al. [12].The second approach to represent the risk proposed by Markowitz [1,2] is the semivariance.In comparison to the mean-variance method, the mean-semivariance approach can better handle the risks in case of asymmetrical security return distributions.Along this line, interested readers can refer to Zhang et al. [13], Huang [14], Najafi and Mushakhian [15], Yan et al. [16], Yang et al. [17], and so forth.Besides the aforementioned two approaches, some effective chance-constrained methods can also be adopted to characterize the risks such as Huang [14], Huang and Zhao [18], and Li et al. [19].
Different from the variance, semivariance and chanceconstrained risk based methods, this paper aims to introduce the reliability term into the portfolio selection problem by using some reliability evaluation indexes, which can be used to optimize the low-risk portfolio strategies according to the real-world applications.We note that the current studies in the literature mainly focus on two specific uncertainties, that is, randomness and fuzziness, in which the returns of the involved securities are often assumed to be independent random variables or fuzzy variables (without correlations).Differently, this paper intends to propose a new representation method for the random return of each security on the basis of sample-based random data framework.This representation can allow for the correlations among returns of different securities in different periods.For instance, in each week (or month), we can collect realistic return data for individual securities, in which the week-specific or month-specific data can be regarded as having correlations among each other, and each week or month can be regarded as a considered stochastic sample with a specified probability.On the basis of these historical data, it is desirable for decision-makers to produce the reliable portfolio strategies to effectively reduce the investigation risk incurred by various uncertain factors.This research will particularly address this issue.To the best of our knowledge, no existing researches paid much attention to the data representation with inherent correlations.
As addressed above, the majority of existing studies always focus on the mean-variance (or semivariance) based approach to decrease the decision risks, in which the objective is often assumed to find the maximal expected return within the given risk threshold.However, in uncertain portfolio selection, the definition of an optimal portfolio may vary, since there are a large number of optimality criteria to measure the existing risks.Through adopting the samplebased random return to capture the randomness in the investigation process, this paper focuses on introducing three risk criteria into the portfolio selection problem to produce the risk-aversion planning, including the min-max criterion, percentile criterion, and expected disutility criterion, following the classical von Neumann and Morgenstern paradigm of decision under risk in economics [20].In particular, these risk-aversion criteria have been successfully applied to some real-world applications, for instance, management, economics, and transportation (see [21][22][23][24]).In particular, with the sample-base random data representation and these risk measures, it is possible for us to transform the formulated models into linear (or mix-integer) programming, which can be expected to solve through either commercial optimization software or the variants of current exiting algorithms.In addition, to effectively solve the proposed models, we in particular design a genetic algorithm with the polishing strategy to search for the near-optimal solutions.Numerical experiments show that the solution quality can be improved greatly in comparison to the traditional genetic algorithm without the polishing strategy, demonstrating the effectiveness of the proposed approaches.
The remainder of this paper is organized as follows.In Section 2, we formulate the problem of interest with three risk evaluation criteria, and the equivalent linear or mixedinteger programming models are also deduced with specific proofs.In Section 3, a genetic algorithm is designed based on the polishing strategy to solve the proposed models when the model cannot be transformed into linear form or solved by commercial optimization software.Section 4 implements different experiments to test performance of the proposed models and algorithms.In Section 5, a conclusion is made finally.

Formulation
Portfolio selection problem involves how to allocate the capital to a given number of securities so that the involved return can be maximized.Because of the uncertainty of the real-world investment environments, the returns of different securities are always set as independent uncertain variables over the entire decision process.This paper aims to handle the portfolio selection problem with a different input data representation, which is termed as sample-based random returns for different securities.Each sample corresponds to a distinct value (vector of values) that random return vector can take.By this method, it can be regarded that the elements in each sample vector are correlated with each other.Particularly, the historical return data over different periods (e.g., month, week) can be used as the input data in producing the least risk investment strategy.In the following, some relevant parameters and notations will be firstly given below.  : the probability of sample  ∈ Ω.

𝜉 𝜔
: the return of the th security over the sample .  : the return of the th security, which is denoted by a random variable.  : the proportion invested for security .
(ii) Constraints.Since   denotes the proportion invested for security , it is required that the sum of all the proportion should be a unity.Thus we need to consider the following constraints; that is, (iii) Objective Functions.In the following discussion, we will specify different types of objective functions in formulating the problem of interest.We have a total of  sets of historical data about the returns of different securities due to the portfolio activities in the past.In each set of historical data, we assume that all the returns of individual securities are correlated with each other.With the current set of historical data which may have a lot of uncertainties among different data set, we aim to find the most reliable decision-making plans to decrease the decision-making risks.
Typically, there are a lot of methods to handle the risks of the securities.Currently, we propose three types of methods to clearly handle the inherent uncertainties in the portfolio process.The first model is referred to as the max-min reliability model, detailed below.

Max-Min Reliable Model of the Portfolio Selection Problem.
According to the a priori information, we have a total of  sample data about the returns of different securities.With this concern, we can actually produce  types of total returns over different samples for each given portfolio selection strategy.Max-min model aims to find the most reliable portfolio plan across different samples such that the risk of decision-making can be decreased as much as possible.In detail, we denote the each sample-based total return by (,   ),  ∈ {1, 2, . . ., }, given below: To integrate these  returns, we first calculate the minimal value over different returns as follows: Then the max-min reliable model can be formulated by max min In this model, it is easy to see that the objective function aims to optimize the lower bound of the  returns over different samples.Typically, this is also a conservative decision-making model for the real-world applications.Figure 1 gives an illustration for the random return for different solutions.As shown, different solutions might correspond to individual return values for different samples.According to the max-min criterion, the least realization value will be used as the evaluation of this solution.In Figure 1, we consider two solutions  and   .Typically, since (  ) ≤ (), solution  is better than   in this criterion.
Clearly, model ( 4) is not a linear programming form for the portfolio selection problem.To handle this model effectively through existing commercial optimization software, we can reformulate this model through introducing an auxiliary variable   .Then, the max-min reliable model can be transformed into the following linear programming model.

Proposition 1. The max-min reliable model is equivalent to the following linear programming model:
≥ 0,  = 1, 2, . . ., ,   ≥ 0. ( Proof.Suppose that  = { 1 ,  2 , . . .,   } is the optimal solution for model (4).Then, let We thus have which implies (, ) is a feasible solution for model (5).On the other hand, let (,   ) be a feasible solution for model ( 5), and we then have which implies (, ) is the optimal solution for model (5).On the other hand, if (, ) is the optimal solution for model (5),  is typically a feasible solution for model (4).Let  be a feasible solution for model ( 4) and set Thus, (, ) is a feasible solution for model (5).We have min Then  is the optimal solution for model (4).The proof is completed.
Remark 2. It is easy to see from the max-min model that all the sample-based total returns have been handled equally, since the objective function has no relationship with the probability of each sample.In other words, this model is suitable for decision-makers with extreme risk-aversion.

Percentile Reliability Model of the Portfolio Selection Problem.
Next, we will introduce a percentile reliability model for the portfolio selection problem, which aims to generate portfolio strategies with varying probability confidence levels.We firstly give a definition of the critical value of random variables as follows.
Definition 3. Let  be a random variable and a probability confidence level  ∈ (0, 1].Then  critical value of  is defined by Remark 4. Assume that  is a discrete random variable with realizations  1 ,  2 , . . .,   and corresponding probabilities where As the return function (, ) is essentially a random variable for each given solution .We can use  critical value of (, ) as the evaluation index of the portfolio selection plan.In detail, we firstly formulate the following model: In this model, we aim to maximize the  critical value of the random return (, ).In practice, the optimal objective implies that the sample-based returns are larger than the optimal objective at least with a probability confidence level .In particular, parameter  is used to denote the riskaversion level in the practice.A risk-aversion decision-maker would like to set a relatively large parameter , while small  leads to the high-risk portfolio selection plans.
Figure 2 is given to show the illustration of random return.In detail, for any  ∈ (0, 1], we have the corresponding critical values of solution  and   , that is, (, ) sup () and (  , ) sup (), respectively.Since (, ) sup () ≤ (  , ) sup (), solution   is better than  in this reliability criterion.
As shown in the percentile reliability model, the form of the first constraint is typically not a linear constraint, which will probably lead to the increase of complexity in the solution process.Next, we aim to transform this constraint into a linear form through introducing two types of decision variables.

Proposition 5. The percentile reliability model is equivalent to the following linear programming model (𝑀 is a sufficiently large number):
Proof.Let   be the optimal solution for models (12).Then we have max    ≥ max   , ∀.
Proof.Since a large parameter  corresponds to a small feasible solution region, this result is obvious.Proposition 7. When  = 1, the percentile reliable model degenerates to the max-min model.

Expected Disutility Model of the Portfolio Selection Problem.
In the following, we focus on minimizing the expected disutility function associated with the total profits to produce the least-risk portfolio strategy, following the classical von Neumann and Morgenstern paradigm of decision under risk in economics [20].For statement simplicity, we next introduce the detailed formulation process of the disutility function in the modeling process.We first give a return target  to denote the upper bound of the historical return data.Then we calculate the following gap functions for the returns of different samples: Here, target  is required to be suitably chosen so as to ensure that the gap functions are all positive, which is the discussion focus of this research.Disutility function will be defined as an increasing (or nondecreasing at least) function over the return gaps that can be either linear or nonlinear.When we take identity mapping function, the return gap itself can be viewed as a special case of the disutility function.With this concern, decisionmakers are required to minimize expected disutility when they select the investment proportion of different securities.In this study, we will use an exponential disutility function () =   to measure the risk-aversion levels.Then we have the following disutility function for each gap function Thus, the expected disutility function can be written as follows: Practically, parameter  in the disutility function represents the level of risk-aversion in the decision-making process.Specifically, if one would like to choose a large parameter , then he/she should be a risk-aversion decision-maker.If parameter  is close to zero, the decision-maker will be riskcompromise.To clearly state the above idea in the portfolio selection process, we give the following illustrations.
Suppose we have two sample-based random returns, that is, 10 and 8 (unit: thousand dollars), with the probabilities 0.6 and 0.4, respectively.Supposing that return target is  = 12, we can obtain the following random gap distribution function (2, 0.6; 4, 0.4).Clearly, this distribution typically corresponds to the expected gap 2.8.We give Table 1 to show the corresponding mapping relationships between parameter  and its corresponding expected disutility.In Table 1, it is obvious that, with the increase of parameter , the corresponding expected disutility increases drastically.Thus, through minimizing the expected disutility, we can obtain the optimal portfolio selection plan with different risk levels.To further show the relationship between parameter  and risk-preference, we here consider a certain gap (, 1.0), which implies that the gap is  with the probability 1.0.We consider the same relationship between the random gap (2, 0.6; 4, 0.4) and certain gap (, 1.0) with respect to different parameter  (i.e., the gaps represented by (2, 0.6; 4, 0.4) and (, 1.0) correspond to the same evaluation value with each parameter ).
Obviously, as shown in Table 2, a larger parameter  will correspond to a poor value  in the deterministic case.In this sense, the risk-aversion decision-makers usually avoid setting the small values of parameter  to optimize a less risky decision plan.
In the following, the expected disutility model for the portfolio selection problem can be formulated as follows: min  [ ( (, ))] Note that the expected disutility function is essentially a nonlinear form in the corresponding optimization model, which might cause potential difficulties by traditional analytical methods in the solution process.In designing the heuristic, we can produce a lower bound model for model ( 22) by Jensen's inequality.Specifically, if function () is a convex function and  is a random variable, we then have [()] ≥ ([]), which leads to the following proposition.
For this model, some existing algorithms can be employed to solve an approximate optimal solution.However, it is difficult for the existing commercial optimization software to solve this model due to the nonlinear form of the objective function.With this concern, we here equivalently transform the nonlinear objective function into a linear formbased equation by introducing the logarithmic function ln(⋅).That is, we have Note that minimizing ([(, )]) is in essence equivalent to minimizing ln(([(, )])) over the region of feasible solutions.Then, we can use ln(([(, )])) as the alternative objective function in the solution process, and the exponential function can be imposed on the optimal objective value to obtain the optimal objective of model (24).Remark 9. Typically, the expected disutility model is a generalization of the expected value model, in which the expected gap can be regarded as an evaluation index of each portfolio selection strategy if we adopt the linear disutility function in formulating the objective function.

Genetic Algorithm
In particular, for the max-min and percentile reliability models, the commercial optimization software (e.g., Lingo and CPLEX) can be used to produce the exact (or near) optimal solutions for these models.For the expected disutility model, we can adopt the genetic algorithm-base heuristics to generate an approximate optimal solution because of the nonlinear objective function.Genetic algorithm is a kind of evolutionary algorithm proposed by Holland [25] in 1975, which can be used for seeking high-quality solutions of mathematical programming models.With the various technical details, the genetic algorithm has been applied to solving a variety of real-world problems, such as vehicle routing problem, transportation problem, and operations management (For more details, please see Aytug et al. [26], Chang and Sim [27], Yang et al. [24], Chung et al. [28], Xu et al. [29], etc.).In the following, we will introduce the detailed techniques in designing the genetic algorithm to solve the problem considered in this paper.
Solution Representation.In the solution algorithm, we use an -dimensional array to denote the chromosome of the proposed model, in which each element is randomly generated in a prespecified positive interval, denoted by [0, ].For clarity, we denote the chromosome by the following form: where   ,  = 1, 2, . . .,  are randomly extracted from interval [0, ].Note that this form is typically infeasible for the portfolio selection process.We then set up the following mapping relationship from a chromosome to a feasible solution: in which   =   / ∑  =1   ,  = 1, 2, . . ., .By this method, it is clear that each chromosome with form  corresponds to its unique feasible solution .Thus, in the following operations, each element of chromosome  should be controlled in interval [0, ], corresponding to a feasible solution.For instance, we have a total of 5 securities, and  is produced in the interval [0, 2] as  = (1, 2, 1, 2, 2).Then the corresponding feasible solution  is  = (0.125, 0.25, 0.125, 0.25, 0.125) . ( In the beginning of the algorithm, we need to produce a total of Pop size chromosomes in the population as the initial input data, denoted by  1 ,  2 , . . .,  Pop size .
Selection Operations.In the solution process, it is required to generate a new population for the selection operations.
In detail, we first rank the chromosomes in a good-to-bad sequence, which is also denoted by  1 ,  2 , . . .,  Pop size for notation simplicity.Then the evaluation value of each chromosome will be given by Eval (  ) =  (1 − ) −1 ,  = 1, 2, . . ., Pop size, (29) where  ∈ (0, 1) is a prespecified parameter in determining evaluation value of each chromosome.With these evaluation values, we select the chromosomes by using the roulette wheel.Specifically, we first produce a sequence by setting a sequence {  }

Pop size 𝑘=0
, where Then, we implement the following procedure for a total of Pop size times: randomly generate a number  in interval (0,  Pop size ); if there exists an index  such that  −1 ≤  ≤   , then chromosome   will be put into the new population for the following crossover and mutation operations.Thus, the newly generated population also contains a total of Pop size chromosomes even if some chromosomes are selected repeatedly.
Crossover Operations.Crossover operations aim to produce the new chromosomes for the population in order to find the approximate optimal solution as soon as possible.For this purpose, we firstly specify the chromosomes that take part in crossover operations.Note that this operation is performed on the basis of a predetermined crossover probability   .That is, the following procedure is employed to determine the crossover chromosomes: for each   , randomly generate a number  in interval (0, 1); if  ≤   , then   will be selected for crossover operations.Typically, a total of   ⋅ Pop size chromosomes can be expectedly chosen for the crossover operations.We perform the crossover operations based on each pair of selected chromosomes.Let   and   be two selected chromosomes.We firstly generate a crossover parameter  ∈ (0, 1), and the newly generated chromosomes   and   are given as follows: It is easy to see that since   and   are generated in interval [0, ],   and   should also lie in this interval, which leads to the feasibility of   and   .We then replace   and   by   and   , respectively, in the initial population.
Mutation Operations.Mutation operations intend to increase the diversity of the chromosomes in the population so as to avoid premature convergence.This operation is also performed under the consideration of the mutation probability   .Specifically, for each   , randomly generate a number  in interval (0, 1); if  ≤   , then   will be selected for mutation operations.In this process, a total of   ⋅ Pop size chromosomes will be expectedly chosen.For each chosen chromosome , we implement the mutation operations according to the following procedure: randomly generate a mutation vector  in interval [−1, 1]; choose a suitable step size  such that the newly produced chromosome   (i.e.,   =  + ) is feasible; replace  by   in the population.
Polishing Chromosomes.Note that, in this solution representation, each security is allowed to own an investment ratio for any feasible solution.Accordingly, in the near-optimal solution, it is possible that some securities have relatively small ratios (probably close to zero), which is typically undesirable in the real-world applications.To avoid this case, we herein give a ratio threshold for the each security in order to improve the solution quality.If the ratio of some security is less than this threshold, we will not consider investing this security in the solution.With this treatment, we can select an approximate optimal solution with guaranteed quality and desirability.
Procedure of Genetic Algorithm.Next, we will give the detailed procedure of the designed algorithm.
Step 1. Initialize Pop size chromosomes in interval [0, ] for the population.
Step 2. Compute the objective value of each chromosome in the population.
Step 3. Compute the evaluation value according to the goodto-bad sequence of chromosomes.
Step 4. Perform the selection operations over the population.
Step 5. Perform the crossover operations.
Step 6. Perform the mutation operations.
Step 7. Repeat Step 2 to Step 6 for a given number of times.
Step 8. Output the best solution encountered as the approximate optimal solution.

Numerical Examples
To test the performance of the proposed models, we will implement different sets of numerical experiments in the following discussion.All the experiments are performed on a personal computer with 4 GB memory and 1.60 GHz processors.
Example 1.In this set of experiments, we assume there are 10 securities.Ten samples are given in Table 3 to show the randomness of the security returns in the proposed models, listed in Table 3.
In this problem, we need to determine the optimal investment ratio of each security.Since the max-min reliability model is a special case of the percentile reliability model, only percentile model will be employed to test the performance of the proposed approaches.Note that the percentile model is equivalent to a linear programming model.In implementations, we can use the CPLEX solver in GAMS commercial software to solve the equivalent mixinteger linear programming models.
In this set of numerical experiments, we test the performance of the optimal objective function with respect to different probability confidence levels, that is,  = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0.In the GAMS software, we set the allowable relative gap parameter as OPTCR = 0.001.That is, when a solution is found within the relative gap 0.1% from the estimated optimal objective value, it will be outputted as a near-optimal solution for the model.With this parameter, all the experiments can output the corresponding optimal solution with zero relative gaps.The variation tendency of the optimal objectives is given in Figure 3, where the optimal objective takes decreasing tendency with the increase of the probability confidence level .This result can be explained suitably because a large probability confident level will lead to a small region of feasible solutions, corresponding to a smaller optimal objective, which also coincides with the results in Proposition 6.Note that when we set  = 1.0, the percentile reliable model outputs the same optimal objective to the max-min reliability model with the same optimal solution.
In the following, we aim to demonstrate the variation of different solutions with respect to varied probability confidence levels, which are given in Figure 4.It is easy to see that, for all the considered probability confidence  levels, the involved nonzero optimal decision variables are  2 ,  3 ,  4 ,  5 ,  6 ,  7 ,  8 , respectively.In this figure, the axis represents the probability confidence levels and -axis denotes the returned optimal solutions for each decision variable.Clearly, in this set of numerical experiments, the most active decision variables are  4 ,  5 , which have a total of five and eight times for investment among these ten experiments.When we set a relative small probability confidence level (e.g.,  ≤ 0.6), at most two securities can be selected for investment.For instance, for  = 0.1 or 0.3, only one security is adopted in the optimal solution, corresponding to  4 or  5 , respectively.When much larger probability confidence levels are adopted, more than three securities need to be invested.For example, when  = 0.9, securities 3, 4, 5, and 6 are needed to be invested with different optimal ratios.
Example 2. Next, the second set of experiments will be performed on the expected disutility model to show the performance of the proposed approaches.Since the expected disutility model is a nonlinear model, it is difficult to use commercial optimization software to solve it.With this concern, we adopt the proposed genetic algorithm to solve the approximate optimal solutions.
In this experiment, we take 20 securities into consideration for possible investment.The sample-based random returns are randomly generated in different intervals given in Table 4.As we need to produce a total of  samples in the problem, we also use the randomly generated sample probabilities as the occurrence chance of each sample.The following procedure is used to determine the probability of each sample.
Step 1. Randomly generate a number in interval [0, 1] for each sample , denoted by   .
In this set of experiments, two solution strategies are adopted to generate approximated optimal solutions, given below.
Strategy 1. Use the chromosomes representation without the polishing strategy.
Strategy 2. Use the chromosomes representation with the polishing strategy.
In this example, we randomly generate ten samples from the given intervals.For these two chromosome representation strategies, different experiments are performed with In these results, we calculate different gaps for individual parameters and strategies.Clearly, the first strategy in genetic algorithm has relatively robust characteristics as the returned best objective varies slightly in interval [1500, 1600].However, this strategy might return a practically undesirable solution since all the securities correspond to their nonzero ratios even with a very small value.This case is practically  undesirable.However, we can improve this situation through the polishing strategy, and the outputted results can be improved to a great extent.As shown, in comparison to Strategy 1, Strategy 2 can reduce the returned objective up to almost 80%.In other words, the returned near-optimal objectives in Strategy 1 are about 4-5 times of the results in Strategy 2, where the best near-optimal objective turns out to be 242.91.These computational results show that the polishing strategy is more effective than the solution Strategy 1 for each set of parameters.For comparison convenience, we also give Figure 5 to show the variation of different Gap 1, in which the largest gap (i.e., 84.70%) occurs in the sixth experiment (the -axis denotes the test indexes and -axis denotes the gaps).

Conclusions
Using the sample-based random data to capture uncertainties of the decision parameters, we developed three different models for the portfolio selection problem with stochastic characteristics of each security return, including max-min reliability model, percentile reliability model and expected disutility model.With our random data representation methods, the max-min and percentile reliability models could be transformed into their linear forms through introducing different auxiliary variables, which could be easily solved by the commercial optimization software.The expected disutility model was formulated based on the disutility function.A lower bound linear programming problem was also deduced for this model based on our random data representation.
To effectively solve the proposed models, we proposed a polishing strategy-based genetic algorithm to produce the approximated optimal solutions.The numerical examples are implemented to specify the detailed characteristics of the proposed models and algorithm.
It is worth mentioning that we propose three models with different decision-making criteria.As for the percentile reliability mode, if the decision-maker is risk-appetite, he/she can take a small parameter ; otherwise, a larger parameter should be considered.In extreme case, if we set  = 1, the percentile reliability model will degenerate to the max-min model reliability model.In this sense, the max-min reliability model is a special case of the percentile reliability model, and it is more suitable for decision-makers with extreme riskaversion.This situation also holds for the expected disutility model with the exponential disutility function.In general, we cannot determine which model is the best in the real-world applications, and the use of different models is closely related to the preferences of the decision-makers.
Further research can focus on the following two aspects.(1) Different uncertainties may occur for the real-world applications with either enough samples or no enough samples.Thus, the reliable portfolio selection model with other uncertainties can be a new topic in the further research.
(2) The proposed models can be easily generalized to the more complicated situations with variance or semivariance threshold constraints.The characteristics of these problems can also be investigated in future study.

Figure 1 :
Figure 1: An illustration of random return and solution comparison.

Figure 2 :
Figure 2: Critical value of random return and solution comparison.

Figure 3 :
Figure 3: Critical value of random return and solution comparison.

Figure 4 :
Figure 4: Variation of different decision variables in the experiments.

Figure 5 :
Figure 5: The variation of Gap 1 in different tests.

Table 1 :
Parameter  and its expected disutility value.

Table 2 :
Parameter  and value .

Table 3 :
Sample data for returns of different securities.

Table 4 :
The intervals for generating sample-based random return data.

Table 5 :
The computational results for different parameters.Test index     Gen. Strategy 1 Strategy 2 Gap 1