Overnight borrowing, interest rates and extreme value theory

Similar documents
Overnight borrowing, interest rates and extreme value theory

Introduction to Algorithmic Trading Strategies Lecture 8

Measuring Financial Risk using Extreme Value Theory: evidence from Pakistan

An Application of Extreme Value Theory for Measuring Financial Risk in the Uruguayan Pension Fund 1

Modelling catastrophic risk in international equity markets: An extreme value approach. JOHN COTTER University College Dublin

Extreme Values Modelling of Nairobi Securities Exchange Index

Analysis of extreme values with random location Abstract Keywords: 1. Introduction and Model

Financial Risk Forecasting Chapter 9 Extreme Value Theory

MEASURING EXTREME RISKS IN THE RWANDA STOCK MARKET

The extreme downside risk of the S P 500 stock index

Mongolia s TOP-20 Index Risk Analysis, Pt. 3

QQ PLOT Yunsi Wang, Tyler Steele, Eva Zhang Spring 2016

Modelling insured catastrophe losses

Relative Error of the Generalized Pareto Approximation. to Value-at-Risk

Risk Analysis for Three Precious Metals: An Application of Extreme Value Theory

An Application of Extreme Value Theory for Measuring Risk

REINSURANCE RATE-MAKING WITH PARAMETRIC AND NON-PARAMETRIC MODELS

Fitting the generalized Pareto distribution to commercial fire loss severity: evidence from Taiwan

Scaling conditional tail probability and quantile estimators

Estimate of Maximum Insurance Loss due to Bushfires

Advanced Extremal Models for Operational Risk

Lecture 6: Non Normal Distributions

Modelling Environmental Extremes

Modelling Environmental Extremes

MODELLING OF INCOME AND WAGE DISTRIBUTION USING THE METHOD OF L-MOMENTS OF PARAMETER ESTIMATION

Comparative Analyses of Expected Shortfall and Value-at-Risk under Market Stress

Application of Conditional Autoregressive Value at Risk Model to Kenyan Stocks: A Comparative Study

MEASURING PORTFOLIO RISKS USING CONDITIONAL COPULA-AR-GARCH MODEL

Characterisation of the tail behaviour of financial returns: studies from India

AN EXTREME VALUE APPROACH TO PRICING CREDIT RISK

A STATISTICAL RISK ASSESSMENT OF BITCOIN AND ITS EXTREME TAIL BEHAVIOR

Random Variables and Probability Distributions

Goran Andjelic, Ivana Milosev, and Vladimir Djakovic*

Heavy Tails in Foreign Exchange Markets: Evidence from Asian Countries

Operational Risk Quantification and Insurance

PARAMETRIC AND NON-PARAMETRIC BOOTSTRAP: A SIMULATION STUDY FOR A LINEAR REGRESSION WITH RESIDUALS FROM A MIXTURE OF LAPLACE DISTRIBUTIONS

Absolute Return Volatility. JOHN COTTER* University College Dublin

Comparing Downside Risk Measures for Heavy Tailed Distributions

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

Value at Risk Estimation Using Extreme Value Theory

Stochastic model of flow duration curves for selected rivers in Bangladesh

ANALYSIS. Stanislav Bozhkov 1. Supervisor: Antoaneta Serguieva, PhD 1,2. Brunel Business School, Brunel University West London, UK

A New Hybrid Estimation Method for the Generalized Pareto Distribution

Modelling of extreme losses in natural disasters

GPD-POT and GEV block maxima

Analysis of the Oil Spills from Tanker Ships. Ringo Ching and T. L. Yip

the display, exploration and transformation of the data are demonstrated and biases typically encountered are highlighted.

Working Paper October Book Review of

2. Copula Methods Background

THRESHOLD PARAMETER OF THE EXPECTED LOSSES

Forecasting Value-at-Risk using GARCH and Extreme-Value-Theory Approaches for Daily Returns

Intraday Volatility Forecast in Australian Equity Market

Amath 546/Econ 589 Univariate GARCH Models: Advanced Topics

Measures of Extreme Loss Risk An Assessment of Performance During the Global Financial Crisis

Portfolio modelling of operational losses John Gavin 1, QRMS, Risk Control, UBS, London. April 2004.

Cross-Sectional Distribution of GARCH Coefficients across S&P 500 Constituents : Time-Variation over the Period

Stock Price Behavior. Stock Price Behavior

VaR versus Expected Shortfall and Expected Value Theory. Saman Aizaz (BSBA 2013) Faculty Advisor: Jim T. Moser Capstone Project 12/03/2012

Fitting financial time series returns distributions: a mixture normality approach

Modelling Joint Distribution of Returns. Dr. Sawsan Hilal space

The tail risks of FX return distributions: a comparison of the returns associated with limit orders and market orders By John Cotter and Kevin Dowd *

Modeling dynamic diurnal patterns in high frequency financial data

Research Article Multiple-Event Catastrophe Bond Pricing Based on CIR-Copula-POT Model

Volatility Clustering of Fine Wine Prices assuming Different Distributions

Modelling Kenyan Foreign Exchange Risk Using Asymmetry Garch Models and Extreme Value Theory Approaches

The Comovements Along the Term Structure of Oil Forwards in Periods of High and Low Volatility: How Tight Are They?

Statistics and Finance

Case Study: Heavy-Tailed Distribution and Reinsurance Rate-making

Intra-Day Seasonality in Foreign Market Transactions

An Improved Saddlepoint Approximation Based on the Negative Binomial Distribution for the General Birth Process

Portfolio Optimization & Risk Management. Haksun Li

Midterm Examination Number 1 February 19, 1996

Analysis of truncated data with application to the operational risk estimation

Window Width Selection for L 2 Adjusted Quantile Regression

