What the hell is a variance matrix?

When I first came to finance, I kept hearing about “risk models”. I wondered, “What the hell is a risk model?” Of course, I didn’t say this out loud — that would have given the game away.  My wife has strict instructions that she is to be the only one to know that I’m an idiot.

Finally enlightenment came: “Oh, you mean a variance matrix.”

So let’s look at what a variance matrix (also called a covariance matrix or a variance-covariance matrix) is, or a risk model if you prefer.

The big idea is that we want something that shows the relationship between each pair of returns. Note “returns” and not “prices”.  You can do the calculations with prices instead of returns, but you end up with useless garbage.

The variance matrix is square with a row and a column — in our case — for each asset.  In practice the number of assets can range from a few to a few thousand.

The diagonal elements of the matrix are the variances of the assets.  If the variance matrix is annualized, then these diagonal elements are the squared volatilities.  It would be nicer if we just had volatilities instead of their squares.  Contrary to popular belief statisticians are not in general sadists — or masochists.  The volatilities are squared because, unpleasant as it is, that is actually the easiest path.

The off-diagonal elements are covariances.  Again, if we have the scaling on an annual basis, then a covariance is the volatility of the first asset times the volatility of the second asset times the correlation of the two assets.  The matrix is symmetric, and hence redundant — the covariance of the first asset with the second asset is the same as the covariance of the second asset with the first asset.

Correlations go from -1 to 1. Here is some R code that will show you what a sample of random normals looks like that have a particular correlation:

> require(MASS)
> mycor <- 0.5
> plot(mvrnorm(100, c(0,0), matrix(c(1,mycor,mycor,1), 2), empirical=TRUE))

Figure 1: 100 random normals with sample correlation = 0.5
100 random normals with sample correlation = 0.5Figure 1 shows an example of the plot command in action.  (The figure cheats slightly by setting xlab and ylab to the empty string.)  You can recall the plot command and re-execute it several times to get a sense of how variable the same sample correlation can look.  You can then change the correlation or the sample size.

Now that we know what a variance matrix is, how do we go about getting one?  We have to estimate it. In the abstract statistical setting we estimate a variance matrix using a matrix of data.  The rows are observations, the columns are variables.  The result is a variance matrix that is number of variables by number of variables.  The assumption is that each observation is independent from the others, and they all have the same distribution.  In particular they all come from the same variance matrix.

We’ll start with a matrix of returns: times (trading days, for instance) in the rows, assets in the columns.  But here is where we and textbook statistics part ways.

The means of returns at different times are (presumably) not the same, like statisticians demand.  But they are close enough that this issue is not worth worrying about.

In market data there is the phenomenon of volatility clustering.  Volatility jumps up during crises and then gradually dies back down again.  The dying down is decidedly non-smooth.

Figure 2: Approximate S&P 500 volatility for 2000 days up to 2010-Aug-23
Figure 2 provides a feel for volatility clustering.  The estimate of volatility was done with a rather crude model, but if we magically had a picture of the true volatility, the overall features would be very similar.

Correlations are dynamic as well.  In general they also increase during crises.

For these reasons using as long of a history as possible to estimate the variance matrix is not a good plan.  We only want to use current data.  What does “current” mean?  Excellent question.

The classical statistics set-up demands that we have more observations than variables (more time periods than assets in our case).  It isn’t that we can’t do the calculation if there are fewer times than assets, it’s that we might not like what we get.  The technical terms for what we get is a singular matrix — one that is not positive-definite.  The practical effect is that the variance matrix says that there exist portfolios with zero volatility. That’s a place not to go — there be dragons.

There are several ways to get around this problem.  We’ll save that discussion for another post or two.

One more caution that you’ll find in the statistics literature is that correlation is only about linear relationships.  There is the possibility of zero correlation, but a strong relationship.  Perhaps there are other opinions, but I don’t think we need to worry about this. Relationships are continuously in flux so we can only ever get an approximation anyway.  A linear approximation is about all we can hope for.

