Journal of Applied Mathematics

Volume 2014 (2014), Article ID 192868, 12 pages

http://dx.doi.org/10.1155/2014/192868

Research Article

Application of Artificial Bee Colony Algorithm to Portfolio Adjustment Problem with Transaction Costs

Wei Chen, Hui Ma, Yiping Yang, and Mengrong Sun

School of Information, Capital University of Economics and Business, Beijing 100070, China

Received 7 February 2014; Accepted 18 May 2014; Published 15 June 2014

Academic Editor: X. Zhang

Copyright © 2014 Wei Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Compared with the conventional probabilistic mean-variance methodology, fuzzy number can better describe an uncertain environment with vagueness and ambiguity. In this paper, we discuss a portfolio adjusting problem under the assumption that the returns of risky assets are fuzzy numbers and there exist general transaction costs in portfolio adjusting process. In the proposed model, we take the first possibilistic moment about zero of a portfolio as the investment return and the second possibilistic moment about the possibilistic mean value of the portfolio as the investment risk. To solve the proposed model, a modified artificial bee colony (ABC) algorithm is developed for calculating the optimal portfolio adjusting strategy. Finally, a numerical example is given to illustrate the effectiveness of the proposed model and approach.

1. Introduction

The mean-variance methodology for the portfolio selection problem, proposed originally by Markowitz [1], has played an important role in the development of modern portfolio selection theory. In Markowitz’s pioneer work, he combined probability theory with optimization tool to study the investment behavior under uncertainty. The key principle of the mean-variance model is to take the expected return of a portfolio as the investment return and to take the variance of the expected return of a portfolio as the investment risk. Following Markowitz’s work, many scholars have studied portfolio selection models, such as Best and Hlouskova [2], Alexander and Baptista [3], Jacobs et al. [4], and Yu and Lee [5]. The basic assumption for using Markowitz’s model is that the future state of the assets can be correctly reflected by asset data in the past. However, this assumption can rarely be satisfied in the ever-changing real asset markets. With the introduction of fuzzy set theory [6] and possibility theory [7], a number of researchers began to employ these theories to study portfolio selection problems in a fuzzy environment. For example, Inuiguchi and Tanino [8] introduced a novel possibilistic programming approach to the portfolio selection problem based on the minimax regret criterion. Carlsson et al. [9] introduced a possibilistic approach to the selection of portfolios with highest utility score. Zhang and Nie [10] introduced the admissible efficient portfolio model under the assumption that the expected returns and risks of assets have admissible errors. Vercher et al. [11] presented a fuzzy downside risk approach for managing portfolio problems in the framework of risk-return trade-off using interval-valued expectations. Gupta et al. [12] applied multicriteria decision making via fuzzy mathematical programming to develop comprehensive models of asset portfolio optimization for the investors’ pursuing either of the aggressive or conservative strategies. Recently, Chen et al. [13] presented a possibilistic portfolio selection model with different interest rates for borrowing and lending. Tsaur [14] developed a fuzzy portfolio model that focuses on different investor risk attitudes. Zhang et al. [15] presented a possibilistic mean-semivariance-entropy model for multiperiod portfolio selection which takes into account four criteria, namely, return, risk, transaction cost, and diversification degree of portfolio. Barak et al. [16] proposed a mean-variance-skewness fuzzy portfolio model with cardinality constraint and by considering the fuzzy chance constraint to measure portfolio liquidity.

Transaction cost is an important factor considered by investors in financial markets. Most of the cases, investors usually start with an existing portfolio and the decisions are how to readjust to the changes in the security market. This adjustment entails both purchases and sales of securities along with transaction costs. Arnott and Wagner [17] found that ignoring transaction costs would result in inefficient portfolio. Recently, several studies have dealt with the portfolio selection problem under consideration of different types of transaction costs. For example, Mao [18], Morton and Pliska [19], Best and Hlouskova [20], and Lobo et al. [21] studied portfolio optimization problem with fixed transaction costs. Konno and Wijayanayake [22] discussed portfolio selection problem under concave transaction costs and minimal transaction unit constraints. Later, they [23] presented a portfolio selection model under nonconvex transaction costs and minimal transaction unit constraints. Zhang et al. [24] dealt with the portfolio selection problem with general transaction costs under the assumption that the returns of assets obey LR-type possibility distributions. However, when complex types of transaction costs are considered, it would be more difficult to obtain the efficient portfolios by using traditional optimization algorithms. Therefore, many scholars applied heuristic algorithm for complex portfolio optimization problems. For example, Chang et al. [25] used three heuristic algorithms based on genetic algorithm (GA), tabu search (TS), and simulated annealing (SA) for the standard M-V model including cardinality and quantity constraints. Crama and Schyns [26] applied SA to the solution of a complex portfolio selection model with realistic constraints. Fernández and Gómez [27] applied artificial neural networks (NN) to trace out the efficient frontier associated with the standard Markowitz mean-variance model with additional cardinality and bounding constraints. Yu et al. [28] proposed an improved radial basis function (RBF) network-based portfolio selection model to realize mean-variance-skewness trade-off. In 2009, Cura [29] utilized particle swarm optimization (PSO) to solve portfolio optimization problems. Later, Chen and Zhang [30] applied PSO algorithm to solve the admissible portfolio selection problem with transaction costs. Krink and Paterlini [31] developed an improved differential evolution (DE) for multiobjective portfolio optimization and compared the proposed algorithm with quadratic programming and the NSGA-II algorithm. Bermúdez et al. [32] introduced a multiobjective GA for cardinality constrained fuzzy portfolio selection.

