Intelligent Black–Litterman Portfolio Optimization Using a Decomposition-Based Multi-Objective DIRECT Algorithm

Abstract

It is agreed that portfolio optimization is of great importance for the financial market. However, input sensitivity and highly-concentrated portfolios have posed a challenge. In this paper, a random forest-based Black–Litterman model is developed, aiming to further enhance the portfolio performance, which adopts a novel method for generating investor views on the basis of random forests. More specifically, the view vector is generated based on the predicted asset returns obtained by random forests, and the confidence matrix which contains the uncertainty of each view is measured by the difference in the predicted values of multiple trees. Furthermore, motivated by decomposition strategy, a novel multi-objective DIRECT algorithm is introduced to effectively resolve the proposed model. Through the construction of a unique indicator, the algorithm possesses the capacity to select potentially-optimal hyperrectangles in all reference directions simultaneously, which will further improve the exploratory nature. Experimental results have demonstrated that the proposed algorithm achieves a better performance over NSGA-II and MOEA/D on the MOP and DTLZ benchmark problems. It is also experimentally verified that the random forest-based Black–Litterman model can obtain higher cumulative returns and Sharpe ratios in the application of Chinese stock markets when compared to the classic MV model.

1. Introduction

Portfolio optimization is the art of asset allocation to achieve the purpose of maximizing investment returns while matching investors’ risk tolerance. As an essential part of practical portfolio management, it has become an attracting issue in both academia and industry during the past few decades. In 1952, Markowitz proposed the Mean-Variance (MV) portfolio optimization model [1], which is regarded as the cornerstone of modern portfolio theory (MPT). Although the MV model has drawn much attention since its introduction, several obvious shortcomings are exposed gradually with the application of the model. The MV model describes the return and the risk of the portfolio by the mean and the variance of the asset returns, respectively [2,3], which means that the forecasting of asset returns will strongly influence portfolio performance. In practice, such forecasts are mainly based on their historical behavior, which may lead to considerable estimation errors. To address this problem, an enormous amount of attempts have been made.

The Black–Litterman (BL) model [4,5] is one of these attempts, which takes views on certain assets into consideration based on a Bayesian approach. Specifically, experts or investors give their subjective opinions on certain assets, and the BL model integrates them with the historical return data to create the future asset returns. Compared with the MV model, it is more stable in the process of forecasting asset returns and covariances, which helps to extend profit advantages and avoid risks, especially in volatile markets [6]. In the past decades, many researchers devoted themselves to the extension of BL model [7,8,9,10,11], in which novel methods to generate views is one of the most promising directions. With the development of artificial intelligence techniques, machine learning algorithms began to be employed for views generating [12,13,14,15] to replace the views given directly through expert experience. However, there is almost no study that places emphasis on the uncertainty introduced by the relevant predictive model itself.

BL portfolio construction is on the basis of the classic return-risk optimization, which is a typical multi-objective optimization problem (MOP). For these kinds of problems, evolutionary multi-objective optimization (EMO) algorithms such as MOEA/D [16] and NSGA-II [17] have fully demonstrated their niche [18]. However, EMO algorithms usually contain random operations in the process of the population initialization and evolution to maintain diversity, which leads to non-reproducibility, and lack a theoretical guarantee of global convergence. From the perspective of practice, multi-objective optimization is expected to be deterministic and theoretically convergent. Inspired by Lipschitzian global optimization [19,20], Jones et al. [21] introduced the DIviding RECTangles (DIRECT) algorithm, which does not require the Lipshitz constant and is more tractable in higher dimensions [22]. Due to the advantages of being easy to implement and having fewer hyperparameters (only one), some researchers explored extensions of DIRECT to make it possible to handle MOPs. However, most of them usually have a high computational complexity. Therefore, some algorithmic changes to the existing algorithms are somehow more practical.

This paper attempts to address a BL portfolio optimization problem in which particular attention is paid to the interesting issue of how to generate investor preferences systematically and precisely. With this goal, a novel random forest-based BL model is developed, and an approach to handle the model mentioned above in a short time, while maintaining high accuracy, is also proposed. The main contribution of this paper could be summarized as follows:

• We develop a new random forest-based BL portfolio optimization model in which a novel method for generating investor views on the basis of random forests is adopted. In this method, a view vector is generated based on the predicted asset returns obtained by random forests, and the confidence matrix which contains the uncertainty of each view is measured by the differences in the predicted values of multiple trees.
• We propose a decomposition-based multi-objective DIRECT algorithm, named multiDecompose, to handle the random forest-based BL model mentioned above in a short time while maintaining high accuracy, in which an indicator is explored to encourage potentially optimal hyperrectangles to be chosen in all directions. With this indicator, the algorithm can provide a better exploration of the search space of an MOP, especially when stuck in a local optimal Pareto set or a part of the global Pareto set.
• We demonstrate the superiority of the proposed algorithm over NSGA-II and MOEA/D on the MOP and DTLZ benchmark problems as well as the effectiveness of solving random forest-based BL model. Moreover, we also test the performance of random forest-based BL model by a comparison with the MV model using the real-world data.

The rest of this paper is organized as follows. We first review some of the past studies on BL models and multi-objective DIRECT algorithms in Section 2. Section 3 is dedicated to relevant preliminaries. Thereafter, in Section 4, a random forest-based BL model is presented, followed by a decomposition-based multi-objective DIRECT algorithm. Then, in Section 5, we demonstrate the superiority of the proposed algorithm over NSGA-II and MOEA/D as well as its effectiveness on the model introduced before. Furthermore, we compare our model with the MV model using the real-world data. Finally, the findings and future perspectives are discussed in Section 6.

2. Related Work

In 1952, Markowitz first applied statistical measures to portfolio selection and proposed the MV model. In 1959, he further introduced the MPT based on the concept of diversification [23] which then became a hot topic in the past few decades. This theory offers an effective tool to obtain the minimized risk under a return target or vice versa. However, several obvious shortcomings have been exposed gradually with the application of the MPT and the MV model, such as depending solely on historical data and being prone to extreme asset allocation. These problems inspired a multitude of researchers to focus on the creative task of model improving.

In 1991, Fischer Black and Robert Litterman proposed the BL model by taking into account investor preferences from the perspective of Bayesian analysis [4]. From then, many researches have proved its advantages in volatile markets [24], which further contributes to its prosperity. For example, Satchell et al. [25] gave a readable mathematical description of the BL model and presented some extensions. Martellini and Ziemann [26] proposed a modified BL model that incorporates positive views of the performance of hedge fund strategies to extend its application to the hedge fund market. Cheung [7] clarified the assumptions and formulation of the BL model and extended it to large portfolio applications using a dimension-reduction technique. Figelman [27] combined a factor structure to BL model and applied it to multi-asset portfolio. Chen and Lim [28] introduced a generalized BL model that accounts for both misspecification and bias in the equilibrium and expert models.

As mentioned earlier, the BL model attaches importance to investor preferences. As a result, particular attention is paid to the interesting issue of how to generate investor preferences systematically and precisely. Alexander and Svetlana [12] made use of random forests to generate investor views and compared the results with other portfolio optimization frameworks. Donthireddy [29] used machine learning (ML) classifiers to acquire investor views and revealed that the enhancement can improve the performance of the BL model significantly. Kara et al. [30] combined GARCH modeling and support vector regression (SVR) to create a novel approach to describe views. In this approach, they forecasted asset returns by SVR with input indicators generated by GARCH and took the returns as the investor views. Min et al. [15] generated quantitative opinions by ML algorithms and demonstrated that random forest maintains the best performance.

BL portfolio construction is on the basis of the classic return-risk optimization, which is a typical MOP. EMO algorithms may be a good choice. However, some difficulties (such as non-reproducibility, the computationally expensive diversity-preservation operator, and the lack of a theoretical guarantee of global convergence) hinder their application. Some researchers change direction to deterministic algorithms. The DIRECT algorithm is an efficient deterministic algorithm that guarantees global convergence. Motivated by its merits, the DIRECT algorithm began to be applied to multi-objective optimization. For example, Wang et al. [31] combined the DIRECT algorithm and multi-objective genetic algorithm (MOGA) to create a novel method called NSDIRECT-GA. In this method, DIRECT is adopted to generate the population used in MOGA. In particular, rank and crowding distance are used to generate potentially optimal hyperrectangles in DIRECT. Al-Dujaili and Suresh [32] extended the DIRECT algorithm and presented the so-called MO-DIRECT. Moreover, they compared MO-DIRECT with other decomposition-based multi-objective techniques and demonstrated its superiority. Wong et al. [33] introduced a multi-objective DIRECT based on the hypervolume indicator, in which the diversification is enhanced by selecting hyperrectangles with higher hypervolume.

