Next Article in Journal
Optimizing Sensor-Controlled Systems with Minimal Intervention: A Fuzzy Relational Calculus Approach
Previous Article in Journal
Development of a Compartment Model to Study the Pharmacokinetics of Medical THC after Oral Administration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling of Mean-Value-at-Risk Investment Portfolio Optimization Considering Liabilities and Risk-Free Assets

1
Department of Mathematics, Faculty of Mathematics and Natural Sciences, Universitas Padjadjaran, Sumedang 45363, Indonesia
2
Faculty of Business and Management, Universiti Sultan Zainal Abidin, Kuala Terengganu 21300, Malaysia
3
Doctoral Program of Mathematics, Faculty of Mathematics and Natural Sciences, Universitas Padjadjaran, Sumedang 45363, Indonesia
4
Faculty of Informatics and Computing, Universiti Sultan Zainal Abidin, Kuala Terengganu 21300, Malaysia
*
Author to whom correspondence should be addressed.
Computation 2024, 12(6), 120; https://doi.org/10.3390/computation12060120
Submission received: 10 May 2024 / Revised: 29 May 2024 / Accepted: 30 May 2024 / Published: 11 June 2024

Abstract

:
This paper aims to design a quadratic optimization model of an investment portfolio based on value-at-risk (VaR) by entering risk-free assets and company liabilities. The designed model develops Markowitz’s investment portfolio optimization model with risk aversion. Model development was carried out using vector and matrix equations. The entry of risk-free assets and liabilities is essential. Risk-free assets reduce the loss risk, while liabilities accommodate a fundamental analysis of the company’s condition. The model can be applied in various sectors of capital markets worldwide. This study applied the model to Indonesia’s mining and energy sector. The application results show that risk aversion negatively correlates with the mean and VaR of the return of investment portfolios. Assuming that risk aversion is in the 5.1% to 8.2% interval, the maximum mean and VaR obtained for the next month are 0.0103316 and 0.0138270, respectively, while the minimum mean and VaR are 0.0102964 and 0.0137975, respectively. The finding of this study is that the vector equation for investment portfolio weights is obtained, which can facilitate calculating investment portfolio weight optimization. This study is expected to help investors control the quality of appropriate investment, especially in some stocks in Indonesia’s mining and energy sector.

1. Introduction

The investing technique of buying several assets with different rates of loss, rates of return, and performance in a portfolio can make fluctuations in returns and losses mutually compensated between assets. It causes the fluctuations to be more controlled overall [1,2,3,4]. It is then called portfolio diversification. In portfolio diversification, allocating investment funds for each asset is critical. The allocation of each asset must be calculated so that the portfolio generates maximum returns and minimum losses [5,6,7]. The problem is determining the allocation weight for each asset.
The problem of determining the weight of asset allocation in a general investment portfolio is solved using the Markowitz asset allocation optimization model [8,9]. For ease of mention, after this, the model is referred to as the Markowitz model. The Markowitz model is based on two simple objective objectives in portfolio preparation: maximizing the mean and minimizing the variance of portfolio returns in a single period. Then, the two objectives are integrated into one as an optimization model for maximizing the subtraction between the mean and variance of portfolio returns in a single period. Because of this, the Markowitz model is often referred to as the mean-variance model.
Previous studies have simulated the Markowitz model. Ivanova and Dospatliev [10] applied the Markowitz model to determine the allocation of asset weights to investment portfolios on stock data in Bulgaria from 2013 to 2016. Kulali [11] and Zavera [12] also applied the same model to stock data in Turkey and Romania in 2015, respectively.
Over time, researchers began to develop the Markowitz model. Kashirina et al. [13] developed a long-term Markowitz model. The long period is divided into several short periods that are mutually exclusive. In more detail, the return earned at the end of each short-term period is not withdrawn but is reinvested in the next period. Another development of the model of Kashirina et al. [13] is that the variance of the portfolio returns in the model is calculated based on the geometric mean. Steinbach [14] also developed the Markowitz model to be multi-period by providing a quantitative treatment of the tradeoff between returns and losses.
The previous developments of the Markowitz model still use the variance of the portfolio return as a minimized objective [15]. The variance measure can be substituted for the maximum possible loss. This measure is referred to as value-at-risk (VaR). The use of the VaR measure has advantages over the use of the variance measure. Using VaR measures can make the size of the maximum loss of an investment portfolio experienced by investors more diverse. In other words, investors can set the maximum loss they may experience [16]. This portfolio asset-allocation optimization model based on VaR is called the mean-VaR model.
Banihashimi et al. [15] applied the mean-VaR model to stock data in Iran. The performance of the weight allocation results obtained is then evaluated for non-parametric efficiency using data envelopment analysis (DEA). Then, Sukono et al. [17] formulated the mean-VaR model using vector and matrix forms. The formulation of the model was based on stock data in Indonesia. They used an additional assumption, where the return of each asset was assumed to have non-constant volatility. This non-constant volatility was accommodated through the generalized autoregressive conditional heteroscedastic (GARCH) model. Sukono et al. [17] used the Lagrangian multiplier method and the Kuhn–Tucker theorem to obtain a solution from the model. Then, Pandiangan et al. [18] developed the Sukono et al. [17] model by adding risk aversion and risk-free assets to its investment portfolio. The analysis showed that the greater the risk aversion, the smaller the mean and VaR of the investment portfolio return.
This paragraph briefly discusses the gaps in previous studies. The existing mean-variance and mean-VaR models do not involve risk-free assets. Including risk-free assets in investment portfolios reduces volatility and enhances stability. Financial reports show that portfolios with risk-free assets, such as government bonds, have higher Sharpe ratios, indicating better returns relative to the risk taken compared to risky assets. Examples are Kashirina et al. [13], Steinbach [14], Banihashemi et al. [15], Miyangaskari et al. [16], and Sukono et al. [17]. Only Pandiangan et al. [18] involved such risk-free assets. Then, the existing mean-variance and mean-VaR models have not considered the assets-owning companies’ liabilities. A liability analysis is crucial in investment portfolio preparation, as it reveals credit risks and potential bankruptcy, which can impact investment performance. Moody’s report indicates that companies with high debt-to-equity ratios are more likely to default, which can negatively affect an investor’s portfolio value if not calculated carefully. Examples of these researchers are Kashirina et al. [13], Steinbach [14], Banihashemi et al. [15], Miyangaskari et al. [16], Sukono et al. [17], and Pandiangan et al. [18].
Based on the explanation of the brief gap, this study aims to formulate a value-at-risk (VaR)-based investment portfolio optimization model involving risk-free assets and company liabilities. The advantage of the model is that the maximum loss of the portfolio can be reduced due to the involvement of risk-free assets, and the maximum loss can be adjusted based on the size of the tolerable quantile. The optimization model designed has a quadratic programming form to accommodate the quadratic form of the VaR measure. In more detail, the solution of the model is calculated using the Lagrange multiplier method and the Khun–Tucker theorem, as calculated by Sukono et al. [17]. Formulating the portfolio optimization model is performed using vector and matrix equations to simplify the form of the equation further. Finally, a numerical illustration is conducted on this model using stock data from the mining and energy sectors as risky assets and bank deposits as risk-free assets. This study is expected to help quality control in forming investment portfolios, especially in the mining and energy sector stocks, for investors both organizationally and individually.

2. Literature Review