The purpose of this paper is to discuss the portfolio adjusting problem with transaction costs under the assumptions that the uncertain returns of assets in financial markets are fuzzy numbers. We propose a possibilistic portfolio model with general transaction costs, in which the first possibilistic moment about zero of a portfolio is used to measure the investment return, and the second possibilistic moment about the possibilistic mean value of the portfolio is used to measure the investment risk. Moreover, we notice that the artificial bee colony (ABC) algorithm proposed by Karaboga [33] is a relatively new algorithm, and numerical comparisons demonstrated that the performance of ABC algorithm is competitive to other population-based algorithms such as genetic algorithm (GA), differential evolution (DE), particle swarm optimization (PSO), and evolutionary algorithm (EA) [34–36]. Due to its simplicity and ease of implementation, ABC algorithm has captured much attention and has been applied to solve many practical optimization problems, such as vehicle routing problem [37], image registration [38], and flowshop scheduling problem [39]. Based on the above discussion, in this paper, a modified ABC algorithm is developed to solve complex portfolio selection problem.

The rest of the paper is organized as follows. In Section 2, we introduce the definitions of the possibilistic moments and some properties. Section 3 presents a portfolio selection model with general transaction costs under the assumption that the returns of assets are trapezoidal fuzzy variables. To solve the proposed portfolio optimization problem, a modified ABC algorithm is described in Section 4. A numeral example is given to illustrate the effectiveness of the proposed models and algorithm in Section 5. Finally, some concluding remarks are given in Section 6.

2. Preliminaries

Let us introduce some definitions which are needed in the following section. A fuzzy number is a fuzzy set of the real line with a normal, fuzzy convex, and continuous membership function of bounder support. The family of fuzzy numbers is denoted by . Moreover, a function is said to be a weighting function if is a nonnegative, monotone increasing and satisfies the normalization condition .

Definition 1 (Fullér and Majlender [40]). Let be a fuzzy number and let be a weighted function. Then -weighted lower possibilistic and upper possibilistic mean values and possibilistic mean of are defined as follows:

Remark 2. The -weighted possibilistic mean of is the arithmetic mean of its -weighted lower and upper possibilistic mean values; that is,

Remark 3. If , then,

That is, the -weighted possibilistic mean value can be considered as a generalization of possibilistic mean value introduced by Carlsson and Fullér [41].

Definition 4 (Saeidifar and Pasha [42]). Let be a fuzzy number and let be a weighted function. Then the th weighted double possibilistic moment of fuzzy number about points , is defined as follows:

Specifically, if , then the th possibilistic moment about the possibilistic mean value of fuzzy number is

Definition 5 (Saeidifar and Pasha [42]). Let be a fuzzy number. Then the th possibilistic moment about zero of a fuzzy number is defined as follows:

Furthermore, Saeidifar and Pasha [42] defined the first possibilistic moment as the crisp possibilistic mean value and the second possibilistic moment as the possibilistic variance of fuzzy number , where

which are consistent with the definitions introduced by Carlsson and Fullér [41].

It is straightforward to show the following formula:

The variance of is defined as the expected value of the squared deviations between the arithmetic mean and the endpoints of its level sets. The variance is always positive and a measure of dispersion or spread of the fuzzy number. In the physical interpretation of the variance, it gives the moment of inertia of the mass distributed about the center of mass; also the variance gives information about the spread of variables around the mean value and it is a very important factor to find out the fluctuations in the observed values (for more see [41, 42]).

3. Portfolio Rebalancing Model

In this paper, we suppose that an investor starts with an existing portfolio and considers reallocate his/her wealth among risky assets. In order to describe conveniently, we use the following notations:

: the proportion invested in asset ; : the random return rate of asset ; : the covariance between assets and ; : the upper bound constraint on asset ; .

Following Markowitz’s idea, we quantify investment return by the expected value of a portfolio and risk by the variance. A rational investor may be interested in obtaining a certain average return from the portfolio at a minimum risk. Thus, a portfolio selection problem in the mean-variance context can be written as

