VOLATILITY. Finance is risk/return trade-off.

Similar documents
VOLATILITY. Time Varying Volatility

Downside Risk: Implications for Financial Management Robert Engle NYU Stern School of Business Carlos III, May 24,2004

Financial Econometrics Jeffrey R. Russell Midterm 2014

Donald Trump's Random Walk Up Wall Street

Chapter 4 Level of Volatility in the Indian Stock Market

Lecture 5a: ARCH Models

Brief Sketch of Solutions: Tutorial 2. 2) graphs. 3) unit root tests

Financial Econometrics: Problem Set # 3 Solutions

WHY IS FINANCIAL MARKET VOLATILITY SO HIGH? Robert Engle Stern School of Business BRIDGES, Dialogues Toward a Culture of Peace

Financial Econometrics Jeffrey R. Russell. Midterm 2014 Suggested Solutions. TA: B. B. Deng

Amath 546/Econ 589 Univariate GARCH Models: Advanced Topics

Appendix. Table A.1 (Part A) The Author(s) 2015 G. Chakrabarti and C. Sen, Green Investing, SpringerBriefs in Finance, DOI /

Research Article The Volatility of the Index of Shanghai Stock Market Research Based on ARCH and Its Extended Forms

Appendixes Appendix 1 Data of Dependent Variables and Independent Variables Period

ARCH and GARCH models

Monetary Economics Measuring Asset Returns. Gerald P. Dwyer Fall 2015

MODELING EXCHANGE RATE VOLATILITY OF UZBEK SUM BY USING ARCH FAMILY MODELS

Brief Sketch of Solutions: Tutorial 1. 2) descriptive statistics and correlogram. Series: LGCSI Sample 12/31/ /11/2009 Observations 2596

Volatility Analysis of Nepalese Stock Market

Lecture 5: Univariate Volatility

Analysis of the Influence of the Annualized Rate of Rentability on the Unit Value of the Net Assets of the Private Administered Pension Fund NN

Graduate School of Business, University of Chicago Business 41202, Spring Quarter 2007, Mr. Ruey S. Tsay. Midterm

LAMPIRAN. Null Hypothesis: LO has a unit root Exogenous: Constant Lag Length: 1 (Automatic based on SIC, MAXLAG=13)

Linda Allen, Jacob Boudoukh and Anthony Saunders, Understanding Market, Credit and Operational Risk: The Value at Risk Approach

Modelling Stock Market Return Volatility: Evidence from India

Financial Econometrics Notes. Kevin Sheppard University of Oxford

Financial Econometrics

1. A test of the theory is the regression, since no arbitrage implies, Under the null: a = 0, b =1, and the error e or u is unpredictable.

Amath 546/Econ 589 Univariate GARCH Models

Conditional Heteroscedasticity

Lecture 1: The Econometrics of Financial Returns

Model Construction & Forecast Based Portfolio Allocation:

Volatility Clustering of Fine Wine Prices assuming Different Distributions

Lecture 6: Non Normal Distributions

Notes on the Treasury Yield Curve Forecasts. October Kara Naccarelli

GARCH Models. Instructor: G. William Schwert

Statistics and Finance

Business Statistics 41000: Probability 3

An Empirical Research on Chinese Stock Market Volatility Based. on Garch

Modeling the volatility of FTSE All Share Index Returns

Monetary Economics Portfolios Risk and Returns Diversification and Risk Factors Gerald P. Dwyer Fall 2015

Graduate School of Business, University of Chicago Business 41202, Spring Quarter 2007, Mr. Ruey S. Tsay. Final Exam

Financial Econometrics (FinMetrics04) Time-series Statistics Concepts Exploratory Data Analysis Testing for Normality Empirical VaR

Econometric Models for the Analysis of Financial Portfolios

LAMPIRAN PERHITUNGAN EVIEWS

Variance clustering. Two motivations, volatility clustering, and implied volatility

BEcon Program, Faculty of Economics, Chulalongkorn University Page 1/7

Per Capita Housing Starts: Forecasting and the Effects of Interest Rate

Lampiran 1 : Grafik Data HIV Asli

MEASURING PORTFOLIO RISKS USING CONDITIONAL COPULA-AR-GARCH MODEL

Financial Time Series Analysis (FTSA)

INFLUENCE OF CONTRIBUTION RATE DYNAMICS ON THE PENSION PILLAR II ON THE

The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2017, Mr. Ruey S. Tsay. Solutions to Final Exam

Assicurazioni Generali: An Option Pricing Case with NAGARCH

V Time Varying Covariance and Correlation. Covariances and Correlations

The Great Moderation Flattens Fat Tails: Disappearing Leptokurtosis

Forecasting Stock Index Futures Price Volatility: Linear vs. Nonlinear Models

9. Appendixes. Page 73 of 95

Implied Volatility v/s Realized Volatility: A Forecasting Dimension

I. Return Calculations (20 pts, 4 points each)