Modeling the extremes of temperature time series. Debbie J. Dupuis Department of Decision Sciences HEC Montréal

Gamma Distribution Fitting

The Application of the Theory of Power Law Distributions to U.S. Wealth Accumulation INTRODUCTION DATA

Paper Series of Risk Management in Financial Institutions

Available online at ScienceDirect. Procedia Economics and Finance 32 ( 2015 ) Andreea Ro oiu a, *

NOTES ON THE BANK OF ENGLAND OPTION IMPLIED PROBABILITY DENSITY FUNCTIONS

1 Volatility Definition and Estimation

Copula-Based Pairs Trading Strategy

Leandro Conte UniSi, Department of Economics and Statistics. Money, Macroeconomic Theory and Historical evidence. SSF_ aa

Extreme Value Theory with an Application to Bank Failures through Contagion

Generalized MLE per Martins and Stedinger

John Cotter and Kevin Dowd

EVA Tutorial #1 BLOCK MAXIMA APPROACH IN HYDROLOGIC/CLIMATE APPLICATIONS. Rick Katz

FAV i R This paper is produced mechanically as part of FAViR. See for more information.

ESTIMATION OF MODIFIED MEASURE OF SKEWNESS. Elsayed Ali Habib *

A market risk model for asymmetric distributed series of return

Modeling Credit Risk of Loan Portfolios in the Presence of Autocorrelation (Part 2)

Folia Oeconomica Stetinensia DOI: /foli A COMPARISON OF TAIL BEHAVIOUR OF STOCK MARKET RETURNS

Bivariate Extreme Value Analysis of Commodity Prices. Matthew Joyce BSc. Economics, University of Victoria, 2011

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

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

Prerequisites for modeling price and return data series for the Bucharest Stock Exchange

Subject CS1 Actuarial Statistics 1 Core Principles. Syllabus. for the 2019 exams. 1 June 2018

A Brief Account of the Turkish Economy,

Key Words: emerging markets, copulas, tail dependence, Value-at-Risk JEL Classification: C51, C52, C14, G17

An Introduction to Statistical Extreme Value Theory

Transcription:

European Economic Review 50 (2006) 547 563 www.elsevier.com/locate/econbase Overnight borrowing, interest rates and extreme value theory Ramazan Genc-ay a,, Faruk Selc-uk b a Department of Economics, Simon Fraser University, 8888 University Drive, Burnaby, British Columbia, V5A 1S6, Canada b Department of Economics, Bilkent University, Bilkent 06800, Ankara, Turkey Received 31 August 2003; accepted 20 October 2004 Available online 23 February 2005 Abstract We examine the dynamics of extreme values of overnight borrowing rates in an inter-bank money market before a financial crisis during which overnight borrowing rates rocketed up to (simple annual) 4000 percent. It is shown that the generalized Pareto distribution fits well to the extreme values of the interest rate distribution. We also provide predictions of extreme overnight borrowing rates using pre-crisis data. The examination of tails (extreme values) provides answers to such issues as to what are the extreme movements to be expected in financial markets; is there a possibility for even larger movements and, are there theoretical processes that can model the type of fat-tails in the observed data? The answers to such questions are essential for proper management of financial exposures and laying ground for regulations. r 2005 Elsevier B.V. All rights reserved. JEL classification: G0; G1; C1 Keywords: Financial crises; Risk management; Extreme value theory; Overnight rate; Federal funds rate; Turkish economy Corresponding author. E-mail addresses: gencay@sfu.ca (R. Genc-ay), faruk@bilkent.edu.tr (F. Selc-uk). URL: http://www.sfu.ca/rgencay, http://www.bilkent.edu.tr/faruk. 0014-2921/$ - see front matter r 2005 Elsevier B.V. All rights reserved. doi:10.1016/j.euroecorev.2004.10.010