where is the expected returns of asset and is the desired mean return of the portfolio.

By solving the above optimization problem continuously with a different each time, a set of efficient points is traced out. This efficient set is called the efficient frontier and is a curve that lies between the global minimum risk portfolio and the maximum return portfolio. In other words, the portfolio selection problem is to find all the efficient portfolios along this frontier.

In order to enrich the model, a weighting parameter is introduced to reflect different investor risk attitudes. With this new parameter, model (9) can be described as follows:

where the parameter can also be interpreted as the risk aversion factor of the investor satisfying . Obviously, the greater the factor is, the more risk aversion the investor has. Especially, when , the investor will be extremely conservative because in this case only the risk of his/her investment is considered and no attention is paid to the returns of his/her investment. Conversely, means that the investor is extremely aggressive to pursue the returns of his/her investment, completely ignoring the risk of investment.

Transaction costs are inevitably present in practical applications of portfolio selection and can be used to capture a number of costs such as brokerage fees, bid-ask spreads, taxes, or even fund loads [21]. The transaction cost associated with a portfolio is usually defined as the sum of individual transaction costs on each asset:

where is the individual cost on the th asset.

Therefore, the portfolio adjusting model with general transaction cost can be expressed as

As discussed in Section 1, the returns of risky assets are in a fuzzy uncertain economic environment and vary from time to time; the future states of returns and risks of risky assets cannot be predicted accurately. Meanwhile, investor can take into account and complete his knowledge with other pieces of information, such as economical and financial behaviors of the companies, government policies, and business strategies. This information can be estimated quantitatively by the expert perceptions. Based on these factors, we consider the portfolio selection problem (12) under the assumption that the returns of assets are trapezoid fuzzy numbers.

Let be a trapezoidal fuzzy number with tolerance interval , left width , and right width , ; that is, can be described with the following membership function:

Then, a  -level sets of can be computed as

Using Definitions 1 and 4 introduced in Section 2, we easily obtain the crisp possibilistic mean value and the possibilistic variance of as follows:

respectively.

Furthermore, the crisp possibilistic mean value of the return associated with the portfolio after paying transaction costs is easily obtained by

and the possibilistic variance of return associated with the portfolio is given by

Analogous to Markowitz’s mean-variance methodology for the portfolio selection problem, the crisp possibilistic mean value corresponds to the return, while the possibilistic variance corresponds to the risk. Thus, the possibilistic portfolio adjusting problem with transaction costs can be formulated as

It should be noted that, in this paper, we consider four kinds of transaction costs, respectively, that is, concave transaction cost, linear transaction cost, convex transaction cost, and a general type of transaction cost, which are shown as follows.

(i)The concave transaction cost function is

(ii)The fixed proportional (linear) transaction cost function is

(iii)The convex transaction cost function is

(iv)The no-convex-no-concave transaction cost function is

where and .

Remark 6. In the general type of transaction cost function , the unit transaction cost is relatively larger when the amount of transaction is smaller and it decreases as the amount increases. Hence it is a nondecreasing concave function up to point. Because of the market impact effect, the transaction cost becomes convex when the amount exceeds point. Between point and point, the transaction cost is fixed proportional to the amount. It can be used to describe the practical situation for the transaction cost precisely.

4. Modified Artificial Bee Colony Algorithm

4.1. Artificial Bee Colony Algorithm

Artificial bee colony (ABC) algorithm is a relatively new algorithm developed by Karaboga [33], which is based on simulating the intelligent forging behavior of honey bee swarm. In ABC algorithm, the colony of artificial bees consists of three groups of bees: employed bees, onlooker bees, and scout bees. Employed bees are responsible for exploiting the nectar sources explored before and they give information to the other waiting bees in the hive about the quality of the food source which they are exploiting. Onlooker bees wait in the hive and establish food source to exploit depending on the information shared by the employed bees. Scouts search environment in order to find a new food source.

In ABC algorithm, each food source is exploited by only one employed bee. In other words, the number of employed bees is equal to the number of food sources around the hive. Moreover, the position of a food source represents a possible solution to the optimization problem and the nectar amount of a food source corresponds to the quality (fitness) of the associated solution. The number of the employed bees or the onlooker bees is equal to the number of solutions in the population.

Initial Population. At the first step, the ABC generates a randomly distributed initial population of solutions (food source positions), where denotes the size of population. Each solution is a -dimensional vector. Here, is the number of optimization parameters.

After initialization, the population of the positions (solutions) is subject to repeated cycles, , where is the maximum cycle number. The search processes of employed bees, onlooker bees, and scout bees are described as follows.

