Interest Rate Volatility : A Consol Rate Approach

In this paper, we propose a new methodology to estimate the volatility of interest rates in the euro area money market. In particular, our approach aims at avoiding the limitations of market implied volatilities, i.e. the dependency on arbitrary choices in terms of maturity and frequencies and/or of other factors like credit and liquidity risks. The measure is constructed as the implied instantaneous volatility of a consol bond that would be priced on the EONIA swap curve over the sample period from 4 January 1999 to 21 November 2013. Our findings show that this measure tracks well the historical volatility since, by dividing the consol excess returns by our volatility measure. This removes nearly entirely excess of kurtosis and volatility clustering, bringing the excess returns close to an ordinary Gaussian white noise.


Introduction
It is common that the logarithm of asset prices is not well depicted by a simple Brownian motion.The first feature of a Brownian motion, which is a Gaussian white noise, is the absence of excess of kurtosis, and the square or the absolute value of this first difference exhibits no autocorrelation.By contrast, the returns of an asset exhibit noticeable excess of kurtosis, and their square or their absolute returns exhibit positive autocorrelations.Attempts to find a better modelling for the log-prices may include the adjunction of jumps or the specification of a variable, generally random, volatility.A variable volatility is sufficient to produce the excess of kurtosis and the positive autocorrelation of squared and absolute values-known as "volatility clustering".Therefore, it is legitimate to wonder whether the recourse to variable volatilities is sufficient to capture the observed features of kurtosis excess and volatility clustering.
To answer this question, one may start with the following thought experiment.Consider the observed time series of a liquidly traded asset price, for instance the euro-dollar exchange rate.Ignoring the short-term interest rates that prevail on those two currencies, we assume for simplicity that two subsequent observations are never identical.The corresponding return is thus the difference of the log-price at time t + 1 and at time t.Then, we define the proxy of its instantaneous volatility, at time t, by the absolute value of that return.The normalised return is therefore defined as the ratio of the return and of that proxy.By construction, this normalized return is a random sequence of +1 and −1.It will exhibit no autocorrelation, and its distribution will be even less leptokurtic than a Gaussian one, which means that it has a negative excess of kurtosis.This can be the basis of a procedure for constructing the instantaneous volatility of the exchange rate, if and only if the value of that rate at time t + 1 is known at time t, which is obviously not the case.So one is led to reformulate the initial question as follows: is it possible to construct a proxy of the instantaneous volatility at time t with data available at time t, in such a way that the return normalized by this proxy exhibits neither excess of kurtosis nor autocorrelation of its square or of its absolute value?
Here, there are two possible ways forward.One might construct that proxy on the basis of the history of the series available at time t, or, on the other hand, one might take recourse to exogenous data available at time t, namely the corresponding market-implied volatilities.The second way should be favoured because it does not assume anything about the statistical properties of the price series.In the example of the euro-dollar exchange rate, the exercise could be described by the following steps: a) to reconstruct an instantaneous volatility from the available market implied volatilities; b) to calculate the excess return of the exchange rate taking into account short-term interest rates of the two currencies; c) to calculate the normalised excess return by dividing the excess return by the instantaneous volatility; and d) to compare the statistical properties of the raw excess returns and of the normalised one.If the raw excess returns exhibit excess of kurtosis and volatility clustering whereas the normalised retruns do not, then it means that it is possible to model the exchange rate as an Ito process with variable volatility.
This paper follows the same approach with an application to the money market yield curve.The case of the yield curve is substantially more complex than the case of the exchange rate since it has a structure in contrast with a single interest rate which is reduced to a number.Furthermore, it is impossible to ignore the short-term interest rate in the case of the yield curve as it is a constitutive element of the curve under study.Finally, there is not one traded asset that could be representative of the yield curve.
The contribution of our analysis to the literature is to suggest a simple way to circumvent the above-mentioned technical difficulties.In particular, it uses as the representative asset a perpetuity paying a constant rate of dividend, called the consol rent, priced from the curve.It defines its excess return, constructs an instantaneous volatility from market-implied ones, constructs the corresponding normalised excess return and compares the statistical features of raw and normalized excess returns.Our results show that the consol rate can be described as an Ito process with variable volatility, hence avoiding the inclusion of jumps.
The paper is organised as follows.After this introduction, Section 2 summarizes the purpose and the main steps of our approach before Section 3 presents the useful definitions of key parameters, namely the consol excess return, the consol volatility and the consol normalized excess returns.Section 4 then compares the statistical features of consol excess return and consol normalized excess return when the process of the curve is simulated according to two arbitrage-free theoretical models.Section 5 compares the statistical features of consol excess return and consol normalized excess return for the EONIA swap curve.Finally, Section 6 concludes while the reconstruction of the consol normalized excess returns and volatility for that empirical EONIA swap curve is described in the annex.

