Selecting Socially Responsible Portfolios : A Fuzzy Multicriteria Approach

We propose a multi-objective approach for portfolio selection, which allows investors to consider not only return and downside risk criteria but also to include environmental, social and governance (ESG) scores in the investment decision-making process. Owing to the uncertain environment of portfolio selection, the return and ESG score of each asset are considered as independent L-R power fuzzy variables. To make the model more realistic, we take budget, floor ceiling and cardinality constraints into account. In order to select the optimal portfolio along the efficient frontier, we apply the Sortino ratio in a credibilistic environment. The subsequent empirical application uses a data set from Bloomberg’s ESG Data in combination with US Dow Jones Industrial Average data. The experimental results show that the proposed model offers promising results for socially responsible investors seeking ethical and sustainability goals beyond the return-risk trade-off and its ability to beat the benchmark.


Introduction
Multiple criteria decision making is a method which pursues making choices in the presence of multiple criteria, goals or objectives [1]. It has been applied in several areas [2][3][4] and in investment decisions are not an exception [5][6][7][8].
In the last decades, financial, social and environmental crises have led to considerable growth of social responsibility (SR) consciousness in society and amongst investors [9][10][11]. Socially responsible investment (SRI), also known as ethical, green, impact, responsible or sustainable investing, is defined as an investment decision-making process that considers not only financial performance but also environmental, social and governance (ESG) criteria [12,13]. In other words, socially responsible investors (SRIs) consider non-financial criteria (e.g., social responsibility performance) in their investment decisions in addition to conventional criteria such as return, risk or liquidity [14,15].
Nowadays, the main investment instrument of SRIs are socially responsible mutual funds [16]. The first modern SR mutual fund, the Pax World Fund, was created in 1971 in the United States and appealed to investors objecting to the Vietnam War [17]. Since then, the SRI industry has expanded by drawing attention of academics and professionals in various cultural contexts and industries around the world [18][19][20][21][22][23][24][25][26][27].
In order to gauge a fuzzy event, Zadeh [58] proposed a non-self-dual possibility measure. However, since the self-dual property is intuitive and important both in theory and practice, Liu and Liu [61] developed a self-dual credibility measure. Different researchers have used credibility distributions to approximate the uncertainty on returns [62][63][64][65][66][67][68]. Nevertheless, to the best of our knowledge, there is no SR portfolio selection model using a credibility-based approach to deal with fuzzy events. Hence, this paper seeks to make a contribution to the SR portfolio selection model approach by assuming that both the return and ESG performance of each asset are an L-R power fuzzy variable whose moments are assessed employing their credibility distributions.
Most previous studies use the variance to measure portfolio risk despite its undesirable properties such as the lack of coherence [48]. However, in real-life decisions investors are more sensitive to downside losses than upside gains but variance does not reflect investors' concern about underperformance [69]. In order to solve this limitation, several downside risk measures have been proposed. The mean-absolute semi-deviation (MASD) is a popular downside risk measure originally introduced by Speranza [70] and it is used in mean-absolute semi-deviation portfolio selection models. MASD was introduced as a downside risk measure in stochastic environments. However, it has only been used in a single study in a credibilistic context [64].
Modern portfolio theory assumes that investors are rational and always try to minimize risk while striving for the highest return possible [48]. Nowadays, the so-called Sortino ratio introduced by Sortino and Van Der Meer [71] is a well-known measure to evaluate portfolio performance. It is a modification of the Sharpe ratio, using downside measure rather than dispersion measure as the measure of risk. Therefore, to select the optimal portfolio along the efficient frontier, this paper defines for the first time the credibilistic Sortino ratio as the ratio between the credibilistic risk premium and the credibilistic mean-absolute semi-deviation.
On the basis of the SRI literature review, the aim and contribution of this paper is two-fold: First, we present a multi-objective approach for the portfolio selection model, which allows investors to consider not only return and downside risk criteria but also ESG performance as part of the investment decision-making process. The return and the ESG score of each asset are quantified by the credibilistic mean values of return and the ESG score, respectively; the downside risk is determined by the credibilistic mean-absolute semi-deviation of return. The return and the ESG score are considered as L-R power fuzzy numbers. In order to make this model more applicable to reality, we consider budget, floor-ceiling and cardinality constraints. To solve the constrained portfolio optimisation problem, we apply the Non-dominated Sorting Genetic Algorithm II (NSGA-II). Second, we present a real-world empirical application using a data set from Bloomberg's ESG Data in combination with US Dow Jones Industrial Average data for the period from February 2014 to December 2017.
The experimental results show that the proposed model offers promising results for socially responsible investors seeking ethical and sustainability goals beyond the return-risk trade-off and its ability to beat the benchmark.
Furthermore, in Section 2, the fundamental concepts relevant to the proposed study are introduced. Thereafter, Section 3 describes the mean-downside risk-socially responsibility model. The solution methodology to solve the above model is presented in Section 4, followed by an empirical study in Section 5, which illustrates the usefulness of the proposed model. The paper closes with the key conclusions.