This section discusses references in portfolio model development. It includes a brief explanation of the general optimization model of the investment portfolio, the quadratic optimization model of the mean-VaR investment portfolio, and the asset-liability model.

2.1. The Explanation of Mathematical Notations

The essential mathematical notations used in this study are as follows:
(a)
t 0 = 0 represents the time an investor starts investing.
(b)
t 1 = 1 represents the time an investor ends investing.
(c)
There is one number of risk-free assets in the portfolio.
(d)
N 1 with N Z represents many risky assets.
(e)
V t 0 represents the total initial investment assets.
(f)
E V t 1 represents the expected total initial investment assets at time t 1 .
(g)
h 0 represents the amount of capital allocated to risk-free assets. We assume that the value of h 0 is determined in advance.
(h)
h k with k = 1 , 2 , , N represents the amount of capital allocated to the k -th risky asset.
(i)
γ k with k = 1 , 2 , , N represents the company’s liability to the k -th risky asset.
(j)
σ k 1 , k 2 with k 1 = 1 , 2 , , N and k 2 = 1 , 2 , , N represents the covariance between the k 1 -th and k 2 -th risky assets.
(k)
w 0 represents the weight of capital allocation on risk-free assets. It can be formulated mathematically as w 0 = h 0 V t 0 . Since h 0 is assumed to be determined in advance, the value of w 0 is also known.
(l)
w k , k = 1 , 2 , , N represents the weight of the capital allocation on the k -th risky asset. It can be formulated mathematically as w k = h k V t 0 .
(m)
μ 0 represents the return from the risk-free asset at time t 1 . We assume the value has been determined in advance.
(n)
r k represents the return from the k -th risky asset at t 1 , assumed to be distributed normal random variables with a mean of E r k = μ k and a variance of V a r r k = σ k 2 .
(o)
R represents the total return of the investment portfolio at time t 1 .
(p)
R represents the total return of the investment portfolio without the return of a risk-free asset at time t 1 .
(q)
B t i , i = 0 , 2 represents the price of the risk-free asset at time t i .
(r)
S t 0 k , k = 1 , 2 , , N represents the initial spot price of the k -th risk asset. The spot price is the asset price when the buyer and seller carry out the transaction [19].
(s)
E S t 1 k , k = 1 , 2 , , N represents the expected spot price of the k -th asset at time t 1 .
(t)
σ k , L represents the covariance between the k -th risky asset and liability.

2.2. The General Optimization Model of the Investment Portfolio

The total assets at time t 0 at time t 1 from the portfolios are, respectively, as follows [20,21,22]:
V t 0 = j = 0 N h j .
Suppose that V t 0 = 1 . Equation (1) can be formulated in the form of weights by substituting h j = w j V t 0 , j = 0 ,   1 ,   2 ,   ,   N . This substitution results in the following equations:
j = 0 N w j = 1 .
Equation (2) shows that the total weight of the assets in the portfolio is 1. In this study, w 0 has been determined first so that we are left with determining w k with k = 1 ,   2 ,   ,   N . Therefore, the total amount of asset weights and returns from a portfolio without risk-free assets is as follows:
k = 1 N w k = 1 w 0 .
The mean return of the portfolio without risk-free assets is formulated as follows:
E R = k = 1 N w k E ( r k ) = k = 1 N w k μ k .
Meanwhile, the variance return of the portfolio without risk-free assets can be briefly stated as follows [23,24]:
V R = k 1 = 1 N k 2 = 1 N w k 1 w k 2 C o v ( r k 1 , r k 2 ) = k 1 = 1 N k 2 = 1 N w k 1 w k 2 σ k 1 , k 2 .
The forms of Equations (3)–(5) can, respectively, be expressed in the form of the vector and matrix multiplication as follows:
w T e = 1 w 0 ,   E R = w T μ , a n d   V R = w T Σ w ,
where e T = 1 1 1 , w T = w 1 w 2 w N , μ T = μ 1 μ 2 μ N , and Σ = σ k 1 , k 2 ; k 1 , k 2 = 1 ,   2 , ,   N , which is a negative semi-definite matrix [25].
In the face of loss from the investment portfolio, investors have different risk aversion. Risk aversion is investors’ tendency to establish certainty in investments. In other words, the smaller the investor’s risk aversion, the more courageous the investor is to face uncertainty and vice versa. Mathematically, risk aversion is expressed as a comparison between the relative constant and twice the initial capital V t 0 = 1 as follows [4]:
ρ = c 2 V t 0 = c 2 ,
where ρ represents the risk aversion and c represents the relative constant [26].
The mean-variance model for asset allocation in the investment portfolio is mathematically expressed as follows:
m a x . E R ρ V R = w T μ c 2 w T Σ w , s . t . w T e = 1 w 0 .
The Lagrange multiplier method can determine the solution of the model in Equation (7). This method results in the Lagrange function for Equation (7) expressed as follows [18]:
L w , λ 1 = w T μ c 2 w T Σ w + λ 1 w T e 1 + w 0 ,
where λ 1 represents the Lagrange multiplier. The solution w in Equation (7) is the one that maximizes (8). Briefly, the solution is as follows [6]:
w * = 2 c Σ 1 μ + λ 1 e ,
where
λ 1 * = c 2 1 w 0 e T Σ 1 μ e T Σ 1 e .
The Hessian matrix of Equation (8) must satisfy Equation (10) as follows to guarantee the solution w * in Equation (9) as a local maximum solution [27]:
x R N + 1 0 ,   x T H w * , λ 1 * x < 0
where
H w , λ 1 = 2 L w , λ 1 w 2 2 L w , λ 1 w λ 1 2 L w , λ 1 λ 1 w L w , λ 1 λ 1 2 .

2.3. Quadratic Optimization Model of the Mean-VaR Investment Portfolio

This section briefly describes the quadratic optimization model of the mean-VaR investment portfolio. Value-at-risk (VaR) measures the maximum potential loss that may occur from a portfolio at a certain level of confidence under normal market conditions [28]. Sukono et al. [29] stated that the VaR of the portfolio without risk-free assets is expressed as follows:
V a R α ( R ) = w T μ + z α w T Σ w 1 2 ,
where z α represents the ( 1 α )-th percentile of the standard normal distribution. The value of α can be adjusted according to the investor’s expectation. The greater the value of α , the smaller the maximum loss from the investment portfolio desired by the investor and vice versa.
The quadratic optimization model of the mean-VaR investment portfolio formulates the objective function as a maximization problem of subtraction returns and multiplying risk aversion with VaR. Mathematically, this is written as follows [22]:
m a x . E R ρ V a R α R = w T μ c 2 w T μ + z α w T Σ w 1 2 , = 1 c 2 w T μ c 2 z α w T Σ w 1 2 , s . t . w T e = 1 w 0 .
The Lagrange multiplier method can determine the solution of the model in Equation (12). This method results in the Lagrange function for Equation (12) expressed as follows [22]:
L w , λ 2 = 1 c 2 w T μ c 2 z α w T Σ w 1 2 + λ 2 ( w T e 1 + w 0 ) ,
where λ 2 represents the Lagrange multiplier. The solution w of Equation (12) is the one that maximizes (13). Briefly, the solution is as follows [22]:
w * = 1 w 0 Σ 1 1 c 2 μ + λ 2 e e T Σ 1 1 c 2 μ + λ 2 e .
In detail,
λ 2 * = b ± b 2 4 a c 2 a ,
where a = e T Σ 1 e , b = 2 1 c 2 μ T Σ 1 e , and c = 1 c 2 2 μ T Σ 1 μ c z α 2 2 . The Hessian matrix of Equation (13) must satisfy Equation (15) as follows to guarantee the solution w * in Equation (14) as a local maximum solution [22]:
x R N + 1 0 ,   x T H w * , λ 2 * x < 0
where
H w , λ 2 = 2 L w , λ 2 w 2 2 L w , λ 2 w λ 2 2 L w , λ 2 λ 2 w L w , λ 2 λ 2 2 .