Presentation of the Exercise
The proposed exercise requires two preliminary steps: first, an accurate definition of the excess return relevant for the consol rate; second, a construction of the consol rate, excess return and volatility.Then, a comparison between the excess return and the normalized excess return of excess of kurtosis and of ACF of squared and absolute values for an empirical yield curve will allow to evaluate the quality of the measure.The empirical yield curve chosen in the paper is the one implied by the EONIA swaps, which is considered as the riskless curve for the euro currency.
The central idea is to claim that if some computed volatility does lead to normalized excess returns free of excess of kurtosis and of autocorrelation in squared and absolute values, then it can be regarded as the volatility of the consol rent 1 .In order to support the previous conclusion, it is therefore useful to see what happens in cases where the true volatility is effectively known.It is effectively known in the case of some arbitrage-free yield curve models.So a way to judge the validity of the proposed criterion is to apply it in the case of simulated trajectories of the yield curve, when the simulation follows an arbitrage-free model for which the volatility can be effectively computed, and for which the length of the sample is the same as the length of the available sample of empirical data.

Consol Rate, Consol Excess Return, Consol Volatility
Let us first recall some key mathematical notions related to the concepts of yield curve and of consol rate, before discussing the basic properties of a consol volatility.From those basic properties, it follows that the consol volatility reduces the consol excess returns to a Gaussian white noise, within the mathematical frameworkof a continuous-time arbitrage-free model.This hints that, in the real world, a correct measure of volatility should also be able to reduce the excess returns to a Gaussian white noise.This can be tested by examining whether dividing the excess returns by the volatility actually reduces the leptokurticity and the volatility clustering.This test will be illustrated with simulated data, for which the true volatility can be known ex ante, before being presented for the actual euro data.In this case, a success of the test indictes that the true volatility can be, and has effectively been, recovered.

Notational Convention
The purpose of this paragraph is to recall and define the key notions of this paper within the continuous-time framework, as well as their basic mathematical properties.For convenience, the following convention is adopted: Latin letters denote dates while Greek letters refer to delay/duration between two dates 2 .

Zero-Coupon and Forward Rates
Yield curves are formally defined as functions of a continuous time parameter, which associate an interest rate to a maturity within a (theoretically unbounded) maturity set.Yield curves are also usually expressed in mainly two ways: a) as zero-coupon interest rate curves; or b) as instantaneous forward interest rate curves.These two ways are equivalent and convey the same quantity of information.With ( ) z τ the zero-coupon interest rate at maturity τ , ( ) f τ the forward interest rate at maturity τ (whereby both interest rates are continuously com- pounded) and ( ) P τ the spot price of the zero coupon of maturity τ , i.e. the present value of one currency unit to be paid over τ , one can write that: and ( ) ( ) ( ), which implies the following mathematical relationship between zero-coupon interest rate and the forward interest rate: or, conversely, as the change of variables is duly invertible: From Equations ( 3) and ( 4), it follows two basic but important properties follow: 1 See also the discussion on the relevance and the scope of the various volatility measures in [1]. 2 For example, a bond observed at time t and having maturity date T shall have a time to maturity τ satisfying to the relation: f . is constant.In this particular case, it means that the constant value taken by ( ) f . is the same as the constant value taken by ( ) z . ,which would imply a flat yield curve while the constant value taken by both ( ) z .and ( ) f . is called the level of the flat yield curve; b) irrespective of the shape of the yield curve, ( ) are always equal and the corresponding unique value defines the short-term interest rate r , i.e.

( ) ( )
. Note also that one defines as parallel shift in the remainder of the analysis a transformation of the yield curve such as a constant is added to the zero-coupon interest rates, or, equivalently, of the forward rates following from Equation (3).