L-R Fuzzy Numbers and Credibility Theory
This section briefly presents the key definitions of fuzzy set theory essential for the subsequent sections.

L-R Power Fuzzy Numbers
The reference functions L π (.), R ρ (.): The membership function of L-R fuzzy number A [i.e. A = (a, b, c, d) L π R ρ ] has the following form: where (b−a) and (d−c) represent the left and right spreads of A, respectively. Following Jalota et al. [67], this research considers the reference functions of the power family of positive parameters π and ρ, where L π (x)= 1 − x π and R ρ (x)= 1 − x ρ , respectively. All over this study, L-R power fuzzy numbers will be defined as A = (a, b, c, d) π,ρ .

Credibility Theory
The credibility theory was founded by Liu [73] in 2004 and refined by Liu [74] in 2007. It has been fairly well applied to the study of the behaviour of fuzzy phenomena. Definition 1. Credibility measure [61]. Let ξ be a fuzzy variable with membership function µ and x a real number. The credibility measure of a fuzzy event, characterised by ξ ≤ x, is defined by Equation (1).
The credibility measure of L-R power fuzzy numbers is expressed as: Definition 2. Expected value [61]. Let ξ be a fuzzy variable. Then the expected value of ξ is defined by Equation (2) provided that at least one of the two integrals is finite.
The crisp equivalent expression for the credibilistic expected value of an L-R power fuzzy number is formulated by Equation (3) [67].
where, Finally, following Vercher and Bermúdez [64], the crisp equivalent expression for the credibility-measure-based MASD of an L-R power fuzzy number is expressed as:

The Mean-Downside Risk-Socially Responsible Portfolio Selection Model
This section formulates the mathematical model of the socially responsible portfolio investment. The portfolio selection model considers three criteria: risk, return and social responsibility (RRS). The parameters and variables used in this study to formulate this model are: ξ r p : L-R fuzzy return of the portfolio, that is, ξ r p = a r p , b r p , c r p , d r p π rp ,ρ rp , e : expected return of the portfolio, u i : upper bound of the proportion invested in the i-th asset, l i : lower bound of the proportion invested in the i-th asset, k: the number of assets held in the portfolio, ω i : the proportion of wealth invested in the i-th asset, y i : a binary variable indicating whether the i-th asset is contained in the portfolio. It takes value 1, if the i-th asset is included in the portfolio, otherwise it takes value 0.

Return
Stock returns are beset with ambiguity and vagueness. For the above reasons, this research assumes that an investor allocates his/her total wealth among n US risky assets which offer L-R fuzzy returns. Following Vercher and Bermúdez [64,75,76] we build the membership function of the L-R fuzzy return by using the empirical percentiles of their historical returns. Concretely, the support and core of ξ r i = a r i , b r i , c r i , d r i π r i ,ρ r i are represented by intervals [q 3 , q 97 ] and [q 45 , q 55 ], respectively.
The shape parameters are obtained in such a way that the fuzzy and empirical quartiles coincide, Then, the maximisation of the expected return of the portfolio can be expressed as:

ESG score
In recent years, portfolio managers are incorporating ESG criteria during portfolio creation. In this study, the ESG score of each asset is defined by the Bloomberg ESG Disclosure Score. Bloomberg provides insight on ESG data for over 9500 companies across over 83 countries. This is collected from company's public filings such as CSR reports, annual reports and websites and a proprietary Bloomberg survey that requests corporate data directly. The ESG score measures the SR performance of a company on a scale from 0 to 100, with 100 being the best-rated ESG score possible and 0 being lowest ESG score possible. In this work, we divide Bloomberg ESG scores by 100, resulting in ESG scores, which vary within a range from 0 to 1.
As already mentioned, social performance criteria used to measure social responsibility of companies are by their very nature ambiguous, uncertain and imprecise. Therefore, this research assumes that the ESG score of each asset is a L-R fuzzy number. The parameters of the L-R power fuzzy numbers are obtained from the empirical percentiles of the historical ESG score data, as previously explained. Then, the maximisation of the expected ESG score of the portfolio can be expressed as:

Risk
Mean-absolute semi-deviation of returns is a more appropriate risk measure than variance of returns, since it is focuses exclusively on adverse deviations. By taking the above into account, portfolio risk is calculated using expected value based MASD measure. Therefore, the minimisation of MASD of the portfolio can be expressed as:

Real-World Constraints
Portfolio managers incorporate different types of restrictions, which consider investor preferences and market conditions. Based on the above statement, our proposed model includes the following five real-world constraints: (i) Capital budget constraint on the assets is defined as (ii) Maximum number of assets held in the portfolio is defined as (iii) Maximal proportion of the capital that can be invested in a single asset is defined as (iv) Minimal proportion of the capital that can be invested in a single asset is defined as (v) No short selling of assets is defined as On the basis of the aforementioned assumptions, the RRS portfolio selection model is formulated as: Here, an admissible portfolio P e is said to be Pareto-efficient if there is no other admissible portfolio P o such that P o , and, with at least one strict inequality [76].
The efficient solutions set is the Pareto optimal set in the decision space and each of their three objective function values constitute the Pareto optimal frontier in the criterion space. The solutions of the Pareto optimal frontier are non-dominated.

Solution Methodology
Incorporating real-world constraints into the portfolio optimisation problem converts it into a constrained portfolio optimisation problem (CPOP) that is NP-hard. In order to overcome this drawback, various multi-objective evolutionary algorithms (MOEAs) have been proposed in the literature. Among them, the Non-dominated Sorting Genetic Algorithm II (NSGA-II) is the most commonly used MOEA for solving CPOPs [77]. In this study, we propose a methodology based on the algorithm NSGA-II introduced by Deb [78] in 2001. The main objective of this algorithm is to find a set of solutions ordered by fronts under the concept of Pareto dominance. Procedural steps of algorithm NSGA-II are those described by Deb et al. [79] in 2002.
The experimental parameter configuration for testing this algorithm is: population size (400), distribution index for crossover (10); probability of crossover (0.9); distribution index for mutation (50); probability of mutation (0.01); and the maximum number of generations (500). Figure 1 shows the overall structure of the MDRSR modelling approach. Subject to , and, with at least one strict inequality [76]. The efficient solutions set is the Pareto optimal set in the decision space and each of their three objective function values constitute the Pareto optimal frontier in the criterion space. The solutions of the Pareto optimal frontier are non-dominated.

Solution Methodology
Incorporating real-world constraints into the portfolio optimisation problem converts it into a constrained portfolio optimisation problem (CPOP) that is NP-hard. In order to overcome this drawback, various multi-objective evolutionary algorithms (MOEAs) have been proposed in the literature. Among them, the Non-dominated Sorting Genetic Algorithm II (NSGA-II) is the most commonly used MOEA for solving CPOPs [77]. In this study, we propose a methodology based on the algorithm NSGA-II introduced by Deb [78] in 2001. The main objective of this algorithm is to find a set of solutions ordered by fronts under the concept of Pareto dominance. Procedural steps of algorithm NSGA-II are those described by Deb et al. [79] in 2002.
The experimental parameter configuration for testing this algorithm is: population size (400), distribution index for crossover (10); probability of crossover (0.9); distribution index for mutation (50); probability of mutation (0.01); and the maximum number of generations (500). Figure. 1 shows the overall structure of the MDRSR modelling approach.

Experimental Results
In this section, we illustrate the usefulness of the MDRSR model with a real-world empirical study using a data set from Bloomberg's ESG Data in combination with data from companies included in the Dow Jones Industrial Average (DJIA). The data corresponds to weekly closing adjusted prices and ESG scores, observed in t = 203 periods from 10 February 2014 until 25 December 2017. In this study, we only consider 29 companies (n = 29), as one company was not included in the Bloomberg ESG Disclosure Score throughout the period under consideration.
The asset returns sample (r it ) is obtained as r it = p it −p it−1 / p it−1 , i = 1, 2, . . . , 29; t = 1, 2, . . . , 203, where p it is the closing price of the i-th asset on Friday of week t. Then, the membership function of the L-R fuzzy return (ξ r i ) and the L-R fuzzy ESG score (ξ ESG i ) are obtained from the sample percentiles of their historical returns and ESG score, respectively, as explained in Section 3.
In order to analyse the performance of our procedure, this study assumes as diversification parameters: l i = 0.05 and u i = 0.3 for every i = 1, 2, . . . , 29. Following Gupta et al. [80], it is not recommended to have either very few or a very large number of assets in the portfolio so as to achieve diversification. Overall, portfolio diversification by investors lies in the narrow range of 3-10 assets. Furthermore, Skolpadungket et al. [81] state that by limiting the number of assets included in the portfolio, transaction costs are reduced significantly due to the reduction of brokerage fees and taxes. By following the above recommendations, this paper considers setting k = 10 assets for an admissible portfolio. Figure 2 shows a three-dimensional representation of the approximated Pareto front (or efficient front) generated by NSGA-II for the RRS model. The set of points represents the set of non-dominated (or efficient) portfolios for which none of the three objectives (i.e., return, downside risk or ESG score) can be improved without sacrificing any others. As can be observed, NSGA-II provides efficient portfolios, which are widely distributed along the efficient front.