548 R. Genc-ay, F. Selc-uk / European Economic Review 50 (2006) 547 563 1. Introduction The Turkish government began implementing a far reaching restructuring and reform program after the general elections in April 1999. The aim of the program was to reduce inflation from its 60 70 percent level per year to single digits by the end of the year 2002. The program gained further momentum after the country made a stand-by arrangement with the International Monetary Fund (IMF) in December 1999 and announced the technical aspects of the disinflation program. The main tool of the disinflation program was the adoption of a tablita with an exit, that is, the percent change in the value of the Turkish Lira (TRL) against a basket of foreign currencies was fixed beginning January 2000. Meanwhile, the stand-by arrangement determined a ceiling for the net domestic assets of the Central Bank. Accordingly, the Central Bank was able to create TRL liquidity only through net foreign capital inflows. The structure of the stabilization program implied that the interest rate would be market determined in line with the exchange rate depreciation and capital flows, and the volatility of interest rates would be higher than the volatility before the program. A close inspection of Fig. 1 reveals that this was the case. It was well-known by the market participants that one of the commercial banks (Demirbank) had an extremely risky position during the year 2000. The bank (with a 4000 Interest Rate (%) 3000 2000 1000 0 1990 1992 1994 1996 1998 2000 9 Interest Rate (log scale) 8 7 6 5 4 3 1990 1992 1994 1996 1998 2000 Days Fig. 1. Daily average overnight interest rates in Turkey (annualized simple) including the February 2001 crisis. The daily rates are weighted average of intraday rates. The sample period is from January 2, 1990 to February 23, 2001 with 2806 observations. Source: Central Bank of Turkey.

R. Genc-ay, F. Selc-uk / European Economic Review 50 (2006) 547 563 549 paid capital of USD 300 million) was funding its estimated USD 7.5 billion government securities portfolio mostly from the money market with very short-term obligations. On Monday, November 20, 2000, Demirbank was not able to borrow from the money market and the Central Bank stepped in to cover Demirbank s position. In the following days, market makers in the government securities market stopped posting prices and Demirbank was not able to liquidate its positions. As a result, overnight interest rates started to increase and there was a rapid capital outflow, starting on Wednesday, November 22. The heavy capital outflow and continuous decrease in the Central Bank reserves caused further liquidity pressure in the market. The Central Bank started to provide liquidity by violating the rule set by the stand-by agreement. However, the injected liquidity bounced back to the Central Bank in the form of additional demand for foreign currency. Therefore, the Central Bank stopped providing liquidity on Thursday, November 30, 2000. Immediately, the overnight interest rate reached its peak at (simple annual) 873 percent on Friday, December 1, 2000. The IMF rushed in an emergency team to discuss an emergency loan. On Tuesday, December 5, Turkish authorities announced a USD 7.5 billion rescue package with the IMF. The following morning before the markets opened, Demirbank became the 11th Turkish bank to be taken under control of the Saving Deposits Insurance Fund. After the IMF s backing of the country in December 2000, the Central Bank reserves returned to its pre-crisis level. Interest rates decreased, albeit stabilizing at a higher level than the pre-crisis average. Nevertheless, there were concerns about the Treasury s ability to borrow from the domestic market at favorable terms. A scheduled domestic debt auction of the Treasury on February 20, 2001, the day before the maturing USD 7 billion domestic debt, was aimed at borrowing approximately USD 5 billion (around 10 percent of the total domestic debt). Suddenly, the day before the auction, Turkish Prime Minister Bulent Ecevit stormed out of a key meeting of top political and military leaders stating that a dispute had arisen between himself and the country s president. The news hit the market and the stock market plunged 18 percent in 1 day. The same day, the Central Bank sold USD 7.5 billion (approximately one-third of the total official reserves) and refused to provide TRL liquidity to the system. The daily average overnight interest rates shot up to (simple annual) 2000 percent on February 20, and 4000 percent on February 21. The government responded by dropping its exchange-rate control and the TRL/ USD exchange rate went up 40 percent in 1 week. The Turkish financial crisis in February 2001 is a case study for extreme risks and risk management practices. In recent years, the problem of extreme risks in financial markets has become topical following the crises in the Asian and Russian markets, and the unexpected big losses of investment banks such as Barings and Daiwa. The value-at-risk (VaR) measures with Gaussian-type innovations failed to cope with these turmoils mainly because the evaluation of extreme risks in the Gaussian model is directly related to the variance. However, the underlying distribution may not even have a finite variance. We investigate the dynamics of the extreme values of overnight borrowing rates in the inter-bank money market for the TRL before the financial crisis of February 2001.

