# Risk Parity portfolio construction

If I would like to construct a fully invested long only portfolio with two asset classes (Bonds $B$ and Stocks $S$) based on the concept of 'risk parity' the weights $W$ of my portfolio would be the following:

Then the weight of the bonds $W_B = \textrm{Vol}(S)/[\textrm{Vol(S)}+\textrm{Vol(B)}]$ and the weights of the stocks $W_S = 1 - W_B$.

Based on this I am going to overweigh the low volatility asset and underweight the high volatility asset. My question is: How to calculate the weights for a portfolio with multiple asset classes, 5 for example, so that each asset class will have the same volatility and contribute the same amount of risk into my portfolio. From historical data I can extract the volatility of each asset class and the correlation between them.

-
Can you show us what you tried and define some variables and equations? In its current form the question is off topic IMHO, see the FAQ. – Bob Jansen Mar 22 '12 at 21:35
You can offset some of the "diversification" (it's diversification only if the numbers hold during high stress periods) by raising the leverage on the low volatility assets. – bill_080 Mar 22 '12 at 22:24
@Bootvis: I don't think it's OT. But the formatting could certainly be improved. But the subject is non-trivial. – SRKX Mar 23 '12 at 7:10
It certainly is an interesting topic but the question, as it is now, does not seem to be written by a professional quant. I would edit the question if I had the time. – Bob Jansen Mar 23 '12 at 8:05

Risk Parity is not about "having the same volatility", it is about having each asset contributing in the same way to the portfolio overall volatility.

The volatility of the portfolio is defined as:

$$\sigma(w)=\sqrt{w' \Sigma w}$$

The risk contribution of asset $i$ is computed as follows:

$$\sigma_i(w)= w_i \times \partial_{w_i} \sigma(w)$$

You can then show that:

$$\sigma(w)=\sum_{i=1}^n \sigma_i(w)$$

The vector of the marginal contributions ($\partial_{w_i} \sigma(w)$) is computed as follows:

$$c(w)= \frac{\Sigma w}{\sqrt{w' \Sigma w}}$$

You can then find the solution by running the following optimization:

$$\underset{w}{\arg \min} \sum_{i=1}^N [\frac{\sqrt{w^T \Sigma w}}{N} - w_i \cdot c(w)_i]^2$$ This article contains all the developments you require to understand how the formulas above are derived.

-
Can you explain what techniques are needed to run that optimization? – nxstock-trader Mar 28 '12 at 22:09
You can basically run this through fmincon in MATLAB for example. Not sure what you mean by "techniques". Are you looking for a specific optimization algorithm? – SRKX Mar 28 '12 at 22:13
I meant what packages/routines to use if I were doing this in R? – nxstock-trader Mar 29 '12 at 1:01
@nxstock-trader: you should be able to find something on this page. I haven't used R for optimization for a long time. You can ask on Mathematics or Stack Overflow as well. – SRKX Mar 29 '12 at 6:09

I am very happy with the following equivalent formulation for the risk budgeting problem (as presented in Bruder, Roncalli, 2012, Managing Risk Exposures using the Risk Budgeting Apporach):

Let $b_i$, $\Sigma_{i=1}^n b_i =1$ be the risk budgets, $y_i$ the unscaled portfolio weights and $S$ the variance covariance matrix and $c$ arbitrary.

$$y^* = \text{arg min}_y \sqrt{y^T S y}, \quad \text{s.t.} \sum_{i=1}^n b_i \ln y_i \geq c, \quad \sum_{i=1}^ny_i=1, \quad y_i \geq 0$$

Now the good thing about this formulation is: It is a quadratic program with convex constraints (assuming $b_i >0$) which is numerically nice. Further more, for numerical implementation one would like to drop the constraint $\sum_{i=1}^ny_i=1$ and manually rescale afterwards $x_i^* = \frac{y_i^*}{\sum_{i=1}^ny_i^*}$. It works better for me than the solution presented in the other answer.

-
Yes, this is a more efficient numerical approach I think. I did not use it in my answer because I find it less intuitive. I'd just add to set $b_i = \frac{1}{n} ~ \forall i$ if he wants an ERC... – SRKX Jul 21 '14 at 13:59