Consol Bond and Consol Rate
A consol bond is defined as a perpetual (infinite horizon) bond paying continuously a constant rate of money, which is called the coupon flow.By definition, the consol price C is defined as the price of the consol bond divided by the coupon flow 3 , which can be expressed, in terms of the yield curve, as follows: where the sensitivity of the consol rate to a parallel shift of the yield curve, ξ , is given by the following product: yD ξ = (9)   This dimensionless number ξ is equal to 1 in case of a flat curve.Empirical evidence shows that yield curves have usually ξ smaller than, but close to, 1.We say that that a financial price is in constant terms, by opposition to current terms, when it is expressed in currency values of a fixed reference date in the past rather than in currency values of the current date.
Finally, the consol wealth process, denoted A hereafter, is defined as the wealth, expressed in constant terms, of an ideal investor facing no transaction costs or short-selling restrictions who holds a portfolio of consol bonds in which he reinvests automatically the whole coupon flow at the prevailing market price.Therefore, the corresponding (infinitesimal) consol excess return, denoted ( ) whereby the first two terms refer to the nominal gain (or loss) due respectively to the change in the market price of the consol bond ( ) C C and to the coupon flow ( ) t C while the third term ( ) d r t refers to the carry cost of the position, i.e. holding the portfolio of consol bonds.

Volatility of a Consol Bond
The quadratic variation of a process t X is denoted [ ] t X , so that Ito's formula is: The volatility of the consol bond, t σ , is defined as: The differential d t L (with t L denoting the logarithm of the consol wealth process, hereafter referenced as to consol performance) can be obtained by applying the Ito's iteration to Equation (10), i.e.: which leads to the following identity using the Ito's formula: Finally, the normalized excess return is defined as: which, by combining Equations (15), ( 16) and (17), yields: We will then apply the mathematical specification of the consol rate, and the calculation of the corresponding volatility as discussed from Equations (11) to (18), to standardise the volatility measure for interest rates in the money market.

Risk-Neutral Probability
Within the framework of Heath-Jarrow-Morton [3], the dynamics of t y , t L and t A are fully specified by the stochastic process of t σ under risk-neutral probability.t A must be a martingale, so Equations ( 10) and (12) imply: with t W denoting a Wiener process, from which it immediately follows that: where t N itself is also a Wiener process under risk-neutral probability.This implies that: [ ] Equation ( 23) is valid not only under risk-neutral probability but also under any equivalent probability to the risk neutral probability, as it pertains to the quadratic variation only.Finally, Equations ( 18) and ( 22) yield that: Note the left-hand side of Equation ( 24) is defined as the normalized excess return.It follows from Equation (24) that the the normalized excess returns-i.e. the excess returns divided by the true volatility-should be akin a Gaussian white noise under risk-neutral probability.As we know, empirical excess returns of any financial asset usually differ from a Gaussian white noise by two properties: a) their empirical distribution more kurtosis than the Gaussian distribution (the property of leptokurticity); and b) their absolute values (or their square values) have a positive serial correlation (the property of volatility clustering).
To the extent that the Ito-process modelisation is a realistic representation of the consol rate dynamics, one should be able to remove those two properties, by dividing the excess returns by a correct measure of the underlying volatility, and that operation of normalization would recover the underlying Gaussian white noise process.

Data
The dataset contains only TARGET working days; it contains all the TARGET working days from 4 January 1999 to 21 November 2013, which represents 3816 TARGET working days.The financial instruments taken into account are handled in the OTC market.They consist into: short-term unsecured deposit of maturity 1-day (overnight, tom-next and spot-next), EONIA swaps from 1-week to 30-year, 6-month EURIBOR swaps for the corresponding maturities, at-the-money implied volatilities of options on the EURIBOR swaps.
We used the options on 6-month EURIBOR swaps with option maturity 1-month, 3-month, 6-month, 1-year, 2-year, 3-year, 4-year and 5-year, and with underlying EURIBOR swap maturity 1-year, 2-year, 3-year, 4-year, 5-year, 7-year, 10-year, 15-year, 20-year, 25-year and 30-year.Other maturities of options or of underlying swaps are represented in the quotes contributed by brokers, but their history may start at relatively recent dates, which makes preferable not to use them.Besides, the EONIA fixing is included in the dataset.Each instrument or fixing is identified in the Reuters database by a unique RIC (Reuters Instrument Code).
As the financial instruments taken into account are handled on the OTC market, we made use of quoted data, generally given as a bid-ask spread from which we retained only the mid.We gave a preference to quotes issued by the broker ICAP, and when not available, our primary fallback was the generic quote of Reuters, which contains the latest quote issued by a bank or broker at the time of its snapshot or of its contribution.In case of missing data, the data set is completed by a reconstruction of data as described in detail in Annex..

