Next Article in Journal
Fractal Generalization of the Scher–Montroll Model for Anomalous Transit-Time Dispersion in Disordered Solids
Next Article in Special Issue
Applying Heath-Jarrow-Morton Model to Forecasting the US Treasury Daily Yield Curve Rates
Previous Article in Journal
A New Kumaraswamy Generalized Family of Distributions with Properties, Applications, and Bivariate Extension
Previous Article in Special Issue
Self-Similar Models: Relationship between the Diffusion Entropy Analysis, Detrended Fluctuation Analysis and Lévy Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

GO-GJRSK Model with Application to Higher Order Risk-Based Portfolio

1
NOMURA Asset Management Co. Ltd., 2-2-1, Toyosu, Koto-ku, Tokyo 135-0061, Japan
2
MAZIN Inc., 3-29-14 Nishi-Asakusa, Tito city, Tokyo 111-0035, Japan
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(11), 1990; https://doi.org/10.3390/math8111990
Submission received: 8 June 2020 / Revised: 30 October 2020 / Accepted: 3 November 2020 / Published: 7 November 2020

Abstract

:
There are three distinguishing features in the financial time series, such as stock prices, are as follows: (1) Non-normality, (2) serial correlation, and (3) leverage effect. All three points need to be taken into account to model the financial time series. However, multivariate financial time series modeling involves a large number of stocks, with many parameters to be estimated. Therefore, there are few examples of multivariate financial time series modeling that explicitly deal with higher-order moments. Furthermore, there is no multivariate financial time series model that takes all three characteristics above into account. In this study, we propose the generalized orthogonal (GO)-Glosten, Jagannathan, and Runkle GARCH (GJR) model which extends the GO-generalized autoregressive conditional heteroscedasticity (GARCH) model and incorporates the three features of the financial time series. We confirm the effectiveness of the proposed model by comparing the performance of risk-based portfolios with higher-order moments. The results show that the performance with our proposed model is superior to that with baseline methods, and indicate that estimation methods are important in risk-based portfolios with higher moments.

1. Introduction

The financial time series is continually brought to our attention. Daily news reports inform us for instance of the latest stock market index values, currency exchange rates, oil prices, and interest rates. It is often desirable to monitor their price behaviors frequently and to try to understand the probable development of the prices in the future. Institutional and individual investors involved in international trades can all benefit from a deeper understanding of price behaviors. Therefore, modeling and understanding the financial time series is a very important task.
There are three distinguishing features in the financial time series, such as stock prices, as follows: (1) Non-normality, (2) serial correlation, and (3) leverage effect. (1) Non-normality means that the returns of financial time series does not follow the normal distribution because of negative skewness and large kurtosis [1,2]. Therefore, we need to model not only mean and variance but also higher-order moments such as skewness and kurtosis. (2) The serial correlation of return is known as momentum or reversal effect. The momentum [3] (reversal [4]) effect means that stock returns with a large (small) return in the past will continue to increase (decrease) during a certain period thereafter. In the case of volatility, the serial correlation is widely known as the volatility clustering [5]. It is a phenomenon that large changes in returns tend to cluster together, resulting in the persistence of the amplitudes of returns. The volatility clustering has intrigued many researchers and oriented in a major way the development of financial time series models. For example, the autoregressive conditional heteroscedasticity (ARCH) [5] and the generalized ARCH (GARCH) model [6] are intended primarily to model this phenomenon. The authors of [7] demonstrated that skewness and kurtosis also have serial correlations, and they developed the GARCH with skewness and kurtosis (GARCHSK) model to capture these correlations of skewness and kurtosis. (3) The leverage effect corresponds to a negative correlation between returns and volatility. The Glosten, Jagannathan, and Runkle GARCH (GJR) model [8] is a representative model that captures the leverage effect. Recently, Many studies [9,10] demonstrated a negative correlation between returns, and skewness and kurtosis.
We need to take all three points into account to model the financial time series. However, when considering multivariate financial time series modeling involving a large number of stocks, there are many parameters to be estimated. Even when we only model the time-varying covariance matrix, it is very difficult to estimate it completely [11]. The generalized orthogonal GARCH (GO-GARCH) model [12] is proposed to estimate the time-varying covariance matrix using independent component analysis (ICA) to reduce the parameters to be estimated.
Therefore, there are few examples of multivariate financial time series modeling that explicitly deal with higher-order moments [13]. Furthermore, there is no multivariate financial time series that takes all three characteristics above into account.
In this study, we propose the GO-GJRSK model which extends the GO-GARCH and incorporates the three features of the financial time series. To define GO-GJRSK model, we first introduce a univariate GJR with skewness and kurtosis (GJRSK) model which is the GARCHSK model with a leverage effect. Furthermore, we extend the GJRSK model to a multivariate model with ICA.
Our proposed method describes the conditional volatility (skewness and kurtosis) in a straightforward form: The linear sum of the previous time volatility (skewness and kurtosis) and the shock with leverage effect. As with the GO-GARCH model, the estimation of the GO-GJRSK model is also easy when the number of stocks is large enough.
We confirm the effectiveness of the proposed methodology by comparing the performance of risk-based portfolios with higher-order moments. We use a representative benchmark dataset to test whether differences in estimation methods improve the performance of risk-based portfolios.
The remainder of the paper is organized as follows. Section 2 gives a brief explanation of risk-based portfolios and the volatility modeling methods. Section 3 describes our proposed GO-GJRSK model, which consists of an ICA and GJRSK model. Section 4 experiments with various risk-based portfolios with different estimation methods. Lastly, Section 5 concludes the paper.