2.4. Asset Liability Model

Suppose that s t i k , i = 0 ,   1 , represents the surplus of the k -th asset at time t i . It can be mathematically expressed as follows:
s t i k = S t i k L t i k = S t i 1 k 1 + μ k L t i 1 k 1 + μ L k ,
where L t i represents the k -th liabilities and assets at time t i and μ L k represents the return of the k -th liabilities. The return of the surplus from the k -th asset can be expressed in the following equation [30,31,32]:
μ S k = s t i k s t i 1 k S t i 1 k = S t i 1 k μ k L t i 1 k μ L k S t i 1 k = μ k L t i 1 k S t i 1 k μ L k = μ k f t i 1 μ L k ,
where f t i is the ratio between liabilities and capital assets at time t i . Meanwhile, the variance of return surplus is stated as follows:
σ S k 2 = σ k 2 2 f 0 σ k , L + 1 f 0 2 σ L k 2 ,
where σ k , L represents the covariance between the k -th return assets and the liability return and σ L k 2 represents the variance of the liability return of k -th asset [33].

3. Materials and Methods

This section briefly describes the materials and methods, which contain the materials and steps in this study.

3.1. Materials

This study had two activities: model development and numerical illustration. The model development was carried out using the model of Pandiangan et al. [18], which formulated an investment portfolio optimization model in the form of mean-VaR involving risk-free assets. This development was in the form of adding the liability variable from the asset owner company. In addition, the model development was carried out using vector and matrix forms, as performed by Sukono et al. [17]. Then, the material for numerical illustration was stock data in the mining and energy sectors, which are traded on the Indonesian capital market through the Indonesia Stock Exchange (IDX). Note that the model can not only be applied to IDX’s stock data but also to other capital markets worldwide. The data are historical monthly return data from 11 stocks obtained through the website http://finance.yahoo.com (accessed on 1 March 2024). The period considered was 1 January 2017 to 1 December 2020. Of the mining and energy stocks, 11 selected stocks had the best performance. In addition to data on 11 stocks, we sued data on risk-free assets in deposits determined by the Central Bank of Indonesia. In this study, the company’s liability data were also used. However, because the company’s liability data are not open, the liability data used were simulation data [18,22].

3.2. Methods

The steps for developing numerical models and illustrations are described in this section. The stages of model development were carried out by examining them in depth, especially in the research of Pandiangan et al. [18]. This study was conducted to develop a mean-VaR quadratic investment portfolio optimization model based on investor risk aversion considering the risk-free assets and liabilities of the company. The formulation of the portfolio weight model here was carried out using a vector and matrix equation approach.
The stages of numerical illustration were carried out to show the application of the model. The data used were the return of 11 stock assets in the mining and energy sectors, returns on company liabilities, and risk-free assets, as described in Section 3.1. The stages of data analysis were broadly carried out in the following way:
(1)
The descriptive analysis was carried out first on the 11 stock return data. The descriptive statistics analyzed were stock code, mean value, variance, covariance, and standard deviation. The mean of each asset was then used to form the vector of return mean μ . Then, the variance and covariance values formed the Σ covariance matrix.
(2)
Next was to generate monthly liability return data from the company via simulation. The covariance between stock asset returns and liability returns was determined first. The covariance value was then used to form the covariance vector of asset and liability returns.
(3)
Next was the stage of determining the optimum asset weight in the investment portfolio based on the mean vector, covariance matrix, and covariance vector between assets and liabilities obtained. In this stage, investors’ risk aversion values were determined via simulation.
(4)
Based on the results of the optimization of asset weights in the investment portfolio, the formation of an efficient portfolio surface graph was carried out. This graph illustrated the feasible points for a rational investor to invest accordingly.

4. Results and Discussions

4.1. Quadratic Optimization Modeling for Mean-VaR Investment Portfolio with Risk-Free Assets and Liabilities

In this section, the development of a quadratic optimization model of the mean-VaR investment portfolio, which includes risk-free assets and liabilities, is designed. The model is inclusive for a single period, e.g., daily, monthly, quarterly, or yearly. Suppose that the covariance vector between asset and liability returns is expressed as γ T = ( γ 1 ,   γ 2 , ,   γ N ) with γ k = 1 f 0 σ k , L . The quadratic optimization model for the mean-VaR investment portfolio without the risk-free asset, involving risk-free assets and liabilities, is expressed as follows:
m a x . w T μ + w T γ c 2 w T μ + z α w T Σ w 1 2 , = 1 c 2 w T μ + w T γ c 2 z α w T Σ w 1 2 , s . t . w T e = 1 w 0 .
The determination of the solution to the problem in Equation (19) is carried out using the Lagrange multiplier method. The Lagrange function of the problem in Equation (19) is expressed as follows:
L w , λ 3 = 1 c 2 w T μ + w T γ c 2 z α w T Σ w 1 2 + λ 3 w T e 1 + w 0 ,
where λ 3 represents the Lagrange multiplier. Based on the Khun–Tucker theorem, the solution of Equation (20) is carried out based on necessary and sufficient conditions [18,33]. The necessary conditions for the solution to Equation (20) fulfill the following equation:
L w , λ 3 w = 1 c 2 μ + γ c 2 z α Σ w w T Σ w 1 2 + λ 3 e = 0 ,
L w , λ 3 λ 3 = w T e 1 + w 0 = 0 .
First, determine the vector w as follows:
1 c 2 μ + γ c 2 z α Σ w w T Σ w 1 2 + λ 3 e = 0 ,
c 2 z α Σ w w T Σ w 1 2 = 1 c 2 μ + γ + λ 3 e .
Multiply both sides of Equation (21) by 2 c z α Σ 1 so that this results in the following equation:
w w T Σ w 1 2 = 2 c z α Σ 1 1 c 2 μ + γ + λ 3 e ,
Then, multiply both sides of Equation (22) by e T so that this results in the following equation:
1 w 0 w T Σ w 1 2 = 2 c z α e T Σ 1 1 c 2 μ + γ + λ 3 e ,
w T Σ w 1 2 = c z α 2 1 w 0 e T Σ 1 1 c 2 μ + γ + λ 3 e .
Substitute Equation (23) into Equation (22) so that this results in the following equation:
w w T Σ w 1 2 = 2 c z α Σ 1 1 c 2 μ + γ + λ 3 e ,
w = w T Σ w 1 2 2 c z α Σ 1 1 c 2 μ + γ + λ 3 e ,
w = c z α 2 1 w 0 e T Σ 1 1 c 2 μ + γ + λ 3 e 2 c z α Σ 1 1 c 2 μ + γ + λ 3 e ,
w * = ( 1 w 0 ) Σ 1 1 c 2 μ + γ + λ 3 e e T Σ 1 1 c 2 μ + γ + λ 3 e .
Next is the determination of the solution λ 3 . Multiply both sides of Equation (24) by w T Σ so that this results in the following equation:
w T Σ w = ( 1 w 0 ) w T 1 c 2 μ + γ + λ 3 e e T Σ 1 1 c 2 μ + γ + λ 3 e ,
or
w T Σ w = ( 1 w 0 ) 1 c 2 μ T + γ T + λ 3 e T w e T Σ 1 1 c 2 μ + γ + λ 3 e .
From Equation (23), obtain the following equation:
1 w 0 e T Σ 1 1 c 2 μ + γ + λ 3 e = w T Σ w 1 2 2 c z α .
Substitute Equation (26) into (25) to obtain the following equation:
w T Σ w = w T Σ w 1 2 2 c z α 1 c 2 μ T + γ T + λ 3 e T w ,
c z α 2 w T Σ w 1 2 = 1 c 2 μ T + γ T + λ 3 e T w .
Substitute Equations (23) and (24) into Equation (27) so that this results in the following equation:
c z α 2 c z α 2 1 w 0 e T Σ 1 1 c 2 μ + γ + λ 3 e = 1 c 2 μ T + γ T + λ 3 e T 1 w 0 Σ 1 1 c 2 μ + γ + λ 3 e e T Σ 1 1 c 2 μ + γ + λ 3 e ,
c z α 2 2 = 1 c 2 μ T + γ T + λ 3 e T Σ 1 1 c 2 μ + γ + λ 3 e ,
e T Σ 1 e λ 3 2 + 2 1 c 2 μ T + γ T Σ 1 e λ 3 + 1 c 2 μ T + γ T Σ 1 1 c 2 μ + γ c z α 2 2 = 0 .
Equation (28) is a quadratic equation, so the possible values of λ 3 > 0 are as follows:
λ 3 = q ± q 2 4 p r 2 p ,
where p = e T Σ 1 e , q = 2 1 c 2 μ T + γ T Σ 1 e , and r = 1 c 2 μ T + γ T Σ 1 1 c 2 μ + γ c z α 2 2 .
After the weight allocation of funds from each asset is obtained, the next step is to examine the terms of sufficient solutions as the maximum solution. It is calculated by examining the Hessian matrix of the Lagrange function in Equation (20). Equation (19) is converted into a minimization problem via the dual method. Then, the solution w * to Equation (24) is called the global maximum solution of Equation (20) if the Hessian matrix is positive definite. The Hessian matrix of Equation (20) is expressed as follows:
H w , λ 3 = 2 L w , λ 3 w 2 2 L w , λ 3 w λ 3 2 L w , λ 3 λ 3 w L w , λ 3 λ 3 2 = c z α 2 w T Σ w 1 2 Σ Σ w w T Σ w T Σ w e e T 0 .
The H ( w , λ 3 ) matrix is called positive definite if x R N + 1 0 ,   x T H ( w , λ 3 ) x > 0 . Since the values of c and z α are positive, it is evident that the problem in Equation (20) is convex [34,35,36].