3. Preliminaries

3.1. Classical Black–Litterman Model

The MV model can be mathematically defined by (1) [1]

$$\begin{array}{c}\hfill \mathrm{minimize}{W}^T\mathsf{\Sigma }W\\ \hfill \mathrm{maximize}{W}^{T}R\\ \hfill \sum _{i=1}^N{w}_i=1\end{array}$$

(1)

where N is the number of assets available, $R={(r_1,r_2,\dots ,r_N)}^T$ is the expected return vector of assets in the portfolio, $W={(w_1,w_2,\dots ,w_N)}^T$ defines the proportion of capital invested in each asset in the portfolio, and $\mathsf{\Sigma }={\left(\sigma \right)}_{N\times N}$ denotes the covariance matrix of the asset returns.

It is obvious that the MV model uses the variance as a risk measure, targeting at obtaining the weight vector that maximizes the expected return of the portfolio but minimizes the risk. However, from a practical perspective, there exist some drawbacks in the model (1), such as highly-concentrated portfolio [5] and sensitive input [34]. To overcome the above problems, the BL model is proposed by using the investor views of the asset returns to revise R and $\mathsf{\Sigma }$ in the model (1).

The cornerstone of the BL model is to introduce investor preferences from the perspective of Bayesian analysis. With the help of the views, a realistic estimate of expected return is no longer demanded. Instead, what is needed is the prior expected return obtained by historical data, so that a posterior estimate of expected returns can be produced. It is assumed that investors will assign their wealth to N assets and $\widehat{E}\left(R\right)$ denotes the posterior estimate of expected returns, which is composed of two parts. The first part is the market equilibrium return $\mathsf{\Pi }$ obtained by the capital asset pricing model, and the second one comes from views on certain assets. Therefore, $\widehat{E}\left(R\right)$ can be formulated as follows:

$$\widehat{E}\left(R\right)={[{\left(\tau \mathsf{\Sigma }\right)}^{-1}+P^T{\Omega }^{-1}P]}^{-1}[{\left(\tau \mathsf{\Sigma }\right)}^{-1}\mathsf{\Pi }+P^T{\mathsf{\Omega }}^{-1}Q]$$

(2)

where

• Q is a K-dimensional view vector which maintains K subjective returns of certain assets.
• $\mathsf{\Omega }$ refers to a $K\times K$ matrix that reflects the confidence in views.
• P is a $K\times N$ mapping matrix, representing the correspondence between K views and corresponding assets.
• $\mathsf{\Sigma }$ is the $N\times N$ covariance matrix of asset returns.
• $\mathsf{\Pi }$ is the market equilibrium return calculated by the Formula (3), where $w_{mkt}$ is the market-cap weights, and $\delta$ denotes the risk aversion coefficient which is obtained by dividing the market excess return $R-R_f$ by its variance ${\sigma }^2$.

$$\mathsf{\Pi }=\delta \mathsf{\Sigma }{w}_{mkt},\delta =\frac{R-R_f}{{\sigma }^2}$$

(3)

Similarly, the posterior estimate of the covariance matrix is shown in (4).

$$\widehat{\mathsf{\Sigma }}=\mathsf{\Sigma }+{[{\left(\tau \mathsf{\Sigma }\right)}^{-1}+P^T{\mathsf{\Omega }}^{-1}P]}^{-1}$$

(4)

Combining the view vector Q with the mapping matrix P can either express relative or absolute opinions on the expected returns. A relative view usually describes statements about the connection between asset returns, for example, ‘A will outperform B by X%’. On the other hand, an absolute view directly provides specific return forecasts for certain assets, such as ‘C will climb to 8%’. $\mathsf{\Omega }$ is another parameter that is critical to BL model. Based on it, each view can be allocated an uncertainty level. In practice, it is believed to be proportional to the variance of asset returns, and can be obtained by the following expression:

$$\mathsf{\Omega }=\tau P\mathsf{\Sigma }{P}^T$$

(5)

where $\tau$ represents the relative proportion of the variance of asset returns.

3.2. DIRECT Framework