2. Related Work

2.1. Risk Based Portfolio

Portfolio construction is one of the most important tasks in finance. Mean Variance Optimization (MVO) proposed by [14] has been the most used in practice for portfolio construction and asset allocation for a long time. The MVO is still a very important framework for investment decision making today. The MVO provides economical justification because it is not only simple and practically tractable but also because it is consistent with the expected utility maximization principle. The goal of MVO is to find portfolios that optimally diversify risk without reducing expected return and to facilitate portfolio construction.
However, there are some practical difficulties with the MVO. First, it is very difficult to estimate the expected return as inputs of the MVO. In addition, the MVO is referred to as “Error Maximization” [15,16]. The reason is that even small changes in the estimated expected returns can result in huge changes in the whole portfolio structure.
Therefore, risk-based portfolio construction methods, which do not require the expected return, which is difficult to estimate, but are based only on risk, are attracting attention as an alternative to the MVO. Minimum variance portfolio [17], risk parity portfolio [18], and maximum diversification portfolio [19] are known as the typical risk-based portfolio construction methods.
In fact, many empirical studies of equity portfolios and asset allocations [17,18,19,20] showed that the risk-based portfolio construction methods perform better than the MVO and market capitalization-weighted portfolios. In [21], it was demonstrated that risk-based portfolios are not very sensitive to the accuracy of estimating covariance matrices, unlike MVO.
More recently, the covariance matrix has been estimated using Gaussian process latent variable models [22], and Student’s t-process latent variable models [23] for the minimum variance portfolio. And new risk parity portfolios are proposed by using the Hilbert transform to make the covariance matrix complex [24] or extended to quaternions [25].
However, these risk-based portfolio construction methods capture risk only in terms of covariance matrix (second moments) and do not capture the non-normality, i.e., skewness and kurtosis. Therefore, risk-based portfolios with higher-order moments have been proposed and their effectiveness has been confirmed [13,26,27]. We examine whether the estimation method improves the performance of risk-based portfolios with higher moments.

2.2. Volatility Modeling

Volatility modeling is important for risk management and the portfolio selection. In [5], it was proposed that the ARCH model could model the serial correlation of volatility. In [6], the GARCH model was proposed, which is a generalization of the ARCH model. The GARCH model assumes that the volatility at time t is formulated with only variables whose values are already known in time t 1 . The GARCH model has been applied to the empirical analysis of many asset prices because the parameters can easily be estimated by the maximum likelihood method. The GARCH model has been modified and extended in a number of ways. An important aspect of improving the GARCH model is the introduction of a leverage effect.
The GJR model [8] is a representative model that incorporates a leverage effect. In another direction of extension of GARCH model, León et al. [7] proposed the GARCHSK model, which also incorporates conditional skewness and conditional kurtosis variations. The feature of the GARCHSK model is that the conditional skewness and conditional kurtosis variations can be explicitly captured with an easy-to-understand structure equivalent to the GARCH model, and the easy to estimate.
Furthermore, extensions to multivariate models are also proposed. In order to efficiently parameterize multivariate GARCH, some methods have been developed based on univariate GARCH models in the following two directions. One is to construct the conditional covariance matrix by explicitly modeling both the volatility and conditional correlation matrix. For example, Engle [28] proposed the dynamic conditional correlation (DCC) model. The other direction is to exploit the idea of factor models to construct multivariate GARCH models. It is widely believed in finance that there exist some common factors, which may be determined empirically or constructed by some statistical methods, driving the evolution of the return series [22,23,29,30]. This approach has commonly been used in much research to model the conditional covariance of financial time series due to its feasibility in estimating large covariance matrices.
Orthogonal GARCH (OGARCH) was first proposed by [31]. The observed time series can be linearly transformed into a set of uncorrelated time series using a principal component analysis. In [12], the concept of ICA is applied to propose the GO-GARCH model for volatility modeling. It consists of a set of conditionally uncorrelated univariate GARCH and a linear map that allows the linkage between these components and observed data. Furthermore, a modified GO-GARCHSK model [32] that incorporates the GO-GARCH model with conditional skewness and conditional kurtosis variation has been proposed. However, there are no multivariate models that have a conditional skewness and kurtosis variation and leverage effect.
Hence, as a multivariate extension of the GJRSK, we propose GO-GJRSK model that has higher order moment variations and a leverage effect (Table 1).

3. GO-GJRSK Model

3.1. Univariate GJRSK Model

We first define a univariate GJRSK model. This model introduces leverage effects in the GARCHSK model. Here, r t is the asset return, h t is conditional volatility, s t is conditional skewness, and k t is conditional kurtosis at time t. The formulation of the GJRSK model is given as follows.
r t = α 0 + α 1 r t 1 + ε t
h t = β 0 + β 1 ε t 1 2 + β 2 h t 1 + β 3 ε t 1 2 I { η t 1 < 0 }
s t = γ 0 + γ 1 η t 1 3 + γ 2 s t 1 + γ 3 η t 1 3 I { η t 1 < 0 }
k t = δ 0 + δ 1 η t 1 4 + δ 2 k t 1 + δ 3 η t 1 4 I { η t 1 < 0 }
η t = h t 1 2 ε t
η t | I t 1 g ( 0 , 1 , s t , k t )
Here α 0 , α 1 are parameters of the AR model and β i , γ i , δ i , i = 0 , 1 , 2 , 3 are parameters of the GJRSK model. I A is an indicator function that returns 1 if A is true and 0 otherwise. g is a probability density function with mean 0, variance 1, skewness s t , and kurtosis k t .
As with GARCHSK model, the probability density function g ( 0 , 1 , s t , k t ) of the GJRSK model can be given by Gram-Charlier expansion using Chebyshev-Hermite polynomials as follows:
g ( η t | I t 1 ) = ϕ ( η t ) ψ 2 ( η t ) Γ t ,
ϕ ( η t ) = 1 2 π h t exp ( η t 2 h t ) ,
ψ ( η t ) = 1 + s t 3 ! ( η t 3 3 η t ) + k t 3 4 ! ( η t 4 6 η t 2 + 3 ) ,
Γ t = 1 + s t 2 3 ! + ( k t 3 ) 2 4 ! .
The probability density function of ε t is f ( η t | I t 1 ) = h t 1 2 g ( η t | I t 1 ) because η t = h t 1 2 ε t . Therefore, the log-likelihood function l t without the constant term is obtained as:
l t = 1 2 ln h t 1 2 η t 2 + ln ( ψ 2 ( η t ) ) Γ t .
We can estimate each parameter of the GJRSK model by maximizing log-likelihood function l t .

3.2. GO-GJRSK Model

In this section, we introduce the GO-GJRSK model. We define the co-skewness matrix M 3 and co-kurtosis matrix M 4 . Let co-skewness of assets i , j , k be s i j k and S i j k be the matrix that can be made by fixing i of s i j k . The co-skewness matrix M 3 is defined as follows.
S i j k = s i 11 s i 1 k s i j 1 s i j k M 3 = [ S 1 j k , S 2 j k , , S N j k ]
Let co-kurtosis of assets i , j , k , l be k i j k l and K i j k l be the matrix that can be made by fixing i , j of k i j k . The co-kurtosis matrix M 4 is defined as follows:
K i j k l = k i j 11 k i j 1 k k i j k 1 k i j k l M 4 = [ K 11 k l , , K 1 N k l | K 21 k l , , K 2 N k l | , , | K N 1 k l , , K N N k l ] .
The GO-GJRSK model assumes that ε t R n is defined by mixing matrix Z R n × n and n independent factors y t = ( y 1 , t , , y n , t ) T R n . When the n assets returns are r t R n , the GO-GJRSK model can be written as follows:
r t = A 0 + A 1 r t 1 + ε t ,
ε t = Z y t ,
η i , t = m i , t H 1 2 y i , t ,
η i , t | I t 1 g ( 0 , 1 , m i , t S , m i , t K ) ,
m i , t H = β i , 0 + β i , 1 y i , t 1 2 + β i , 2 m i , t 1 H + β i , 3 y i , t 1 2 I { η i , t 1 < 0 } ,
m i , t S = γ i , 0 + γ i , 1 η i , t 1 3 + γ i , 2 m i , t 1 S + γ i , 3 η i , t 1 3 I { η i , t 1 < 0 } ,
m i , t K = δ i , 0 + δ i , 1 η i , t 1 4 + δ i , 2 m i , t 1 K + δ i , 3 η i , t 1 4 I { η i , t 1 < 0 } ,
where A 0 R n and A 1 R n × n are parameters of the VAR model and β i , j , γ i , j , δ i , j , i = 1 , , n , j = 0 , , 3 are parameters of the GJRSK model. The GO-GJRSK model follows a probability density function p ( 0 , H t , S t , K t ) with a conditional covariance matrix H t , a co-skewness matrix S t , and a co-kurtosis matrix K t . Here, the following relationship holds:
H t = Z M t H Z T
S t = Z M t S ( Z T Z T )
K t = Z M t K ( Z T Z T Z T )
M t H = diag H { m 1 , t H , , m n , t H }
M t S = diag S { m 1 , t S , , m n , t S }
M t K = diag K { m 1 , t K , , m n , t K } .
The operator ⊗ represents the Kronecker product. diag H ( e 1 , , e n ) makes a ( n × n ) sparse covariance matrix with elements e i when i = j and zero otherwise. diag S ( e 1 , , e n ) makes a ( n × n 2 ) sparse co-skewness matrix with elements e i when i = j = k and zero otherwise. diag K ( e 1 , , e n ) makes a ( n × n 3 ) sparse co-kurtosis matrix with elements e i when i = j = k = l and zero otherwise. Using this relationship, we can obtain the conditional covariance matrix H t , the co-skewness matrix S t , and the co-kurtosis matrix K t .

3.3. Estimation Procedure

The parameters of GO-GJRSK model are estimated by the following procedure.
  • Step1 VAR Model:
    First, we estimate the vector autoregression model (VAR) model in Equation (15). After estimating the parameter of VAR model, the residuals ε t have been collected for further GO-GJRSK modeling.
  • Step2 Independent Component Analysis:
    ICA is a signal processing method to express a set of random variables as linear combinations of statistically independent component variables [33]. If y i , t and y j , t , ( i j ) are mutually statistically independent, and mixing matrix Z exists, such that:
    ε t = Z y t
    then y t are n independent component of the residuals ε t . The mixing matrix Z can be solved by fast fixed-point algorithm given by [34].
  • Step3 Estimate Univariate GJRSK Model:
    Once the independent components y t are obtained, the conditional variance, conditional skewness, and conditional kurtosis of independent components can be estimated by the univariate GJRSK model. The parameters of the univariate GARCHSK model can be estimated by maximizing the log-likelihood function in Equation (11).

4. Experiment

In this section, we report the results of our empirical studies with well-known benchmarks. We compared the out-of-sample performance among several portfolio strategies to check the effectiveness of our GO-GJRSK model.

4.1. Dataset

In the experiments, we used well-known academic benchmarks called Fama and French (FF) datasets (https://mba.tuck.dartmouth.edu/pages/faculty/ken.french/data_library.html) to ensure the reproducibility of the experiment. This FF dataset is public and is readily available to anyone. The FF datasets have been recognized as standard datasets and heavily adopted in finance research because of their extensive coverage of asset classes and very long historical data series. We used the FF17 dataset and FF25 dataset. For example, the FF17 dataset contains monthly returns of 17 portfolios representing different industrial sectors while the FF25 dataset includes 25 portfolios formed on the basis of size and book-to-market ratio. We used both datasets as monthly data from January 1989 to December 2018.

4.2. Higher Order Risk-Based Portfolio

Here, we formulate the risk-based portfolio construction methodology in the experiment. In this study, two risk-based portfolios with higher-order moments are used to confirm the effectiveness of our proposed model. One is a minimum variance portfolio and the other is a risk parity portfolio.
Let w = ( w 1 , , w n ) T be the weights of each of the n assets, M 2 is the covariance matrix, M 3 is the co-skewness matrix, and M 4 is the co-kurtosis matrix. Then, the portfolio variance σ P 2 , skewness s P 3 , and kurtosis k P 4 can be expressed as follows [35]:
σ P 2 = w T M 2 w
s P 3 = w T M 3 ( w w )
k P 4 = w T M 4 ( w w w )
We then use these portfolio characteristics to extend the minimum variance portfolio and the risk-based portfolio to include higher-order moments.

4.3. Minimum Variance Portfolio with Higher Moments

The minimum variance portfolio is extended with the skewness in Equation (29) and kurtosis in Equation (30) as follows [26]:
min w λ 1 σ P 2 λ 2 s P 3 + λ 3 k P 4 , s . t . i = 1 N w i = 1 , w i > 0 ,
where λ i > 0 , i = 1 , 2 , 3 are weights to account for the trade-offs between the moments of each order. Here, if λ 1 = 1 , λ 2 = 0 , λ 3 = 0 , Equation (31) is the minimum variance portfolio.

4.4. Risk Parity Portfolio with Higher Moments

The extension of risk parity portfolios to higher-order moments is discussed in [27]. First, as in the risk parity portfolio, we define from the second to fourth order marginal risk contribution ( M R C ). M R C i , i = 2 , 3 , 4 are derived by differentiating the portfolio variance of Equation (28), skewness of Equation (29), and the kurtosis of Equation (30) by the weight w, respectively.
M R C 2 = 1 2 σ P 2 w = M 2 w
M R C 3 = 1 3 s P 3 w = M 3 ( w w )
M R C 4 = 1 4 k P 3 w = M 4 ( w w w ) .
Then, we can decompose the portfolio variance, skewness and kurtosis using M M C i as follows:
σ P 2 = i = 1 N w i × M R C 2 , i = w T M R C 2
s P 3 = i = 1 N w i × M R C 3 , i = w T M R C 3
k P 4 = i = 1 N w i × M R C 4 , i = w T M R C 4 .
We define the contribution of the moments ( R M ) by dividing Equations (35)–(37) by each portfolio moment, respectively.
R M 2 , i = w i × M R C 2 , i σ P 2
R M 3 , i = w i × M R C 3 , i s P 3
R M 4 , i = w i × M R C 4 , i k P 4
Based on the above, the risk parity portfolio is extended as follows [27]:
min w λ 1 f 1 ( w ) + λ 2 f 2 ( w ) + λ 3 f 3 ( w ) s . t . i = 1 N w i = 1 , w i > 0 , f 1 ( w ) = i = 1 N j = 1 N M C 2 , i M C 2 , j 2 , f 2 ( w ) = i = 1 N j = 1 N M C 3 , i M C 3 , j 2 , f 3 ( w ) = i = 1 N j = 1 N M C 4 , i M C 4 , j 2 ,
where λ i > 0 , i = 1 , 2 , 3 are weights to account for the trade-offs among the M C i . Here, if λ 1 = 1 , λ 2 = 0 , λ 3 = 0 , Equation (41) is the risk parity portfolio.

4.5. Experimental Settings

In our empirical studies, the tested portfolios have the following meanings:
  • “EW” stands for equally-weighted portfolio [36];
  • “MV” stands for minimum-variance portfolio [14];
  • “RP” stands for risk parity portfolio [18];
  • “HMV” stands for higher order minimum-variance portfolio [26]. Here, λ 1 = 1 , λ 2 = 1 , λ 3 = 1 in Equation (31);
  • “HRP” stands for higher order risk parity portfolio [27]. Here, λ 1 = 1 , λ 2 = 1 , λ 3 = 1 in Equation (41).
The following method is used to estimate the covariance, co-skewness, and co-kurtosis matrix:
  • “MLE” stands for maximum likelihood estimation of covariance, co-skewness and co-kurtosis matrix;
  • “GO” stands for covariance matrix in GO-GARCH model [12];
  • “GO-SK” stands for covariance, co-skewness and co-kurtosis matrix in GO-GARCHSK model [32];
  • “GO-GJRSK” stands for covariance, co-skewness, and co-kurtosis matrix in our proposed GO-GJRSK model.
We used the latest 10 years (120 months) to compute for the parameter estimation, and used each dataset from July 1936 to March 2020, as the out-of-sample period. Each portfolio was updated by sliding one-month-ahead.
We performed our experiment with R language.

4.6. Performance Measures

We compared the out-of-sample performance of the portfolios. In evaluating the portfolio strategy, we used the following measures.
Here, the portfolio return R t at time t is defined as:
R t = i = 1 n r i , t w i , t 1 ,
where r i , t is the return of i asset at time t, w i , t 1 is the weight of i asset in the portfolio at time t 1 , and n is the number of assets.
We evaluated the portfolio strategy using its annualized return (AR), the downside risk as the standard deviation of negative return (DR), and the return/downside risk (R/R) (return divided by downside risk) as the portfolio strategy. As a measure of risk, standard deviation is usually used [14]. However, the standard deviation as a risk measure has several drawbacks. It punishes both upside and downside deviations of return. Moreover, the variance is justified only when the joint distribution of asset returns is symmetric. In this study, we assume that the asset returns are asymmetric and non-normal as several empirical study shows. On the other hand, DR only penalizes the downside deviations of returns. Hence, it is intuitively much more appealing.
AR = 12 T t = 1 T R t
DR = 12 T 1 t = 1 T ( R t × 1 { R t 0 } ) 2
R / R = AR / DR
Here, μ = ( 1 / T ) t = 1 T R t is the average return of the portfolio.

4.7. Results

The results are shown in Table 2. For both FF17 and FF25, the AR and R/R of HMV portfolio using the GO-GJRSK model are the highest in MV and HMV portfolios. The DR of MV portfolio using the GO-GARCH model has the lowest in MV and HMV portfolios.
Furthermore, interestingly, extending MV to HMV using a simple MLE does not improve R/R. Furthermore, when the GO-GARCHSK model is used in HMV, the R/R is subordinated to the HMV using MLE.
On the other hand, for both FF17 and FF25, HRP using the GO-GJRSK model gives the best results in terms of AR, DR, and R/R in RP and HRP portfolios. As in the result of MV portfolio, extending RP to HRP with a simple MLE does not improve R/R. This is a different result from the previous study [27]. Furthermore, when the GO-GARCHSK model is used in HRP, the R/R is inferior to the HMV using MLE. As in [21], risk-based portfolios without higher moments have little impact on performance due to differences in estimation methods.
We can confirm the effectiveness of our proposed GO-GJRSK model for simple MLE and GO-GARCHSK.

5. Conclusions

In this study, we proposed the GO-GJRSK model which extends the GO-GARCH and incorporates the three features of a financial time series: (1) non-normality, (2) serial correlation, and (3) leverage effect.
We confirmed the effectiveness of the proposed method by comparing the performance of risk-based portfolios with higher-order moments. The results showed that the performance of the proposed method was superior to that of the MLE and GO-GARCHSK. We showed that extending MV to HMV and RP to HRP with a simple MLE did not improve R/R. This result indicated that estimation methods are important in risk-based portfolios with higher moments as with MVO.
A future challenge is to extend our model to stochastic volatility, skewness, and kurtosis models because stochastic volatility model might provide a better alternative to the GARCH models [37].

Author Contributions

Conceptualization, K.N.; Data curation, K.N.; Formal analysis, K.N. and Y.U.; Methodology, K.N. and Y.U.; Project administration, K.N. and Y.U.; Software, K.N.; Validation, K.N.; Writing–original draft, K.N.; Writing–review & editing, K.N. and Y.U. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ARCHautoregressive conditional heteroscedasticity
GARCHgeneralized ARCH
GARCHSKGARCH with skewness and kurtosis
GJRGlosten, Jagannathan and Runkle GARCH
GJRSKGJR with skewness and kurtosis
DCCdynamic conditional correlation
O-GARCHorthogonal GARCH
GO-GARCHgeneralized O-GARCH
GO-GARCHSKGO-GARCH with skewness and kurtosis
GO-GJRSKGO-GJR with skewness and kurtosis
ICAindependent component analysis
VARvector autoregression
FF17Fama French dataset for 17 sectors
FF25Fama French dataset for 25 factors
EWequally weight portfolio
MVminimum variance portfolio
HMVminimum variance portfolio with higher moments
RPrisk parity portfolio
HRPrisk parity portfolio with higher moments
MLEmaximum likelihood estimator
ARannualized return
DRannualized downside risk
R/RAR/DR

References

  1. Mandelbrot, B.B. The variation of certain speculative prices. In Fractals and Scaling in Finance; Springer: Berlin/Heidelberg, Germany, 1997; pp. 371–418. [Google Scholar]
  2. Cont, R. Empirical properties of asset returns: Stylized facts and statistical issues. Quant. Financ. 2001, 1, 223–236. [Google Scholar] [CrossRef]
  3. Jegadeesh, N.; Titman, S. Returns to buying winners and selling losers: Implications for stock market efficiency. J. Financ. 1993, 48, 65–91. [Google Scholar] [CrossRef]
  4. Jegadeesh, N. Evidence of predictable behavior of security returns. J. Financ. 1990, 45, 881–898. [Google Scholar] [CrossRef]
  5. Engle, R.F. Autoregressive conditional heteroscedasticity with estimates of the variance of United Kingdom inflation. Econom. J. Econom. Soc. 1982, 50, 987–1007. [Google Scholar] [CrossRef]
  6. Bollerslev, T. Generalized autoregressive conditional heteroskedasticity. J. Econom. 1986, 31, 307–327. [Google Scholar] [CrossRef] [Green Version]
  7. León, Á.; Rubio, G.; Serna, G. Autoregresive conditional volatility, skewness and kurtosis. Q. Rev. Econ. Financ. 2005, 45, 599–618. [Google Scholar] [CrossRef]
  8. Glosten, L.R.; Jagannathan, R.; Runkle, D.E. On the relation between the expected value and the volatility of the nominal excess return on stocks. J. Financ. 1993, 48, 1779–1801. [Google Scholar] [CrossRef]
  9. Boyer, B.; Mitton, T.; Vorkink, K. Expected idiosyncratic skewness. Rev. Financ. Stud. 2010, 23, 169–202. [Google Scholar] [CrossRef]
  10. Blau, B.M.; Whitby, R.J. Idiosyncratic kurtosis and expected returns. J. Investig. 2017, 26, 81–88. [Google Scholar] [CrossRef]
  11. Engle, R.F.; Ledoit, O.; Wolf, M. Large dynamic covariance matrices. J. Bus. Econ. Stat. 2019, 37, 363–375. [Google Scholar] [CrossRef] [Green Version]
  12. Van der Weide, R. GO-GARCH: A multivariate generalized orthogonal GARCH model. J. Appl. Econom. 2002, 17, 549–564. [Google Scholar] [CrossRef]
  13. Boudt, K.; Carl, P.; Peterson, B.G. Asset allocation with conditional value-at-risk budgets. J. Risk 2012, 15, 39–68. [Google Scholar] [CrossRef]
  14. Markowitz, H. Portfolio selection. J. Financ. 1952, 7, 77–91. [Google Scholar]
  15. Michaud, R.O. The Markowitz optimization enigma: Is ‘optimized’optimal? Financ. Anal. J. 1989, 45, 31–42. [Google Scholar] [CrossRef]
  16. Chopra, V.K.; Ziemba, W.T. The effect of errors in means, variances, and covariances on optimal portfolio choice. In Handbook of the Fundamentals of Financial Decision Making: Part I; World Scientific: Singapore, 2013; pp. 365–373. [Google Scholar]
  17. Clarke, R.G.; De Silva, H.; Thorley, S. Minimum-variance portfolios in the US equity market. J. Portf. Manag. 2006, 33, 10–24. [Google Scholar] [CrossRef]
  18. Maillard, S.; Roncalli, T.; Teïletche, J. The properties of equally weighted risk contribution portfolios. J. Portf. Manag. 2010, 36, 60–70. [Google Scholar] [CrossRef]
  19. Choueifaty, Y.; Coignard, Y. Toward maximum diversification. J. Portf. Manag. 2008, 35, 40–51. [Google Scholar] [CrossRef] [Green Version]
  20. Lee, W. Risk-Based Asset Allocation: A New Answer to an Old Question? J. Portf. Manag. 2011, 37, 11. [Google Scholar] [CrossRef]
  21. Nakagawa, K.; Imamura, M.; Yoshida, K. Risk-based portfolios with large dynamic covariance matrices. Int. J. Financ. Stud. 2018, 6, 52. [Google Scholar] [CrossRef] [Green Version]
  22. Nirwan, R.S.; Bertschinger, N. Applications of gaussian process latent variable models in finance. In Proceedings of the SAI Intelligent Systems Conference, London, UK, 5–6 September 2019; Springer: Berlin/Heidelberg, Germany, 2019; pp. 1209–1221. [Google Scholar]
  23. Uchiyama, Y.; Nakagawa, K. TPLVM: Portfolio Construction by Student’s t-Process Latent Variable Model. Mathematics 2020, 8, 449. [Google Scholar] [CrossRef] [Green Version]
  24. Uchiyama, Y.; Kadoya, T.; Nakagawa, K. Complex Valued Risk Diversification. Entropy 2019, 21, 119. [Google Scholar] [CrossRef] [Green Version]
  25. Sugitomo, S.; Maeta, K. Quaternion Valued Risk Diversification. Entropy 2020, 22, 390. [Google Scholar] [CrossRef] [Green Version]
  26. Martellini, L.; Ziemann, V. Improved estimates of higher-order comoments and implications for portfolio selection. Rev. Financ. Stud. 2010, 23, 1467–1502. [Google Scholar] [CrossRef]
  27. Baitinger, E.; Dragosch, A.; Topalova, A. Extending the risk parity approach to higher moments: Is there any value added? J. Portf. Manag. 2017, 43, 24–36. [Google Scholar] [CrossRef]
  28. Engle, R. Dynamic conditional correlation: A simple class of multivariate generalized autoregressive conditional heteroskedasticity models. J. Bus. Econ. Stat. 2002, 20, 339–350. [Google Scholar] [CrossRef]
  29. Fama, E.F.; French, K.R. Common risk factors in the returns on stocks and bonds. J. Financ. Econ. 1993, 33, 3–56. [Google Scholar] [CrossRef]
  30. Fama, E.F.; French, K.R. A five-factor asset pricing model. J. Financ. Econ. 2015, 116, 1–22. [Google Scholar] [CrossRef] [Green Version]
  31. Ding, Z. Time Series Analysis of Speculative Returns; Department of Economics, University of California: San Diego, CA, USA, 1994. [Google Scholar]
  32. Xu, Q.F.; Jiang, C.X.; Kang, P. Dynamic portfolio selection under higher moments. In Proceedings of the 2007 International Conference on Machine Learning and Cybernetics, Hong Kong, China, 19–22 August 2007; IEEE: Piscataway, NJ, USA, 2007; Volume 5, pp. 2488–2493. [Google Scholar]
  33. Comon, P. Independent component analysis, a new concept? Signal Process. 1994, 36, 287–314. [Google Scholar] [CrossRef]
  34. Hyvärinen, A.; Oja, E. A fast fixed-point algorithm for independent component analysis. Neural Comput. 1997, 9, 1483–1492. [Google Scholar] [CrossRef]
  35. Jondeau, E.; Rockinger, M. Optimal portfolio allocation under higher moments. Eur. Financ. Manag. 2006, 12, 29–55. [Google Scholar] [CrossRef] [Green Version]
  36. DeMiguel, V.; Garlappi, L.; Uppal, R. Optimal versus naive diversification: How inefficient is the 1/N portfolio strategy? Rev. Financ. Stud. 2009, 22, 1915–1953. [Google Scholar] [CrossRef] [Green Version]
  37. Chan, J.C.; Grant, A.L. Modeling energy price dynamics: GARCH versus stochastic volatility. Energy Econ. 2016, 54, 182–189. [Google Scholar] [CrossRef] [Green Version]
Table 1. Comparison between our proposed and the models of related work.
Table 1. Comparison between our proposed and the models of related work.
Non-NormalityHigher MomentLeverage Effect
GO-GARCH [12]
GO-GARCHSK [32]
GO-GJRSK
(Our proposed)
Table 2. AR, DR, and R/R for an EW, MV, RP, HMV, and HRP portfolios using three methods. In the rows labeled AR, DR, and R/R, the best results appear in bold for MV, HMV portfolios, and RP, HRP portfolios.
Table 2. AR, DR, and R/R for an EW, MV, RP, HMV, and HRP portfolios using three methods. In the rows labeled AR, DR, and R/R, the best results appear in bold for MV, HMV portfolios, and RP, HRP portfolios.
FF17EWMVHMVRPHRP
MLEGOMLEGO-SKGO-GJRSKMLEGOMLEGO-SKGO-GJRSK
AR [%]11.3110.6910.2611.2411.0212.0311.3611.3611.3511.0112.37
DR [%]12.069.179.0711.5511.819.4711.5211.5211.5311.6410.61
R/R0.941.171.130.970.931.270.990.990.980.951.17
FF25EWMVHMVRPHRP
MLEGOMLEGO-SKGO-GJRSKMLEGOMLEGO-SKGO-GJRSK
AR [%]13.2011.5811.6912.8312.4513.8113.0613.0712.9712.7114.05
DR [%]14.3410.9410.8213.9114.3211.5413.8913.8913.8814.2011.38
R/R0.921.061.080.920.871.200.940.940.930.901.23
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nakagawa, K.; Uchiyama, Y. GO-GJRSK Model with Application to Higher Order Risk-Based Portfolio. Mathematics 2020, 8, 1990. https://doi.org/10.3390/math8111990

AMA Style

Nakagawa K, Uchiyama Y. GO-GJRSK Model with Application to Higher Order Risk-Based Portfolio. Mathematics. 2020; 8(11):1990. https://doi.org/10.3390/math8111990

Chicago/Turabian Style

Nakagawa, Kei, and Yusuke Uchiyama. 2020. "GO-GJRSK Model with Application to Higher Order Risk-Based Portfolio" Mathematics 8, no. 11: 1990. https://doi.org/10.3390/math8111990

APA Style

Nakagawa, K., & Uchiyama, Y. (2020). GO-GJRSK Model with Application to Higher Order Risk-Based Portfolio. Mathematics, 8(11), 1990. https://doi.org/10.3390/math8111990

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