INFORMATION EFFICIENCY HYPOTHESIS THE FINANCIAL VOLATILITY IN THE CZECH REPUBLIC CASE

Booth School of Business, University of Chicago Business 41202, Spring Quarter 2014, Mr. Ruey S. Tsay. Solutions to Midterm

ANALYSIS OF CORRELATION BETWEEN THE EXPENSES OF SOCIAL PROTECTION AND THE ANTICIPATED OLD AGE PENSION

Booth School of Business, University of Chicago Business 41202, Spring Quarter 2012, Mr. Ruey S. Tsay. Solutions to Midterm

Time series: Variance modelling

Lecture Note of Bus 41202, Spring 2008: More Volatility Models. Mr. Ruey Tsay

FINANCIAL ECONOMETRICS PROF. MASSIMO GUIDOLIN

1 Volatility Definition and Estimation

Market Risk Analysis Volume I

A market risk model for asymmetric distributed series of return

THE INFORMATION CONTENT OF IMPLIED VOLATILITY IN AGRICULTURAL COMMODITY MARKETS. Pierre Giot 1

Lecture 6: Univariate Volatility

Research on the GARCH model of the Shanghai Securities Composite Index

Structural GARCH: The Volatility-Leverage Connection

Modeling Exchange Rate Volatility using APARCH Models

Financial Times Series. Lecture 8

Graduate School of Business, University of Chicago Business 41202, Spring Quarter 2007, Mr. Ruey S. Tsay. Solutions to Final Exam

IMPACT OF MACROECONOMIC VARIABLE ON STOCK MARKET RETURN AND ITS VOLATILITY

Statistical Models and Methods for Financial Markets

Booth School of Business, University of Chicago Business 41202, Spring Quarter 2016, Mr. Ruey S. Tsay. Solutions to Midterm

ANALYSIS OF THE RETURNS AND VOLATILITY OF THE ENVIRONMENTAL STOCK LEADERS

Time series analysis on return of spot gold price

Financial Time Series Analysis: Part II

starting on 5/1/1953 up until 2/1/2017.

Booth School of Business, University of Chicago Business 41202, Spring Quarter 2012, Mr. Ruey S. Tsay. Midterm

Methods for A Time Series Approach to Estimating Excess Mortality Rates in Puerto Rico, Post Maria 1 Menzie Chinn 2 August 10, 2018 Procedure:

Market Risk Analysis Volume II. Practical Financial Econometrics

Jaime Frade Dr. Niu Interest rate modeling

Fin285a:Computer Simulations and Risk Assessment Section 7.1 Modeling Volatility: basic models Daníelson, ,

The Effects of Oil Price Volatility on Some Macroeconomic Variables in Nigeria: Application of Garch and Var Models

Web Appendix. Are the effects of monetary policy shocks big or small? Olivier Coibion

Booth School of Business, University of Chicago Business 41202, Spring Quarter 2013, Mr. Ruey S. Tsay. Midterm

Lampiran 1. Data Penelitian

Properties of financail time series GARCH(p,q) models Risk premium and ARCH-M models Leverage effects and asymmetric GARCH models.

The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2009, Mr. Ruey S. Tsay. Solutions to Final Exam

BUSM 411: Derivatives and Fixed Income

A Note on the Oil Price Trend and GARCH Shocks

MODELING ROMANIAN EXCHANGE RATE EVOLUTION WITH GARCH, TGARCH, GARCH- IN MEAN MODELS

MODELING VOLATILITY OF BSE SECTORAL INDICES

CHAPTER 5 MARKET LEVEL INDUSTRY LEVEL AND FIRM LEVEL VOLATILITY

Transcription:

VOLATILITY RISK Finance is risk/return trade-off. Volatility is risk. Advance knowledge of risks allows us to avoid them. But what would we have to do to avoid them altogether??? Imagine! How much should I get paid to take a given risk? Which risks are not worth taking? 1

DEFINE VARIANCE E(x) is the expected value of a random variable x x-e(x) is the unexpected part of x V(x) is the variance of x V(x)=E(x-E(x)) 2 Standard Deviation is the square root of Variance Time Varying Volatility VOLATILITY IS THE STANDARD DEVIATION OF RETURNS - the unpredictable part of asset prices. The unpredictable part is the surprise part, or the error in the prediction. Consider the AR(1) model: rt 0 1rt 1t t t terpast t Expected part given the past Unexepected part We define the conditional variance as the variance of the surprise part or t. 2

More generally for any model we can write it as rt t t Expected part Unexepected part The variance is then defined as the variance of the unexpected part or the deviation of the return from its prediction. 2 2 2 t E t past E rt t past Before we modeled returns where the variance of the surprise or error was constant. If the errors are Normal then the conditional distribution is given by: 2 t t, f r past N Now we have acknowledged that the variance may be time varying as well as 2 the mean! t t, t f r past N 3