Testing the Robustness of the Benchmark Rule
To assess the correctness of a measure of consol volatility ("benchmark rule" hereafter), we will test whether the excess returns of the consol bond, when normalized by our volatility measure, resemble a Gaussian white noise process, i.e. that both leptokurticity and volatility clustering are essentially reduced.Removing the volatility clustering is not sufficient to assess the correctness of the measure.It is easy to see that for an indicator X that oscillates rapidly enough, the absolute value, and the square, of the ratio dL X have small autocorrelations: the rapid oscillations may remove the volatility clustering but only at the expense of an increase of the leptokurticity.This makes necessary to require the reduction of both elements (volatility clustering and leptokurticity) in order to assess the quality of the volatility indicator.
Two formal tests are presented in this section.First, we conduct simulation of affine factor models (which allow the knowledge of the true volatility ex ante) and apply the test on the result of those simulations: this aims at checking that our implementation of the test actually behaves as it is supposed to.Second, we apply the test to the actual history of the EONIA curve, and to our reconstruction of the consol volatility, and assess by this way the quality of this reconstruction of the consol volatility.

Test Based on Simulation
The consol performance L is also well defined in the case of simulated data.In simulating the evolution of the yield curve under an arbitrage-free affine factor model, one can check whether a) d t L exhibits (or not) lepto- kurticity and volatility clustering; and b) d d 2 does exhibit less (or no) leptokurticity and less less (or no) volatility clustering.This type of exercise is of particular interest as it would provide a benchmark for the results of our test and hence a natural comparison point for the results based on the empirical data set.

General Setting of the Simulation Exercise
Denote with 1 t and 2 t two consecutive TARGET working days.Excess returns are constructed by setting the cost of carry equal to the EONIA observed at the close of business of 1 t .Normalized excess returns are con- structed using the volatility computed at the close of business of 1 t , (and not 2 t ).To assess the existence of leptokurticity, we examine whether the excess of kurtosis differs from zero (contrary to a normal distribution where its value is zero).To test the existence of volatility clustering, we use the correlation of the absolute values of two consecutive returns, and the correlation of the square of two consecutive returns.
Both simulations are run on 3816 TARGET working days.