4.2. Numerical Illustration Results

4.2.1. Descriptive Statistics of Stock Data

This section analyzes descriptive statistics of the stock data used in the numerical illustration. Table 1 gives the stock code, mean, variance, and standard deviation of the 11 stock data.
Table 1 shows that of the 11 stocks, the most significant mean is in PTBA, while the most minor mean is obtained from RUIS. The values from Table 1 are then used for the mean vector, vector one, and the covariance matrix. The vector of the mean return of 11 stocks is expressed as follows:
μ = 0.0199 0.0314 0.0374 0.0214 0.0272 0.0469 0.0354 0.0135 0.0855 0.0340 0.0135
Furthermore, the covariance between 11 stock asset returns must first be determined to form a covariance matrix. We formed it with the help of Microsoft Excel software, and the results are given in Table A1 (See Appendix A). The next step is determining the covariance vector between stock asset returns and company liability returns. Monthly liability data are calculated with a simulation. To generate monthly liability data, we assume that it is normally distributed based on the parameters of each stock asset liability return. Furthermore, the covariance value between stock asset returns and liability returns is determined, and the results are given as follows:
γ = 0.0016 0.1449 0.0009 0.0020 0.0268 0.1198 0.6424 0.0024 0.5944 0.0127 0.0025
In this study, risk-free assets are deposits with a monthly interest rate set by the Central Bank of Indonesia, namely μ 0 = 0.07 per year 0.07 12 per month. Then, the weight allocation for risk-free assets w 0 is assumed to be 0.5.

4.2.2. Determination of Optimum Weight Allocation

In this section, the determination of the optimum weight allocation for the investment portfolio uses Equation (24). Risk aversion is determined with a simulation [37,38] whose value is 5.1 < c 8.2 with an increase of Δ c = 0.05 . In other words, c is in the set of 5.1 ,   5.15 ,   , 8.15 ,   8.2 . This interval is chosen because the weight allocation of each stock obtained is positive. Calculations are carried out using the Matlab version 2007a and the Microsoft Excel version 2403. The optimum weighting results for each of these shares are given in Table A2 (See Appendix B).
Table A2 shows that with the value of c = 5.1 , there is a weight of w 9 = 0.0000004 . It shows the existence of short sales. This analysis assumes that the capital market authority does not permit short sales, so this is not analyzed further. Table A2 shows that the greater the value of risk aversion, the smaller the mean and VaR of the portfolio return. Using the mean and VaR of the portfolio returns in Table A2, the efficient surface graph given in Figure 1 is presented.
Figure 1 informs investors that investments in 11 stocks worth making are at points along the efficient surface line. In this study, an efficient portfolio surface graph is formed for the constant level of investor risk aversion in the interval 5.1 < c 8.2 . In that interval, the minimum values of VaR and mean of the portfolio returns are 0.013797 and 0.0102941, respectively. The VaR and the mean return are obtained when the risk aversion investor is c = 8.2 . Meanwhile, the maximum values of VaR and the mean of the portfolio return are 0.013827 and 0.0103316, respectively. The VaR and mean values are obtained when the risk aversion investor is c = 5.2 .

4.3. Discussion

4.3.1. Relationship Analyses

The relationship between risk aversion and mean of the portfolio return is the first thing to be analyzed in this section. Using the values of risk aversion and mean of the portfolio return in Table A2 (see Appendix B), the relationship between the two can be visualized in a Cartesian diagram. The visualization of the relationship between risk aversion and mean of the portfolio return is presented in Figure 2.
Figure 2 shows that risk aversion and mean of the portfolio return have a negative relationship. It means that the greater the risk aversion of the investor, the smaller the mean of the portfolio return becomes and vice versa. It is logical because if investors avoid uncertainty (this means they have high-risk aversion), the mean of the portfolio return that may be obtained is also very small [39,40].
Next, the relationship between risk aversion and the VaR of the portfolio is analyzed. Using the values of risk aversion and VaR portfolio in Table A2 (see Appendix B), the relationship between the two can be visualized in a Cartesian diagram. The visualization of the relationship between risk aversion and portfolio VaR is presented in Figure 3.
Figure 3 also shows that risk aversion has a relationship that is not in line with the VaR of the portfolio. In other words, this is similar to the relationship between risk aversion and mean of the portfolio return, where the greater the risk aversion of the investor, the smaller the VaR of the portfolio return and vice versa. Practically speaking, if investors are highly uncertainty averse (this means the risk aversion is large), their VaR of portfolio returns will also be small.
Next is a discussion on choosing the optimum portfolio. Optimum portfolio selection can be reviewed by comparing the largest comparison of return and VaR. After this, the comparison is referred to as the ratio. The smallest ratio value indicates that the portfolio has the lowest efficiency level, while the largest ratio value indicates that the portfolio has a high efficiency level (most optimum). To further illustrate this point, we present a visualization of the relationship between portfolio VaR and its ratio in Figure 4.
Figure 4 shows that the smallest ratio value is 0.7448. This ratio is obtained from the mean and VaR of portfolio returns 0.0102964 and 0.0137975, respectively. The portfolio is then said to be the portfolio with the lowest efficiency level. Then, the largest ratio value is 0.7486182. This ratio is obtained from the mean and VaR of the portfolio returns 0.0103316 and 0.0138270, respectively. This portfolio is then said to be the optimum portfolio. Coincidentally, this optimum portfolio occurs at the maximum portfolio for the value of investors’ risk aversion in the interval 5.1 < c 8.2 .
Based on the discussion in the previous sections, a summary of the risk aversion intervals, minimum portfolio, maximum portfolio, optimum portfolio, and composition of weights on the optimum portfolio are given in Table 2.
The summary in Table 2 is expected to be a reference for investors, individuals, and organizations in making investment decisions, especially in the 11 mining and energy sector stocks, including risk-free assets and company liabilities in forming investment portfolios.