Experimental Results
In this section, we illustrate the usefulness of the MDRSR model with a real-world empirical study using a data set from Bloomberg's ESG Data in combination with data from companies included in the Dow Jones Industrial Average (DJIA). The data corresponds to weekly closing adjusted prices and ESG scores, observed in t = 203 periods from 10 February 2014 until 25 December 2017. In this study, we only consider 29 companies (n=29), as one company was not included in the Bloomberg ESG Disclosure Score throughout the period under consideration.
The asset returns sample (rit) is obtained as r it = p it -p it-1 p it-1 , i = 1, 2, …, 29; t = 1, 2, …, 203, where pit is the closing price of the i-th asset on Friday of week t. Then, the membership function of the L-R fuzzy return ( ξ r i ) and the L-R fuzzy ESG score (ξ ESG i ) are obtained from the sample percentiles of their historical returns and ESG score, respectively, as explained in Section 3.
In order to analyse the performance of our procedure, this study assumes as diversification parameters: li= 0.05 and ui= 0.3 for every i = 1,2,…, 29. Following Gupta et al. [80], it is not recommended to have either very few or a very large number of assets in the portfolio so as to achieve diversification. Overall, portfolio diversification by investors lies in the narrow range of 3-10 assets. Furthermore, Skolpadungket et al. [81] state that by limiting the number of assets included in the portfolio, transaction costs are reduced significantly due to the reduction of brokerage fees and taxes. By following the above recommendations, this paper considers setting k = 10 assets for an admissible portfolio. Figure. 2 shows a three-dimensional representation of the approximated Pareto front (or efficient front) generated by NSGA-II for the RRS model. The set of points represents the set of non-dominated (or efficient) portfolios for which none of the three objectives (i.e., return, downside risk or ESG score) can be improved without sacrificing any others. As can be observed, NSGA-II provides efficient portfolios, which are widely distributed along the efficient front.
The influence of the inclusion of the ESG score criterion into the Markowitz base model is examined in Figure. 2: non-dominated portfolios with the highest expected return and downside risk achieved the lowest expected ESG score. Furthermore, non-dominated portfolios with the highest expected ESG score obtain both a medium expected return and downside risk. As previously mentioned, the table reveals that there is no portfolio that dominates all others in terms of achieving both a high-expected return and ESG score and low risk simultaneously.  The influence of the inclusion of the ESG score criterion into the Markowitz base model is examined in Figure 2: non-dominated portfolios with the highest expected return and downside risk achieved the lowest expected ESG score. Furthermore, non-dominated portfolios with the highest expected ESG score obtain both a medium expected return and downside risk. As previously mentioned, the table reveals that there is no portfolio that dominates all others in terms of achieving both a high-expected return and ESG score and low risk simultaneously. In addition, we applied the Spearman's rank partial correlation to analyse the relationship between the three selected criteria (i.e., return, ESG score and downside risk). Thanks to this measure, we can check the coherence of the relationship between return and downside risk, removing the influence of the ESG score on both variables. It was not possible to identify the relationship of both variables with the ESG score. As shown in Table 1, there is a positive correlation between return and downside risk, which means that more profitable portfolios are riskier. The relationship between ESG Scores and downside risk could be accounted for by the fact that stocks of companies with a high ESG scores are traded more often than firms with low ESG scores and hence are subject to more volatility, hence increasing their risk. Finally, the relationship between ESG Scores and return suggest that stocks of companies with high ESG scores are not as profitable as their peers with low ESG scores. This paper has proposed a fuzzy multi-objective approach that optimises the expected return, the expected ESG score and the downside risk of a given portfolio, taking budget, floor-ceiling and cardinality constraints into consideration. Each portfolio on the efficient frontier ( Figure 2) offers the maximum possible expected return and expected ESG score for a given level of downside risk. However, investors and portfolio managers want to select the optimal portfolio along the efficient frontier that meets their preferences and expectations.
Modern portfolio theory assumes that investors are rational and always try to minimise risk while striving for the highest return possible [48]. The so-called Sortino ratio introduced by Sortino and Van Der Meer [71] in 1991 is a well-known measure to evaluate portfolio performance. It is a modification of the Sharpe ratio [82], using a downside measure rather than a dispersion measure to quantify risk. To select the optimal portfolio along the efficient frontier, we will define the credibilistic Sortino ratio as the ratio between the credibilistic risk premium and the credibilistic mean-absolute semi-deviation.
According to Sharpe [82], an optimal portfolio is one that seeks to maximise expected risk-adjusted returns. In order to select the SR optimal portfolio, we maximise the credibilistic Sortino ratio; it gauges the risk-adjusted return of an asset or portfolio in a downside risk framework. The Sortino ratio is expressed as: where, E ξ p is the expected fuzzy return of the portfolio, MASD(ξ p ) is the fuzzy downside risk of the portfolio and E(ξ Rf ) is the target or required rate of return (i.e., US 1-Year Treasury Constant Maturity Rate). Finally, we compare the performance of the SR optimal portfolio selected using the proposed Sortino ratio with a similar investment alternative in the financial markets. Table 2 compares the SR optimal portfolio with its benchmark, the SPDR Dow Jones Industrial Average ETF (DIA), which replicates the performance of DJIA Index, during the period from January 2018 to December 2018. Note that the selected SR portfolio beats the benchmark regarding return, value at risk and ESG score by an average of 1.007, 729 y 25 bps, respectively. On the basis of the above knowledge, it can be concluded that the applied model in this research offers promising results for SR-investors seeking ethical and sustainability goals beyond the return-risk trade-off.