Simulation Models
For the sake of simplicity, we will use the special case of arbitrage-free models known as affine models with constant parameters, with continuous time setting, and continuous trajectories.We will perform our simulations in the case of the generic 1-factor5 affine model and in the case of the generic 2-factor affine model.The 1-factor model is in effect the simplest choice, and the 2-factor model is its more natural generalisation.
In an arbitrage-free model, the zero-coupon bond price takes necessarily the form of the expectation: ( ) whereby the expectation refers to the so-called risk-neutral probability.Furthermore, the probability under which the short-term rate s r actually diffuses should be equivalent, in the probabilistic sense6 , to the risk- neutral one.We will refer to that second probability as to the data-generating probability.
We have performed several attempts with different choices of parameters.The results are always that the normalised excess returns behave close to a gaussian white noise, and that the raw excess returns behave less close to a gaussian white noise.Yet the contrast between the normalized excess return's and the raw excess return's behaviours may be more or less pronounced.Typically, we find a excess of kurtosis ranging between zero and six for the raw case, and close to zero for the normalized case.The simulations that we present here will roughly correspond to a median case.
Simulation with 1-factor model-The generic 1-factor affine model is coincident with the Duffie and Kan [4] 1-factor model (hereafter DK1).In the DK1 model, the risk-neutral probability is the solution 7 of the stochastic differential equation (SDE): ( ) with a and b positive constants, c a real constant, and ν a positive or zero constant, c and ν not being both equal to zero at the same time.By construction, this model has thus four parameters ( a , b , c and ν ) and one unique factor which can be identified to the short-term interest rate, s r .It evolves between 2 c ν − and +∞ .
The data-generating probability should be equivalent (in probabilistic terms) to the risk-neutral probability.Since we are here in a continuous-time setting, the equivalence condition implies that the Brownian part of the data-generating probability is also provided by the expression 2 d s s c r W ν + with the same c and v whereas it has no implication whatsoever as regards the form of the drift terms of the data-generating probability.Nevertheless, again for the sake of simplicity, one focuses on the specific case for which the stochastic equation defining the data-generating probability takes a form similar to Equation (26), only allowing for different values for the drift parameters a * and b * .The historical, denoted with stars, determine the computation of the diffusion of factors.The risk neutral, denoted without stars, determine the computation of the curve at each step.
Note that the DK1 model is the generic case of the 1-factor affine model with constant coefficients.It is also the simplest possible arbitrage-free model of the yield curve, which includes two specific cases: a) the original Vasicek model [5], obtained when ν is set to zero; b) the original Cox-Ingersoll-Ross model [6] (the CIR model hereafter), obtained when c is set to zero 8 .
The functional form resulting from the DK1 model happens to be realistic (see [7]).Yet, the DK1 model appears in practice relatively far from the real motion of actual yield curves.
To perform the simulation, we need to choose values for the parameters and for the initial value of the factor.There is no compelling reason to choose one set of parameters rather than another.We have adopted values producing yields which are realistic for the euro, but this, strictly speaking, is not a constraint for the current purpose.We perform the simulation on the basis of the values for the parameters and initial values of the factors listed in Table 1.
Based on this simulation model, the following Table 2 reports the results for the (raw) consol excess returns and the normalized consol excess returns: Simulation with 2-factor model-The generic 2-factor affine model is coincident with the model presented in [8] (hereafter GS2), but we will rephrase it with another choice of parameters and factors, in order to ensure formal consistency with the previous discussion.
In the GS2 model, excluding again the case where the short-term rate is bounded away from zero, the riskneutral probability can be seen as the solution of SDE described in Equation (26): with 1 a , 2 a , 11 b , 12 b , 22 b , c , ν , s p and s r satisfying certain constraints, namely: ( ) and either: ( ) ( ) The square root sign over the matrix has to be interpreted as the matrix square root operator, as opposed to an operator acting component by component.The model has seven parameters, 1 a , 2 a , 11 b , 12 b , 22 b , c and ν , and two factors of which the first one is identified to the short-term rate s r .As was the case for the 1-factor affine model, s r evolves between 2 c ν − and +∞ .The second factor, denoted with p , has no particular economic interpretation.It evolves between −∞ and +∞ , and its physical dimension is the same one as for a volatility, or equivalently, as the square root of a rate.
From the Equation ( 27), it appears that if 12 b is set to zero, the GS2 model is reduced to the DK1 model.The factor s r in this case is not influenced by the dynamics of the other factor p and follows simply the solution of Equation (26).
Again, while the mathematic structure of the model only obliges us to have the same Brownian part for the risk-neutral and data-generating probabilities, we focus nevertheless on data-generating probabilities sharing the same algebraic form with the risk-neutral one.The data-generating probability is then given by an equation similar to (27), in which parameters 1 a , 2 a , 11 b , 12 b , 22 b are replaced by counterparts denoted with an asterisk, following similar constraints than the parameters without asterisk.
We perform the simulation on the basis of the following values for the parameters and initial values listed in Table 3.
The results for the consol excess returns on the left side and for the normalized excess return on the right side of Equation (24) become those listed in Table 4.
We obtain again the expected results, regarding the fact that leptokurticity and volatility clustering are present in the excess returns and removed from the normalized excess returns.Leptokurticity and volatility clustering reach values comparable to those of the 1-factor model.Yet, as we will see, they still cannot be compared with what is observed on empirical excess returns.

Test Based on Empirical Data
Similarly to the general specifications recalled in the previous sections, the consol rate and the corresponding volatility for the EONIA is calculated over the period from 4 January 1999 to 21 November 2013, i.e. 3816 TARGET working days.Table 5 reports the results based on actual consol volatility calculatd from the empirical data.A graphical representation of the volatility measure based on the consol rate specification together with the level of the consol interest rate is reported in Figure 1.
As shown by Table 5, the excess of kurtosis and the volatility clustering exhibited by the normalised consol excess returns is substantially lower than those exhibited by (raw) excess returns.It is also interesting to underline that the excess of kurtosis and the volatility of the normalized consol excess returns based on empirical data appears even lower that the value they take for the raw excess returns (i.e.before nomalisation) in the case of the simulations.