Let”s review:

  • The true variance matrix is constantly changing.
  • In many cases we will have fewer time periods than assets.

The first point suggests that what we want (most often) is a prediction for some future time period rather than an estimation of the past.  The second suggests that that prediction is going to be noisy.  Very noisy to a statistician’s eye.

Using a financial perspective, the prediction error of the variance matrix is rather modest though.  Predicting returns is close to impossible, so relatively speaking the variance matrix is predicted with brilliant accuracy.

Subscribe to the Portfolio Probe blog by Email

This entry was posted in R language, Risk, Statistics and tagged , , . Bookmark the permalink.

16 Responses to What the hell is a variance matrix?

  1. Pat says:

    It would have been good to have had a specific example of a variance matrix in this post. There is such an example in Implied alpha — almost wordless.

  2. Pingback: Anomalies meet volatility | Portfolio Probe | Generate random portfolios. Fund management software by Burns Statistics

  3. Nice post.

    Simply two suggestions about your conclusions:

    - “The true variance matrix is constantly changing”. True. But there exist recent estimation methods that allow for time-varying volatilities and correlations. See for instance the dynamic conditional correlation (DCC) model of Engle (2002).

    - “In many cases we will have fewer time periods than assets”. True. In these cases, and even with more time periods than assets, the best estimation method I know is the shrinkage covariance matrix by Ledoit and Wolf.

    Finally, I totally agree with your last sentence. In Finance, the error incurred when estimating mean asset returns is much larger than that incurred when estimating the covariance matrix.

    I found your blog recently, it is very interesting. I really like it.

    • Pat says:


      Thanks for the comment. I agree with both of your points. The plan is for future posts to cover those issues.

      A quick comment on the Ledoit-Wolf shrinkage estimate: I think that we need more evidence, but I’d be a bit surprised if Ledoit-Wolf doesn’t turn out to be very good for a lot of applications. An R implementation of Ledoit-Wolf is available in the BurStFin package (which still is not yet on CRAN as of this writing, but available via burns-stat.com).

  4. Pat,

    I totally agree with you regarding the good performance of the Ledoit-Wolf method. I have used it in many portfolio strategies and under many different datasets, and the results are consistently acceptable, even with a very large number of assets.

    Thanks for the R recommendation, although I use the Matlab version from Wolf’s web page.

    Looking forward to reading your future posts.

  5. Pat says:

    A subsequent post explains (well, at least tries to) what goes wrong when you use prices instead of returns in the estimation process. The post is The number 1 novice quant mistake.

  6. Joseph says:

    Hi Pat,

    I am a student using R and am just starting out. I am trying to graph my dcc-garch estimations but cannot figure out how to do it. I was just wondering what command do I use to plot figure 2 above?


  7. Pingback: Information flows like water | Portfolio Probe | Generate random portfolios. Fund management software by Burns Statistics

  8. Pingback: Unproxying weight constraints | Portfolio Probe | Generate random portfolios. Fund management software by Burns Statistics

  9. Pingback: Factor models of variance in finance | Portfolio Probe | Generate random portfolios. Fund management software by Burns Statistics

  10. Pingback: The number 1 novice quant mistake | Portfolio Probe | Generate random portfolios. Fund management software by Burns Statistics

  11. Pingback: The basics of Value at Risk and Expected Shortfall | Portfolio Probe | Generate random portfolios. Fund management software by Burns Statistics

  12. Pingback: Blog year 2010 in review | Portfolio Probe | Generate random portfolios. Fund management software by Burns Statistics

  13. Pingback: Variance matrix differences | Portfolio Probe | Generate random portfolios. Fund management software by Burns Statistics

  14. Pingback: Blog year 2013 in review | Portfolio Probe | Generate random portfolios. Fund management software by Burns Statistics

Leave a Reply

Your email address will not be published. Required fields are marked *

You may use these HTML tags and attributes: <a href="" title=""> <abbr title=""> <acronym title=""> <b> <blockquote cite=""> <cite> <code> <del datetime=""> <em> <i> <q cite=""> <strike> <strong>