Conclusions
This paper contributes to the SRI literature by proposing a novel fuzzy multi-objective approach that optimises the expected return, the expected ESG score and the downside risk of a given portfolio, subject to real-world constraints such as budget, floor-ceiling and cardinality. Owing to the uncertain environment of portfolio selection, we considered the return and the ESG score of each as L-R power fuzzy numbers, where all their parameters were obtained from the sample percentiles of the historical data set of the returns and ESG score, respectively. Furthermore, in order to quantify the return and ESG score of the portfolio we used the credibilistic mean and to quantify the downside risk of the portfolio we used credibilistic mean-absolute semi-deviation. A real-world empirical study was presented using a data set from Bloomberg's ESG Data in combination with Dow Jones Industrial Average (DJIA). Additionally, by maximising the expected risk-adjusted returns (via credibilistic Sortino ratio), we selected the SR optimal portfolio and evaluated its performance in relation to another similar investment alternative in financial markets. The results of the numerical experiment establish that the proposed model offers promising results for socially responsible investors seeking ethical and sustainability goals beyond the return-risk trade-off.
Our proposed approach has several advantages over previous work on SR portfolio selection as we use, for the first time: i) the credibilistic fuzzy framework to deal with uncertainty, vagueness and imprecision environment of the SR portfolio selection; ii) the credibilistic mean-absolute semi-deviation of returns to gauge the portfolio risk; iii) the Bloomberg ESG Disclosure Score to measure the level of social responsibility of each asset; iv) the credibilistic Sortino ratio as the ratio between the credibilistic risk premium and the credibilistic mean-absolute semi-deviation to select the optimal portfolio along efficient frontier.
Finally, there are several future research possibilities for overcoming the limitations of our study. First, this study considers a specific period of time of 47 months. However, stock prices are sensitive to news to changes in macroeconomic and financial conditions. Therefore, a future interesting research would be to consider a wider time window in order to observe and identify stock market behaviour (i.e., bear and bull markets). On the other hand, this research uses a data set from Bloomberg's ESG Data in combination with US Dow Jones Industrial Average. Nevertheless, a future research objective would be to consider other ESG rating scores, such as the Thomson Reuters ESG Score, the MSCI ESG Score and the Vigeo-Eiris ESG Score, and compare their results with the proposed model. Additionally, this study considers the below-mean absolute semi-deviation as a measure of portfolio risk. However, in the literature has been proposed other coherent risk measures such as expected shortfall, conditional value at risk and lower partial moment. Therefore, it would be interesting to consider some of these risk measures in order to quantify portfolio risk. Moreover, in this research, we consider five real-world constraints. Nevertheless, it is still far from reflecting the real-world trading scenarios. Thus, a future research opportunity would be to extend the proposed model and incorporate other realistic constraints, such as transaction costs, class, trading and turnover and so forth. Funding: This research was not funded.

Conflicts of Interest:
The authors declare no conflict of interest.