In 1993, Jones et al. [21] took inspiration from Lipschitzian optimization and introduced the DIRECT algorithm, which does not require prior knowledge of the Lipschitz constant. The DIRECT algorithm regards the entire solution space as a hyperrectangle, and continuously divides it to obtain new sub-hyperrectangles. In each iteration, it selects the potentially optimal hyperrectangles according to certain rules which are then subdivided into equal parts along the longest side, and calculates the corresponding function at the center points. It is obvious that one of the keys to determining the algorithm’s efficiency is the strategy for selecting potentially optimal hyperrectangles. Assuming that there are currently m hyperrectangles in the solution space, the jth hyperrectangle is considered as the potentially optimal hyperrectangle if there is a constant $\tilde{K}$ satisfies the Formula (6)

$$\begin{array}{ccc}\hfill f\left(c_j\right)-\tilde{K}{d}_j& \le & f\left(c_i\right)-\tilde{K}{d}_i,\mathrm{for}\mathrm{all}i=1,2,\dots ,m\hfill \\ \hfill f\left(c_j\right)-\tilde{K}{d}_j& \le & f_{min}-ϵ\left|f_{min}\right|\hfill \end{array}$$

(6)

where $ϵ$ is a positive constant, $f\left(c_j\right)$ is the function value at the center point of hyperrectangle j, and $d_i$ refers to the center-vertex distance of hyperrectangle i, and $f_{min}$ denotes the current optimal function value. Figure 1 shows the properties of potentially optimal hyperrectangles more intuitively.

After the potentially optimal hyperrectangles are determined, they will be sub-divided along the longest side to form new sub-hyperrectangles. For an n-dimensional solution space, there may be multiple longest sides of the same length. In such a situation, it is necessary to specify the partition strategy because different strategies form different sub-hyperrectangles. For example, the DIRECT algorithm could successively select the coordinate direction with the smallest value of $v_i=min\{f(c+\delta e_i),f(c-\delta e_i)\}$. The detailed partition procedure is provided in Algorithm 1.

Algorithm 1 Procedure for hyperrectangle partition

Input: Function to be minimized f, the current potentially optimal hyperrectangle i Output: Newly acquired hyperrectangles from i  1: Find the dimension set I with the largest side length and set $\delta$ to be $1/3$ of the largest side length.  2: Calculate the function value at the points $c\pm \delta e_i$ for all $i\in I$, where c refers to the center of a hyperrectangle and $e_i$ denotes the unit vector of corresponding direction.  3: Divide all the longest side dimensions of the current hyperrectangle into 3 equal parts, starting from the dimension of minimum $v_i=min\{f(c+\delta e_i),f(c-\delta e_i)\}$ continuing to the dimension of maximum $v_i$.

Then, the DIRECT algorithm can be expressed in Algorithm 2.

Algorithm 2 DIRECT algorithm

Input: Function to be minimized f, evaluation budget t Output: Approximate minimum of function f  1: Normalize the solution space to a hyperrectangle with $c_1$ to be the center point.  2: Calculate the function value $f\left(c_1\right)$ for $c_1$ and set $f_{min}=f\left(c_1\right),m=1,t=0$.  3: while evaluation budget is not exhausted do  4:     Identify the current set S of potentially optimal hyperrectangles.  5:     while $S\ne \varnothing$ do  6:         Choose a potentially optimal hyperrectangle $j\in S$.  7:         Divide the hyperrectangle j to generate new hyperrectangles by Algorithm 1.  8:         Update $f_{min}$ and $m=m+\Delta m$.      ▹$\Delta m$ is the number of new hyperrectangles generated during the partition procedure.  9:         $S\leftarrow S-\left\{j\right\}$  10:     $t=t+1$

4. Materials and Methods

In this section, we propose a novel model for BL portfolio optimization problems in which investor views are generated on the basis of random forests. Then, a decomposition-based multi-objective DIRECT algorithm is developed to handle the proposed model.

4.1. Random Forest-Based Black–Litterman Portfolio Optimization Model

4.1.1. Generating Views Using Random Forests

As mentioned before, Q and $\mathsf{\Omega }$ revise the priori returns, which is extremely important to the BL model. Nowadays, many attempts are made to generate views by ML algorithms which try to train a predictor based on historical data to form specific return forecasts for certain assets, i.e., Q. However, this kind of methods may not provide the view uncertainty. Therefore, we employ the random forests to form view vector and generate the view uncertainty at the same time.