550 R. Genc-ay, F. Selc-uk / European Economic Review 50 (2006) 547 563 It is shown that the extreme value theory (EVT) is a useful tool in assessing the financial fragility. Estimated parameters imply that the Turkish interest rates have unbounded support and the interest rates may come from a distribution with infinite variance. We also provide estimates of overnight borrowing rates at 0.999 percentile (1 day in every 4 years) on pre-crisis data. The results indicate that every 4 years one can expect to see a day with overnight interest rates as high as 1000 percent (simple annual) in Turkey. In other words, the extraordinary levels observed during the February 2001 crisis were in the nature of the economy before they actually materialized. We also show that such extraordinary levels are a statistical impossibility for the United States. Estimated parameters of the daily US Effective Federal Funds Rate distribution imply that the distribution of US interest rates has a bounded support and all moments exist. As a result, estimated interest rates even at the 0.9999th quantile (1 day in every 10,000 days, approximately 30 years) for the US are less than the historical high of 22.4 percent. We postulate that the difference between the US and Turkey stems from the fact that the two countries have a different degree of development of the financial markets and institutions. This paper is structured as follows. In Section 2, the methodological framework is introduced. Section 3 reports the descriptive statistics and estimation results for the Turkish and the US overnight interest rates. We conclude afterwards. 2. Extreme value theory Extreme value theory is a powerful and fairly robust framework to study the tail behavior of a distribution. Embrechts et al. (1997) is a comprehensive source of the EVT for the finance and insurance literature. Reiss and Thomas (1997) and Beirlant et al. (1996) also provide extensive coverage. 1 In the following section, we present the parametric framework for our study. 2.1. Fisher Tippett theorem The theorem of Fisher and Tippett (1928) is the core of the EVT. The theory deals with the convergence of maxima. Suppose that X 1 ; X 2 ;...; X n is a sequence of independently and identically distributed random variables from an unknown distribution function F ðxþ: Denote the maximum of the first mon observations of X by M m ¼ maxðx 1 ; X 2 ;...; X n Þ: Given a sequence of a m 40 and b m such that ðm m b m Þ=a m ; the sequence of normalized maxima converges in distribution to the following so-called Generalized Extreme Value (GEV) distribution: ( e ð1þxxþ 1=x if xa0; H x ðxþ ¼ (1) e e x if x ¼ 0; 1 There have been a limited number of extreme value studies in finance literature. See Danielsson and de Vries (1997), Embrechts (2000) and Genc-ay and Selc-uk (2004) and references therein for EVT applications in finance.

R. Genc-ay, F. Selc-uk / European Economic Review 50 (2006) 547 563 551 where X is such that 1 þ xx40 and x is the shape parameter. When x40; the distribution is known as the Fre chet distribution and it has a fat-tail. The larger the shape parameter, the more fat-tailed the distribution. If xo0; the distribution is known as the Weibull distribution. Finally, if x ¼ 0; it is the Gumbel distribution. 2 In general, we are not only interested in the maxima of observations, but also in the behavior of large observations which exceed a high threshold. Given a high threshold u, the distribution of excess values of X over threshold u is defined by F u ðyþ ¼PfX upyjx4ug ¼ Fðy þ uþ FðuÞ 1 FðuÞ (2) which represents the probability that the value of X exceeds the threshold u by at most an amount yx0 given that X exceeds over the threshold u. A theorem by Balkema and de Haan (1974) and Pickands (1975) shows that for a sufficiently high threshold u, the distribution function of the excess can be approximated by the generalized Pareto distribution (GPD), i.e., the excess distribution F u ðyþ converges to the GPD below as the threshold gets large: 8 G x;b ðxþ ¼ 1 1 þ x x 1=x >< if xa0; b >: 1 e x=b if x ¼ 0; where x is the shape parameter. The GPD embeds a number of other distributions. When x40; it takes the form of the ordinary Pareto distribution and E½X k Š is infinite for kx1=x: For the security returns or high frequency foreign exchange returns, the estimates of x are usually less than 0.5, implying that the returns have finite variance. See, for instance, Longin (1996) and Dacorogna et al. (2001). The importance of the Balkema and de Haan (1974) and Pickands (1975) results is that the distribution of excesses may be approximated by the GPD by estimating x and b as a function of a high threshold u. The parameters of the GPD can be estimated with various methods such as the method of probability weighted moments or the maximum likelihood method. For x4 0:5 which corresponds to heavy-tails, Hosking and Wallis (1987) presents evidence that maximum likelihood regularity conditions are fulfilled and the maximum likelihood estimates are asymptotically normally distributed. Therefore, the approximate standard errors for the estimators of b and x can be obtained through maximum likelihood estimation. Notice that even if the identical and independent distribution condition fails, the EVT may still be an accurate approximation of the actual distribution function of maxima (Reiss and Thomas, 1997, pp. 10 and 172). If there is dependency in the data, the estimation result should be modified with an estimated extremal index. In this case, the quantile estimate would be different than the one estimated with weak dependency or independency assumption. We refer the reader to Longin (2000), Embrechts et al. (1997, Chapter 8) and Reiss and Thomas (1997, 2 Extensive coverage can be found in Gumbel (1958). (3)

552 R. Genc-ay, F. Selc-uk / European Economic Review 50 (2006) 547 563 Chapter 6) for further discussion on dependency, extremal index and its implications in practice. 2.2. The tail estimation Since F u ðyþ in Eq. (2) converges to the GPD for sufficiently large u, we have the following representation: FðxÞ ½1 FðuÞŠG x;b;u ðx uþþfðuþ, (4) where x ¼ y þ u and yx0: The last term on the right-hand side can be determined by the empirical estimator ðn N u Þ=n where N u is the number of exceedances of u and n is the sample size. The tail estimator, therefore, is given by dfðxþ ¼1 N u n 1 þ! 1=^x x u ^x. (5) ^b For a given probability level q, a percentile ð ^x q Þ at the tail is estimated by inverting the tail estimator in (5), ^x q ¼ u þ ^b! n ^x ð1 qþ 1. (6) ^x N u In statistics, this is the quantile estimation and it can be utilized in VaR estimations in finance applications. 2.3. Preliminary data analysis In statistics, a QQ-plot (quantile quantile plot) is a convenient visual tool to examine whether a sample comes from a specific distribution. Specifically, the quantiles of a hypothesized distribution are plotted against the quantiles of an empirical distribution. If the sample comes from the hypothesized distribution, the QQ-plot is linear. In EVT and its applications, the QQ-plot is typically plotted against the Gumbel distribution, i.e., a distribution with a thin-sized tail. If the data come from a distribution in the maximum domain of attraction of the Gumbel, the points on the graph would lie along a positively sloped straight line. If there is a concave presence, this would indicate a fat-tailed distribution, whereas a convex departure is an indication of short-tailed distribution. A second tool is the sample mean excess function (MEF) which is the sum of the excesses over the threshold u divided by the number of data points which exceed the threshold u. It is an estimate of the MEF which describes the expected overshoot of a threshold once an exceedance occurs. If the empirical MEF is a positively sloped straight line above a certain threshold u, it is an indication that the data follows the GPD with a positive shape parameter x: On the other hand, if the data comes from a distribution which is in the maximum domain of attraction of the Gumbel, one would observe a horizontal MEF while short-tailed data would result

R. Genc-ay, F. Selc-uk / European Economic Review 50 (2006) 547 563 553 Table 1 Descriptive statistics of the (simple annual) overnight interest rates before and after the February 2001 crisis n Mean Std. Ku. Sk. Med. Min. Max. Low High Pre-crisis 2801 73.0 49.4 82.62 7.5 67.7 13.6 873.1 34.0 107.7 Full sample 2806 75.7 99.6 943.2 26.4 67.7 13.6 4018.6 34.4 109.6 n: Sample size; Mean: Sample mean; Std: Standard deviation; Ku: Kurtosis; Sk: Skewness; Min: Minimum observed rate; Max: Maximum observed rate; Low: Rate corresponding to the 5th percentile; High: Rate corresponding to the 95th percentile. The full sample period is from January 2, 1990 to February 23, 2001 with 2806 observations. Pre-crisis period excludes the last five business days. Source: Central Bank of Turkey. Overnight interest rates are daily weighted average of inter-bank interest rates as calculated by the Central Bank of Turkey in a negative slope. See Reiss and Thomas (1997, Chapter 2) and Embrechts et al. (1997, Chapter 6) for these and other diagnostic tools and exploratory data analysis for extremes. 3. Empirical results The data source for the daily overnight interest rates (simple annual) is the Central Bank of Turkey. The daily rates are calculated by the Central Bank as a weighted average of intraday transactions in the inter-bank money market. 3 The descriptive statistics of daily average simple annual overnight interest rates before and after the February 2001 crisis are given in Table 1. The full sample period is January 2, 1990 February 23, 2001 with 2806 observations. It includes all available daily interest rate data from the inter-bank money market in Turkey. In this paper, all calculations and predictions are carried out with simple annual interest rates. In Table 1, the sample means of 73.0 and 75.7 percent correspond to compound annual interest rates of approximately 107 and 113 percent, respectively. Although it is high by the standards of a developed market, it reflects the high inflation levels and associated risk in the economy. The annual average percent increase in consumer prices (inflation) during the sample period was 75.4 percent, implying an average annual real interest rate of 18 21 percent. Both kurtosis and skewness estimates show that the interest rates are far from being normally distributed. The estimated kurtosis 82.62 before the crisis shows that the interest rate distribution has a fat-tail. The estimated skewness of 7.5 before the crisis points out that the distribution is skewed. After the crisis, both skewness (26.4) and kurtosis (943.2) estimates indicate that fat-tailness and skewness of the distribution have substantially increased. 3 The observed maximum intraday interest rate is (simple annual) 6200 percent. This means borrowing TRL at an interest rate of 17.2 percent for 1 day ð17:2 360 6200Þ: We use weighted average of intraday rates as provided by the Central Bank.

554 R. Genc-ay, F. Selc-uk / European Economic Review 50 (2006) 547 563 3.1. Tail estimation of excess interest rates It is necessary to determine a threshold interest rate to estimate the parameters of the GPD. As presented in Section 2.4, the QQ-plot and the MEF are two empirical tools for this task. In a QQ-plot, the quantiles of the empirical distribution function on the x-axis are plotted against the quantiles of the Gumbel distribution function on the y-axis. The points should lie approximately along a straight line if the data is from the assumed Gumbel distribution. Since the Gumbel distribution has a medium-sized tail, a concave relationship between the quantiles of the empirical and the Gumbel distributions indicate a heavy-tailed distribution for the time series under study. The top panel of Fig. 2 indicates that the sample points start deviating from linear behavior at around 80 percent and form a concave pattern. The sample MEF is another diagnostic tool used to determine a threshold. If the points of the MEF exhibit an upward trend, this indicates heavy-tailed behavior. Short-tailed data exhibit negatively sloped behavior whereas an exponential distribution has a flat MEF. The bottom panel of Fig. 2 demonstrates that the sample MEF is approximately linear and positively sloped after the 80 percent interest rate threshold. The examination of both panels in Fig. 2 indicate that the approximate threshold value corresponds to 80 percent. For a given threshold level, a tail estimation involves the estimation of the parameters of the GPD. The maximum likelihood estimation is a convenient way to obtain standard errors of the parameter estimates. 4 For the threshold u ¼ 80; the threshold exceedances are 389 data points which constitute the upper 13.9 percent tail of the original sample of 2801 sample points. Notice that an 80 percent overnight interest rate (simple annual) is equivalent to an annualized compound interest rate of 122 percent. During the sample period, the average annual percent increase in TRL/ USD exchange rate is 75 percent. Therefore, the threshold implies a 23 percent annualized dollar interest rate. The maximum likelihood estimates of x and b at threshold u ¼ 80 are (standard errors are in parenthesis) 0.73 (0.086) and 22 (2.03), respectively. The estimated parameters are statistically significant at the 1 percent level. As discussed in Section 2.1, when x40; the distribution is known as the Fre chet distribution and it has a fattail with tail index 1=x: The larger the shape parameter, the more fat-tailed the distribution. The value of ^x ¼ 0:73 indicates that the overnight interest rate series may come from a distribution with infinite variance. The estimated GPD model with a threshold of 80 percent interest rate (u ¼ 80) before the February 2001 crisis is presented in the top panel of Fig. 3. The estimated model is plotted in a solid curve while the interest rates above the threshold are shown in circles (on logarithmic scale). The estimated model successfully captures 4 EVIM (Extreme Value Analysis in Matlab, Genc-ay et al. (2001)) is a comprehensive toolbox for EVT analysis. This toolbox is available at http://www.bilkent.edu.tr/faruk. EVIS (Extreme Values in S-Plus) is another free suit of S-Plus functions for EVT analysis. The package may be obtained from http:// www.math.ethz.ch/mcneil. We used both packages in our calculations and estimations.

R. Genc-ay, F. Selc-uk / European Economic Review 50 (2006) 547 563 555 Gumbel Distribution 9 8 7 6 5 4 3 2 1 0 0 100 200 300 400 500 600 700 800 900 Ordered Data 160 140 120 Mean Excess 100 80 60 40 20 0 20 40 60 80 100 120 140 Threshold Fig. 2. Top: QQ-plot of daily overnight interest rates (simple annual) before the February 2001 crisis against the Gumble distribution. The quantiles of the empirical distribution function on the x-axis are plotted against the quantiles of the Gumbel distribution function on the y-axis. The points should lie along the straight line if the data is from the assumed Gumbel distribution. A concave presence indicates a fattailed distribution. Source: Central Bank of Turkey. Bottom: Sample mean excesses of daily overnight interest rates (simple annual) before the February 2001 crisis over increasing thresholds. A straight line with positive slope above a given threshold u is a sign of the GPD in tail. Notice that the plot is approximately linear and positively sloped after 80 percent, indicating the Pareto behavior in the tail.

556 R. Genc-ay, F. Selc-uk / European Economic Review 50 (2006) 547 563 1.0 0.8 Fu(x-u) 0.6 0.4 0.2 0.0 100 500 1000 x (log scale) Threshold 1.4 80.0 81.0 82.4 84.7 86.0 88.6 89.9 94.1 95.1 99.0 100.0 Shape Parameter (xi) and CI=95% 1.2 1.0 0.8 0.6 0.4 400 385 369 353 336 320 303 287 271 254 238 221 205 189 172 Exceedances Fig. 3. Top: The estimated GPD model for the excess interest rates before the February 2001 crisis. The estimated model is plotted as a curve while the actual daily overnight interest rates (simple annual) above the threshold are shown in circles (in logarithmic scale). The threshold is 80 percent interest rate, u ¼ 80: The estimated parameters are x ¼ 0:73 and ^b ¼ 22: The number of exceedances is 389. The sample period is from January 2, 1990 to February 16, 2001 which ends 1 week before the February 2001 crisis. Bottom: Estimates of shape parameter, x; at different thresholds (upper x-axis) or alternatively with different number of exceedances (lower x-axis) before the February 2001 crisis. The parameter estimates are from 100 GPD models with thresholds ranging between 80 percent and 100 percent interest rates (simple annual). The sample period is from January 2, 1990 to February 16, 2001. The sample period ends at the week before the February 2001 crisis.

R. Genc-ay, F. Selc-uk / European Economic Review 50 (2006) 547 563 557 the underlying extreme values, and the tail behavior of excess interest rates is successfully approximated by a Fre chet-type distribution. To examine the robustness of the results to the choice of the threshold value, the maximum likelihood estimation of the shape parameter, x; is carried out for GPD models with a range of thresholds between 80 and 100 percent. The estimates of x; as well as their asymptotic confidence intervals, are reported in the bottom panel of Fig. 3. On the lower x-axis the number of data points exceeding the threshold is plotted and on the upper x-axis the threshold is located. The estimates of the shape parameter, x; are plotted on the y-axis. The results indicate that the shape parameter is fairly stable within the range of 80 90 percent interest rate. Of course, there are fewer observations at higher thresholds and less statistical precision. This is reflected in the wider confidence intervals. 3.2. Point predictions at the tails It is possible to estimate a percentile value at the tail from the estimated parameters of the GPD. In Fig. 4, the tail estimates are reported with the pre- February 2001 crisis data. Since EVT is used to calculate unconditional probabilities, we do not predict the timing of the crisis but report that certain extreme interest rate levels are expected to happen every so often. On a log log plot as in Fig. 4, GPD becomes linear with negative ð1=xþ slope and it is easier to interpret results visually. In this figure, the left y-axis indicates the tail probabilities, 1 FðxÞ: The vertical dotted line starting from the x-axis is the estimated overnight interest rate which corresponds to the 1 percent tail probability on the y-axis. This value is obtained from the intersection of the 1 percent tail probability with the estimated tail and traced down to the x-axis. Therefore, the 99 percent quantile (1 day in every 100 days) interest rate corresponds to 257 percent. The dotted curve is the confidence interval calculated by the profile likelihood method. 5 The 95 percent confidence level corresponds to the intersection of the horizontal dotted line starting from the right y-axis. The range of the 95 percent confidence interval is 221 percent for the lower end and 313 percent for the higher end. Fig. 4 has a straightforward message. When the interest rate corresponding to the 99 percent quantile (1 day in every 100 days) is evaluated at the 95 percent confidence level, the maximum simple annual interest rate would go as high as 313 percent. This is a substantially high overnight interest rate but it is not in the order of ten thousands or millions percent compound interest rate as it was observed during the February crisis. Let us now go further in the tail to evaluate further extreme possibilities and investigate the 99.9 percent quantile estimate (1 day in every 4 years). This is presented in Fig. 5 where the interest rate estimates at the 99.9 percent quantile (1 day in every 4 years) are plotted as a function of the threshold (upper x-axis) or alternatively as a function of exceedances (lower x-axis) with the 5 The usual Wald standard errors are computed from the inversion of the Hessian of the log-likelihood function. Confidence limits from the inversion of the likelihood ratio statistic are called profile likelihood confidence limits.

558 R. Genc-ay, F. Selc-uk / European Economic Review 50 (2006) 547 563 0.0500 95% CI (lower) 95% CI (upper) Tail Probability (log scale) 0.0050 Estimated 99% quantile 95 99 Level of Confidence 0.0005 x-value corresponding to estimated 99% quantile Confidence Interval (CI) 100 257 500 1000 x (log scale) Fig. 4. Tail estimate for the excess interest rates before the February 2001 crisis. The left y-axis indicates the tail probabilities 1 FðxÞ: The vertical dotted line from x-axis is the estimated interest rate (at 257 percent) and intersects with the 1 percent tail probability on the left y-axis. Therefore, it is the estimated interest rate at the 99 percent quantile. The dotted curve is the confidence interval around this estimate, calculated by the profile likelihood method. It intersects with the 95 percent horizontal confidence line from the right y-axis and the corresponding 95 percent confidence interval is ð221; 313Þ in percent. Notice that the confidence interval becomes larger when the confidence level is increased (a down movement in the level of confidence on the right y-axis). The sample period is from January 2, 1990 to February 16, 2001 which ends 1 week before the February 2001 crisis. pre-february crisis data. The estimated interest rate settles above 1000 percent (simple annual) for most thresholds although variations are higher at higher thresholds. The upper confidence interval indicates that before the crisis, every 4 years one could expect to see a day with an overnight interest rate as high as 4800 percent. The results show that overnight interest rates observed during the financial crisis were in the nature of the economy before they materialized. The policy makers and the financial authorities in Turkey, along with the Turkish commercial banks, simply ignored the possibility that such extreme tail events might occur. 3.3. A comparison with the US overnight interest rates Although the Turkish daily overnight rate is an excellent case study with high volatility and a thick-tailed distribution, this data set has not been studied widely in the literature and is not well-known. Hence, we have repeated the extreme value analysis with the daily US Effective Federal Funds Rate as a comparison. The Federal Funds Rate (FFR) is the interest rate that banks with excess reserves at a

R. Genc-ay, F. Selc-uk / European Economic Review 50 (2006) 547 563 559 5000 Threshold 80.0 81.0 82.4 84.7 86.0 88.6 89.9 94.1 95.1 99.0 100.0 4000 0.999 Quantile and CI=0.95 3000 2000 1000 0 400 385 369 353 336 320 303 287 271 254 238 221 205 189 172 Exceedances Fig. 5. Interest rate estimates of quantile 0.999 as a function of the threshold (upper x-axis) or alternatively as a function of exceedances (lower x-axis) before the February crisis. These estimates are from 100 GPD models with thresholds ranging between 80 percent and 100 percent interest rates (simple annual). The estimated interest rate at 0.999 quantile settles above 1000 percent at most thresholds although variations are higher at the higher thresholds. The Wald confidence intervals indicate that interest rates as high as 4000 percent is a possibility. The sample period is from January 2, 1990 to February 16, 2001 which ends 1 week before the February 2001 crisis. Federal Reserve District Bank charge other banks that need overnight loans. The sample period is from July 1, 1954 to December 31, 2000. The sample size is 16,986 daily observations. The data source is the Federal Reserve Board of Governors. 6 In Fig. 6 (bottom, left), the quantiles of the empirical distribution function of the daily US Effective Federal Funds Rate on the x-axis are plotted against the quantiles of the Gumbel distribution function on the y-axis. The convex relationship between the quantiles of the empirical and the Gumbel distributions indicate a thin-tailed distribution for the US interest rates. The GPD estimations with different numbers of exceedances in this data set indicate that the estimated shape parameter ^x is 6 We do not intend to model the FFR dynamics in the US: It is far beyond the scope of this study. There is ongoing discussion as to whether interest rates are stationary or not. Conventional economic and finance theory often assumes that interest rates are stationary and utilize a model which is a mean-reverting process. However, some empirical studies provide evidence that interest rates are nonstationary processes.

560 R. Genc-ay, F. Selc-uk / European Economic Review 50 (2006) 547 563 25 25 24 24 0.9999 Quantile (CI = 0.95) 23 22 21 20 0.999 Quantile (CI = 0.95) 23 22 21 20 19 100 150 200 250 300 350 400 450 500 Exceedances 19 100 150 200 250 300 350 400 450 500 Exceedances 7 0.25 13.51 14.04 Threshold 14.60 15.15 16.01 17.53 6 0.3 Gumbel Distribution 5 4 3 2 1 0 13 14 15 16 17 18 19 20 21 22 Ordered Data Shape(xi) (CI = 0.95) 0.35 0.4 0.45 0.5 0.55 0.6 0.65 700 600 500 400 Exceedances 300 200 Fig. 6. Top: Effective Federal Funds Rate (daily) estimates of 0.9999th (top, left) and 0.999 (top, right) quantiles as a function of exceedances. The estimated interest rate even at 0.9999 quantile (1 day in every 30 years) is less than the historical high of 22.36 percent (July 22, 1981). Bottom left: QQ-plot of Effective Federal Funds Rate (daily, simple annual) against the Gumbel distribution. The points should lie along the straight line if the data are from a Gumbel distribution. A convex presence indicates a thin-tailed distribution. Bottom right: Estimates of shape parameter, x; at different thresholds (upper x-axis) or alternatively with different number of exceedances (lower x-axis) for the US. The parameter estimates are from 100 GPD models with thresholds ranging between 13 percent and 18 percent interest rates (simple annual). The sample period is from July 1, 1954 to December 31, 2000. Data source: Federal Reserve Board of Governors, H.15 Release. remarkably stable around 0:40 (see Fig. 6). This implies that the maximum daily US Effective Federal Funds Rate distribution has a bounded support and all moments exist. As a result, the extraordinary overnight rates observed in Turkey are a statistical impossibility for the United states. Indeed, Fig. 6 plots estimated interest rates at the 0.9999th quantile (1 day in every 10,000 days, approximately 30 years), and at the 0.999th quantile (1 day in every 1000 days, approximately 3 years) as a function of the number of exceedances. The estimated interest rate is very stable around 20.5 percent at the 0.1 percent tail and at around 22 percent at the 0.01 percent tail regardless of the number of exceedances. Moving in the tail of the distribution from the 0.1 percent to the 0.01 percent region increases the estimated interest rate by only around 150 bases points (1.5 percentage point). Notice that it is not the observed maximum that would determine the possible extreme values at

R. Genc-ay, F. Selc-uk / European Economic Review 50 (2006) 547 563 561 certain quantiles. It is the shape of the extreme value distribution at which maximum interest rates converge. Clearly, the maximum interest rates in Turkey converge to a distribution which may not even exhibit a finite variance while the maximum interest rates in the U.S. has bounded support. 7 We postulate that the difference between the U.S. and Turkey stems from the fact that the two countries have a different degree of development of the financial markets and institutions. First, financial markets in Turkey are very shallow and the financial structure is very weak. Total assets of all commercial banks in Turkey (around USD 150 billion) are comparable to the assets of a medium size international bank. Second, the Turkish economy has a long history of high inflation with accommodative monetary policies conducted by the Central Bank which has been under the strong influence of the government. The FED, on the other hand, may be considered as a truly independent monetary authority. It is reasonable to expect that the distributional characteristics of interest rates in Turkey are changing in view of the fact that Turkey now has an independent central bank and the financial markets are developing very fast. 4. Lessons from the Turkish crises and conclusions Financial crises in emerging markets in general and the Turkish crises in November 2000 and February 2001 in particular, provide several lessons for investors both in developing and developed countries. Since a significant portion of total savings in developed economies are invested in emerging markets by hedge funds, mutual funds and other institutions in the form of portfolio investment, the costs of financial crises are not confined to the residents of emerging market countries. Therefore, a careful investigation of the market dynamics and the causes of crises in these economies would benefit investors at large by increasing the investor awareness. Fundamental macroeconomic indicators such as growth rate, current account balance, real exchange rate, budget deficit, export import ratio and debt income ratio are the main sources for assessing the current and future status of an economy. Therefore, they play a significant role in the decision making process of the IMF, credit rating agencies and multinational fund managers. One of the lessons from the Turkish crisis is that even if there is no deterioration in fundamental indicators, the balance sheet issues in the financial sector may create an environment in which even a small shock can lead to a total collapse of the system. In particular, a balance sheet mismatch situation (funding long-term illiquid assets with short-term obligations) combined with slack supervision and regulation is an invitation for a liquidity and currency crisis. Among others, Eichengreen (2001) investigates both the Argentinian and Turkish crises in detail. He points to the vulnerability of the banking sector as the main 7 Given estimated values ^x ¼ 0:40 and ^b ¼ 3:75 with a threshold interest rate 13.5 percent, the estimated endpoint of the maximum interest rate distribution in the US is 22.67 percent.

562 R. Genc-ay, F. Selc-uk / European Economic Review 50 (2006) 547 563 source of crisis in Turkey. Although it was expected that the lira was going to appreciate during the program, the government could not commit to an upfront devaluation as the Turkish banks and the private sector had large unhedged foreign exchange exposures. It was hoped by the program designers that the banking sector would strengthen before the economy would move into a floating exchange rate regime. As Eichengreen (2001) points out, this strategy created a moral hazard in the system and an incentive to strengthen both balance sheets and supervision diminished. 8 As a result, short foreign exchange positions of the banks doubled during the initial phase of the program. The direct consequence of this was that the fear of destabilizing the economy forced the authorities to resist any correction in the exchange rate, even if it meant extraordinary increases in overnight interest rates. Given the government s strong commitment to exchange rates, the relevant question from the investors point of view was what extraordinary interest rates may be observed under extreme situations?. Our estimation results from the pre-crisis data indicate that every 4 years one could expect to see a day with overnight interest rates over 1000 percent (simple annual). In other words, the extraordinary levels observed during the crisis were in the nature of the economy before they actually materialized. The developments during and after the crisis indicate that both policy designers and the financial sector failed to take this characteristic of the economy into account. Another important implication of our results for risk management is that the degree of development of the financial markets and institutions in different countries implies different levels of risk and associated risk premium. Therefore, risk managers should abstain from applying the same risk model and the same value-at-risk methodology to different economies. A statistical impossibility for a developed economy might easily be 1 day in every 100 days in an emerging market. Acknowledgements Ramazan Genc- ay gratefully acknowledges financial support from the Swiss National Science Foundation under NCCR-FINRISK, Natural Sciences and Engineering Research Council of Canada and the Social Sciences and Humanities Research Council of Canada. We are very grateful to Alexander J. McNeil for providing his EVIS software. Faruk Selc- uk gratefully acknowledges the financial support from the Research Development Grant Program of Bilkent University. References Balkema, A.A., de Haan, L., 1974. Residual lifetime at great age. Annals of Probability 2, 792 804. Beirlant, J., Teugels, J., Vynckier, P., 1996. Practical Analysis of Extreme Values. Leuven University Press, Leuven. 8 According to Eichengreen (2001), the moral hazard was created because exchange risk was socialized, that is to say, its strategy committed the government to preventing the exchange rate from moving and therefore to compensating the banks for their losses if the policy failed.

R. Genc-ay, F. Selc-uk / European Economic Review 50 (2006) 547 563 563 Dacorogna, M.M., Genc-ay, R., Mu ller, U.A., Olsen, R.B., Pictet, O.V., 2001. An Introduction to High Frequency Finance. Academic Press, San Diego. Danielsson, J., de Vries, C.G., 1997. Tail index and quantile estimation with very high frequency data. Journal of Empirical Finance 4, 241 257. Eichengreen, B., 2001. Crisis prevention and management: Any new lessons from Argentina and Turkey. Background paper written for the World Bank s Global Development Finance 2002. Economics Department, University of California, Berkeley, CA. Embrechts, P., 2000. Extremes and Integrated Risk Management. Risk Books and UBS Warburg, London. Embrechts, P., Klu ppelberg, C., Mikosch, T., 1997. Modeling Extremal Events for Insurance and Finance. Springer, Berlin. Fisher, R.A., Tippett, L.H.C., 1928. Limiting forms of the frequency distribution of the largest or smallest member of a sample. Proceedings of Cambridge Philosophical Society 24, 180 190. Genc-ay, R., Selc-uk, F., 2004. Extreme value theory and value-at-risk: Relative performance in emerging markets. International Journal of Forecasting 20, 287 303. Genc-ay, R., Selc-uk, F., Ulugülyaǧcı, A., 2001. EVIM: A software package for extreme value analysis in Matlab. Studies in Nonlinear Dynamics and Econometrics 5, 213 239. Gumbel, E., 1958. Statistics of Extremes. Colombia University Press, New York. Hosking, J.R.M., Wallis, J.R., 1987. Parameter and quantile estimation for the generalized Pareto distribution. Technometrics 29, 339 349. Longin, F.M., 1996. The asymptotic distribution of extreme stock market returns. Journal of Business 63, 383 408. Pickands, J., 1975. Statistical inference using extreme order statistics. Annals of Statistics 3, 119 131. Reiss, R., Thomas, M., 1997. Statistical Analysis of Extreme Values. Birkhäuser, Basel.