So the conditional variance of a return is the variance of the surprise. IF the mean t =0 then we have rt 0 t t In other words, the surprise is the return. For simplicity, we will assume the mean is zero for now. We do this for convenience. IF the mean were not zero we would have to first subtract t from r t and then model the variance of the demeaned returns. Simplification for now If markets are nearly efficient (in Fama sense) then returns are nearly unpredictable. If the drift is small then the conditional variance is given by: 2 1 1 2 2 t E rt t Ft E rt Ft 4

ANNUALIZED VOLATILITY Returns are essentially unpredictable if the efficient market hypothesis is a good approximation. Hence the variance over multiple days is the sum of the variances over each of the days (no covariance terms). If the variance were the same every day of the year then the annualized variance is given by n*(daily variance) where n is the number of trading days. In the US n~~252. Annualized Volatility is sqr(252)*daily volatility or sqr(12)*monthly volatility Let s first take a look at historical volatility. Then we will consider models that capture the time varying features observed in the data. 5

HISTORY OF THE US EQUITY MARKET: S&P500 S&P500 daily returns 1990 to present.12 R.08.04.00 -.04 -.08 -.12 90 92 94 96 98 00 02 04 06 08 10 12 14 16 6

.15.10.05 3,000 2,500 2,000 1,500.00 -.05 -.10 1,000 500 0 90 92 94 96 98 00 02 04 06 08 10 12 14 16 ADJ_CLOSE R Let s look more closely at the characteristics of financial returns. ALMOST UNPREDICTABLE EFFICIENT MARKET HYPOTHESIS SURPRISINGLY LARGE NUMBER OF EXTREMES FAT TAIL DISTRIBUTIONS PERIODS OF HIGH AND LOW VOLATILITIES VOLATILITY CLUSTERING 7

CHECK IT OUT! HOW TO CHECK FOR EXCESSIVE EXTREMES HOW TO CHECK FOR VOLATILITY CLUSTERING? HISTORICAL VOLATILITY Estimate the standard deviation of a random variable T 2 ˆ 252 rj / jtk What assumptions do we need? Choose K small so that the variance is constant Choose K large to make the estimate as accurate as possible Funny boxcars and shadow volatility movements!! K 8

EXPONENTIAL SMOOTHING Volatility Estimator used by RISKMETRICS Updating 1 r 2 2 2 t t1 t1 AN EXAMPLE WEAKNESSES How to choose lambda No mean reversion II ARCH/GARCH MODELS GARCH VOLATILITY FORECASTING WITH GARCH ESTIMATING AND TESTING GARCH MANY MODELS 9

Decisions in Finance often require measures of the volatility Risk can be thought of as the variance. Expected returns of an asset are related to their variance. Optimal portfolio selection Markowitz (1952) and Tobin (1958) (Nobel Prize). The value of an option BLACK-SCHOLES AND MERTON Options can be used as insurance policies. For a fee we can eliminate financial risk for a period. What is the right fee? Black and Scholes(1972) and Merton(1973) (Nobel Prize) developed an option pricing formula from a dynamic hedging argument. 10

IMPLEMENTING THESE MODELS PRACTITIONERS REQUIRED ESTIMATES OF VARIANCES AND COVARIANCES EQUIVALENTLY WE SAY VOLATILITIES (Standard deviation which is the square root of the variance) AND CORRELATIONS (which is the covariance divided by the product of the standard deviations) ESTIMATES DIFFER FOR DIFFERENT TIME PERIODS Volatility is apparently varying over time What is the volatility NOW! What is it likely to be in the future? How can we forecast something we never observe? 11

Along comes ARCH 2003 Nobel prize to Engle. WHAT IS ARCH? Autoregressive Conditional Heteroskedasticity THE SIMPLEST PROBLEM WHAT IS VOLATILITY NOW? One answer is the standard deviation over the last 5 years But this will include lots of old information that may not be relevant for short term forecasting Another answer is the standard deviation over the last 5 days But this will be highly variable because there is so little information 12

THE ARCH ANSWER Use a weighted average of the volatility over a long period with higher weights on the recent past and small but non-zero weights on the distant past. Choose these weights by looking at the past data; what forecasting model would have been best historically? This is a statistical estimation problem. The ARCH Model The ARCH model of Engle(1982) is a family of specifications for the conditional variance. The q th order ARCH or ARCH(q) model is q ht jrt j1 Where in the ARCH notation 2 2 h E r F is the conditional variance. t t t t1 2 j 13

FROM THE SIMPLE ARCH GREW: GENERALIZED ARCH (Bollerslev) a most important extension Tomorrow s variance is predicted to be a weighted average of the Long run average variance Today s variance forecast The news (today s squared return) GARCH h r h 2 t t1 t1 Generalization of Exponential Smoothing Generalization of ARCH Generalization of constant volatility 14

UPDATING Suppose the model is: h t =.00001+.05r 2 t-1 +.9h t-1 And today annualized volatility is 20% and the market return is -3%, what is my estimate of tomorrow s volatility from this model?.000198=.00001+.05(-.03 2 )+.9*(.2 2 /252) Or 22.3% annualized. REPEAT STARTING AT T=1 IF WE KNOW THE PARAMETERS AND SOME STARTING VALUE FOR h 1, WE CAN CALCULATE THE ENTIRE HISTORY OF VOLATILITY FORECASTS OFTEN WE USE A SAMPLE VARIANCE FOR h 1. 15