Random forest [35] is an ensemble method which constructs a set of base classification models, and synthesizes the final output by averaging the predicted value of each base classification model. The main training steps of random forests are shown in Algorithm 3.

Algorithm 3 Main training steps of random forests

Input:S data samples Output: M decision trees 1: for$i\leftarrow 1toM$do 2:     Sample randomly from the original training data set to form S new data samples. 3:     Construct the decision tree based on the new data set and the corresponding node split rules.

As a decision tree based ensemble method, random forest is appealing due to the advantages it provides, including the ability to maintain high prediction accuracy and not being prone to overfitting. Moreover, it can be directly applied in high-dimensional problems. In 2011, Hutter et al. [36] used random forests to replace the Gaussian process (GP) to ease the updating difficulty of GP models under a large-scale known point set in Bayesian optimization. Furthermore, the variance of the predicted values of multiple trees is used as an uncertainty measure. The core idea of this work is that consistent predictions based on different models and different training data indicate more reliable results. In this paper, we follow this idea to form the absolute view vector based on the random forests and measure the view uncertainty by the difference in the predicted values of multiple trees. Then, the procedure of obtaining view vector $Q^{rf}$ and the view uncertainty ${\mathsf{\Omega }}^{rf}$ are summarized as follows.

• Construct a random forest model based on the historical return data to predict asset returns, and use the model output $y_{pred}$ as an absolute view.
• Calculate the maximum absolute deviation between the output of M decision trees and the model output $\mathrm{MAE}=max|y_i-y_{pred}|,\mathrm{for}i\in [1,2,\dots ,M]$.
• Calculate the proportion of subtrees whose absolute deviation from $y_{pred}$ is less than $\frac{\mathrm{MAE}}{10}$ and take it as the uncertainty of the corresponding view.
• Repeat the above process until the K views and corresponding uncertainty are constructed.

Based on the view vector and the view uncertainty generated from random forests concurrently, we can calculate the posterior estimate of expected returns and covariance matrix further. We believe that the adopting of random forests will significantly weaken the conservatism and greatly enhance the generalization of the BL model.

4.1.2. Random Forest-Based Black–Litterman Model

Based on the view vector and the view uncertainty generating from random forests, the posterior estimate of expected returns and covariance matrix can be reformed as follows:

$$\begin{array}{c}\hfill \widehat{E}{\left(R\right)}^{rf}={[{\left(\tau \mathsf{\Sigma }\right)}^{-1}+P^T{\left({\mathsf{\Omega }}^{rf}\right)}^{-1}P]}^{-1}[{\left(\tau \mathsf{\Sigma }\right)}^{-1}\mathsf{\Pi }+P^T{\left({\mathsf{\Omega }}^{rf}\right)}^{-1}{Q}^{rf}]\end{array}$$

(7)

$${\widehat{\mathsf{\Sigma }}}^{rf}=\mathsf{\Sigma }+{[{\left(\tau \mathsf{\Sigma }\right)}^{-1}+P^T{\left({\mathsf{\Omega }}^{rf}\right)}^{-1}P]}^{-1}$$

(8)

Then, we denote the random forest-based BL model with the following formulation:

$$\begin{array}{c}\hfill \mathrm{minimize}{W}^T{\widehat{\mathsf{\Sigma }}}^{rf}W\\ \hfill \mathrm{maximize}{W}^T\widehat{E}{\left(R\right)}^{rf}\\ \hfill \sum _{i=1}^N{w}_i=1\end{array}$$

(9)

which is a typical MOP.

4.2. DIRECT Algorithm for Multi-Objective Optimization

The random forest-based BL model is a typical MOP which is much more complicated than the single-objective ones. The current mainstream algorithms for MOPs are basically based on evolutionary computing, which contains random initialization and various random evolution operators to explore the solution space. However, there are some inherent problems (such as non-reproducibility, computationally expensive diversity-preservation operators, and a lack of a theoretical guarantee of global convergence) limit their application. The DIRECT algorithm is an efficient deterministic algorithm that guarantees global convergence. Motivated by its merits, we extend it to MOPs and propose a new multi-objective DIRECT algorithm to construct Pareto-optimal portfolios of the aforementioned model.