Concluding Remarks
This paper proposes a new measure of volatility derived from the specification of the consol rate for the EONIA swap curve in order to have an accurate estimation of volatility free from any model-based specifications and relaxed from maturity and frequency constraints.We demonstrate that this volatility measure is very close to the true (unobserved implicit) instantaneous volatility as it allows the excess returns of the consol rate to display a Gaussian white noise process (under risk-neutral probability or any similar probability) once normalised by this measure.This finding is quite powerful for several reasons.
First, it legitimates the use of yield curve dynamics being free of jumps.It is thus more parsimonious.The discrepancy between the statistical features of the excess returns and those of a Gaussian white noise can be brought back to the mere variability of the volatility and do not require the intervention of jumps, at least for what regards the particular EONIA swap curve.Second, our findings allow an homogeneisation of volatility measure (with a forward-looking feature), providing information for the entire market without being restricted to one particular maturity or to the frequency of coupons/cash flows idiosyncratic to a particular benchmark instrument.
The restrictive nature of standard volatility measures (due to the strong link to a certain maturity and/or frequency) usually limit the use of volatility measure in times series regressions.The consol volatility provides new research avenues as regards volatility transmission and/or assessment of market stress.should be such that: • it re-prices all the instruments contained in the above described yield curves to their exact original observed price; • it can be entirely described through a finite (albeit not a priori specified) number of rates and dates.
The instruments are sorted by ascending maturity.One constructs the curve by recurrence up to each maturity.Each of those maturities is termed a "knot point".Rates at intermediate points on the curve can be estimated by assuming a shape for the curve either in zero-coupon price or rate space.The choice of that interpolation rule constitutes the signature of the bootstrapping method.This choice is not conditioned by any theoretical reason, but by several practical reasons.The chosen rule should be such that: 1) The curve is smooth.
2)The curve does not have strong oscillations.
To those requirements, we add the supplementary one that: 3) The integral of the zero-coupon price between two knots can be computed in closed analytical form.The same holds for the integral of the zero-coupon price multiplied by the time to maturity.
The two first conditions are antagonist: an interpolation rule that favours one requirement will generally disfavour the other one.
The third requirement reflects the necessity of performing consol-related calculations, to be described in paragraph 2.3.2..We have tested four different rules, satisfying to those three criteria, among which the "unsmoothed Fama-Bliss" bootstrapping, whose interpolation rule results into stepwise constant forward rates.All of the methods performed close in term of the realistic aspect of the constructed curve.On the basis of some minute differences, we have chosen as default bootstrapping method one of the three other methods, namely, the one used in [9].

Extrapolation
It order to handle swaptions of maturity 5-year on 30-year swaps, one needs a yield curve covering a range of 35 years.Yet, quotes for the EONIA swaps stop at the 30-year tenor, at least for the price source that we have chosen to privilege.We have therefore extrapolated the EONIA curve.We have chosen to also extrapolate the EURIBOR swap curve to 35-year, to ensure the similarity of their treatment.The extrapolation of a curve has been achieved by adding the zero-coupon rate z(35) defined as 2 × z(30) − z(25).

Converted Volatilities
The primary input for the calculation of the consol volatility is a set of implied volatilities quoted in the market.The implied volatilities that we have used as raw data are those of EURIBOR swaptions, which are standardized options on 6 month EURIBOR swaps.They cannot be directly used, and this, for two reasons.
Firstly, the swaptions volatilities are implied by a Black and Scholes model in which the logarithm of the swap rate is assumed to be a Brownian motion.By contrast, while the consol volatility, resulting from the equations presented in the text, has to be implied by a Black and Scholes model.A conversion will thus be necessary, changing the raw swaptions volatilities into other ones implied by the second model.
Secondly, the swaptions volatilities pertain to EURIBOR-linked instruments, whereby the consol volatility pertains to the EONIA curve.
The change of model cannot be done by an exact calculation (or that exact calculation would be too complicated).However, as we already mentioned, we know that the motion of the empirical yield curves are primarily composed of parallel shifts.We then make an approximation and assume that those movements consist purely of parallel shifts.Thus, we only need to do the calculation at the first order, i.e. to multiply the raw swaptions volatility by the sensitivity of the log zero coupon price w.r.t. the log swap rate.
The two conversions are simultaneously achieved as follows: The underlying of the option-which is a forward EURIBOR swap-is priced from the EURIBOR curve.Then, one computes the quantity of parallel shift to apply to the EONIA curve to let it price the forward swap at its present market price.That quantity is called the "z-spread", we denote it with h.
The sensitivity of the log zero coupon price (of the EONIA curve) w.r.t. the log swap rate (non compounded, with day-count actual/360) is then given by: 100 365 log 100 1 36000 where CONSOL(t), is the consol rate observed at the close of business of t, continuously compounded, with a day count actual/365, and expressed in percentage points.The normalized excess return is the quotient of the excess return between 1 t and 2 t and of the consol volatility observed at the close of business of 1 t .

Table 1 .
Parameters and initial values of the 1-factor model simulation.

Table 2 .
Results of the 1-factor model simulation.

Table 3 .
Parameters and initial values of the 2-factor model simulation.

Table 4 .
Results of the 2-factor model simulation.

Table 5 .
Results from the empirical data.