GARCH(p,q) The Generalized ARCH model of Bollerslev(1986) is an ARMA version of this model. GARCH(p,q) is q p 2 t j tj j tj j1 j1 h r h Asymmetric Volatility Often negative shocks have a bigger effect on volatility than positive shocks Nelson(1987) introduced the EGARCH model to incorporate this effect. I will use a Threshold GARCH or TARCH q q p 2 2 h r r I 0 h rt j t j tj j tj j tj j1 j1 j1, 16

NEW ARCH MODELS GJR-GARCH TARCH STARCH AARCH NARCH MARCH SWARCH SNPARCH APARCH TAYLOR-SCHWERT FIGARCH FIEGARCH Component Asymmetric Component SQGARCH CESGARCH Student t GED SPARCH Autoregressive Conditional Density Autoregressive Conditional Skewness ROLLING WINDOW VOLATILITIES NUMBER OF DAYS=5,260,1300.6.4.2.0 64 66 68 70 72 74 76 78 80 82 84 86 88 90 92 94 96 98 00 02 V5 V260 V1300 17

1.4 ARCH/GARCH VOLATILITIES 1.2 1.0 0.8 0.6 0.4 0.2 0.0 65 70 75 80 85 90 95 00 GARCHVOL CONFIDENCE INTERVALS.10.05.00 -.05 -.10 1990 1992 1994 1996 1998 2000 2002 3*GARCHSTD SPRETURNS -3*GARCHSTD 18

UNCONDITIONAL, OR LONG RUN, OR AVERAGE VARIANCE WHAT IS E(r 2 )? Eh 2 2 E r 2 E h E r E h t t1 t1 Eh 2 2 2 1 UNCONDITIONAL VARIANCE Suppose the model is: h t =.00001+.05r t-1 2 +.9h t-1 What is the unconditional annualized volatility?.0002 or 22.44% annualized. 19

The GARCH Model Again r h r h t t 2 t t 1 t 1 2 2 1 rt 1 ht 1 The variance of r t is a weighted average of three components a constant or unconditional variance yesterday s forecast yesterday s news ERRORS r t t t WE MODEL SO THE CONDITIONAL VARIANCE IS THE VARIANCE OF THE s. IF t =0 then 2 2 E rt t F t1 E t F t1 2 2 E rt F t1 E t F t1 20

We saw that the errors have fat tails and may be skewed. This is inconsistent with the assumption that errors are normally distributed. We now introduce a way of thinking about GARCH models where we allow the error to have a different distribution than the normal. The key is we still need the error to be mean zero and we need h t to still be the conditional variance how do we do that? An interesting alternative way to think about the GARCH model WE CAN WRITE t hz t t WHERE z t is iid (not necessarily normal) with mean zero and variance one. 2 NOW, 2 2 1 E t Ft E htzt past hte zt h t h t is known given past and z t is iid with variance 1. 21

THE ERRORS z t MUST HAVE VARIANCE 1 THEY COULD BE NORMAL THEY MIGHT HAVE FATTER TAILS LIKE THE STUDENT T OR GENERALIZED EXPONENTIAL DIAGNOSTIC CHECKING r h z t t t so r t h t standardized series When we divide r t by its standard deviation, it should have a variance of 1 for all t. We can check this by creating the standardized series and looking for volatility clustering in them. z t 22

Intuition 2 Suppose r t has variance. 2 Then rt 1 Var Var 2 rt 2 1 So if we divide a return by it s standard deviation, the new return should have variance 1. Now suppose that returns have different variances. Let Var r 2 t t So the r t Var 1 t If the GARCH model is the right model then the GARCH series should be the right variances for the r t s and and h t should be the right standard deviation so Var r t 1 h t 23

This suggests that we can standardize each return by dividing by its conditional standard deviation and the resulting series should no longer have time varying volatility. The variance should be constant. We have a test for time varying volatility so we can check and see if the GARCH model is right by testing to see if the standardized returns no longer have time varying volatility. If we divide by the wrong standard deviation then the resulting series will not have constant variance, but rather will have time varying variance. Time varying volatility is revealed by volatility clustering. These are measured by the Ljung Box statistic on squared returns. A failure to reject the null of the standardized series constant volatility suggests the GARCH model is a good one. 24

Time varying volatility is revealed by volatility clusters These are measured by the Ljung Box statistic on squared returns The standardized returns zt rt / ht no longer should show significant volatility clustering FORECASTING VOLATILITY 25

FORECASTING FOR GARCH ONE STEP AHEAD FORECAST TWO STEP FORECAST h r h 2 t2 t1 t1 E h F h t2 t t1 2 2 t2 t t1 E h F h MEAN REVERTING VOLATILITY Forecasts converge to the same value no matter what the current volatility 2 k 1 2 tk t1 E h past h E h 2 tk past if + <1 LITTLE UPDATING FOR LONG HORIZON VOLATILITY 26