4.3.2. Comparation with Other Models

This subsection compares the accuracy of the mean of the optimum portfolio returns of the introduced model (we call it model I) and previous models. The previous models were explained in Section 2.2 (we call it model II) and 2.3 (we call it model III). To measure this accuracy, we first conducted repeated calculations of the mean of the portfolio return from each model for different observation periods. The initial observation period is 1 January 2017 to 1 December 2020, and the subsequent observation period was an additional one month from the previous one. The number of observation periods is thirteen. Then, we determine the optimum portfolio for each observation period. Each mean of the return from the optimum portfolio is used as a forecast of the mean of the portfolio return from the portfolio one month ahead. After that, we compare each forecast with the actual return from each model. We then measure the accuracy of the forecast using the absolute error. The accuracy measurement results show that model I tends to have a smaller absolute error than models II and III in each observation period. It is also shown by the mean absolute error (MAE) of thirteen predictions from the three models, where the MAE from model I is smaller than models II and III. Although this is not analytical evidence, it can show that the model introduced in this study is reliable practically. In summary, the results of the accuracy comparison can be seen in Table 3.

5. Conclusions

This study designs the development of a quadratic optimization model for the mean-VaR investment portfolio, which involves risk-free assets and liabilities. The model can determine the combination of capital allocation weights on each stock asset in forming an investment portfolio. The optimum portfolio weight determination model is expressed using vector and matrix equations, which can simplify and speed up computing with specific software, i.e., the Matlab version 2007a and the Microsoft Excel version 2403.
The numerical illustration is completed to illustrate the use of the model. Note that the model can be illustrated in various capital markets worldwide. In this study, the illustration is based on data on 11 stock assets of Indonesia’s mining and energy sectors. Then, risk-free assets are time deposits whose interest rates are adjusted to Bank Indonesia’s policy. Then, the company’s liabilities are used as simulation data. The numerical illustration shows that risk aversion and the mean of the portfolio return have a negative relationship. Therefore, it means that the greater the risk aversion of the investor, the smaller the mean of the portfolio return becomes and vice versa. This relationship is similar to the relationship between risk aversion and the VaR of the portfolio return. It is logical because if investors avoid uncertainty (this means that their risk aversion is significant), their mean and VaR of portfolio returns will be small and vice versa. We also compare the accuracy of the mean of the optimum portfolio returns of the introduced model and previous models. This comparation shows that the mean of the return from the model introduced in this study has the best accuracy based on its MAPE. Therefore, the model introduced in this study is reliable practically.
This study has its shortcomings. The weaknesses include the fact that the model development is carried out for only a single period, so it cannot accommodate the wishes of investors who invest for multiple periods. Therefore, future research needs further development for models with multi-period time indices. In addition, weaknesses occur in the analysis of numerical illustrations, where the monthly liability data used are the result of a simulation.

Author Contributions

Conceptualization, S. and R.A.I.; methodology, P.L.B.G.; software, R.A.I.; validation, M.D.J., S. and M.M.; formal analysis, A.S.; investigation, S., M.M. and R.A.I.; resources, S.; data curation, R.; writing—original draft preparation, S. and P.L.B.G.; writing—review and editing, R.A.I.; visualization, M.D.J. and R.; supervision, S.; project administration, A.S.; funding acquisition, S. All authors have read and agreed to the published version of the manuscript.

Funding

DRPM Universitas Padjadjaran, number: 1549/UN6.3.1/PT.00/2023 and UniSZA/2023/FELLOWPENELIDIKAN(011).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

Thank you to Padjadjaran University for the Padjadjaran University Internal Research Grant, Fiscal Year 2023, the “ACADEMIC LEADERSHIP GRANT (ALG)” program under Sukono, MM., M.Si., with contract number: 1549/UN6.3.1/PT.00/2023. Also, thanks to Universiti Sultan Zainal Abidin for their cooperation in the research fellow scheme with reference code: UniSZA/2023/FELLOWPENELIDIKAN(011).

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Table A1. Covariance matrix of 11 Stocks.
Table A1. Covariance matrix of 11 Stocks.
BSSRBYANCITAHRUMMBAPMDKAMEDCPSABPTBAPTRORUIS
BSSR0.01860.0028−0.00030.00680.0031−0.00030.0114−0.00550.00370.00300.0011
BYAN0.00280.02420.00560.00340.01300.00010.0018−0.00050.00320.0032−0.0017
CITA−0.00030.00560.0288−0.00690.00000.0101−0.00750.00020.0016−0.0061−0.0035
HRUM0.00680.0034−0.00690.02830.0066−0.00130.01460.0053−0.00220.00980.0055
MBAP0.00310.01300.00000.00660.01990.00020.0046−0.0015−0.00770.00600.0022
MDKA−0.00030.00010.0101−0.00130.00020.01410.00440.0053−0.00030.0032−0.0012
MEDC0.01140.0018−0.00750.01460.00460.00440.04870.00320.02510.01620.0053
PSAB−0.0055−0.00050.00020.0053−0.00150.00530.00320.0270−0.00350.0105−0.0007
PTBA0.00370.00320.0016−0.0022−0.0077−0.00030.0251−0.00350.35040.0341−0.0004
PTRO0.00300.0032−0.00610.00980.00600.00320.01620.01050.03410.02390.0003
RUIS0.0011−0.0017−0.00350.00550.0022−0.00120.0053−0.0007−0.00040.00030.0184

Appendix B

