Calculating Value at Risk (VaR) of a stock portfolio using Python
What is Value at risk (VaR)?
Value at risk (VaR) is a statistic used to try and quantify the level of financial risk within a firm or portfolio over a specified time frame. VaR provides an estimate of the maximum loss from a given position or portfolio over a period of time, and you can calculate it across various confidence levels.
Estimating the risk of a portfolio is important to longterm capital growth and risk management, particularly within larger firms or institutions. VaR is typically framed as something like this:
 "We have a portfolio VaR of 250,000 USD over the next month at 95% confidence"
 This means that, with 95% confidence, we can say that the portfolio's loss will not exceed 250,000 USD in a month
In this post I'll walk you through the steps to calculate this metric across a portfolio of stocks.
How is VaR calculated?
There are two main ways to calculate VaR:
 Using Monte Carlo simulation
 Using the variancecovariance method
In this post, we'll focus on using method (2) (variancecovariance). In short, the variancecovariance method looks at historical price movements (standard deviation, mean price) of a given equity or portfolio of equities over a specified lookback period, and then uses probability theory to calculate the maximum loss within your specified confidence interval. You can read more detail here, but we'll calculate it step by step below using Python.
Before we get started, note that the standard VaR calculation assumes the following:

Normal distribution of returns  VaR assumes the returns of the portfolio are normally distributed. This is of course not realistic for most assets, but allows us to develop a baseline using a much more simplistic calculation.

(Modifications can be made to VaR to account for different distributions, but here we'll focus on the standard VaR calculation)

Standard market conditions  Like many financial instruments, VaR is best used for considering loss in standard markets, and is not wellsuited for extreme/outlier events.
Steps to calculate the VaR of a portfolio
In order to calculate the VaR of a portfolio, you can follow the steps below:
 Calculate periodic returns of the stocks in the portfolio
 Create a covariance matrix based on the returns
 Calculate the portfolio mean and standard deviation
(weighted based on investment levels of each stock in portfolio)
 Calculate the inverse of the normal cumulative distribution (PPF) with a specified confidence interval, standard deviation, and mean
 Estimate the value at risk (VaR) for the portfolio by subtracting the initial investment from the calculation in step (4)
1) Calculate periodic returns of the stocks in the portfolio
import pandas as pd
from pandas_datareader import data as pdr
import fix_yahoo_finance as yf
import numpy as np
import datetime as dt
# Create our portfolio of equities
tickers = ['AAPL','FB', 'C', 'DIS']
# Set the investment weights (I arbitrarily picked for example)
weights = np.array([.25, .3, .15, .3])
# Set an initial investment level
initial_investment = 1000000
# Download closing prices
data = pdr.get_data_yahoo(tickers, start="20180101", end=dt.date.today())['Close']
#From the closing prices, calculate periodic returns
returns = data.pct_change()
returns.tail()
2) Create a covariance matrix based on the returns
# Generate VarCov matrix
cov_matrix = returns.cov()
cov_matrix
This will allow us to calculate the standard deviation and mean of returns across the entire portfolio.
3) Calculate the portfolio mean and standard deviation
# Calculate mean returns for each stock
avg_rets = returns.mean()
# Calculate mean returns for portfolio overall,
# using dot product to
# normalize individual means against investment weights
# https://en.wikipedia.org/wiki/Dot_product#:~:targetText=In%20mathematics%2C%20the%20dot%20product,and%20returns%20a%20single%20number.
port_mean = avg_rets.dot(weights)
# Calculate portfolio standard deviation
port_stdev = np.sqrt(weights.T.dot(cov_matrix).dot(weights))
# Calculate mean of investment
mean_investment = (1+port_mean) * initial_investment
# Calculate standard deviation of investmnet
stdev_investment = initial_investment * port_stdev
Next, we can plug these variables into our percentage point function (PPF) below.
4) Calculate the inverse of the normal cumulative distribution (PPF) with a specified confidence interval, standard deviation, and mean
# Select our confidence interval (I'll choose 95% here)
conf_level1 = 0.05
# Using SciPy ppf method to generate values for the
# inverse cumulative distribution function to a normal distribution
# Plugging in the mean, standard deviation of our portfolio
# as calculated above
# https://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.norm.html
from scipy.stats import norm
cutoff1 = norm.ppf(conf_level1, mean_investment, stdev_investment)
5) Estimate the value at risk (VaR) for the portfolio by subtracting the initial investment from the calculation in step 4
#Finally, we can calculate the VaR at our confidence interval
var_1d1 = initial_investment  cutoff1
var_1d1
#output
#22347.7792230231
Here we are saying with 95% confidence that our portfolio of 1M USD will not exceed losses greater than 22.3k USD over a one day period.
Value at risk over nday time period
What if we wanted to calculate this over a larger window of time? Below we can easily do that by just taking our 1 day VaR and multiplying it by the square root of the time period (this is due to the fact that the standard deviation of stock returns tends to increase with the square root of time).
# Calculate n Day VaR
var_array = []
num_days = int(15)
for x in range(1, num_days+1):
var_array.append(np.round(var_1d1 * np.sqrt(x),2))
print(str(x) + " day VaR @ 95% confidence: " + str(np.round(var_1d1 * np.sqrt(x),2)))
# Build plot
plt.xlabel("Day #")
plt.ylabel("Max portfolio loss (USD)")
plt.title("Max portfolio loss (VaR) over 15day period")
plt.plot(var_array, "r")
1 day VaR @ 95% confidence: 20695.24
2 day VaR @ 95% confidence: 29267.49
3 day VaR @ 95% confidence: 35845.21
4 day VaR @ 95% confidence: 41390.49
5 day VaR @ 95% confidence: 46275.97
6 day VaR @ 95% confidence: 50692.79
7 day VaR @ 95% confidence: 54754.47
8 day VaR @ 95% confidence: 58534.99
9 day VaR @ 95% confidence: 62085.73
10 day VaR @ 95% confidence: 65444.11
11 day VaR @ 95% confidence: 68638.36
12 day VaR @ 95% confidence: 71690.43
13 day VaR @ 95% confidence: 74617.76
14 day VaR @ 95% confidence: 77434.51
15 day VaR @ 95% confidence: 80152.33
(Extra) Checking distributions of our equities against normal distribution
As mentioned in the calculation section, we are assuming that the returns of the equities in our portfolio are normally distributed when calculating VaR. Of course, we can't predict that moving forward, but we can at least check how the historical returns have been distributed to help us assess whether VaR is suitable to use for our portfolio.
import matplotlib.mlab as mlab
import matplotlib.pyplot as plt
# Repeat for each equity in portfolio
returns['AAPL'].hist(bins=40, normed=True,histtype="stepfilled",alpha=0.5)
x = np.linspace(port_mean  3*port_stdev, port_mean+3*port_stdev,100)
plt.plot(x, scipy.stats.norm.pdf(x, port_mean, port_stdev), "r")
plt.title("AAPL returns (binned) vs. normal distribution")
plt.show()
AAPl returns vs. normal distribution
FB returns vs. normal distribution
C returns vs. normal distribution
DIS returns vs. normal distribution
From the above we can see the returns have all been fairly normally distributed for our chosen stocks since 2018.
//You can view the code in Google Colab here  to run it just make a copy of the notebook.