Monotonic Term Structure of Volatility FORECAST PERIOD DOW JONES SINCE 1990 Dependent Variable: DJRET Method: ML - ARCH (Marquardt) - Normal distribution Date: 01/10/08 Time: 13:42 Sample: 1/02/1990 1/04/2008 Included observations: 4541 Convergence achieved after 15 iterations GARCH = C(2) + C(3)*RESID(-1)^2 + C(4)*GARCH(-1) Coefficient Std. Error z-statistic Prob. C 0.000527 0.000119 4.414772 0.0000 Variance Equation C 1.00E-06 1.37E-07 7.290125 0.0000 RESID(-1)^20.064459 0.004082 15.79053 0.0000 GARCH(-1) 0.925645 0.005025 184.2160 0.0000 R-squared -0.000371 Mean dependent var 0.000338 Adjusted R-squared -0.001032 S.D. dependent var 0.009795 S.E. of regression 0.009800 Akaike info criterion -6.640830 Sum squared resid 0.435759 Schwarz criterion -6.635174 Log likelihood 15082.00 Hannan-Quinn criter. -6.638838 Durbin-Watson stat 2.001439 27

VOLGARCH.45.40.35.30.25.20.15.10.05 90 92 94 96 98 00 02 04 06 FORECAST.184 VOLFORECAST.180.176.172.168.164.160 08M01 08M02 08M03 08M04 08M05 08M06 08M07 08M08 08M09 08M10 08M11 08M12 28

.028.024.020.016.012.008.004 1998 1999 2000 2001 2002 2003 2004 2005 2006 DJSD DJSD1 DJSD2 DJSD3 DJSD4 DJSD5 DJSD0 DJSDEND MULTI-STEP FORECASTS tk t tk t tk1 t 2 k 1 2 ht 1 2 2 2 E h F E r F E h F Forecast over long horizons is the sum of the forecasts from t+1 to t+k. V R V r E h F K K K t t t t j t j t j1 j1 29

EXOGENOUS VARIABLES IN A GARCH MODEL h t Include predetermined variables into the variance equation Easy to estimate and forecast one step Multi-step forecasting is difficult Timing may not be right 2 rt 1 ht 1 z t 1 EXAMPLES Non-linear effects Deterministic Effects News from other markets Heat waves vs. Meteor Showers Other assets Implied Volatilities Index volatility MacroVariables or Events 30

PARAMETER ESTIMATION Historical data reveals when volatilities were large Pick parameters to match the historical volatility episodes Maximum Likelihood is the estimation procedure used by EVIEWS. PLAUSIBLE ANSWERS WE EXPECT ALL THREE PARAMETERS OF A GARCH(1,1) TO BE POSITIVE. WE EXPECT THE SUM OF ALPHA AND BETA TO BE VERY CLOSE TO ONE BUT LESS THAN ONE. WE EXPECT THE UNCONDITIONAL VARIANCE TO BE CLOSE TO THE DATA VARIANCE. 31

DID THE ESTIMATION ALGORITHM CONVERGE? Generally the software will reliably find the maximum of the likelihood function and will report it. Sometimes it does not. You may get silly values. What then? Check with other starting values Check with other iterations Scale the data so the numbers are not so small Often the problem is the data. Look for outliers or peculiar features. Use longer data set NORMALITY THIS ESTIMATION METHOD IS OPTIMAL IF THE ERRORS ARE NORMAL AND IF THE SAMPLE IS LARGE AND THE MODEL IS CORRECT. IT IS STILL GOOD WITHOUT NORMALITY BUT OTHER ESTIMATORS COULD BE BETTER SUCH AS STUDENT-T. 32

COMPARE MODELS MODELS WHICH ACHIEVE THE HIGHEST VALUE OF THE LOG LIKELIHOOD ARE PREFERRED. IF THEY HAVE DIFFERENT NUMBERS OF PARAMETERS THIS IS NOT A FAIR COMPARISON. USE AIC OR BIC (SCHWARZ) INSTEAD. THE SMALLEST VALUE IS BEST. WHAT IS THE BEST MODEL? The most reliable and robust is GARCH(1,1) A student-t error assumption gives better estimates of tails. For equities asymmetry is almost always important. See asymmetric volatility section. For long term forecasts, a component model is often needed. Even better is a model which incorporates economic variables 33

III VALUE AT RISK ESTIMATION VALUE AT RISK GARCH ASYMMETRIC VOLATILITY DOWNSIDE RISK BUBBLES AND CRASHES DOWNSIDE RISK THE RISK OF A PORTFOLIO IS THAT ITS VALUE WILL DECLINE, NOT THAT IT WILL INCREASE HENCE DOWNSIDE RISK IS NATURAL. MANY THEORIES AND MODELS ASSUME SYMMETRY: c.f. MARKOWITZ, TOBIN, SHARP AND VOLATILITY BASED RISK MANAGEMENT SYSTEMS. DO WE MISS ANYTHING IMPORTANT? 34