Table A2. Optimization Process with Risk-Free Assets w 0 = 0.5 and Liabilities.
Table A2. Optimization Process with Risk-Free Assets w 0 = 0.5 and Liabilities.
c w 1 w 2 w 3 w 4 w 5 w 6 w 7 w 8 w 9 w 10 w 11 w T e E R V a R 0.05 R Ratio
5.100.00600.00060.05520.16200.00490.00280.00020.01080.00000.03690.22060.50.01030.01380.7486
5.150.00600.00060.05520.16210.00480.00270.00020.01080.00000.03700.22060.50.01030.01380.7485
5.200.00600.00060.05510.16210.00480.00270.00020.01080.00000.03700.22060.50.01030.01380.7484
5.250.00610.00060.05500.16220.00470.00270.00020.01080.00000.03710.22060.50.01030.01380.7483
5.300.00610.00060.05500.16220.00470.00270.00020.01080.00000.03710.22060.50.01030.01380.7483
5.350.00610.00060.05490.16230.00470.00270.00020.01080.00000.03710.22060.50.01030.01380.7481
5.400.00610.00060.05480.16230.00460.00270.00020.01080.00000.03720.22060.50.01030.01380.7481
5.450.00610.00060.05480.16240.00460.00270.00020.01080.00000.03720.22070.50.01030.01380.7480
5.500.00610.00060.05470.16240.00460.00270.00020.01080.00000.03720.22070.50.01030.01380.7479
5.550.00610.00060.05460.16250.00450.00270.00020.01080.00000.03720.22070.50.01030.01380.7477
5.600.00610.00060.05460.16250.00450.00270.00020.01080.00000.03730.22070.50.01030.01380.7477
5.650.00620.00060.05450.16260.00450.00270.00020.01080.00000.03730.22070.50.01030.01380.7476
5.700.00620.00060.05450.16260.00440.00270.00020.01080.00000.03730.22070.50.01030.01380.7476
5.750.00620.00060.05440.16270.00440.00270.00020.01080.00000.03740.22070.50.01030.01380.7475
5.800.00620.00060.05440.16270.00440.00270.00020.01080.00000.03740.22070.50.01030.01380.7475
5.850.00620.00060.05430.16280.00430.00270.00020.01070.00000.03740.22080.50.01030.01380.7473
5.900.00620.00050.05430.16280.00430.00270.00020.01070.00000.03740.22080.50.01030.01380.7472
5.950.00620.00050.05420.16280.00430.00270.00020.01070.00000.03750.22080.50.01030.01380.7472
6.000.00620.00050.05420.16290.00430.00270.00020.01070.00000.03750.22080.50.01030.01380.7471
6.050.00620.00050.05410.16290.00420.00270.00020.01070.00010.03750.22080.50.01030.01380.7470
6.100.00620.00050.05410.16300.00420.00270.00020.01070.00010.03750.22080.50.01030.01380.7469
6.150.00630.00050.05400.16300.00420.00270.00020.01070.00010.03760.22080.50.01030.01380.7469
6.200.00630.00050.05400.16300.00410.00270.00020.01070.00010.03760.22080.50.01030.01380.7469
6.250.00630.00050.05390.16310.00410.00270.00020.01070.00010.03760.22080.50.01030.01380.7467
6.300.00630.00050.05390.16310.00410.00270.00020.01070.00010.03760.22090.50.01030.01380.7467
6.350.00630.00050.05380.16310.00410.00270.00020.01070.00010.03770.22090.50.01030.01380.7466
6.400.00630.00050.05380.16320.00400.00270.00020.01070.00010.03770.22090.50.01030.01380.7465
6.450.00630.00050.05370.16320.00400.00270.00020.01070.00010.03770.22090.50.01030.01380.7465
6.500.00630.00050.05370.16330.00400.00270.00020.01070.00010.03770.22090.50.01030.01380.7464
6.550.00630.00050.05370.16330.00400.00270.00020.01070.00010.03770.22090.50.01030.01380.7464
6.600.00630.00050.05360.16330.00400.00270.00020.01070.00010.03780.22090.50.01030.01380.7464
6.650.00630.00050.05360.16340.00390.00260.00020.01070.00010.03780.22090.50.01030.01380.7462
6.700.00630.00050.05350.16340.00390.00260.00020.01070.00010.03780.22090.50.01030.01380.7461
6.750.00640.00050.05350.16340.00390.00260.00020.01070.00010.03780.22090.50.01030.01380.7461
6.800.00640.00050.05350.16350.00390.00260.00020.01070.00010.03780.22090.50.01030.01380.7461
6.850.00640.00050.05340.16350.00380.00260.00020.01070.00010.03790.22100.50.01030.01380.7459
6.900.00640.00050.05340.16350.00380.00260.00020.01060.00010.03790.22100.50.01030.01380.7459
6.950.00640.00050.05340.16350.00380.00260.00020.01060.00010.03790.22100.50.01030.01380.7459
7.000.00640.00050.05330.16360.00380.00260.00020.01060.00010.03790.22100.50.01030.01380.7458
7.050.00640.00050.05330.16360.00380.00260.00020.01060.00010.03790.22100.50.01030.01380.7458
7.100.00640.00050.05320.16360.00370.00260.00020.01060.00010.03800.22100.50.01030.01380.7457
7.150.00640.00050.05320.16370.00370.00260.00020.01060.00010.03800.22100.50.01030.01380.7457
7.200.00640.00050.05320.16370.00370.00260.00020.01060.00010.03800.22100.50.01030.01380.7457
7.250.00640.00050.05310.16370.00370.00260.00020.01060.00010.03800.22100.50.01030.01380.7456
7.300.00640.00050.05310.16370.00370.00260.00020.01060.00010.03800.22100.50.01030.01380.7456
7.350.00640.00050.05310.16380.00370.00260.00020.01060.00010.03800.22100.50.01030.01380.7455
7.400.00640.00050.05300.16380.00360.00260.00020.01060.00010.03810.22100.50.01030.01380.7454
7.450.00640.00050.05300.16380.00360.00260.00020.01060.00010.03810.22100.50.01030.01380.7454
7.500.00650.00050.05300.16380.00360.00260.00020.01060.00010.03810.22110.50.01030.01380.7454
7.550.00650.00050.05300.16390.00360.00260.00020.01060.00010.03810.22110.50.01030.01380.7453
7.600.00650.00050.05290.16390.00360.00260.00020.01060.00010.03810.22110.50.01030.01380.7453
7.650.00650.00050.05290.16390.00360.00260.00020.01060.00010.03810.22110.50.01030.01380.7453
7.700.00650.00050.05290.16390.00350.00260.00020.01060.00010.03820.22110.50.01030.01380.7453
7.750.00650.00050.05280.16400.00350.00260.00020.01060.00010.03820.22110.50.01030.01380.7451
7.800.00650.00050.05280.16400.00350.00260.00020.01060.00010.03820.22110.50.01030.01380.7451
7.850.00650.00050.05280.16400.00350.00260.00020.01060.00010.03820.22110.50.01030.01380.7451
7.900.00650.00050.05280.16400.00350.00260.00010.01060.00010.03820.22110.50.01030.01380.7451
7.950.00650.00050.05270.16410.00350.00260.00010.01060.00010.03820.22110.50.01030.01380.7450
8.000.00650.00050.05270.16410.00340.00260.00010.01060.00010.03820.22110.50.01030.01380.7450
8.050.00650.00050.05270.16410.00340.00260.00010.01060.00010.03830.22110.50.01030.01380.7450
8.100.00650.00050.05260.16410.00340.00260.00010.01060.00010.03830.22110.50.01030.01380.7449
8.150.00650.00050.05260.16410.00340.00260.00010.01060.00010.03830.22110.50.01030.01380.7449
8.200.00650.00050.05260.16420.00340.00260.00010.01060.00010.03830.22110.50.01030.01380.7448