4.2.1. Decomposition Based Strategy for Potentially Optimal Hyperrectangles Selecting

As mentioned earlier, the selecting strategy for potentially optimal hyperrectangles is essential to the DIRECT algorithm, which means a possible modification or alternative that makes it suitable for multiple objectives is a straightforward manner for the DIRECT algorithm’s extension. Decomposition is a basic idea which has been applied in many popular multi-objective algorithms, such as MOEA/D. Based on decomposition, the original problem is transformed into multiple scalar sub-problems in the direction of weight vectors, and then achieve optimization simultaneously. The penalty-based boundary intersection (PBI) approach is a commonly used decomposition method [16]. Through PBI, a MOP is converted into the following scalar optimization subproblem using the weight vector $\lambda$ and the reference point $z^{*}$ (see Figure 2).

$$\begin{array}{cc}\hfill \mathrm{minimize}g\left(x\right|\lambda ,z^{*})=d_1+\theta d_2& \\ \hfill \mathrm{where}{d}_1=\frac{\left|\right|{(z^{*}-F\left(x\right))}^T\lambda \left|\right|}{\left|\right|\lambda \left|\right|}& \\ \hfill \mathrm{and}{d}_2=\left|\right|F\left(x\right)-(z^{*}-d_1\lambda )\left|\right|.\end{array}$$

(10)

$\theta >0$ is the pre-configured penalty parameter. Motivated by this idea, we propose a novel strategy to encourage to choose potentially optimal hyperrectangles in all directions.

To make the strategy extend to multiple objectives, we first introduce a new indicator based on decomposition to quantify the quality of hyperrectangles (more specifically, the solutions at their centers). Suppose there are H weight vectors ${\lambda }_i,i\in [1,2,\dots ,H]$ and N hyperrectangles, the objective function value at the center of hyperrectangle j is $f_j,j\in [1,2,\dots ,N]$, and the reference point is denoted by $z_{ref}$. Then, the procedure of generating the indicator is outlined in Algorithm 4. Following the idea of decomposition, the indicator maintains the ability to consider sub-problems in multiple reference directions at the same time.

Algorithm 4 Generation procedure of the decomposition based indicator

Input:H weight vectors ${\lambda }_i,i\in [1,2,\dots ,H]$, N objective function vectors at the center of N hyperrectangles $f_j,j\in [1,2,\dots ,N]$ and reference point $z_{ref}$ Output: $f_j^s,j\in [1,2,\dots ,N]$ for N hyperrectangles 1: for$j\leftarrow 1toN$do 2:     Calculate the projection of the $(f_j-z_{ref})$ vector in each weight direction using $d_1$ in the Formula (10), select the direction with the largest projection to be the direction h of reference point and denote the corresponding projection value by $p_h^j$. 3:     Update the maximum and minimum projection $p_h^{max},p_h^{min}$ in the direction h by $p_h^j$. 4: for$j\leftarrow 1toN$do 5:     $f_j^s\leftarrow \frac{p_h^j-p_h^{min}}{p_h^{max}+p_h^{min}}$ 6: return$f_j^s,j\in [1,2,\dots ,N]$

As listed in Algorithm 4, the direction with largest $d_1$ (maximum projection value) among all H sub-problems is used as the direction of the proposed indicator of every hyperrectangle, and then the largest $d_1$ of each hyperrectangle will be normalized on the corresponding direction to obtain the value of the indicator (i.e., their contribution on the non-dominated front).

On the basis of the proposed indicator, we can define a strategy for selecting potentially optimal hyperrectangles. That is, the hyperrectangle j is considered potentially optimal if there exists $\tilde{K}$ such that

$$\begin{array}{ccc}\hfill f_j^s-\tilde{K}{d}_j& \le & f_i^s-\tilde{K}{d}_i,\mathrm{for}\mathrm{all}i=1,2,\dots ,m\hfill \\ \hfill f_j^s-\tilde{K}{d}_j& \le & f_{min}^s-ϵ\left|f_{min}^s\right|\hfill \end{array}$$

(11)

where $ϵ$ is a positive constant, m is the number of hyperrectangles to be compared, $f_j^s$ is the indicator value at the center point of hyperrectangle j, and $d_j$ refers to the center-vertex distance of hyperrectangle j, where $f_{min}^s$ denotes the current optimal indicator value.