The Employed Bees Phase. An employed bee produces a modification on the position (solution) in her memory depending on the local information (visual information) and tests the nectar amount (fitness value) of the new source (new solution). Provided that the nectar amount of the new one is higher than that of the previous one, the bee memorizes the new position and forgets the old one. Otherwise she keeps the position of the previous one in her memory. The expression, which is used by an employed bee to produce a modification, is given as follows:

where and are randomly chosen indexes. Although is determined randomly, it has to be different from . is a random number between .

The Onlooker Bees Phase. After all employed bees complete the search process, they share the nectar information of the food sources and their position information with the onlooker bees on the dance area. Then, an onlooker bee evaluates the nectar information taken from all employed bees and chooses a food source depending on the probability, , associated with that food source:

where is the fitness value of the solution which is proportional to the nectar amount of the food source in the position and is the number of food sources which is equal to the number of employed bees.

The Scout Bees Phase. Any food source position that does not improve the fitness value will be abandoned and replaced by a new position that is randomly determined by a scout bee. This helps avoid suboptimal solutions. The new random position chosen by the scout bee will be calculated from

where is determined randomly. It should be noticed that it has to be different from .

Based on the above discussions, we formally describe the main steps of the standard ABC algorithm, given in Algorithm 1.

(1) Initialize the population of solutions , ,

(2) Evaluate the population

(3) cycle = 1

(4) Repeat

(5) Produce new solutions for the employed bees by using (24) and evaluate them

(6)  Apply the greedy selection process for the employed bees

(7)  Calculate the probability values for the solutions by (25)

(8)  Produce the new solutions for the onlookers from the solutions selected depending on and evaluate them

(9)  Apply the greedy selection process for the onlooker bees

(10)  Determine the abandoned solution for the scout, if exists, and replace it with a new

        randomly produced solution by (26)

(11)   Memorize the best solution achieved so far

(12)  cycle = cycle + 1

(13) Until cycle = MCN

Algorithm 1: Standard ABC algorithm.

4.2. Modified Artificial Bee Colony Algorithm

Although some researchers [34–36] have demonstrated that the performance of the ABC algorithm is better than or similar to those of other population-based algorithms with the advantage of employing fewer control parameters, ABC also faces up some insufficiencies. For example, ABC can get trapped in local optima when solving complex multimodel function optimization problems [36]. Recently, a few modified or improved algorithms based on the classical ABC algorithm are proposed. For example, Akay and Karaboga [43] introduced a modified ABC algorithm and applied it for efficiently solving real-parameter optimization problems. Alatas [44] proposed chaotic bee colony algorithms for global numerical optimization. Xiang and An [45] proposed an efficient and robust artificial bee colony algorithm for numerical optimization. Kalayci and Gupta [46] proposed an ABC algorithm for a sequence-dependent disassembly line balancing problem with multiple objectives and compared it with GA, PSO, SA, TS, river formation dynamics (RFD), and ant colony optimization (ACO). Tsai [47] integrated ABC algorithm and bees' algorithms (BA) into a hybrid ABC-BA algorithm and conducted experiments on six constrained optimization problems with equality or inequality constraints. These modified or improved ABC algorithms have shown a better performance than the classical ABC algorithm. In view of the above, considering the complexity of the portfolio optimization proposed, we present an improved ABC algorithm based on chaotic theory to solve the problem.

Chaotic Initialization. Population initialization is a crucial task in evolutionary algorithms. If no information about the solution is available, then random initialization is the most commonly used method to initialize the population. However, it may affect the convergence speed and the quality of the final solution. Having the characteristics of the certainty, ergodicity, and randomicity, the chaotic map is suitable to initialize the population for the purpose of increasing the population diversity and achieving high-quality solution. The logistic map, which is well known, is given as follows:

where is a control parameter, is a chaotic variable, and . Obviously, under the conditions that the initial . Especially, behaves chaotic dynamics when and . The pseudocode of proposed chaotic initialization is given in Algorithm 2.

(1) Set population size SN and dimension of solution

     //Randomly generate the initial value of chaotic variable

(2) for   to   do

(3)    

(4)    while     do

(5)      

(6)    end while

(7) end for

     // The iterative process of chaotic initialization

(8) for   to SN  do

(9)    for   to   do

(10)      

(11)      

(12)   end for

(13) end for

Algorithm 2: Chaotic initialization.

A Modified Search Equation. The employed bee and onlooker bee produce a modification on the position in their memory by using (24). However, the convergence performance is poor in some cases. Therefore, here we proposed a modified search equation, in which a forgetting factor and a neighborhood factor are considered. It is given as follows:

where is the forgetting factor, which expresses the memory strength for current food source when searching the next food source. It dynamically decreases with the increase of iterations. In addition, the differences of current food source and neighbourhood food source also affect the convergence speed. Thus, the neighborhood factor is introduced to accelerate the convergence speed by adjusting the radius of the search for new candidates. The parameters and dynamically change as follows:

where the value of and represent the maximum and minimum percentage of the position adjustment for employed bees and onlooker bees. Therefore, with these selected values, the value of linearly decreases and the value of linearly increases.

Chaotic Search for Scout Bees. To avoid the phenomenon of stagnation and premature convergence in the standard ABC algorithm, we employed chaotic search technique to get out of local optimum and get better results. Therefore, the chaotic search technique for scout bees is specifically illustrated in Algorithm 3.

(1) Set the maximum number of chaotic iteration and iterative variable

    // The variable denotes the th chaotic variable in the th iteration

(2) for   to   do

(3)    

(4) end for

     // The iterative process of chaotic search

(5) while     do

(6)    for   to   do

(7)       

            // represents a candidate solution

(8)      

(9)    end for

         //Greedy selection is applied between

(10)   if  

(11)      Replace solution with candidate solution

(12)      Set

(13)      Break

(14)   else

(15)      Set

(16)   end if

(17)   Set

(18) end while

Algorithm 3: Chaotic search for scout bees.

The Modified Algorithm. According to the analysis of the standard ABC algorithm and the modification mentioned above, it is clear that the modified ABC algorithm can be well balanced between the exploration and exploitation. Naturally, the modified ABC algorithm is given in Algorithm 4.

(1) Set population size SN, the number of maximum cycles MCN and the control parameter

(2) Perform Algorithm 2 to fulfill the chaotic initialization of population , and calculate

     their objective function values and fitness value values. Set and

(3) while  iteration   do

        // The employed bees phase

(4)   for   to SN  do

(5)       Produce a new candidate food source corresponding to food source using (24)

(6)       if   ()   Then  

(7)       else  

(8)       end if

(9)   end for

(10)  Calculate the fitness values of all food sources and the probability values by using (15)

        // The onlooker bees phase

(11)  Set

(12)   while     do

(13)       if     then

(14)         Set

(15)         Produce a new candidate food source for the onlooker bee corresponding to food source using (24)

(16)         if     then

(17)             

(18)         else  

(19)         end if

(20)       end if

(21)      Set

(22)      if     then  

(23)      end if

(24)   end while

        //The scout bees phase

(25)   for   to SN  do

(26)       if     then

(27)        Perform Algorithm 3 to implement chaotic search

(28)       end if

(29)   end for

(30)  Set

(31) end while

Algorithm 4: Modified artificial bee colony algorithm.

5. Numerical Example

In order to illustrate our proposed effective means and approaches for the portfolio selection problem in this paper, we give a numerical example. Consider a 4-security problem with the following possibility distributions.

Now we use the possibilistic portfolio selection model (19) proposed in this study to determine the future optimal portfolio selection strategy. The rate of transaction costs for 4 securities is , the upper bound , and , , . Moreover, we let the maximum cycle , the size of the population , and the control parameter (Table 1).

Table 1: Possibility distributions of 4 securities.

0.05
0.06
0.01
0.01

0.06
0.065
0.02
0.005

0.068
0.075
0.02
0.005

0.065
0.07
0.015
0.03

To demonstrate the effects of four kinds of transaction costs on the portfolio selection, we let the risk tolerance parameter , , and , respectively. Some optimal results obtained by using modified ABC algorithm for model (19) are shown in Tables 2, 3, and 4. It should be noticed that is the objective value and is the total transaction costs.

Table 2: Comparison results with different transaction costs with .

0
0.2
0.381
0.298

0.2
0.8
0.261
0.298

0.8
0
0.223
0.203

0
0
0.135
0.201

0.0309
0.0225
0.0168
0.0239

TC
0.0328
0.0340
0.0430
0.0357

Table 3: Comparison results with different transaction costs with .

0.04
0.2
0.350
0.291

0.18
0.8
0.305
0.308

0.78
0
0.223
0.201

0
0
0.122
0.200

0.0171
0.0125
0.0093
0.0133

TC
0.0328
0.0340
0.0429
0.0357

Table 4: Comparison results with different transaction costs with .

0.200
0.193
0.372
0.272

0
0.800
0.304
0.326

0.8
0
0.201
0.203

0
0.007
0.123
0.199

0.0034
0.0025
0.0018
0.0026

TC
0.0319
0.0341
0.0429
0.0358

As is shown in Tables 2, 3, and 4, the total transaction cost under concave function is the smallest among all the transaction cost functions and the biggest under convex function . This is because that concave transaction cost has a decreasing slope, which means that the more the amount of asset is allocated, the less transaction cost increases, while the increasing derivative of convex transaction cost will induce much more transaction cost if the holding of asset becomes more. Moreover, for concave transaction cost function and the linear transaction cost function , the portfolio proportions are concentrated on certain securities. Oppositely, the convex transaction cost function and the no-convex-no-concave transaction cost function can decentralize portfolio proportions, and do better than . This means that the general type of transaction cost can better avoid the amount of each asset being too high or too low. For example, Table 3 shows that the holding of asset 4 for transaction cost reaches to its lower bound with and the total transaction cost is the smallest. What is more, for linear transaction cost , the holding of assets 3 and 4 reaches their lower bound with and the total transaction cost is bigger than that with . Comparing with and , the portfolio proportions are decentralized and its total transaction cost is the biggest. Furthermore, for , the portfolio proportions are further decentralized with proportions , , , and .

To show the effectiveness of the proposed algorithm, modified ABC algorithm, ABC algorithm, and GA are compared for . Results are listed in Tables 5, 6, 7, and 8. We can see that under the same transaction cost function, the return obtained by modified ABC algorithm is the biggest among three algorithms and the risk obtained by modified ABC algorithm is the smallest. It illustrates that the modified ABC algorithm is competitive to ABC algorithm and GA.

Table 5: Comparison results of modified ABC, ABC, and GA with .

Return
Risk

Modified ABC
0.040
0.180
0.780
0
0.06721

0.0154

ABC
0.002
0.190
0.800
0.008
0.06719

0.0153

GA
0.023
0.172
0.800
0.005
0.06712

0.0153

Table 6: Comparison results of modified ABC, ABC, and GA with .

Return
Risk

Modified ABC
0.200
0.800
0
0
0.059

0.0125

ABC
0.200
0.800
0
0
0.059

0.0125

GA
0.274
0.725
0
0.001
0.058

0.0112

Table 7: Comparison results of modified ABC, ABC, and GA with .

Return
Risk

Modified ABC
0.339
0.306
0.221
0.134
0.06163

0.0084

ABC
0.342
0.300
0.223
0.135
0.06154

0.0083

GA
0.343
0.300
0.217
0.139
0.06160

0.0083

Table 8: Comparison results of modified ABC, ABC, and GA with .

Return
Risk

Modified ABC
0.274
0.326
0.201
0.199
0.06242

0.0120

ABC
0.254
0.346
0.200
0.200
0.06222

0.0119

GA
0.307
0.291
0.202
0.200
0.06228

0.0119

Furthermore, in order to test the robustness of the proposed algorithm, we solve the problem by setting the different parameters in the modified ABC algorithm. We employ the relative error which is defined by (maximal objective − actual objective)/maximal objective , where the maximal objective is the maximum of all the computational results obtained and the actual objective is . Here we choose and the no-convex-no-concave transaction cost function . The detailed results are shown in Table 9. The results show that all the relative errors are less than , which means that the modified ABC algorithm has strong robustness to parameter variations and resistance disturbances.

Table 9: Relative error with .

SN
MCN
Limit

Relative error

10
1000
100
0.0159515
0.0357

20
1000
50
0.0159450
0.0765

20
1500
100
0.0159567
0.0031

20
2000
150
0.0159552
0.0125

30
1500
100
0.0159472
0.0627

30
1000
150
0.0159482
0.0564

40
1500
200
0.0159537
0.0219

50
2000
100
0.0159556
0.0100

60
1000
50
0.0159436
0.0852

70
1500
100
0.0159572
0

Finally, in order to show that the modified ABC algorithm is more effective and has a faster convergence speed than ABC algorithm and GA, the comparison of convergence characteristic is given in Figure 1. It shows that the modified ABC algorithm converges within 500 iterations, while ABC algorithm and GA converge within 1000 iterations. What is more, the global optimal value obtained by modified ABC algorithm is the biggest. That is to say, the modified ABC algorithm converges to the global optimal solution fastest and is the most effective algorithm among the three algorithms to solve complex portfolio problems.

Figure 1: Performance comparison of the modified ABC algorithm, ABC algorithm, and GA.

6. Conclusion

Based on the possibilistic moments, this paper deals with a fuzzy portfolio selection model under four kinds of transaction costs. Then modified ABC algorithm is applied to solve the optimal problem. Furthermore, we get the portfolio proportions for 4 securities under four kinds of transaction costs. The results show that the portfolio proportions change with different types of transaction costs. Finally, we compared the results of modified ABC algorithm with that of ABC algorithm and GA, which showed that the modified ABC algorithm is better than ABC algorithm and GA in solving fuzzy portfolio selection problem.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors thank Professor Zhijun Wu for his valuable comments and suggestions. Moreover, this research was supported by the Humanity and Social Science Youth Foundation of the Ministry of Education of China (no. 13YJC630012) and the Training Programme Foundation for the Beijing Municipal Excellent Talents (no. 2013D005019000007). The second author acknowledges the support of Beijing Natural Science Foundation (no. 9122003) and Beijing Philosophical Social Science Project (no. 11JGB077). The third author acknowledges the support of the National Natural Science Foundation (no. 71240002). The fourth author acknowledges the support of Beijing Natural Science Foundation (no. 9142003).

References

H. Markowitz, “Portfolio selection,” Journal of Finance, vol. 7, pp. 77–91, 1952.

M. J. Best and J. Hlouskova, “The efficient frontier for bounded assets,” Mathematical Methods of Operations Research, vol. 52, no. 2, pp. 195–212, 2000.

G. J. Alexander and A. M. Baptista, “A comparison of VaR and CVaR constraints on portfolio selection with the mean-variance model,” Management Science, vol. 50, no. 9, pp. 1261–1273, 2004.

B. I. Jacobs, K. N. Levy, and H. M. Markowitz, “Portfolio optimization with factors, scenarios, and realistic short positions,” Operations Research, vol. 53, no. 4, pp. 586–599, 2005.

J. R. Yu and W. Y. Lee, “Portfolio rebalancing model using multiple criteria,” European Journal of Operational Research, vol. 209, no. 2, pp. 166–175, 2011.

L. A. Zadeh, “Fuzzy sets,” Information and Control, vol. 8, no. 3, pp. 338–353, 1965.

L. A. Zadeh, “Fuzzy sets as a basis for a theory of possibility,” Fuzzy Sets and Systems, vol. 1, no. 1, pp. 3–28, 1978.

M. Inuiguchi and T. Tanino, “Portfolio selection under independent possibilistic information,” Fuzzy Sets and Systems, vol. 115, no. 1, pp. 83–92, 2000.

C. Carlsson, R. Fullér, and P. Majlender, “A possibilistic approach to selecting portfolios with highest utility score,” Fuzzy Sets and Systems, vol. 131, no. 1, pp. 13–21, 2002.

W. G. Zhang and Z. K. Nie, “On admissible efficient portfolio selection problem,” Applied Mathematics and Computation, vol. 159, no. 2, pp. 357–371, 2004.

E. Vercher, J. D. Bermúdez, and J. V. Segura, “Fuzzy portfolio optimization under downside risk measures,” Fuzzy Sets and Systems, vol. 158, no. 7, pp. 769–782, 2007.

P. Gupta, M. K. Mehlawat, and A. Saxena, “Asset portfolio optimization using fuzzy mathematical programming,” Information Sciences, vol. 178, no. 6, pp. 1734–1755, 2008.

W. Chen, Y. Yang, and H. Ma, “Fuzzy portfolio selection problem with different borrowing and lending rates,” Mathematical Problems in Engineering, vol. 2011, Article ID 263240, 15 pages, 2011.

R. C. Tsaur, “Fuzzy portfolio model with different investor risk attitudes,” European Journal of Operational Research, vol. 227, no. 2, pp. 385–390, 2013.

W. G. Zhang, Y. J. Liu, and W. J. Xu, “A possibilistic mean-semivariance-entropy model for multi-period portfolio selection with transaction costs,” European Journal of Operational Research, vol. 222, no. 2, pp. 341–349, 2012.

S. Barak, M. Abessi, and M. Modarres, “Fuzzy turnover rate chance constraints portfolio model,” European Journal of Operational Research, vol. 228, no. 1, pp. 141–147, 2013.

R. D. Arnott and W. H. Wanger, “The measurement and control of trading costs,” Financial Analysts Journal, vol. 46, no. 6, pp. 73–80, 1990.

J. C. T. Mao, “Essentials of portfolio diversification strategy,” Journal of Finance, vol. 25, pp. 1109–1121, 1970.

A. J. Morton and S. R. Pliska, “Optimal portfolio management with transaction costs,” Mathematical Finance, vol. 5, pp. 337–356, 1995.

M. J. Best and J. Hlouskova, “Portfolio selection and transactions costs,” Computational Optimization and Applications, vol. 24, no. 1, pp. 95–116, 2003.

M. S. Lobo, M. Fazel, and S. Boyd, “Portfolio optimization with linear and fixed transaction costs,” Annals of Operations Research, vol. 152, no. 1, pp. 341–365, 2007.

H. Konno and A. Wijayanayake, “Portfolio optimization problem under concave transaction costs and minimal transaction unit constraints,” Mathematical Programming B, vol. 89, no. 2, pp. 233–250, 2001.

H. Konno and A. Wijayanayake, “Portfolio optimization under D.C. transaction costs and minimal transaction unit constraints,” Journal of Global Optimization, vol. 22, no. 2, pp. 137–154, 2002.

X. L. Zhang, W. G. Zhang, W. J. Xu, and W. L. Xiao, “Possibilistic approaches to portfolio selection problem with general transaction costs and a CLPSO algorithm,” Computational Economics, vol. 36, no. 3, pp. 191–200, 2010.