References

  1. Siregar, B.; Pangruruk, F.A. Portfolio Optimization Based on Clustering of Indonesia Stock Exchange: A Case Study of Index LQ45. Indones. J. Bus. Anal. 2021, 1, 59–70. [Google Scholar] [CrossRef]
  2. Wang, Y.; Chen, Y.; Liu, Y. Modeling Portfolio Optimization Problem by Probability-Credibility Equilibrium Risk Criterion. Math. Probl. Eng. 2016, 2016, 9461021. [Google Scholar] [CrossRef]
  3. Baumann, P.; Trautmann, N. Portfolio-Optimization Models for Small Investors. Math. Methods Oper. Res. 2013, 77, 345–356. [Google Scholar] [CrossRef]
  4. Sukono; Rosadi, D.; Maruddani, D.A.I.; Ibrahim, R.A.; Johansyah, M.D. Mechanisms of Stock Selection and Its Capital Weighing in the Portfolio Design Based on the MACD-K-Means-Mean-VaR Model. Mathematics 2024, 12, 174. [Google Scholar] [CrossRef]
  5. Hidayat, Y.; Purwandari, T.; Sukono; Prihanto, I.G.; Hidayana, R.A.; Ibrahim, R.A. Mean-Value-at-Risk Portfolio Optimization Based on Risk Tolerance Preferences and Asymmetric Volatility. Mathematics 2023, 11, 4761. [Google Scholar] [CrossRef]
  6. Purwandari, T.; Riaman; Hidayat, Y.; Sukono; Ibrahim, R.A.; Hidayana, R.A. Selecting and Weighting Mechanisms in Stock Portfolio Design Based on Clustering Algorithm and Price Movement Analysis. Mathematics 2023, 11, 4151. [Google Scholar] [CrossRef]
  7. Novais, R.G.; Wanke, P.; Antunes, J.; Tan, Y. Portfolio Optimization with a Mean-Entropy-Mutual Information Model. Entropy 2022, 24, 369. [Google Scholar] [CrossRef] [PubMed]
  8. Markowitz, H. The Utility of Wealth. J. Political Econ. 1952, 60, 151–158. [Google Scholar] [CrossRef]
  9. Sharpe, W.F. A Simplified Model for Portfolio Analysis. Manag. Sci. 1963, 9, 277–293. [Google Scholar] [CrossRef]
  10. Ivanova, M.; Dospatliev, L. Application of Markowitz Portfolio Optimization on Bulgarian Stock Market from 2013 to 2016. Int. J. Pure Appl. Math. 2018, 117. [Google Scholar] [CrossRef]
  11. Kulali, I. Portfolio Optimization Analysis with Markowitz Quadratic Mean-Variance Model. Eur. J. Bus. Manag. 2016, 8, 73–79. [Google Scholar]
  12. Zavera, I.C. Application of Markowitz Model on Romanian Stock Market. HOLISTICA J. Bus. Public Adm. 2017, 8, 97–103. [Google Scholar] [CrossRef]
  13. Kashirina, I.L.; Azarnova, T.V.; Bondarenko, Y.V.; Shchepina, I.N. Modeling and Optimization of Assets Portfolio with Consideration of Profits Reinvestment. Glob. J. Pure Appl. Math. 2016, 12, 2023–2033. [Google Scholar]
  14. Steinbach, M.C. Markowitz Revisited: Mean-Variance Models in Financial Portfolio Analysis. Soc. Ind. Appl. Math. 2001, 43, 31–85. [Google Scholar] [CrossRef]
  15. Banihashemi, S.; Azarpour, A.M.; Navvabpour, H. Portfolio Optimization by Mean-Value at Risk Framework. Appl. Math. Inf. Sci. 2016, 10, 1935–1948. [Google Scholar] [CrossRef]
  16. Ghaoui, L.E.; Oks, M.; Oustry, F. Worst-Case Value-At-Risk and Robust Portfolio Optimization: A Conic Programming Approach. Oper. Res. 2003, 51, 543–556. [Google Scholar] [CrossRef]
  17. Sukono; Lesmana, E.; Napitupulu, H.; Hidayat, Y.; Saputra, J.; Ghazali, P.L.B. Mean-VAR Portfolio Optimisations: An Application of Multiple Index Models with Non-Constant Volatility and Long Memory Effects. Int. J. Innov. Creat. Chang. 2019, 9, 364–381. [Google Scholar]
  18. Pandiangan, N.; Sukono; Hasbullah, E.S. Comparison of Quadratic Investment Portfolio on Five Stocks of Mining Companies with Risk Free Assets and without Risk Free Assets. J. Phys. Conf. Ser. 2021, 1722, 012069. [Google Scholar] [CrossRef]
  19. Sun, J. The Application of Markowitz Model and Index Model on Portfolio Optimization. In 2022 7th International Conference on Social Sciences and Economic Development (ICSSED 2022); Atlantis Press: Amsterdam, The Netherlands, 2022. [Google Scholar]
  20. Chaweewanchon, A.; Chaysiri, R. Markowitz Mean-Variance Portfolio Optimization with Predictive Stock Selection Using Machine Learning. Int. J. Financ. Stud. 2022, 10, 64. [Google Scholar] [CrossRef]
  21. Chen, P.; Lezmi, E.; Roncalli, T.; Xu, J. A Note on Portfolio Optimization with Quadratic Transaction Costs. SSRN Electron. J. 2019. [Google Scholar] [CrossRef]
  22. Pandiangan, N.; Sukono, S.; Hasbullah, E.S. Quadratic Investment Portfolio Based on Value-At-Risk with Risk-Free Assets: For Stocks of the Mining and Energy Sector. Int. J. Energy Econ. Policy 2021, 11, 175–184. [Google Scholar] [CrossRef]
  23. Stowe, D.L. Portfolio Mathematics with General Linear and Quadratic Constraints. J. Math. Financ. 2019, 09, 675–690. [Google Scholar] [CrossRef]
  24. Bai, Z.; Liu, H.; Wong, W.-K. Making Markowitz’s Portfolio Optimization Theory Practically Useful. SSRN Electron. J. 2010. [Google Scholar] [CrossRef]
  25. Hasbullah, E.S.; Halim, N.B.A.S.; Putra, A.S.; Bon, A.T. Mean-Variance Portfolio Optimization on Islamic Stocks by Using Non Constant Mean and Volatility Models and Genetic Algorithm. Int. J. Eng. Technol. 2018, 7, 366. [Google Scholar] [CrossRef]
  26. Van Mieghem, J.A. Capacity Management, Investment, and Hedging: Review and Recent Developments. Manuf. Serv. Oper. Manag. 2003, 5, 269–302. [Google Scholar] [CrossRef]
  27. Purwani, S.; Ibrahim, R.A. Using Simple Fixed-Point Iterations to Estimate Generalized Pareto Distribution Parameters. IAENG Int. J. Appl. Math. 2024, 54, 194–204. [Google Scholar]
  28. Kumari, S.K.; Kumar, P.; Priya, J.; Surya, S.; Bhurjee, A.K. Mean-Value at Risk Portfolio Selection Problem Using Clustering Technique: A Case Study. In AIP Conference Proceedings; AIP Publishing: College Park, MD, USA, 2019; p. 020178. [Google Scholar]
  29. Sukono; Lesmana, E.; Johansyah, M.D.; Napitupulu, H.; Hidayat, Y.; Purwandari, T. Value-at-Risk and Minimum Variance in the Investment Portfolio with Non Constant Volatility. Int. J. Recent Technol. Eng. 2019, 8, 197–202. [Google Scholar] [CrossRef]
  30. Zhu, H.-N.; Zhang, C.-K.; Jin, Z. Continuous-Time Mean-Variance Asset-Liability Management with Stochastic Interest Rates and Inflation Risks. J. Ind. Manag. Optim. 2020, 16, 813–834. [Google Scholar] [CrossRef]
  31. Zhou, Z.; Zeng, X.; Xiao, H.; Ren, T.; Liu, W. Multiperiod Portfolio Optimization for Asset-Liability Management with Quadratic Transaction Costs. J. Ind. Manag. Optim. 2018, 15, 1493–1515. [Google Scholar] [CrossRef]
  32. Sheng, D.-L.; Shen, P. Portfolio Optimization with Asset-Liability Ratio Regulation Constraints. Complexity 2020, 2020, 1435356. [Google Scholar] [CrossRef]
  33. Platanakis, E.; Sutcliffe, C. Asset–Liability Modelling and Pension Schemes: The Application of Robust Optimization to USS. Eur. J. Financ. 2017, 23, 324–352. [Google Scholar] [CrossRef]
  34. Pan, J.; Zhang, Z.; Zhou, X. Optimal Dynamic Mean-Variance Asset-Liability Management under the Heston Model. Adv. Differ. Equations 2018, 2018, 258. [Google Scholar] [CrossRef]
  35. Trabelsi, N.; Tiwari, A.K. Market-Risk Optimization among the Developed and Emerging Markets with CVaR Measure and Copula Simulation. Risks 2019, 7, 78. [Google Scholar] [CrossRef]
  36. Han, C.-H.; Wang, K. Stressed Portfolio Optimization with Semiparametric Method. Financ. Innov. 2022, 8, 27. [Google Scholar] [CrossRef] [PubMed]
  37. Consiglio, A.; Cocco, F.; Zenios, S.A. Scenario Optimization Asset and Liability Modelling for Individual Investors. Ann. Oper. Res. 2007, 152, 167–191. [Google Scholar] [CrossRef]
  38. Daulay, S.N.R.; Halim, N.A.; Hidayana, R.A. Investment Portfolio Optimization with a Mean-Variance Model Without Risk-Free Assets. Int. J. Quant. Res. Model. 2022, 3, 113–117. [Google Scholar] [CrossRef]
  39. Gusliana, S.A.; Salih, Y. Mean-Variance Investment Portfolio Optimization Model Without Risk-Free Assets in Jii70 Share. Oper. Res. Int. Conf. Ser. 2022, 3, 101–106. [Google Scholar] [CrossRef]
  40. Sirait, E.P.; Salih, Y.; Hidayana, R.A. Investment Portfolio Optimization Model Using The Markowitz Model. Int. J. Quant. Res. Model. 2022, 3, 124–132. [Google Scholar] [CrossRef]