MEASURING DOWNSIDE RISK Many measures have been proposed. Let r be the one period continuously compounded return with distribution f(r) and mean zero. Let x be a threshold. Skewness= E r / 3 2 Er 3/2 Probability of loss = P r x, Expected loss = E r r x x is the Value at risk if P rx Value at Risk For a portfolio the future value is uncertain VaR is a number of $ that you can be 99% sure, is worse than what will happen. It is the 99% of the loss distribution (or the 1% quantile of the gain distribution) Simple idea, but how to calculate this? 35

PREDICTIVE DISTRIBUTION OF PORTFOLIO GAINS 1% $ GAINS ON PORTFOLIO CONCERNS This single number (a quantile) is used to represent a full distribution. It can be misleading. It fails to satisfy some basic principles in rather pathological cases. It assumes that you do not change your portfolio over the next day Bad assumption for day traders and proprietary trading desks Does not recognize role of dynamic hedging But it can be calculated on a higher frequency too. 36

PUZZLER $1,000,000 portfolio of SP500. Find one day 99% VaR. First, just take a guess. HISTORICAL VaR If History repeats, look at worst outcomes in the past For example, S&P500 over the last year. On a $1,000,000 portfolio, the 99% VaR is? 37

HISTOGRAM OF S&P500 GAINS 1% quantile = -0.0445 60 50 40 30 20 10 0-0.06-0.04-0.02 0.00 0.02 0.04 Series: RET Sample 11/29/2010 10/28/2011 Observations 233 Mean 0.000427 Median 0.000947 Maximum 0.047407 Minimum -0.066634 Std. Dev. 0.013729 Skewness -0.639278 Kurtosis 6.977393 Jarque-Bera 169.4528 Probability 0.000000 HISTORICAL S&P500 VaR If I use 2 years of data, it is $34,411 With 3 years, it is $44,594 And with 30 years it is $30,108 Which is more accurate? 38

Dow Jones 99% VaR USING ONE YEAR HISTORICAL QUANTILES 120000 100000 80000 60000 40000 20000 0 30 35 40 45 50 55 60 65 70 75 80 85 90 95 00 05 VAR 45000 40000 35000 30000 25000 20000 15000 2000 2001 2002 2003 2004 VAR 39

HOWEVER The extreme events of more than one year ago are not considered All outcomes of the last year are considered to be equally likely How should you pick the window? Look at a picture of VaR over time it has a boxcar shape implying that risk goes down exactly one year after a big market decline. VOLATILITY BASED VaR With a good volatility forecast, predict the standard deviation of tomorrows return. Assume a Normal Distribution. Then VaR is 2.33* t But what do we use for the volatility? GARCH forecasts! Other volatility estimates? 40

GARCH MODEL FOR DJ USE FOR EXAMPLE DATA FOR 10 YEARS (95-05) FORECAST OUT OF SAMPLE AND RECORD THE DAILY STANDARD DEVIATION MULTIPLY BY 2.33 WE GET RESULTS GARCH MODEL C 1.30E-06 2.90E-07 4.474868 0.0000 RESID(-1)^2 0.085773 0.006714 12.77441 0.0000 GARCH(-1) 0.907123 0.007644 118.6731 0.0000 DATE RETURN DAILY SD VaR 2005-01-07-0.001783 NA NA 2005-01-10 NA 0.006166 14367.77 41

VOLATILITY BASED VaR WITH STUDENT-T ERRORS Assume that: Because Then let z * ~ Student t, t * * z Vz v/ v2, V 1 v/ v 2 r h z h z * / v / v 2, t t t t t And estimate volatility and the shape of the error distribution jointly. In EViews =@qtdist(.01,v)/sqr[v/(v-2)] STUDENT-T RESULTS GARCH WITH STUDENT T ERRORS C 1.01E-06 3.20E-07 3.147330 0.0016 RESID(-1)^2 0.063884 0.009483 6.736429 0.0000 GARCH(-1) 0.929008 0.009920 93.65090 0.0000 T-DIST. DOF 8.839721 1.240570 7.125529 0.0000.01 QUANTILE OF UNIT STUDENT T DISTRIBUTION(8.8DF) IS 2.49 DATE RETURN DAILY SD VaR 2005-01-07-0.001783 NA NA 2005-01-10 NA 0.006260 155874 42

VOLATILITY BASED VaR WITHOUT NORMALITY or t What is the right multiplier for the true distribution? Maybe neither the normal nor the student t are correct! If: r h z, z ~ iid... t t t t Then 1% quantile of the standardized residuals should be used. This is the bootstrap estimator or Hull and White s volatility adjustment. HISTOGRAM OF STANDARDIZED RESIDUALS.01 QUANTILE = -2.55 300 250 200 150 100 50 0-5.0-2.5 0.0 2.5 Series: GARCHRESID Sample 1/09/1995 1/20/2005 Observations 2520 Mean -0.040388 Median -0.021369 Maximum 3.036714 Minimum -6.058763 Std. Dev. 0.999691 Skewness -0.383229 Kurtosis 4.474130 Jarque-Bera 289.8543 Probability 0.000000 43