With this unique selecting strategy, the DIRECT algorithm possesses the capacity to select potentially optimal hyperrectangles in all reference directions simultaneously, which will further improve the exploratory nature when the algorithm is stuck in a locally optimal Pareto set or a part of the global Pareto front.

4.2.2. Decomposition Based Partition Procedure

The original definition of the partition procedure also does not fit the MOPs. Therefore, based on the indicator mentioned earlier, we devised a novel idea to split the solution space into smaller hyperrectangles. Specifically, after the set of potentially optimal hyperrectangles is defined, we evaluate all those hyperrectangles at the points of $c\pm \frac{2}{3}e{L}_{max}$ in every longest side dimension where c and $L_{max}$ are their centers and the maximum side length, respectively, and e denotes the unit vector of corresponding direction. Followed by generating the indicator value $f^s(c+\frac{2}{3}e{L}_{max})$ and $f^s(c-\frac{2}{3}e{L}_{max})$ for all these hyperrectangles, we divide each hyperrectangle into thirds starting from the dimension of minimum $v=min\{f^s(c+\frac{2}{3}e{L}_{max}),f^s(c-\frac{2}{3}e{L}_{max})\}$ continuing to the dimension of maximum v.

4.2.3. Decomposition Based Multi-Objective DIRECT Algorithm

Lying on the solid foundation discussed in Section 4.2.1 and Section 4.2.2, we extend the DIRECT algorithm to MOPs, and introduce a decomposition based multi-objective DIRECT (multiDecompose) algorithm. Figure 3 provides a description of multiDecompose summarizing the adaptations to MOPs presented earlier. Despite the modification, multiDecompose will not come with a increased computational complexity due to the the complexity of Algorithm 4 is $O\left(MN\right)$ which is of the same nature of the original one, where M denotes the number of weight vectors and N refers to the quantity of the hyperrectangles remaining to be compared.

5. Results

In this section, we provide the results comparing multiDecompose with two widely adopted multi-objective algorithms, and discuss them from different dimensions. Then, after the forecasting performance of the proposed technique random forest model is verify, we address the random forest-based BL model with multiDecompose and compare the results with the classic MV model.

5.1. Performance of multiDecompose on Benchmark Problems

We use 8 widely used MOP and DTLZ test instances (listed in

Table 1. Details of all test instances.

Test Instance	Number of Objectives	Solution Space Dimension
MOP1	2	1
MOP2	2	3
MOP3	2	2
MOP4	2	3
DTLZ1	3	7
DTLZ2	3	10
DTLZ3	3	10
DTLZ4	3	10

) to compare multiDecompose with NSGA-II and MOEA/D, two representative multi-objective algorithms. NSGA-II and MOEA/D are implemented by the pymoo library [37] with default parameters, except that the population size is set to be 200 for all test instances. Additionally, in multiDecompose, we set $H=50,ϵ={10}^{12},tol=0.005$. Then, a set of experiments is carried out to provide the comparison between these three algorithms in terms of Pareto front ($PF$) and the most popular quality indicator, hypervolume ($HV$) [38]. Both these algorithms have been run 30 times independently for each test instance.

The $PFs$ of the three algorithms on the MOP and DTLZ test instances are summarized in Figure 4 and Figure 5, respectively. It is clear from the two figures that, as to the quantity and uniformness of final solutions, multiDecompose is superior to the other two algorithms for all test instances, and even illustrates an overwhelming advantage on the DTLZ instances. NSGA-II shows a mediocre performance while MOEA/D finds few non-dominated solutions under a specific number of function evaluations in most test instances. Furthermore, from the perspective of the number of function evaluations, multiDecompose converges much faster than NSGA-II and MOEA/D on all test instances.

Table 2. Average $HV$ value of the solutions obtained by multiDecompose, NSGA-II and MOEA/D.

	multiDecompose	NSGA-II	MOEA/D
MOP1	13494.20024	12157.93554	8265.419011
MOP2	0.304549570	0.302570066	0.296286100
MOP3	360.3988959	360.3714609	354.1858850
MOP4	27.37042078	27.35540827	26.62617847
DTLZ1	4.947239949	4.937894772	4.071856060
DTLZ2	0.506487643	0.447724955	0.123985535
DTLZ3	31.63621196	30.60244867	23.59450020
DTLZ4	0.547895448	0.506647045	0.073129642