Figure 1. Efficient Portfolio Surface Graph.
Figure 1. Efficient Portfolio Surface Graph.
Computation 12 00120 g001
Figure 2. The Relationship between Risk Aversion and Mean of the Portfolio Return.
Figure 2. The Relationship between Risk Aversion and Mean of the Portfolio Return.
Computation 12 00120 g002
Figure 3. The Relationship between Risk Aversion and the VaR of the Portfolio Return.
Figure 3. The Relationship between Risk Aversion and the VaR of the Portfolio Return.
Computation 12 00120 g003
Figure 4. The Relationship between the VaR of the Portfolio Return and the Ratio.
Figure 4. The Relationship between the VaR of the Portfolio Return and the Ratio.
Computation 12 00120 g004
Table 1. Descriptive Statistics of Stock Data.
Table 1. Descriptive Statistics of Stock Data.
Stock CodeMeanVarianceStandard Deviation
BSSR0.01990.01850.1362
BYAN0.03140.02420.1556
CITA0.03740.02880.1698
HRUM0.02140.02830.1683
MBAP0.02720.01990.1412
MDKA0.04690.01410.1187
MEDC0.03540.04870.2207
PSAB0.01350.02700.1642
PTBA0.08550.35040.5919
PTRO0.03400.02390.1546
RUIS0.01350.01840.1357
Table 2. Summary of Stock Weight Optimization Results on Investment Portfolios.
Table 2. Summary of Stock Weight Optimization Results on Investment Portfolios.
VariableValue(s)
Risk Aversion 5.1 < c 8.2 ;   c = 0.05
Minimum Portfolio E R = 0.0102964 ; V a R 0.05 R = 0.0137975
Maximum Portfolio E R = 0.0103316 ; V a R 0.05 R = 0.0138270
Optimum Portfolio E R = 0.0103316 ; V a R 0.05 R = 0.0138270
Optimum Capital Weight Allocation of Each Stocks w 1 = 0.0497
w 2 = 0.0545
w 3 = 0.0510
w 4 = 0.0333
w 5 = 0.0239
w 6 = 0.1275
w 7 = 0.0027
w 8 = 0.0248
w 9 = 0.0085
w 10 = 0.0294
w 11 = 0.0947
Table 3. Accuracy comparison between the mean of the optimum portfolio returns from models I, II, and III.
Table 3. Accuracy comparison between the mean of the optimum portfolio returns from models I, II, and III.
Observation PeriodForecasting DateAbsolute Error of Portfolio Return Mean
Model IModel IIModel III
1 January 2017 to 1 December 20201 January 20210.07560.01340.0535
1 January 2017 to 1 January 20211 February 20210.03010.02810.0277
1 January 2017 to 1 February 20211 March 20210.03820.11500.0653
1 January 2017 to 1 March 20211 April 20210.00510.07990.0161
1 January 2017 to 1 April 20211 May 20210.00300.03900.0141
1 January 2017 to 1 May 20211 June 20210.01900.01490.0105
1 January 2017 to 1 June 20211 July 20210.00380.01280.0114
1 January 2017 to 1 July 20211 August 20210.02040.02070.0218
1 January 2017 to 1 August 20211 September 20210.05600.11020.0750
1 January 2017 to 1 September 20211 October 20210.00630.03940.0275
1 January 2017 to 1 October 20211 November 20210.01400.07150.0389
1 January 2017 to 1 November 20211 December 20210.01770.00890.0147
1 January 2017 to 1 December 20211 January 20220.05070.02700.0407
Mean Absolute Error (MAE)0.02610.04470.0321
The smallest absolute error for each forecasting date and the smallest MAE are boldface.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sukono; Ghazali, P.L.B.; Johansyah, M.D.; Riaman; Ibrahim, R.A.; Mamat, M.; Sambas, A. Modeling of Mean-Value-at-Risk Investment Portfolio Optimization Considering Liabilities and Risk-Free Assets. Computation 2024, 12, 120. https://doi.org/10.3390/computation12060120

AMA Style

Sukono, Ghazali PLB, Johansyah MD, Riaman, Ibrahim RA, Mamat M, Sambas A. Modeling of Mean-Value-at-Risk Investment Portfolio Optimization Considering Liabilities and Risk-Free Assets. Computation. 2024; 12(6):120. https://doi.org/10.3390/computation12060120

Chicago/Turabian Style

Sukono, Puspa Liza Binti Ghazali, Muhamad Deni Johansyah, Riaman, Riza Andrian Ibrahim, Mustafa Mamat, and Aceng Sambas. 2024. "Modeling of Mean-Value-at-Risk Investment Portfolio Optimization Considering Liabilities and Risk-Free Assets" Computation 12, no. 6: 120. https://doi.org/10.3390/computation12060120

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop