Next Article in Journal
Measures of Efficiency of Agricultural Insurance in Italy, Economic Evaluations
Next Article in Special Issue
On Identifying the Systemically Important Tunisian Banks: An Empirical Approach Based on the △CoVaR Measures
Previous Article in Journal
Alpha Beta Risk and Stock Returns—A Decomposition Analysis of Idiosyncratic Volatility with Conditional Models
Previous Article in Special Issue
A General Framework for Portfolio Theory. Part II: Drawdown Risk Measures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Target Matrix Estimators in Risk-Based Portfolios

Department of Economics and Business Studies, University of Genova, Via Vivaldi, 5, 16126 Genova, Italy
Risks 2018, 6(4), 125; https://doi.org/10.3390/risks6040125
Submission received: 14 October 2018 / Revised: 29 October 2018 / Accepted: 2 November 2018 / Published: 5 November 2018
(This article belongs to the Special Issue Computational Methods for Risk Management in Economics and Finance)

Abstract

:
Portfolio weights solely based on risk avoid estimation errors from the sample mean, but they are still affected from the misspecification in the sample covariance matrix. To solve this problem, we shrink the covariance matrix towards the Identity, the Variance Identity, the Single-index model, the Common Covariance, the Constant Correlation, and the Exponential Weighted Moving Average target matrices. Using an extensive Monte Carlo simulation, we offer a comparative study of these target estimators, testing their ability in reproducing the true portfolio weights. We control for the dataset dimensionality and the shrinkage intensity in the Minimum Variance (MV), Inverse Volatility (IV), Equal-Risk-Contribution (ERC), and Maximum Diversification (MD) portfolios. We find out that the Identity and Variance Identity have very good statistical properties, also being well conditioned in high-dimensional datasets. In addition, these two models are the best target towards which to shrink: they minimise the misspecification in risk-based portfolio weights, generating estimates very close to the population values. Overall, shrinking the sample covariance matrix helps to reduce weight misspecification, especially in the Minimum Variance and the Maximum Diversification portfolios. The Inverse Volatility and the Equal-Risk-Contribution portfolios are less sensitive to covariance misspecification and so benefit less from shrinkage.

1. Introduction

The seminal contributions of Markowitz (1952, 1956) laid the foundations for his well-known portfolio building technique. Albeit elegant in its formulation and easy to be implemented in real-world applications, the Markowitz model relies on securities returns sample mean and sample covariance as inputs to estimate the optimal allocation. However, there is a large consensus on the fact that sample estimators perpetuate large estimation errors; this directly affects portfolio weights that often exhibit extreme values, fluctuating over time with very poor performance out-of-sample (De Miguel et al. 2009).
This problem has been tackled from different perspectives: Jorion (1986) and Michaud (1989) suggested Bayesian alternatives to the sample estimators; Jagannathan and Ma (2003) added constraints to the Markowitz model limiting the estimation error; Black and Litterman (1992) derived an alternative portfolio construction technique exclusively based on the covariance matrix among asset returns, avoiding estimating the mean value for each security and converging to the Markowitz Minimum Variance portfolio with no short-sales. This latter technique is supported by results in Merton (1980) and Chopra and Ziemba (1993), who clearly demonstrated how the mean estimation process can lead to more severe distortions than those in the case of the covariance matrix.
Following this perspective, estimation error can be reduced by considering risk-based portfolios: findings suggest they have good out-of-sample performance without much turnover (De Miguel et al. 2009). There is a recent research strand focused on deriving risk-based portfolios other than the Minimum Variance one. In this context, Qian (2006) designed a way to select assets by assigning to each of them the same contribution to the overall portfolio risk; Choueifaty and Coignard (2008) proposed a portfolio where diversification is the key criterion in asset selection; Maillard et al. (2010) offered a novel portfolio construction technique where weights perpetuate an equal risk contribution while maximising diversification. These portfolios are largely popular among practitioners1: they highlight the importance of diversification, risk budgeting; moreover, they put risk management in a central role, offering a low computational burden to estimate weights. They are perceived as “robust” models since they do not require the explicit estimation of the mean. Unfortunately, limiting the estimation error in this way poses additional problems related to the ill-conditioning of the covariance matrix that occurs when the number of securities becomes sensitively greater than the number of observations. In this case, the sample eigenvalues become more dispersed than the population ones (Marčenko and Pastur 1967), and the sample covariance matrix directly affects weight estimation. This means that for a high-dimensional dataset, the sample covariance matrix is not a reliable estimator.
To reduce misspecification effects on portfolio weights, more sophisticated estimators than the sample covariance have been proposed, for example the Bayes-Stein shrinkage technique (James and Stein 1961), henceforth shrinkage stems for its practical implementation and related portfolio performance. This technique reduces the misspecification in the sample covariance matrix by shrinking it towards an alternative estimator. Here, the problem is to select a convenient target estimator as well as to find the optimal intensity at which to shrink towards the sample covariance matrix. The latter is usually derived by minimising a predefined loss function to obtain the minimum distance between the true and the shrunk covariance matrices (Ledoit and Wolf 2003). A comprehensive overview on shrinkage intensity parameters can be found in (De Miguel et al. 2013), where the authors proposed an alternative way of deriving the optimal intensity based on the smoothed bootstrap approach. On the other hand, the target matrix is often selected among the class of structured covariance estimators (Briner and Connor 2008), especially because the matrix which shrinks is the sample one. As noted in (Candelon et al. 2012), the sample covariance matrix is the Maximum Likelihood Estimator (MLE) under the Normality of asset returns, hence it lets the data speak without imposing any structure. This naturally suggests it might be pulled towards a more structured alternative. Dealing with financial data, the shrinkage literature proposes six different models for the target matrix: the Single-Index market model (Ledoit and Wolf 2003; Briner and Connor 2008; Candelon et al. 2012; Ardia et al. 2017); the Identity matrix (Ledoit and Wolf 2004a; Candelon et al. 2012); the Variance Identity matrix (Ledoit and Wolf 2004a); the Scaled Identity matrix (De Miguel et al. 2013); the Constant Correlation model (Ledoit and Wolf 2004b) and (Pantaleo et al. 2011); the Common Covariance (Pantaleo et al. 2011). All these targets belong to the class of more structured covariance estimators than the sample one, thus implying the latter is the matrix to shrink.
Despite the great improvements in portfolio weight estimation under the Markowitz portfolio building framework, the shrinkage technique has only been applied in one work involving risk-based portfolios: that of Ardia et al. (2017), who comprehensively described the impacts of variance and covariance misspecifications in risk-portfolio weights. Ardia et al. (2017) tested four alternative covariance estimators to reduce weight misspecification; among those, only one refers to shrinkage as in (Ledoit and Wolf 2003), leaving room open for further research. In our work, we contribute to the existing literature, filling this gap and offering a comprehensive overview of shrinkage in risk-based portfolios. In particular, we study the effect of six target matrix estimators on the weights of four risk-based portfolios. To achieve this goal, we provide an extensive Monte Carlo simulation aimed at (1) assessing estimators’ statistical properties and similarity with the true target matrix; (2) addressing the problem of how the selection of a specific target estimator affects the portfolio weights. We find out that the Identity and Variance Identity hold the best statistical properties, being well conditioned even in a high-dimensional dataset. These two estimators also represent the more efficient target matrices towards which to shrink the sample one. In fact, portfolio weight derived shrinking towards the Identity and Variance Identity minimise the distance from their true counterparts, especially in the case of Minimum Variance and Maximum Diversification portfolios.
The rest of the paper is organised as follows. Section 2 introduces the risk-based portfolios employed in the study. Section 3 illustrates the shrinkage estimator, the moves to the six target matrix estimators and provides useful insights into misspecification when shrinkage is applied to risk-based portfolios. In Section 4, we run an extensive Monte Carlo analysis for describing how changes in the target matrix affect risk-based portfolio weights. Section 5 concludes.

2. Risk-Based Portfolios

Risk-based portfolios are particularly appealing since they rely only on the estimation of a proper measure of risk, i.e., the covariance matrix between asset returns. Assume an investment universe made by p assets:
X =   ( x 1 , , x p )
is a n × p containing a history of n log-returns for the i-th asset, where i = 1 , , p . The covariance matrix among asset log-returns is the symmetric square matrix Σ 2 of dimension p × p , and the unknown optimal weights form the vector ω of dimension p × 1 . Our working framework assumes to consider four risk-based portfolios: the Minimum Variance (MV), the Inverse Volatility (IV), the Equal-Risk-Contribution (ERC), and the Maximum Diversification (MD) upon two constraints; no short-selling ( ω + p ) and full allocation of the available wealth ( ω 1 p = 1 , where 1 p is the vector of ones of length p ).
The Minimum Variance portfolio (Markowitz 1952) derives the optimal portfolio weights by solving this minimization problem with respect to ω :
ω M V   argmin ω { ω Σ ω   |   ω + p , ω 1 p = 1 }
where ω Σ ω is the portfolio variance.
In the Inverse Volatility portfolio, also known as the equal-risk-budget portfolio (Leote de Carvalho et al. 2012), a closed form solution is available. Each element of the vector ω is given by the inverse of the i-th asset variance (denoted by Σ i , i 1 ) divided by the inverse of the sum of all asset variances:
ω I V   ( 1 , 1 1 i = 1 p i , i 1 , , p , p 1 i = 1 p i , i 1 )
In the Equal-Risk-Contribution portfolio, as the name suggests, the optimal weights are calculated by assigning to each asset the same contribution to the whole portfolio volatility, thus originating a minimization procedure to be solved with respect to ω :
ω E R C   argmin ω { i = 1 p ( % R C i 1 p ) 2 | ω + p , ω 1 p = 1 }
Here, % R C i ω i c o v i , π ω Σ ω is the percentage risk contribution for the i-th asset, ω Σ ω is the portfolio volatility as earlier defined, and ω i c o v i , π provides a measure of the covariance of the i-th exposure to the total portfolio π , weighted by the corresponding ω i .
Turning to the Maximum Diversification, as in Choueifaty and Coignard (2008) we preliminary define D R ( ω ) as the portfolio’s diversification ratio:
D R ( ω ) ω d i a g ( Σ ) ω Σ ω
where d i a g ( Σ ) is a p × 1 vector which takes all the asset variances Σ i , i and ω d i a g ( Σ ) is the weighted average volatility. By construction, it is D R ( ω ) 1 , since the portfolio volatility is sub-additive (Ardia et al. 2017). Hence, the optimal allocation is the one with the highest DR:
ω M D   argmax ω { D R ( ω ) |   ω + p , ω 1 p = 1 } .

3. Shrinkage Estimator

The shrinkage technique relies upon three ingredients: the starting covariance matrix to shrink, the target matrix towards which the former is shrunk, and the shrinkage intensity, or roughly speaking the strength at which the starting matrix must be shrunk.
In financial applications, the starting matrix which is to shrink is always the sample covariance matrix. This is a very convenient choice that helps in the selection of a proper shrinkage target: being the sample covariance a model-free estimator that completely reflects the relationships among data3, it becomes natural to select a target in the class of more structured covariance estimators (Briner and Connor 2008). In addition, this strategy allows direct control over the trade-off between estimation error and model error in the resulting shrinkage estimates. In fact, the sample covariance matrix is usually affected by a large amount of estimation error. This is reduced when shrinking towards a structured target which minimizes the sampling error at the cost of adding some misspecification by imposing a specific model. At this point, the shrinkage intensity is crucial because it must be set in such a way to minimize both errors.
To define the shrinkage estimator, we start from the definition of sample covariance matrix S . Recalling Equation (1), S is given by
S =   1 n 1 X ( I n 1 n 1 n 1 n ) X ,
where I n denotes the n × n identity matrix and 1 n is the ones column vector of length n . The shrinkage methodology enhances the sample covariance matrix estimation by shrinking S towards a specific target matrix T :
Σ s =   δ T + ( 1 δ ) S
where Σ s is the shrinkage estimator; δ the shrinkage parameter and T the target matrix. In this work, we focus on the problem of selecting the target matrix. After a review of the literature on target matrices, in the following rows we present the target estimators considered in this study and we assess through a numerical illustration the impact of misspecification in the target matrix for the considered risk-based portfolios.

3.1. Target Matrix Literature Review

The target matrix should fit a desirable number of requirements: First, it should be structured much enough to lower the estimation error of the sample covariance matrix while not bringing too much error from model selection. Second, it should reflect the important features of the true covariance matrix (Ledoit and Wolf 2004b). The crucial question is: how much structure should we impose to fill in the requirements? Table 1 shows the target matrices employed so far in the literature, summarising information about the formula for the shrinkage intensity, the wealth allocation rule, and the addressed research question. Not surprisingly, all the papers shrink the sample covariance matrix. What surprises is that only six target matrices have been examined: the one relying on the Single-Index market model, the Identity matrix, the Scaled Identity, and the Variance Identity, the Constant Correlation model and the Common Covariance model. Previously, four were proposed by Ledoit and Wolf in separate works (Ledoit and Wolf 2003, 2004a, 2004b) and were again proposed in subsequent works, while the Common Covariance appears only in (Pantaleo et al. 2011) and the Scaled Identity only in (De Miguel et al. 2013).
In Table 1, papers have been listed taking into account their contribution to the literature as regards the adoption of a novel target matrix estimator, the re-examination of a previously proposed target, and the comparison among different estimators. Ledoit and Wolf popularised the shrinkage methodology in portfolio selection: in Ledoit and Wolf (2003), they were also the first to compare the effects of shrinking towards different targets in portfolio performance. Shrinking towards the Variance Identity and shrinking towards the Market Model are two out of the eight estimators for the covariance matrix compared with respect to the reduction of estimation error in portfolio weights. They found significant improvements in portfolio performance when shrinking towards the Market Model. Briner and Connor (2008) well described the importance of selecting the target matrix among the class of structured covariance estimators, hence proposing to shrink the asset covariance matrix of demeaned returns towards the Market model as in Ledoit and Wolf (2003). Candelon et al. (2012) compared the effect of double shrinking the sample covariance either towards the Market Model and the Identity, finding that both estimators carry on similar out-of-sample performances. De Miguel et al. (2013) compared the effects of different shrinkage estimators on portfolio performance, highlighting the importance of the shrinkage intensity parameter and proposing a scaled version of the Identity Matrix as a target. Another important comparison among target matrices is that of Pantaleo et al. (2011), who compared the Market and Constant Correlation models as in Ledoit and Wolf (2003, 2004b) with the Common Covariance of Schäfer and Strimmer (2005) implemented as target matrix for the first time in finance. The authors assessed the effects on portfolio performances while controlling for the dimensionality of the dataset, finding that the Common Covariance should not be used when the number of observations is less than the number of assets. Lastly, Ardia et al. (2017) is the only work to implement shrinkage in risk-based portfolios. They shrunk the sample covariance matrix as in Ledoit and Wolf (2003), finding that the Minimum Variance and the Maximum Diversification portfolios are the most affected from covariance misspecification, hence they benefit the most from the shrinkage technique.

3.2. Estimators for the Target Matrix

We consider six estimators for the target matrix: the Identity and the Variance Identity matrix, the Single-index, the Common Covariance, the Constant Correlation and the Exponential Weighted Moving Average (EWMA) models. They are all structured estimators, in the sense that the number of parameters to be estimated is far less the 1 2 p ( p + 1 ) required in the sample covariance case. Compared with the literature, we take into account all the previous target estimators,4 adding to the analysis the EWMA: this estimator well addresses the problem of serial correlation and heteroskedasticity in asset returns.
The identity is a matrix with ones on the diagonal and zeros elsewhere. Choosing the Identity as the target is justified by the fact that is shows good statistical properties: it is always well conditioned and hence invertible (Ledoit and Wolf 2003). Besides the identity, we also consider a multiple of the identity, named the Identity Variance. This is given by:
T V I d I p d i a g ( S ) I p ,
where d i a g ( S ) is the main diagonal of the sample covariance matrix (hence the assets variances) and I p the identity matrix of dimension p.
The Single Index Model (Sharpe 1963) assumes that the returns r t can be described by a one-factor model, resembling the impact of the whole market:
r t = α + β r m k t + ε t ,   with   t = 1 , , n ,
where r m k t is the overall market returns; β is the vector of factor estimates for each asset; α is the market mispricing, and ε t the model error. The Single-Index market model represents a practical way of reducing the dimension of the problem, measuring how much each asset is affected by the market factor. The model implies the covariance structure among asset returns is given by:
T s i s m k t 2 β β + Ω ,
where s m k t 2 is the sample variance of asset returns; β is the vector of beta estimates and Ω contains the residual variance estimates.
The Common Covariance model is aimed at minimizing the heterogeneity of assets variances and covariances by averaging both of them (Pantaleo et al. 2011). Let var i j , i = j and covar i j ,   i j being the variances and covariances of the sample covariance matrix, respectively, their averages are given by:
var ¯ = 1 p k = 1 p var k , i = j ;
covar ¯ = 1 p ( p 1 ) / 2 k = 1 p ( p 1 ) / 2 covar k ,   i j ;
where p is the number of securities. The resulting target matrix T c v has its diagonal elements all equal to the average of the sample covariance, while non-diagonal elements are all equal to the average of sample covariances.
In the Constant Correlation model the main diagonal is filled with sample variances, and elsewhere a constant covariance parameter which is equal for all assets. The matrix can be written according to the following decomposition:
T c c P   d i a g ( S ) P ,
where P is the lower triangular matrix filled with the constant correlation parameter ρ ¯ = 1 p ( p 1 ) / 2 i = 1 p ρ i j for i < j and ones in the main diagonal. Here, d i a g ( S ) represents the main diagonal of the sample covariance matrix.
The EWMA model (J. P. Morgan and Reuters Ltd. 1996) was introduced by JP Morgan’s research team to provide an easy but consistent way to assess portfolio covariance. RiskMetrics EWMA considers the variances and covariance driven by an Integrated GARCH process:
T EWMA , t X X + λ T EWMA , t 1 ,
with T EWMA , 0 =   I p T EWMA , t 1 is the target matrix at time t 1 and λ is the smoothing parameter: the higher λ , the higher the persistence in the variance.

3.3. The Impact of Misspecification in the Target Matrix

We are now going to show to which extent risk-based portfolios can be affected by misspecification in the target matrix. To do so, we provide a numerical illustration, merely inspired by the one in Ardia et al. (2017). Assume an investment universe made by three securities: a sovereign bond (Asset-1), a corporate bond (Asset-2), and equity (Asset-3), we are able to impose an arbitrary structure to the related 3   ×   3 true covariance matrix5. We preliminary recall that Σ can be written according to the following decomposition:
Σ   ( d i a g ( Σ ) ) 1 / 2 P Σ ( diag ( Σ ) ) 1 / 2
where ( d i a g ( Σ ) ) 1 / 2 is a diagonal matrix with volatilities on the diagonal and zeros elsewhere and P Σ is the related correlation matrix, with ones on the diagonal and correlations symmetrically displaced elsewhere. We impose
( Σ 1 , 1 1 / 2 ,   Σ 2 , 2 1 / 2 ,   Σ 3 , 3 1 / 2 ) = ( 0.1 ,   0.1 ,   0.2 )
and
( P Σ ; 1 , 2 ,   P Σ ; 1 , 3 ,   P Σ ; 2 , 3 ) =   ( 0.1 ,   0.2 ,   0.7 )  
hence, the true covariance matrix is:
Σ     [ 0.010 0.001 0.004 0.001 0.010 0.014 0.004 0.014 0.040 ]
Now assume that the true covariance matrix Σ is equal to its shrunk counterpart when δ =   1 2 :
Σ   Σ s = 1 2 S + 1 2 T
that is both the sample covariance matrix S and the target matrix T must be equal to 1 2 Σ and the true target matrix is:
S   T     [ 0.005 0.0005 0.002 0.0005 0.005 0.007 0.002 0.007 0.020 ]
with few algebraic computations, we can obtain the volatilities and correlations simply by applying the covariance decomposition, ending up with
( T 1 , 1 1 / 2 ,   T 2 , 2 1 / 2 ,   T 3 , 3 1 / 2 )   = ( 0.0707 ,   0.0707 ,   0.1414 ) ;
( P T ; 1 , 2 ,   P T ; 1 , 3 ,   P T ; 2 , 3 ) = ( 0.1 ,   0.2 ,   0.7 ) .
In this case, we can conclude that the target matrix T is undervaluing all the covariance and correlation values.
At this point, some remarks are needed. First, as summarised in Table 2, we work out the true risk-based portfolio weights. Weights are differently spread out: the Minimum Variance equally allocates wealth to the first two assets, excluding equities. This because it mainly relies upon the asset variance, limiting the diversification of the resulting portfolio. The remaining portfolios allocate wealth without excluding any asset; however, the Maximum Diversification overvalues Asset-1 assigning to it more than 50% of the total wealth. The Inverse Volatility and Equal-Risk-Contribution seem to maximise diversification under a risk-parity concept, similarly allocating wealth among the investment universe.
Second, assuming Σ as the true covariance matrix allows us to simulate misspecification both in the volatility and in the correlation components of the target matrix T by simply increasing or decreasing the imposed true values. Since we are interested in investigating misspecification impact on the true risk-based portfolio weights, we measure its effects after each shift with the Frobenius norm between the true weights and the misspecified ones:
ω ˜ F 2 =   i = 1 p ω ˜ i 2
where ω ˜ = ω ω ^ .
Third, turning the discussion onto the working aspects of this toy example, we will separately shift the volatility and the correlation of Asset-3, as done in Ardia et al. (2017). The difference with them is that we modify the values in the true target matrix T . Moreover, in order to also gauge how shrinkage intensity affects the portfolio weights, we perform this analysis for 11 values of δ , spanning from 0 to 1 (with step 0.1). This allows us to understand both extreme cases, i.e., when the true covariance matrix is only estimated with the sample estimator ( δ = 0 ) and only with the target matrix ( δ = 1 ). Remember that the true shrinkage intensity is set at δ = 1 2 .
Moving to the core of this numerical illustration, we proceed as follows. First, for what is concerning the volatility, we let T 3 , 3 1 / 2 vary between 0 and 0.5, ceteris paribus. Results are summarised in Figure 1, row 1. As expected, there is no misspecification in all the risk-based portfolio at the initial state T 3 , 3 1 / 2 = 0.1414 , i.e., the true value. All the portfolio weights are misspecified in the range [ 0 ; 0.1414 ) , with the Minimum Variance portfolio showing the greatest departure from the true portfolio weights when the Asset-3 volatility is undervalued below 0.12. The absence of misspecification effects in its weights is due to the initial high-risk attributed to Asset-3; in fact, it is already excluded from the optimal allocation at the initial non-perturbated state. Regarding the other portfolios, their weights show a similar behaviour to the one just described: the Maximum Diversification weights depart from the non-misspecified state to reach the maximum distance from the true weights of 0.4; however, this effect dissipates as soon as the shrinkage intensity grows. The same applies for the Inverse Volatility and the Equal-Risk-Contribution. On the contrary, when volatility is overvalued in the range ( 0.1414 ; 0.5 ] , the Minimum Variance is not misspecified, since Asset-3 is always excluded from the allocation. This fact helps to maintain the stability of its weights: this portfolio is not affected by shifts in the shrinkage intensity when there is over-misspecification. All the remaining portfolios show low levels of misspecification due to diversification purposes. In particular, they react in the same way to shrinkage intensity misspecification, showing an increase in the Frobenius norm especially for low values of Asset-3 variance. A common trait shared by all the considered portfolios is that when weights are estimated with the sample covariance, only the distance from true portfolios is at the maximum.
Second, we assess the correlation misspecification impact. We let the correlation between Asset-3 and Asset-2 ( P T ; 2 , 3 ) vary from 0 to 1, ceteris paribus. In this case, the greatest signs of perturbation are in the Minimum Variance and in the Maximum Diversification portfolios, while the Equal-Risk-Contribution shows far less distortions, as presented in Figure 1, row 2. The Minimum Variance portfolio is again misspecified in one direction: when the correlation parameter is undervalued and the sample covariance matrix dominates the target matrix in the shrinkage. On the other hand, the Maximum Diversification shows the highest departure from no-misspecification levels in both senses. However, as the Maximum Diversification it seems to benefit from high values of shrinkage intensity. The Equal-Risk-Contribution reacts similarly to the Maximum Diversification, but with a far lower level of misspecification, the Inverse Volatility is not affected at all by misspecification in the correlation structure of the target matrix T. This is due to the specific characteristics of Asset-3 and the way in which the Inverse Volatility selects to allocate weights under a risk-parity scheme.
In conclusion, with this numerical illustration, we assess the effects of target matrix misspecification in risk-based portfolios: the four risk-based portfolios react similarly to what previously found in Ardia et al. (2017), even if in our case, shifts originated in the target matrix. The Minimum Variance and the Maximum Diversification portfolios are the most impacted: the weights of the former are severely affected by volatility and covariance shifts undervaluing the true values; the latter shows perturbations in weights when shifts are more extreme. Both portfolios benefit from a higher level of shrinkage intensity. On the other hand, the Inverse Volatility and the Equal-Risk-Contribution weights suffer less from both sources of misspecification. Overall, weights are affected by shifts in the shrinkage intensity: when sample covariance is the estimator ( δ = 0 ), the distance from the true weights stands at the maximum level in all the considered portfolios.

4. Case Study—Monte Carlo Analysis

This section offers a comprehensive comparison of the six target matrix estimators by means of an extensive Monte Carlo study. The aim of this analysis is twofold: (1) assessing estimators’ statistical properties and similarity with the true target matrix; (2) addressing the problem of how selecting a specific target estimator impacts on the portfolio weights. This investigation is aimed at giving a very broad overview about (1) and (2) since we monitor both the p / n ratio and the whole spectrum of shrinkage intensity. We run simulations for 15 combinations of p and n , and for 11 different shrinkage intensities spanning in the interval [ 0   ; 1 ] , for an overall number of 165 scenarios.
The Monte Carlo study is designed as follows. Returns are simulated assuming a factor model is the data generating process, as in MacKinlay and Pastor (2000). In detail, we impose a one-factor structure for the returns generating process:
r t = ξ f t + ε t ; with   t = 1 , , n ,
where f t is the k × 1 vector of returns on the factor, ξ is the p × 1 vector of factor loadings, and ε t the vector of residuals of p length. Under this framework, returns are simulated implying multivariate normality and absence of serial correlation. The asset factor loadings are drawn from a uniform distribution and equally spread, while returns on the single factor are generated from a Normal distribution. The bounds for the uniform distribution and the mean and the variance for the Normal one are calibrated on real market data, specifically on the empirical dataset “49-Industry portfolios” with monthly frequency, available on the Kennet French website6. Residuals are drawn from a uniform distribution in the range [ 0.10 ; 0.30 ] so that the related covariance matrix is diagonal with an average annual volatility of 20%.
For each of the 165 scenarios, we apply the same strategy. First, we simulate the n × p matrix of asset log-returns, then we estimate the six target matrices and their corresponding shrunk matrices Σ ^ s . Finally, we estimate the weights of the four risk-based portfolios. Some remarks are needed. First, we consider the number of assets as p = { 10 ,   50 ,   100 } and number of observations as n = { 60 ,   120 ,   180 ,   3000 ,   6000 } months, which correspond to 5, 10, 15, 250 and 500 years. Moreover, the shrinkage intensity is allowed to vary between their lower and upper bounds as δ = { 0 ,   0.1 ,   0.2 ,   0.3 ,   0.4 ,   0.5 ,   0.6 ,   0.7 ,   0.8 ,   0.9 ,   1 } . For each of the 165 scenarios we run 100 Monte Carlo trials7, giving robustness to the results.
We stress again the importance of Monte Carlo simulations, which allow us to impose the true covariance Σ and hence the true portfolio weights ω . This is crucial because we can compare the true quantities with their estimated counterparts.
With respect to the point (1), we use two criteria to assess and compare the statistical properties of target matrices: the reciprocal 1-norm condition number and the Frobenius Norm. Being the 1-norm condition number defined as:
C N ( A ) =   κ ( A ) = A 1
for a given A. It measures the matrix sensitivity to changes in the data: when it is large, it indicates that a small shift causes important changes, offering a measure of the ill-conditioning of A. Since C N ( A ) takes value in the interval [ 0   ; + ) , it is more convenient to use its scaled version, the R C N ( A ) :
R C N ( A ) = 1 /   κ ( A )
It is defined in the range [ 0   ; 1 ] : the matrix is well-conditioned if the reciprocal condition number is close to 1 and ill-conditioned vice-versa. Under the Monte Carlo framework, we will study its Monte Carlo estimator:
E [ C N ( A m ) ] =   1 M m = 1 M C N ( A m )
where M is the number of Monte Carlo simulations. On the other hand, the Frobenius norm is employed to gauge the similarity between the estimated target matrix and the true one. We define it for the p × p symmetric matrix Z as:
F N ( Z ) = Z F 2 =   i = 1 p j = 1 p z i j 2
In our case, Z = Σ ^ s Σ . Its Monte Carlo estimator is given by the following:
E [ F N ( A m ) ] = 1 M m = 1 M F N ( A m )
Regarding (2), we assess the discrepancy between true and estimated weights again with the Frobenius norm. In addition, we report the values at which the Frobenius norm attains its best results, i.e., when the shrinkage intensity is optimal.

4.1. Main Results

Figure 2 summarises the statistical properties of the various target matrices.
Figure 2 illustrates the reciprocal 1-norm condition number: the matrix is well-conditioned when the value is closer to 1, vice-versa is ill-conditioned the more it tends to zero. Overall, the Identity and the Variance Identity stem for being always well-conditioned: across all the combinations of p and n their reciprocal condition number is always one. Therefore, we focus our analysis on the remaining target matrices. In the case where p = 10 , the Common Covariance dominates the other three alternatives, who perform poorly. As the number of assets increases, the reciprocal condition number deteriorates, especially for the EWMA, which now performs worse than the others, and for the Common Covariance, which is now aligned to the Single-Index and the Constant Correlation model. In conclusion, excluding the Identity and the Variance Identity, the considered targets show poor statistical properties.
Then, we turn to the study of similarity among true and estimated target matrices. Figure 3 represents the Monte Carlo Frobenius norm between the true and the estimated target matrices. The surfaces give a clear overview about the relation among the Frobenius norm itself, the p / n ratio and the shrinkage intensity. Overall, the Frobenius norm is minimised by the Single-Index and the Common Covariance: in these cases, the target matrices are not particularly affected by the shrinkage intensity, while their reactions to increases in the p / n ratio are controversial. In fact, quite surprisingly the distance between true and estimated weights diminishes as both p and n increases. For p = 50 and p = 100 , there is a hump for small p / n values; however, the Frobenius norm increases when p n 1 . Despite the low condition number, the EWMA shows a similar behaviour to the Single-Index and the Constant Correlation target matrices, especially with respect to p / n values. On the other hand, it is more affected by shifts in the shrinkage parameters; the distance from the true weights increases moving towards the target matrix. Lastly, the Identity and the Variance Identity show a similar behaviour: their distances from the true target matrix increase for higher values of δ and p / n . Lastly, the Common Covariance is the most far away from the true target matrix, being very sensitive both to high shrinkage intensity and p / n values.
To conclude, the Identity and the Variance Identity are the most well-conditioned matrices, being stable across all the examined p / n combinations. Nevertheless, the Single-Index and the Common Covariance target matrices show the greatest similarity with the true target matrix minimizing Frobenius norm, while both the Identity and the Variance Identity seem less similar to the true target.

4.1.1. Results on Portfolio Weights

Table 3 and Table 4 present the main results of the Monte Carlo study: for each combination of p and n , we report the Monte Carlo estimator of the Frobenius norm between the true and estimated weights. In particular, Table 3 reports averaged Frobenius norm along with the shrinkage intensity (excluding the case δ = 0 , which corresponds to the sample covariance matrix), while Table 4 lists the minimum values for the optimal shrinkage intensity.
In both tables, we compare the six target matrices by examining one risk-based portfolio at a time and the effect of increasing p for fixed n . Special attention is devoted to the cases when p > n : the high-dimensional sample. We have this scenario only when p = 100 and n = 60 . Here, the sample covariance matrix becomes ill-conditioned (Marčenko and Pastur 1967), thus it is interesting to evaluate gains obtained with shrinkage. The averaged Frobenius norm values in Table 3 give us a general overview about how target matrices perform across the whole shrinkage intensity spectrum in one goal. We aim to understand if, in average terms, shrinking the covariance matrix benefits risk-portfolio weights. On the other hand, the minimum Frobenius norm values help us understanding to what extent the various target matrices can help reproducing the true portfolio weights: the more intensity we need, the better the target is. In both tables, sample values are listed in the first row of each Panel.
Starting from Table 3, Panel A, the Minimum Variance allocation seems better described by the Identity and the Variance Identity regardless of the number of assets p . In particular, we look at the difference between the weights calculated entirely on the sample covariance matrix and those of the targets: the Identity and the Variance Identity are the only estimator to perform better. In fact, shrinking towards the sample is not as bad as shrinking towards the Common Covariance. By increasing n and moving to Panel B, similar results are obtained. This trend is confirmed in Panel C, while in the cases of n = 3000 and n = 6000 , all the estimators perform similarly. Hence, for the Minimum Variance portfolio the Identity matrix works best at reproducing portfolio weights very similar to the true ones. The same conclusions apply for the Maximum Diversification portfolio: when p and n are small, the Identity and the Variance Identity outperform other alternatives. On the other hand, we get very different results for the Inverse Volatility and Equal-Risk-Contribution. Both portfolios seem not gaining benefits from the shrinkage procedure, as the Frobenius norm is very similar to that of the sample covariance matrix for all the target matrices under consideration. This is true for all pairs of p and n . In the high-dimensional case ( p = 100 ; n = 60 ), the Identity matrix works best in reducing the distance between true and estimated portfolio weights, both for the Minimum Variance and Maximum Diversification portfolios. On average, shrinkage does not help too much when alternative target matrices are used; only in the case of Common Covariance is shrinking worse than using the sample covariance matrix. All these effects vanish when we look at the Inverse Volatility and Equal-Risk-Contribution portfolios: here, shrinkage does not help too much, whatever the target is.
Overall, results are in line with the conclusions of the numerical illustrations in Section 3. Indeed, the Minimum Variance portfolio shows the highest distance between true and estimated weights, similar to the Maximum Diversification. Both portfolios are affected by the dimensionality of the sample: shrinkage always help in reducing weights misspecification; it improves in high-dimensional cases. On contrary, estimated weights for the remaining portfolios are close to the true ones by construction, hence, shrinkage does not help too much.
Switching to Table 4, the results illustrate again that the Identity and the Variance Identity attain the best reduction of the Frobenius norm for the Minimum Variance and Maximum Diversification portfolios. If results are similar to those of Table 3 for the former, results for the latter show an improvement in using the shrinkage estimators. The Identity, Variance Identity, Common Covariance, and Constant Correlation target matrices outperform all the alternatives, including the sample estimator, minimising the Frobenius norm in a similar fashion. This is true also for the high-dimensional case. On the contrary, the other two portfolios do not benefit from shrinking the sample covariance matrix, even in high-dimensional samples, confirming the insights from Table 3. Lastly, we look at the shrinkage intensity at which target matrices attain the highest Frobenius norm reduction. Those values are displayed in Figure 4. The intensity is composed of the interval [ 0 ; 1 ] : the more it is close to 1, the more the target matrix helps in reducing the estimation error of the sample covariance matrix. Interestingly, the Identity and the Variance Identity show shrinkage intensities always close to 1, meaning that shrinking towards them is highly beneficial, as they are fairly better than the sample covariance matrix. This is verified either for the high-dimensional case and for those risk portfolios (Inverse Volatility and Equal-Risk-Contribution), who do not show great improvements when shrinkage is adopted.

4.1.2. Sensitivity to Shrinkage Intensity

To have a view on the whole shrinkage intensity spectrum (i.e., the interval ( 0 ; 1 ) ) we refer to Figure 5, where we report the Frobenius Norms for the weights (y-axis) with regard to the shrinkage intensity (x-axis). Each column corresponds to a specific risk-based portfolio: from left to right, the Minimum Variance, the Inverse Volatility, the Equal-Risk-Contribution, and the Maximum Diversification, respectively. Each row corresponds to the p / n ratio in n ascending order. For each subfigure, the Identity is blue circle-shaped, the Variance Identity is green square-shaped, the Single-Index is red hexagram-shaped, the Common Covariance is black star-shaped, the Constant Correlation is cyan plus-shaped, and the EWMA is magenta diamond-shaped.
Figure 5 illustrates the case p = 100 , as to include the high-dimensional scenario. Starting from the latter (first row, n = 60 ), the Variance Identity is the only target matrix to always reduce weight misspecification for all the considered portfolios, for all shrinkage levels. The Identity do the same, excluding the ERC case where it performs worse than the sample covariance matrix. The remaining targets behave very differently across the four risk-based portfolios: the Common Covariance is the worst in both the Minimum Variance and Maximum Diversification and the EWMA is the worst in both remaining portfolios. The Market Model and the Constant Correlation do not improve much from the sample estimator across all portfolios.
Looking at the second row ( n = 120 ), the Identity is the most efficient target, reducing the distance between estimated and true portfolio weights in all the considered portfolios. The Variance Identity is also very efficient in Minimum Variance and Maximum Diversification portfolios, while the remaining targets show similar results as in the previous case. The same conclusions apply for the case n = 180 .
When the number of observations is equal to or higher than n = 3000 , results do not change much. The Identity, the Variance Identity, the Market model, and the Constant Correlation are the most efficient target matrices towards which to shrink, while the EWMA is the worst for both Inverse Volatility and Equal-Risk-Contribution portfolios and the Common Covariance is the worst for the Minimum Variance and Maximum Diversification ones.
In conclusion, for the Minimum Variance portfolio, the Common Covariance should not be used, since it always produces weights very unstable and distant from the true ones. At the same time, the EWMA should not be used to shrink the covariance matrix in the Inverse Volatility and Equal-Risk-Contribution portfolios. The most convenient matrices towards which to shrink are the Identity and the Variance Identity. Overall, the Minimum Variance and the Maximum Diversification portfolio weights gain more from shrinkage than those of the Inverse Volatility and Equal-Risk-Contribution allocations.

5. Conclusions

In this article, we provide a comprehensive overview of shrinkage in risk-based portfolios. Portfolios solely based on the asset returns covariance matrix are usually perceived as “robust” since they avoid estimating the asset returns mean. However, they still suffer from estimation error when the sample estimator is used, causing misspecification in the portfolio weights. Shrinkage estimators have been proved to reduce the estimation error by pulling the sample covariance towards a more structured target.
By means of an extensive Monte Carlo study, we compare six different target matrices: the Identity, the Variance Identity, the Single-index model, the Common Covariance, the Constant Correlation, and the Exponential Weighted Moving Average, respectively. We do so considering their effects on weights for the Minimum Variance, Inverse Volatility, Equal-risk-contribution, and Maximum diversification portfolios. Moreover, we control for the whole shrinkage intensity spectrum and for dataset size, changing observation length and number of assets. Therefore, we are able to (1) assess estimators’ statistical properties and similarity with the true target matrix; (2) address the problem of how selecting a specific target estimator affects the portfolio weights.
Regarding point (1), the findings suggest the Identity and the Variance Identity matrices hold the best statistical properties, being well conditioned across all the combinations of observations/assets, especially for high-dimensional datasets. Nevertheless, these targets are both not very similar to the true target matrix. The Single-Index and the Constant Correlation target matrices show the greater similarity with the true target matrix, minimizing the Frobenius norm, albeit being poorly conditioned when observations and assets share similar sizes. Turning to point (2), the Identity attains the best results in terms of distance reduction between the true and estimated portfolio weights for both the Minimum Variance and Maximum Diversification portfolio construction techniques. The Variance Identity shows a similar performance. Both estimators are also stable against shifts in the shrinkage intensity.
Overall, selecting the target matrix is very important, since we verified that there are large shifts in the distance between true and estimated portfolio weights when shrinking towards different targets. In risk-based portfolio allocations the Identity and the Variance Identity matrices represent the best target among the six considered in this study, especially in the case of Minimum Variance and Maximum Diversification portfolios. In fact, they are always well conditioned and outperform their competitor in deriving the most similar weights to the true ones.
Lastly, the findings confirm that the Minimum Variance and Maximum Diversification portfolios are more sensitive to misspecification in the covariance matrix, therefore, they benefit the most when the sample covariance matrix is shrunk. These findings are in line with what was previously found in Ardia et al. (2017): the Inverse Volatility and the Equal-Risk-Contribution are more robust to covariance misspecification; hence, allocations do not improve significantly when shrinkage is used.

Funding

This research received no external funding.

Acknowledgments

I would thank my Ph.D. supervisor Marina Resta for having spent time in reading the manuscript and giving precious suggestions. I also thank three anonymous referees for helpful comments and effort into improving this paper.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Ardia, David, Guido Bolliger, Kris Boudt, and Jean Philippe Gagnon-Fleury. 2017. The Impact of Covariance Misspecification in Risk-Based Portfolios. Annals of Operations Research 254: 1–16. [Google Scholar] [CrossRef]
  2. Black, Fischer, and Robert Litterman. 1992. Global Portfolio Optimization. Financial Analysts Journal 48: 28–43. [Google Scholar] [CrossRef] [Green Version]
  3. Briner, Beat G., and Gregory Connor. 2008. How Much Structure Is Best? A Comparison of Market Model, Factor Model and Unstructured Equity Covariance Matrices. Journal of Risk 10: 3–30. [Google Scholar] [CrossRef]
  4. Candelon, Bertrand, Christophe Hurlin, and Sessi Tokpavi. 2012. Sampling Error and Double Shrinkage Estimation of Minimum Variance Portfolios. Journal of Empirical Finance 19: 511–27. [Google Scholar] [CrossRef]
  5. Chopra, Vijay Kumar, and William T. Ziemba. 1993. The Effect of Errors in Means, Variances, and Covariances on Optimal Portfolio Choice. The Journal of Portfolio Management 19: 6–11. [Google Scholar] [CrossRef]
  6. Choueifaty, Yves, and Yves Coignard. 2008. Toward Maximum Diversification. The Journal of Portfolio Management 35: 40–51. [Google Scholar] [CrossRef]
  7. De Miguel, Victor, Lorenzo Garlappi, and Raman Uppal. 2009. Optimal versus Naive Diversification: How Inefficient Is the 1/N Portfolio Strategy? Review of Financial Studies 22: 1915–53. [Google Scholar] [CrossRef]
  8. De Miguel, Victor, Alberto Martin-Utrera, and Francisco J. Nogales. 2013. Size Matters: Optimal Calibration of Shrinkage Estimators for Portfolio Selection. Journal of Banking and Finance 37: 3018–34. [Google Scholar] [CrossRef]
  9. J., P. Morgan, and Reuters Ltd. 1996. Risk Metrics—Technical Document, 4th ed. New York: Morgan Guaranty Trust. [Google Scholar]
  10. Jagannathan, Ravi, and Tongshu Ma. 2003. Risk Reduction in Large Portfolios: Why Imposing the Wrong Constraints Helps. Journal of Finance 58: 1651–84. [Google Scholar] [CrossRef]
  11. James, W., and Charles Stein. 1961. Estimation with Quadratic Loss. In Proceedings of the 4th Berkeley Symposium on Probability and Statistics, Volume 1. Paper presented at Fourth Berkeley Symposium on Mathematical Statistics and Probability, Berkeley, CA, USA, June 20–July 30. [Google Scholar]
  12. Jorion, Philippe. 1986. Bayes-Stein Estimation for Portfolio Analysis. The Journal of Financial and Quantitative Analysis 21: 279–92. [Google Scholar] [CrossRef]
  13. Ledoit, Olivier, and Michael Wolf. 2003. Improved Estimation of the Covariance Matrix of Stock Returns with an Application to Portfolio Selection. Journal of Empirical Finance 10: 603–21. [Google Scholar] [CrossRef]
  14. Ledoit, Olivier, and Michael Wolf. 2004a. A Well-Conditioned Estimator for Large-Dimensional Covariance Matrices. Journal of Multivariate Analysis 88: 365–411. [Google Scholar] [CrossRef]
  15. Ledoit, Olivier, and Michael Wolf. 2004b. Honey, I Shrunk the Sample Covariance Matrix—Problems in Mean-Variance Optimization. Journal of Portfolio Management 30: 110–19. [Google Scholar] [CrossRef]
  16. Leote de Carvalho, Raul, Xiao Lu, and Pierre Moulin. 2012. Demystifying Equity Risk–Based Strategies: A Simple Alpha plus Beta Description. The Journal of Portfolio Management 38: 56–70. [Google Scholar] [CrossRef]
  17. MacKinlay, A. Craig, and Lubos Pastor. 2000. Asset Pricing Models: Implications for Expected Returns and Portfolio Selection. Review of Financial Studies 13: 883–916. [Google Scholar] [CrossRef]
  18. Maillard, Sébastien, Thierry Roncalli, and Jérôme Teïletche. 2010. The Properties of Equally Weighted Risk Contribution Portfolios. The Journal of Portfolio Management 36: 60–70. [Google Scholar] [CrossRef]
  19. Marčenko, Vladimir Alexandrovich, and Leonid Andreevich Pastur. 1967. Distribution of Eigenvalues for Some Sets of Random Matrices. Mathematics of the USSR-Sbornik 1: 507–36. [Google Scholar] [CrossRef]
  20. Markowitz, Harry. 1952. Portfolio Selection. The Journal of Finance 7: 77–91. [Google Scholar] [CrossRef]
  21. Markowitz, Harry. 1956. The Optimization of a Quadratic Function Subject to Linear Constraints. Naval Research Logistics Quarterly 3: 111–33. [Google Scholar] [CrossRef]
  22. Merton, Robert C. 1980. On Estimating the Expected Return on the Market: An Exploratory Investigation. Topics in Catalysis 8: 323–61. [Google Scholar]
  23. Michaud, Richard O. 1989. The Markowitz Optimization Enigma: Is Optimized Optimal? Financial Analysts Journal 45: 31–42. [Google Scholar] [CrossRef]
  24. Pantaleo, Ester, Michele Tumminello, Fabrizio Lillo, and Rosario N. Mantegna. 2011. When Do Improved Covariance Matrix Estimators Enhance Portfolio Optimization? An Empirical Comparative Study of Nine Estimators. Quantitative Finance 11: 1067–80. [Google Scholar] [CrossRef]
  25. Qian, Edward. 2006. On the Financial Interpretation of Risk Contribution: Risk Budgets Do Add Up. Journal of Investment Management 4: 41–51. [Google Scholar] [CrossRef]
  26. Schäfer, Juliane, and Korbinian Strimmer. 2005. A Shrinkage approach to large-scale covariance matrix estimation and implications for functional genomics. Applications in Genetics and Molecular Biology 4: 1175–89. [Google Scholar] [CrossRef] [PubMed]
  27. Sharpe, William F. 1963. A simplified model for portfolio analysis. Management Science 9: 277–93. [Google Scholar] [CrossRef]
1
The majority of papers on risk-based portfolios are published in journal aimed at practitioners, as the Journal of Portfolio Management.
2
With this we refer to the population covariance matrix, which by definition is not observable and then unfeasible. Hence, Σ is estimated taking into account the observations stored in X: we will deeply treat this in the next section.
3
The sample covariance matrix is the Maximum Likelihood Estimator (MLE) under Normality, therefore it lets data speaks without imposing any structure.
4
In reality, we exclude the Scaled Identity of De Miguel et al. (2013) because of its great similarity with the Identity and Variance Identity implemented in our study.
5
Ardia et al. (2017) imposes Asset-1 and Asset-2 to have 10% annual volatility; Asset-3 to have 20% annual volatility; correlations between Asset-1/Asset-2 and Asset-1/Asset-3 are set as negative and correlation between corporate bonds and equities (Asset-2/Asset-3) is set as positive. However, to better resemble real data, specifically the S&P500, the US corporate index and the US Treasury Index total returns, we assume all three correlation parameters to be positive.
6
7
Simulations were done in MATLAB setting the random seed generator at its default value, thus ensuring the full reproducibility of the analysis. Related code available at the GitHub page of the author: https://github.com/marconeffelli/Risk-Based-Portfolios.
Figure 1. Frobenius norm between true and estimated weights; first row reports misspecification in volatility, while second row in correlation. The surfaces’ three dimensions are: the shrinkage intensity in y axis (from 0 to 1); the misspecification in the volatility (from 0 to 0.5) or in the correlation (from 0 to 1) in x axis and the Frobenius norm in z axis. Each column refers to a specific risk-based portfolio. From the left to the right: Minimum Variance (MV), Inverse Volatility (IV), Equal-Risk-Contribution (ERC), Maximum Diversification (MD), respectively.
Figure 1. Frobenius norm between true and estimated weights; first row reports misspecification in volatility, while second row in correlation. The surfaces’ three dimensions are: the shrinkage intensity in y axis (from 0 to 1); the misspecification in the volatility (from 0 to 0.5) or in the correlation (from 0 to 1) in x axis and the Frobenius norm in z axis. Each column refers to a specific risk-based portfolio. From the left to the right: Minimum Variance (MV), Inverse Volatility (IV), Equal-Risk-Contribution (ERC), Maximum Diversification (MD), respectively.
Risks 06 00125 g001
Figure 2. The reciprocal 1-norm condition number (y-axis) as the p / n ratio moves from p 60 to p 6000 (x-axis). Each column corresponds to a specific target matrix: from left to right, the Identity (Id): blue circle-shaped; the Variance Identity (VId): green square-shaped; the Single-Index (SI): red hexagram-shaped; the Common Covariance (CV): black star-shaped; the Constant Correlation (CC): cyan plus-shaped; and the Exponential Weighted Moving Average (EWMA): magenta diamond-shaped, respectively. Each row corresponds to a different p : in ascendant order from 10 (first row) to 100 (third row).
Figure 2. The reciprocal 1-norm condition number (y-axis) as the p / n ratio moves from p 60 to p 6000 (x-axis). Each column corresponds to a specific target matrix: from left to right, the Identity (Id): blue circle-shaped; the Variance Identity (VId): green square-shaped; the Single-Index (SI): red hexagram-shaped; the Common Covariance (CV): black star-shaped; the Constant Correlation (CC): cyan plus-shaped; and the Exponential Weighted Moving Average (EWMA): magenta diamond-shaped, respectively. Each row corresponds to a different p : in ascendant order from 10 (first row) to 100 (third row).
Risks 06 00125 g002
Figure 3. Surfaces representing the Frobenius norm (z-axis) between the true and the estimated target matrices, considering the shrinkage intensity (y-axis) and the p / n ratio (x-axis). Each column corresponds to a specific target matrix: from left to right, the Identity (Id), the Variance Identity (VId), the Single-Index (SI), the Common Covariance (CV), the Constant Correlation (CC), and the EWMA, respectively. Each row corresponds to a different p : in ascendant order from p = 10 (first row) to p = 100 (third row).
Figure 3. Surfaces representing the Frobenius norm (z-axis) between the true and the estimated target matrices, considering the shrinkage intensity (y-axis) and the p / n ratio (x-axis). Each column corresponds to a specific target matrix: from left to right, the Identity (Id), the Variance Identity (VId), the Single-Index (SI), the Common Covariance (CV), the Constant Correlation (CC), and the EWMA, respectively. Each row corresponds to a different p : in ascendant order from p = 10 (first row) to p = 100 (third row).
Risks 06 00125 g003
Figure 4. Optimal shrinkage intensity parameters for which the Frobenius norm is minimized. Each column corresponds to a specific target matrix: from left to right, the Identity, the Variance Identity, the Single-Index, the Common Covariance, the Constant Correlation, and the EWMA, respectively. Each row corresponds to a different n : in ascendant order from n = 60 (first row) to p = 6000 (fifth row). In each subplot, the MV portfolio is blue circle-shaped; the IV is green-square shaped; the ERC is red-triangle shaped; and the MD is black-cross shaped.
Figure 4. Optimal shrinkage intensity parameters for which the Frobenius norm is minimized. Each column corresponds to a specific target matrix: from left to right, the Identity, the Variance Identity, the Single-Index, the Common Covariance, the Constant Correlation, and the EWMA, respectively. Each row corresponds to a different n : in ascendant order from n = 60 (first row) to p = 6000 (fifth row). In each subplot, the MV portfolio is blue circle-shaped; the IV is green-square shaped; the ERC is red-triangle shaped; and the MD is black-cross shaped.
Risks 06 00125 g004
Figure 5. Frobenius norm for portfolio weights with regard to the shrinkage intensity parameter, when p = 100 .
Figure 5. Frobenius norm for portfolio weights with regard to the shrinkage intensity parameter, when p = 100 .
Risks 06 00125 g005
Table 1. Literature review of target matrices. “SCVm” = sample covariance matrix. “N.A.” = not available. “GMVP” = Global Minimum Variance Portfolio.
Table 1. Literature review of target matrices. “SCVm” = sample covariance matrix. “N.A.” = not available. “GMVP” = Global Minimum Variance Portfolio.
ReferenceMatrix to ShrinkTarget MatrixShrinkage IntensityPortfolio Selection RuleResearch Question
(Ledoit and Wolf 2003)SCVmMarket Model and Variance IdentityRisk-function minimisationClassical Markowitz problemPortfolio Performance comparison
(Ledoit and Wolf 2004a)SCVmIdentityRisk-function minimisationN.A.Theoretical paper to gauge the shrinkage asymptotic properties
(Ledoit and Wolf 2004b)SCVmConstant Correlation ModelOptimal shrinkage constantClassical Markowitz problemPortfolio Performance comparison
(Briner and Connor 2008)Demeaned
SCVm
Market ModelSame as (Ledoit and Wolf 2004b)N.A.Analysis of the trade-off estimation error and model specification error
(Pantaleo et al. 2011)SCVmMarket Model, Common Covariance and Constant Correlation ModelUnbiased estimator of (Schäfer and Strimmer 2005)Classical Markowitz problemPortfolio Performance comparison
(Candelon et al. 2012)SCVmMarket Model and IdentitySame as (Ledoit and Wolf 2003, 2004b)Black-Litterman GMVPPortfolio Performance comparison
(De Miguel et al. 2013)SCVmScaled IdentityExpected quadratic loss and bootstrapping approachClassical Markowitz problemComprehensive investigation of shrinkage estimators
(Ardia et al. 2017)SCVmMarket ModelSame as (Ledoit and Wolf 2003)Risk-based portfoliosTheoretical paper to assess effect on risk-based weights
Table 2. True weights of the four risk-based portfolios.
Table 2. True weights of the four risk-based portfolios.
AssetMinimum Variance (MV)Inverse Volatility (IV)Equal-Risk-Contribution (ERC)Maximum Diversification (MD)
Asset-10.5000.4000.4480.506
Asset-20.5000.4000.3740.385
Asset-30.0000.2000.1770.108
Table 3. Frobenius norm for the portfolio weights. Values are averaged along with the shrinkage intensity (excluding the case δ = 0 ). For each n , the first line reports the Frobenius norm for the sample covariance matrix. Abbreviations in use are: S for sample covariance; Id for identity matrix; VId for Variance Identity; SI for Single-Index; CV for Common Covariance; CC for Constant Correlation and EWMA for Exponentially Weighted Moving Average.
Table 3. Frobenius norm for the portfolio weights. Values are averaged along with the shrinkage intensity (excluding the case δ = 0 ). For each n , the first line reports the Frobenius norm for the sample covariance matrix. Abbreviations in use are: S for sample covariance; Id for identity matrix; VId for Variance Identity; SI for Single-Index; CV for Common Covariance; CC for Constant Correlation and EWMA for Exponentially Weighted Moving Average.
P = 10P = 50P = 100
MVIVERCMDMVIVERCMDMVIVERCMD
Panel A: n = 60
S0.8340.15850.17360.58420.77210.05730.06370.49330.75550.04090.04470.4565
Id0.68630.14250.15280.50450.62150.05590.06310.38730.49670.04040.04510.3652
VId0.69350.15830.17320.51760.59990.05670.06340.40920.59010.04040.04450.3686
SI0.8380.15850.17360.56780.76850.05730.06370.47090.750.04090.04470.4288
CV1.24380.15830.17311.0111.14840.05670.06280.93811.13860.04040.04380.9185
CC0.83530.15850.17330.53610.78080.05730.06350.43280.76630.04090.04450.3922
EWMA0.84730.15930.17450.5950.78110.05750.0640.51420.73250.04110.0450.4431
Panel B: n = 120
S0.90640.08770.09890.46490.78140.0590.06560.50650.65190.04240.04720.4332
Id0.81570.0870.09830.42560.62590.06130.06880.43540.63070.03890.04310.328
VId0.82350.08710.09850.42840.62590.06130.06880.43540.4890.04210.04710.3712
SI0.90970.08770.09890.45630.77770.0590.06560.49250.64580.04240.04720.419
CV1.32690.08710.09820.96671.18060.05870.06511.01381.09740.04210.04670.8951
CC0.9050.08770.09880.43570.78220.0590.06550.46360.65660.04240.04710.3856
EWMA0.92810.08830.09960.48590.79940.05920.06580.52460.67880.04270.04750.4601
Panel C: n = 180
S0.79890.13110.14230.50070.79320.05640.06270.46310.69050.04040.0440.4065
Id0.72060.13080.1420.47360.67050.05620.06250.4050.54770.03750.03990.3748
VId0.72730.13080.14210.47570.68380.05620.06260.41270.57540.04020.0440.3556
SI0.80010.13110.14230.49540.79040.05640.06270.45450.68730.04040.0440.3982
CV1.27150.13080.14190.99611.20730.05620.06240.99881.14220.04020.04370.8705
CC0.79570.13110.14230.48030.7920.05640.06260.42590.6920.04040.0440.3672
EWMA0.84150.13220.14350.5260.82840.05670.06310.50050.72060.04080.04450.4429
Panel D: n = 3000
S0.75040.14760.15960.39570.7340.0490.05390.39880.5130.03840.04280.3259
Id0.74410.14770.15970.39460.70090.0490.05390.38720.46150.03840.04280.3096
VId0.74370.14770.15960.39450.70430.0490.05390.38860.46730.03840.04280.312
SI0.75160.14760.15960.39550.73390.0490.05390.39840.51230.03840.04280.3252
CV1.28640.14770.15970.9631.22810.0490.05380.99541.10410.03840.04280.6822
CC0.74880.14760.15960.39490.73160.0490.05390.39040.50960.03840.04280.3143
EWMA0.85630.14890.16110.44520.81610.04970.05470.46520.62440.03890.04350.4076
Panel E: n = 6000
S0.96720.13020.14090.48210.57370.05390.05890.34810.57720.04020.04370.3436
Id0.94960.13010.14080.48130.60950.05750.06390.40760.54490.04020.04370.3342
VId0.9510.13010.14090.48150.54190.0540.05890.34010.54830.04020.04370.3354
SI0.96880.13020.14090.4820.5740.05390.05890.34790.57720.04020.04370.3434
CV1.41420.13010.14081.00341.14360.0540.05890.97061.14220.04020.04370.7031
CC0.96560.13020.14090.48140.57090.05390.05890.34150.5750.04020.04370.3368
EWMA1.04320.13120.14220.52320.69460.05470.05990.43190.6810.04070.04440.4229
Table 4. Frobenius norm for the portfolio weights. Values corresponds to the optimal shrinkage intensity, listed after the Frobenius norm for each portfolio. We report values for the sample covariance matrix ( δ = 0 ) separately in the first row of each panel. For each n , the first line reports the Frobenius norm for the sample covariance matrix. Abbreviations used are: S for sample covariance; Id for identity matrix; VId for Variance Identity; SI for Single-Index; CV for Common Covariance; CC for Constant Correlation and EWMA for Exponentially Weighted Moving Average.
Table 4. Frobenius norm for the portfolio weights. Values corresponds to the optimal shrinkage intensity, listed after the Frobenius norm for each portfolio. We report values for the sample covariance matrix ( δ = 0 ) separately in the first row of each panel. For each n , the first line reports the Frobenius norm for the sample covariance matrix. Abbreviations used are: S for sample covariance; Id for identity matrix; VId for Variance Identity; SI for Single-Index; CV for Common Covariance; CC for Constant Correlation and EWMA for Exponentially Weighted Moving Average.
P = 10P = 50P = 100
MVIVERCMDMVIVERCMDMVIVERCMD
Panel A: n = 60
S0.83400.15850.17360.58420.77210.05730.06370.49330.75550.04090.04470.4565
Id0.67780.14240.15250.5010.59970.05580.06240.37040.4710.04030.04460.3462
VId0.66890.15810.1730.50840.55390.05650.06270.37950.54280.04020.04370.3331
SI0.83450.15850.17350.5580.76660.05730.06370.46330.74790.04090.04470.4195
CV1.23920.15810.17290.5091.1170.05650.06270.37951.10680.04020.04370.3331
CC0.83350.15850.17310.50810.77330.05730.06340.37950.7570.04090.04440.3332
EWMA0.83310.15860.17370.58520.77060.05730.06370.49530.72130.04090.04470.4395
Panel B: n = 120
S0.90640.08770.09890.46490.78140.0590.06560.50650.65190.04240.04720.4332
Id0.81210.0870.09810.42410.61190.06130.06850.42550.6130.03880.04280.3111
VId0.81210.0870.09820.42420.61190.06130.06850.42550.44250.0420.04670.3445
SI0.9070.08770.09890.45260.7760.0590.06560.48720.64310.04240.04720.414
CV1.32690.0870.09810.42451.17560.05860.06510.43021.09160.0420.04670.3445
CC0.90430.08770.09870.42410.7810.0590.06540.43020.65270.04240.04710.3446
EWMA0.90520.08760.09880.46510.77970.05890.06550.50560.65540.04240.04720.4331
Panel C: n = 180
S0.79890.13110.14230.50070.79320.05640.06270.46310.69050.04040.0440.4065
Id0.71770.13070.14190.47240.66130.05620.06240.39770.5340.03750.03980.3645
VId0.7180.13070.14190.47240.66140.05620.06240.39790.54280.04020.04370.3331
SI0.7990.13110.14230.49290.78970.05640.06270.45150.68630.04040.0440.3955
CV1.27150.13070.14180.47241.20730.05620.06240.39791.14220.04020.04370.3331
CC0.79420.13110.14220.47250.79120.05640.06260.39770.69040.04040.04390.3331
EWMA0.80350.13120.14240.50080.79510.05640.06260.46530.69380.04040.0440.4074
Panel D: n = 3000
S0.75040.14760.15960.39570.7340.0490.05390.39880.5130.03840.04280.3259
Id0.74250.14770.15960.39410.69880.0490.05380.38590.45730.03840.04280.3072
VId0.74260.14760.15960.39410.69880.0490.05380.38590.45730.03840.04280.3072
SI0.75060.14760.15960.39530.73390.0490.05390.39830.5120.03840.04280.325
CV1.28640.14760.15960.39511.22810.0490.05380.38591.10410.03840.04280.3072
CC0.74770.14760.15960.39460.72990.0490.05390.3860.50730.03840.04280.3072
EWMA0.76150.14770.15970.39810.74390.04910.05390.40430.52630.03840.04290.3346
Panel E: n = 6000
S0.96720.13020.14090.48210.57370.05390.05890.34810.57720.04020.04370.3436
Id0.94860.130.14080.48110.60850.05750.06390.40720.54280.04020.04370.3331
VId0.94860.130.14080.48110.53650.0540.05890.33810.54280.04020.04370.3331
SI0.96750.13020.14090.4820.57380.05390.05890.34780.57720.04020.04370.3433
CV1.41420.130.14080.48111.14360.0540.05890.33811.14220.04020.04370.3331
CC0.96440.13020.14090.48120.56870.05390.05890.33810.57330.04020.04370.3331
EWMA0.97650.13020.14090.48320.59010.0540.0590.35610.590.04020.04380.3524

Share and Cite

MDPI and ACS Style

Neffelli, M. Target Matrix Estimators in Risk-Based Portfolios. Risks 2018, 6, 125. https://doi.org/10.3390/risks6040125

AMA Style

Neffelli M. Target Matrix Estimators in Risk-Based Portfolios. Risks. 2018; 6(4):125. https://doi.org/10.3390/risks6040125

Chicago/Turabian Style

Neffelli, Marco. 2018. "Target Matrix Estimators in Risk-Based Portfolios" Risks 6, no. 4: 125. https://doi.org/10.3390/risks6040125

APA Style

Neffelli, M. (2018). Target Matrix Estimators in Risk-Based Portfolios. Risks, 6(4), 125. https://doi.org/10.3390/risks6040125

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