BOOTSTRAP VaR DATE RETURN DAILY SD VaR 2005-01-07-0.001783 NA NA 2005-01-10 NA 0.006166 15724.38 OVERVIEW AND REVIEW HISTORICAL QUANTILES RESULT IS SENSITIVE TO SAMPLE INCLUDED VOLATILITY BASED RESULT IS SENSITIVE TO THE ERROR DISTRIBUTION NORMAL UNDERSTATES EXTREME RISK T AND BOOTSTRAP ARE BETTER. RESULTS ARE NOT SENSITIVE TO THE SAMPLE INCLUDED 44

ASYMMETRIC VOLATILTY MODELS ASYMMETRIC VOLATILITY Positive and negative returns might have different weights. For example: h r I r I h 2 2 t 1 t1 r 0 2 t1 r 0 t1 t1 t1 h r r I h 2 2 t t1 t1 r 0 t1 We typically find for equities that or equivalently >0 2 1 t1 45

NEWS IMPACT CURVE TOMORROWS VARIANCE TODAY S NEWS = RETURNS Other Asymmetric Models EGARCH: NELSON(1989) rt 1 log( ht ) log( ht 1) ht 1 r h t1 t1 NGARCH: ENGLE(1990) 2 h r h t ( ) t1 t1 46

PARTIALLY NON-PARAMETRIC ENGLE AND NG(1993) VOLATILITY NEWS WHERE DOES ASYMMETRIC VOLATILITY COME FROM? LEVERAGE - As equity prices fall the leverage of a firm increases so that the next shock has higher volatility on stock prices. This effect is usually too small to explain what we see. RISK AVERSION News of a future volatility event will lead to stock sales and price declines. Subsequently, the volatility event will occur. Since events are clustered, any news event will predict higher volatility in the future. This effect is more plausible on broad market indices since these have systematic risk. 47

BACK TO VALUE AT RISK FIND QUANTILE OF FUTURE RETURNS One day in advance Many days in advance REGULATORY STANDARD IS 10 DAY 1% VaR. MULTI-DAY VaR What is the risk over 10 days if you do no more trading? Clearly this is greater than for one day. Now we need the distribution of multiday returns. 48

10 Day VaR If volatility were constant, then the multi-day volatility would simply require multiplying by the square root of the days. With normality and constant variance this becomes 7.36 or sqr(10)*2.33 VaR is 7.36 * sigma What is sigma? MULTI-DAY HORIZONS Because volatility is dynamic and asymmetric, the lower tail is more extreme and the VaR should be greater. 49

TWO PERIOD RETURNS Two period return is the sum of two one period continuously compounded returns Look at binomial tree version Asymmetry gives negative skewness Low variance High variance MULTIPLIER FOR 10 DAYS For a 10 day 99% value at risk, conventional practice multiplies the daily standard deviation by 7.36 For the same multiplier with asymmetric GARCH it is simulated from the example to be 7.88 Bootstrapping from the residuals the multiplier becomes 8.52 50

CALCULATION BY SIMULATION EVALUATE ANY MEASURE BY REPEATEDLY SIMULATING FROM THE ONE PERIOD CONDITIONAL DISTRIBUTION: METHOD: f t Draw r t+1 Update density and draw observation t+2 Continue until T returns are computed. Repeat many times Compute measure of downside risk rt 1 ESTIMATE TARCH MODEL VARIABLE COEF STERR T-STAT P-VALUE C 1.68E-06 2.58E-07 6.519983 0.0000 RESID(-1)^2 0.005405 0.008963 0.603066 0.5465 RESID(-1)^2*(RESID(-1)<0) 0.123800 0.010668 11.60488 0.0000 GARCH(-1) 0.918895 0.008211 111.9126 0.0000 DATE CONDITIONAL VARIANCE 2005-01-07 0.006835 2005-01-10 0.006726 51

TARCH STANDARD DEVIATIONS.035.030.025.020.015.010.005.000 95 96 97 98 99 00 01 02 03 04 DJSDGARCH DJSDTARCH TARCH STANDARD DEVIATIONS.035.030.025 DJSDTARCH.020.015.010.005.000.000.005.010.015.020.025.030 DJSDGARCH 52

DOWNSIDE RISK DOWNSIDE RISK With Asymmetric Volatility, the multiperiod returns are asymmetric with a longer left tail. For long horizons, the central limit theorem will reduce this effect and returns will be approximately normal. This is observed in data too. 53

1 DAY RETURNS ON D.J. 600 500 400 300 200 100 0-0.075-0.050-0.025 0.000 0.025 0.050 Series: DJRET Sample 1/03/1995 1/20/2005 Observations 2524 Mean 0.000403 Median 0.000557 Maximum 0.061547 Minimum -0.074549 Std. Dev. 0.011179 Skewness -0.250246 Kurtosis 7.024450 Jarque-Bera 1729.644 Probability 0.000000 10 DAY RETURNS ON D.J. 360 320 280 240 200 160 120 80 40 0-0.1 0.0 0.1 Series: @MOVSUM(DJRET,10) Sample 1/03/1995 1/20/2005 Observations 2524 Mean 0.004091 Median 0.006289 Maximum 0.153400 Minimum -0.189056 Std. Dev. 0.033474 Skewness -0.619562 Kurtosis 6.014598 Jarque-Bera 1117.209 Probability 0.000000 54

0.2 0 SKEWNESS OF MULTIPERIOD RETURNS 0 25 50 75 100 125 150 175 200 225-0.2-0.4-0.6 SKEW_ALL SKEW_TRIM SKEW_PRE SKEW_POST -0.8-1 EVIDENCE FROM DERIVATIVES THE HIGH PRICE OF OUT-OF-THE- MONEY EQUITY PUT OPTIONS IS WELL DOCUMENTED THIS IMPLIES SKEWNESS IN THE RISK NEUTRAL DISTRIBUTION MUCH OF THIS IS PROBABLY DUE TO SKEWNESS IN THE EMPIRICAL DISTRIBUTION OF RETURNS. DATA MATCHES EVIDENCE THAT THE OPTION SKEW IS ONLY POST 1987. 55

MATCHING THE STYLIZED FACTS ESTIMATE DAILY MODEL SIMULATE 250 CUMULATIVE RETURNS 10,000 TIMES WITH SEVERAL DATA GENERATING PROCESSES CALCULATE SKEWNESS AT EACH HORIZON SKEWS FOR SYMMETRIC AND ASYMMETRIC MODELS 0.2 0.0-0.2-0.4-0.6-0.8-1.0-1.2-1.4 25 50 75 100 125 150 175 200 225 250 SKEW_EX SKEW_BOOT_EX SKEW_EXS SKEW_BOOT_EXS 56

IMPLICATIONS Multi-period empirical returns are more skewed than one period returns (omitting 1987 crash) Asymmetric volatility is needed to explain this. Skewness has increased since 1987, particularly for longer horizons. Simulated skewness is noisy because higher moments do not exist when the persistence is so close to one. Presumably this is true for the data too. Many other asymmetric models could be compared on this basis. SPIKES AND CRASHES RARE BUT EXTREME MOVES DOWN OR UP WHAT ARE THE CHARACTERISTICS OF SPIKES AND CRASHES? CAN WE PREDICT THEM? 57

.08.06.04.02.00 -.02 -.04 -.06 LOOK AT 3%MARKET DROPS 10 8 6 4 2 -.08 0 1990 1992 1994 1996 1998 2000 2002 2004 2006 DJRET DOWN3 CORRELATIONS Find the correlations between the indicator variable for a 3% crash and several independent variables 22 day average of down events 22 day average of up events Up within 3 days Exponential smoothed.06 Exponential smoothed.02 Garch volatility estimated from 1928-1990 and forecast thereafter 58

CORRELATIONS 0.12 0.1 0.08 0.06 0.04 0.02 0-0.02 DOWN3(1) 3-D Column 1 DOWNMA22 UP3 UPMA22 VOL06 VOL02 VOL22 VOL252 TARCH is estimated from 1928-1990 and forecast to 2006. Downma22 is a moving average of past down days, and similar is Upma22. VOLGARCH PREDICT CRASHES Dependent Variable: DOWN3 Method: ML - Binary Probit (Quadratic hill climbing) Sample (adjusted): 1/02/1990 12/27/2005 Convergence achieved after 8 iterations Variable Coefficient Std. Error z-statistic Prob. C -3.599880 0.311942-11.54024 0.0000 VOLGARPRE90 3.506755 1.125058 3.116955 0.0018 DJRET(-1) -13.95141 5.379376-2.593499 0.0095 DJRET(-2) -7.714360 5.270780-1.463609 0.1433 VOL252 2.577129 1.983686 1.299162 0.1939 Mean dependent var 0.006434 S.D. dependent var 0.079964 S.E. of regression 0.080011 Akaike info criterion 0.069986 Sum squared resid 25.83771 Schwarz criterion 0.077786 Log likelihood -136.4069 Hannan-Quinn criter. 0.072750 Restr. log likelihood -157.1161 Avg. log likelihood -0.033756 LR statistic (4 df) 41.41845 McFadden R-squared 0.131808 Probability(LR stat) 2.20E-08 59

CONCLUSION GARCH vols are the best predictors of an impending crash. Market declines also have short run predictability. Many other variables are also significant in 30 s. But results from 60-now are very similar. CRASHES AND SPIKES The higher volatility, the more likely the market will crash or soar. Looking just at the crashes, they are likely preceded by high volatility and market declines Spikes also are more likely with high volatility and market declines. 60

Markowitz Optimal Portfolios 61