T. J. Chang, N. Meade, J. E. Beasley, and Y. M. Sharaiha, “Heuristics for cardinality constrained portfolio optimisation,” Computers and Operations Research, vol. 27, no. 13, pp. 1271–1302, 2000.

Y. Crama and M. Schyns, “Simulated annealing for complex portfolio selection problems,” European Journal of Operational Research, vol. 150, no. 3, pp. 546–571, 2003.

A. Fernández and S. Gómez, “Portfolio selection using neural networks,” Computers and Operations Research, vol. 34, pp. 1177–1191, 2007.

L. Yu, S. Wang, and K. K. Lai, “Neural network-based mean-variance-skewness model for portfolio selection,” Computers and Operations Research, vol. 35, no. 1, pp. 34–46, 2008.

T. Cura, “Particle swarm optimization approach to portfolio optimization,” Nonlinear Analysis: Real World Applications, vol. 10, no. 4, pp. 2396–2406, 2009.

W. Chen and W. G. Zhang, “The admissible portfolio selection problem with transaction costs and an improved PSO algorithm,” Physica A: Statistical Mechanics and its Applications, vol. 389, no. 10, pp. 2070–2076, 2010.

T. Krink and S. Paterlini, “Multiobjective optimization using differential evolution for real-world portfolio optimization,” Computational Management Science, vol. 8, no. 1-2, pp. 157–179, 2011.

J. D. Bermúdez, J. V. Segura, and E. Vercher, “A multi-objective genetic algorithm for cardinality constrained fuzzy portfolio selection,” Fuzzy Sets and Systems, vol. 188, no. 1, pp. 16–26, 2012.

D. Karaboga, “An idea based on honeybee swarm for numerical optimization,” Tech. Rep. TR06, Erciyes University, Engineering Faculty, Computer Engineering Department, 2005.

D. Karaboga and B. Basturk, “A powerful and efficient algorithm for numerical function optimization: artificial bee colony (ABC) algorithm,” Journal of Global Optimization, vol. 39, no. 3, pp. 459–471, 2007.

D. Karaboga and B. Basturk, “On the performance of artificial bee colony (ABC) algorithm,” Applied Soft Computing Journal, vol. 8, no. 1, pp. 687–697, 2008.

D. Karaboga and B. Akay, “A comparative study of Artificial Bee Colony algorithm,” Applied Mathematics and Computation, vol. 214, no. 1, pp. 108–132, 2009.

W. Y. Szeto, Y. Wu, and S. C. Ho, “An artificial bee colony algorithm for the capacitated vehicle routing problem,” European Journal of Operational Research, vol. 215, no. 1, pp. 126–135, 2011.

A. Banharnsakun, T. Achalakul, and B. Sirinaovakul, “The best-so-far selection in Artificial Bee Colony algorithm,” Applied Soft Computing Journal, vol. 11, no. 2, pp. 2888–2901, 2011.

X. Wang, X. Xie, and T. C. E. Cheng, “A modified artificial bee colony algorithm for order acceptance in two-machine flow shops,” International Journal of Production Economics, vol. 141, no. 1, pp. 14–23, 2013.

R. Fullér and P. Majlender, “On weighted possibilistic mean and variance of fuzzy numbers,” Fuzzy Sets and Systems, vol. 136, no. 3, pp. 363–374, 2003.

C. Carlsson and R. Fullér, “On possibilistic mean value and variance of fuzzy numbers,” Fuzzy Sets and Systems, vol. 122, no. 2, pp. 315–326, 2001.

A. Saeidifar and E. Pasha, “The possibilistic moments of fuzzy numbers and their applications,” Journal of Computational and Applied Mathematics, vol. 223, no. 2, pp. 1028–1042, 2009.

B. Akay and D. Karaboga, “A modified Artificial Bee Colony algorithm for real-parameter optimization,” Information Sciences, vol. 192, pp. 120–142, 2012.

B. Alatas, “Chaotic bee colony algorithms for global numerical optimization,” Expert Systems with Applications, vol. 37, no. 8, pp. 5682–5687, 2010.

W. L. Xiang and M. Q. An, “An efficient and robust artificial bee colony algorithm for numerical optimization,” Computers and Operations Research, vol. 40, no. 5, pp. 1256–1265, 2013.

C. B. Kalayci and S. M. Gupta, “Artificial bee colony algorithm for solving sequence-dependent disassembly line balancing problem,” Expert Systems with Applications, vol. 40, no. 18, pp. 7231–7241, 2013.

H. C. Tsai, “Integrating the artificial bee colony and bees algorithm to face constrained optimization problems,” Information Sciences, vol. 258, pp. 80–93, 2014.