presents the mean of $HV$-metric value of final solutions obtained by each algorithm on each test instance. It is obvious that our algorithm also performs splendidly on each test instance in terms of $HV$-metric, i.e., our algorithm maintains the maximum $HV$-metric value for all test instances. NSGA-II is close behind, while MEAD/D performs the worst. It can be inferred from the aforementioned facts that multiDecompose can achieve a better performance over NSGA-II and MOEA/D in terms of the solution quality and convergence speed on the 2-objective and 3-objective optimization problems.

5.2. Portfolio Performance of the Random Forest-Based BL Model

Based on the data set from the 10 representative stocks (shown in the

Table 3. Stock list used in the experiments.

	Code	Name	Location
1	600519	Kweichow Moutai Co., Ltd.	Renhuai, China
2	601318	Ping An Insurance (Group) Co., Ltd.	Shenzhen, China
3	600036	China Merchants Bank	Shenzhen, China
4	601398	Industrial and Commercial Bank of China	Beijing, China
5	602276	Jiangsu Hengrui Medicine Co., Ltd.	Lianyungang, China
6	601186	China Railway Construction	Beijing, China
7	601288	Agricultural Bank of China	Beijing, China
8	603288	Foshan Haitian Flavouring and Food Co., Ltd.	Foshan, China
9	601012	Longi Green Energy Technology Co., ltd.	Xi’an, China
10	600031	Sany Heavy Industry Co., ltd.	Beijing, China

) of Shanghai Stock Exchange 50 Index, which covers the monthly volume and price data over a period of 5 years from January 2014 to December 2018, we apply multiDecompose to solve the random forest-based BL model and compare the obtained portfolio with the classic MV portfolio.

More specifically, we first train a random forest model on the basis of the dataset mentioned earlier to predict monthly stock returns. Then, we summarize the squared prediction errors for all 10 stocks in 2019 and display the distribution in Figure 6. It can be seen from the histogram that most of the squared prediction errors are maintained below 0.005, which verifies that the forecast values obtained from the random forest model are qualified to generate a view vector and the confidence matrix for the subsequent BL model.

Furthermore, we set up the monthly warehouse transfer rule following the convention in the financial industry, which means that we redistribute wealth monthly based on the solutions obtained by the portfolio optimization models. In this case, the stock returns predicted by the trained random forest model are used to form a view vector and the view uncertainty monthly, followed by the generation of a posterior estimate of expected returns, and a covariance matrix for the proposed model. After that, the proposed models are solved based on the multiDecompose algorithm and the solution with the largest Sharpe ratio is selected as the portfolio in the following month.

Figure 7 reveals the cumulative returns of the proposed model and classic MV model in 2019. Obviously, the proposed model significantly outperforms the classic MV model during the whole period of backtesting, and its annual return reaches $70\%$ compared to $15\%$ for the MV model, which is consistent with our analysis before. From the perspective of the Sharpe ratio, the proposed model has a better Sharpe ratio of $0.903$ while the ratio is $0.359$ for the MV model. We can conclude from the above results that the proposed model is superior to the MV model in terms of portfolio performance.

6. Discussion

In this paper, a random forest-based BL portfolio optimization model was developed to handle the portfolio selection problem. In this model, a novel method for generating investor views on the basis of random forests was adopted. More specifically, the view vector was generated based on the predicted asset returns obtained by random forests and the confidence matrix which contains the uncertainty of each view was measured by the difference in the predicted values of multiple trees. Aiming to effectively resolve the proposed model while maintaining high accuracy, a novel decomposition based multi-objective DIRECT algorithm was introduced, which possesses the capacity to select potentially optimal hyperrectangles in all reference directions simultaneously. It was demonstrated that the proposed algorithm achieves a better performance over NSGA-II and MOEA/D on the MOP and DTLZ benchmark problems in terms of the quality of the Pareto front. Moreover, the experiments certified that the random forest-based BL portfolio optimization model can obtain higher cumulative returns and Sharpe ratio in the application of Chinese stock market when compared to the classic MV model.

For future work, the proposed model can be extended with more trading constraints. Furthermore, the proposed algorithm can be tested on larger historical financial datasets and generalized to fuzzy portfolio selection problems.