Next Article in Journal
Measuring the Distance between Sets of ARMA Models
Next Article in Special Issue
Jump Variation Estimation with Noisy High Frequency Financial Data via Wavelets
Previous Article in Journal
Estimation of Gini Index within Pre-Specified Error Bound
Previous Article in Special Issue
Continuous and Jump Betas: Implications for Portfolio Diversification
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Market Microstructure Effects on Firm Default Risk Evaluation

by
Flavia Barsotti
1,* and
Simona Sanfelici
2
1
Risk Methodologies, Group Financial Risks, Group Risk Management, UniCredit Spa, Piazza Gae Aulenti, Tower A, Floor 20, Milano 20154, Italy
2
Department of Economics, University of Parma, Parma 43125, Italy
*
Author to whom correspondence should be addressed.
Econometrics 2016, 4(3), 31; https://doi.org/10.3390/econometrics4030031
Submission received: 27 February 2016 / Revised: 1 June 2016 / Accepted: 22 June 2016 / Published: 8 July 2016
(This article belongs to the Special Issue Financial High-Frequency Data)

Abstract

:
Default probability is a fundamental variable determining the credit worthiness of a firm and equity volatility estimation plays a key role in its evaluation. Assuming a structural credit risk modeling approach, we study the impact of choosing different non parametric equity volatility estimators on default probability evaluation, when market microstructure noise is considered. A general stochastic volatility framework with jumps for the underlying asset dynamics is defined inside a Merton-like structural model. To estimate the volatility risk component of a firm we use high-frequency equity data: market microstructure noise is introduced as a direct effect of observing noisy high-frequency equity prices. A Monte Carlo simulation analysis is conducted to (i) test the performance of alternative non-parametric equity volatility estimators in their capability of filtering out the microstructure noise and backing out the true unobservable asset volatility; (ii) study the effects of different non-parametric estimation techniques on default probability evaluation. The impact of the non-parametric volatility estimators on risk evaluation is not negligible: a sensitivity analysis defined for alternative values of the leverage parameter and average jumps size reveals that the characteristics of the dataset are crucial to determine which is the proper estimator to consider from a credit risk perspective.

1. Introduction

Many empirical works in the financial literature highlight the poor performance of structural credit risk models for defaultable bonds in predicting default probabilities and credit spreads, especially when considering very short maturities. The classical approach to firm’s value models is the one proposed by Merton [1]: a firm’s equity and debt are studied as contingent claims partitioning the total value of the firm; the underlying assets value dynamic is defined by a geometric Brownian motion. The empirical evidence indicates that classical standard structural models still have some difficulties in accurately explaining default rates and spreads simultaneously. Estimation or calibration methods provide evidence that predicted credit spreads are far below observed ones [2], the structural variables explain little of the credit spread variation [3], pricing error is large for corporate bonds [4]. From an econometric point of view, the calibration of these models on real data is non trivial due to key variables which are not directly observable: the underlying firm’s asset value and the volatility of firm’s asset value returns. This critical issue concerning the calibration of the underlying asset dynamics becomes even more problematic when one wants to consider also market microstructure effects into the structural modeling framework. The market microstructure literature strongly indicates that trading noises can affect equity prices so that the estimation of equity volatility and other related quantities may become a very difficult task. Observed equity prices can diverge from their equilibrium value due to illiquidity, asymmetric information, price discreteness and other measurement errors. Among many others, the works by [5,6] study the effects of trading noise on the optimal sampling frequency of observed equity prices since sparse sampling may reduce the bias in volatility measures due to such microstructure effects. The aim is to use as much as possible the information content of high-frequency noisy equity prices and define a procedure to back out the firm’s asset volatility.
In the specific context of structural credit risk models, the relationship between the unobservable asset volatility and the observed equity value predicted by the pricing model is masked by trading noise; ignoring microstructure effects could non-trivially inflate estimates for the “true” equity volatility, and this would lead to poor estimates for spreads and default rates. This issue has been analyzed in [7], where the authors extend [8] to explicitly account for trading noise contamination on equity prices: they devise a particle filter-based maximum likelihood method based solely on the time series of observed equity values, showing its robustness with respect to market microstructure effects. The importance of using high frequency data to back out parameters involved in firm’s value dynamics has been highlighted in [9]: the authors propose a novel approach to identify the volatility and jump risks components of individual firms from high-frequency observed equity prices, arguing that these new measures can allow to capture the effects of stochastic volatility and jumps in the firm’s underlying asset dynamics ([3,10]). Their analysis suggests that high-frequency-based volatility measures can help to better explain credit spreads, above and beyond what is already captured by the true leverage ratio. However, the highest frequency considered in this paper is the 5-min conservative sampling frequency which allows to filter microstructure effects. Recently, several non-parametric estimators of daily stock volatility have been proposed, allowing to exploit the information content of intra-day high-frequency data without being affected by microstructure effects [5,6,11,12,13,14,15].
In this paper, we carry on the analysis in [16] considering a structural approach à la Merton (see [1]) for defaultable corporate bonds in the presence of market microstructure noise and study a class of models for the underlying firm’s assets value dynamics allowing for: (i) stochastic asset return volatility; (ii) jump component. This extension to stochastic volatility and jump component for the underlying assets is deemed to be particularly significant from a credit risk modeling perspective: indeed, pricing and hedging issues related to equity markets suggest that stochastic volatility and jumps are crucial in order to improve credit spreads predictions. In addition to this, a negative leverage effect between stock prices and volatility is a typical feature observed in real data: stock prices tend to decrease when volatility rises.
The mathematical framework for a firm’s assets value process is defined through, alternatively, Heston [17] and Bates [18] models; as particular cases, the classical Merton model [1] and Merton model with jumps in the underlying asset dynamics [19] can also be derived. We propose a particular econometric approach to structural models calibration based on a non-parametric estimation of equity volatility from high-frequency intra-day equity prices. Market microstructure effect is introduced by assuming equity prices contaminated by trading noise under two different scenarios: (a) trading noise independent log-Gaussian distributed; (b) trading noise correlated with intra-day equity log-returns. From both an econometric and credit risk point of view, the interest of the analysis is concentrated on three key variables: equity volatility, underlying (unobservable) asset volatility, default probability. The aim of the paper is to exploit the information content of intra-day high-frequency prices and compare the performance of alternative non-parametric equity estimators in filtering out microstructure effects, retrieve the underlying asset volatility and thus evaluate the corresponding predicted default probabilities under the historical measure. The estimated underlying asset volatility is obtained as result of a specific calibration procedure matching the actual default probability (for a given maturity) and considering equity volatility estimation as input: this technique allows us to overcome the problem of its non-observability even when dealing with a structural approach. The results of our Monte Carlo simulation analysis highlight that asset volatility and default probabilities are deeply affected by the choice of the non-parametric estimator for equity volatility. From an econometric point of view, the commonly used 5-min Realized Volatility estimator is unable to provide reliable estimates for equity volatility in the presence of market microstructure noise, leading to a significant underestimation of both asset volatility and default probabilities and causing a wrong prediction of the real firm’s risk profile. Obviously, the Realized Volatility estimator based on noisy 1-min returns would provide even worse estimates, thus confirming results already discussed in [16]. For this reason, in this paper we discuss full results for the 5-min Realized Volatility estimator providing a more fair benchmark for all the other estimates. The Monte Carlo analysis has been developed for firms belonging to A rating class, revealing that the choice of the non-parametric volatility estimator has a strong influence on risk evaluation. This confirms the results of the analysis in [16] and extends their validity to the case of stochastic volatility models with jump components for the firm’s assets.
The paper also provides empirical evidence highlighting the relevance of this study from a practical credit risk perspective: the idea underlying the empirical analysis is to consider high frequency equity prices of a US company and the corresponding time series of credit default swap (CDS) premium as direct measure of a firm’s default probability (see also [9]). A regression analysis between daily CDS quotes and the corresponding volatility measures based on high-frequency intra-day equity prices is performed to show that alternative non parametric equity volatility measures can have a different explanatory power and thus a different impact for a firm default probability evaluation.
The paper is organized as follows: Section 2 describes the alternative non-parametric estimation techniques for equity volatility considered in our paper. Section 3 provides empirical evidence related to the role of equity volatility estimates based on high-frequency equity prices for default probability evaluation. Section 4 defines the structural credit risk model. Section 5 reports the results of our Monte Carlo analysis for backing out assets volatility and evaluating default probabilities. Section 6 analyzes the sensitivity of the asset volatility calibration and default probability estimation w.r.t. increasing levels of the leverage effect and of the jump component. Finally, Section 7 concludes.

2. Non-Parametric Volatility Estimation under Microstructure Noise

This section describes the approach we follow to estimate the equity volatility Σ t s from high-frequency equity data and gives details about the implementation of the specific volatility measures considered in our analysis. We set p ˜ t : = log S ˜ t , the noisy equity log-price. Time is measured in daily units. We build daily measure of volatility by considering daily windows of n intra-day equity data p ˜ t , j , j = 1 , , n . Besides the well known Realized Volatility estimator Σ t R V : = j = 1 n δ j ( p ˜ ) 2 , where δ j ( p ˜ ) : = p ˜ t , j - p ˜ t , j - 1 is the j-th within-day equity log-return on day t, we consider the following estimators of the volatility process Σ t s : the bias corrected estimator by Hansen and Lunde [12]
Σ t H L : = Σ t R V + 2 n n - 1 j = 1 n - 1 δ j ( p ˜ ) δ j + 1 ( p ˜ ) ;
the flat-top realized kernels by [6,20]
Σ t K : = h = - H H k h H + 1 j = | h | + 1 n δ j ( p ˜ ) δ j - | h | ( p ˜ ) .
with kernels of T H 2 type k ( x ) = sin 2 π 2 ( 1 - x ) 2 . The realized kernels may be considered as unbiased corrections of the Realized Volatility by means of the first H autocovariances of the returns. In particular, when H is selected to be zero the realized kernels become the Realized Volatility. Our analysis includes also the two-scale estimator by [15]
Σ t T S : = S S - 1 1 S s = 1 S Σ t G ( s ) - 1 S Σ t R V .
The two-scale (subsampling) estimator is a bias-adjusted average of lower frequency realized volatility computed on S non-overlapping observation sub-grids G ( s ) containing n S observations. Recently, [13] proposed a pre-averaging technique as an alternative to sub-sampling in order to reduce the microstructure effects. The idea is that if one averages a number of observed log-prices, one is closer to the latent process p ( t ) . This approach, when well implemented, gives rise to rate optimal estimators of power variations. In particular, a consistent estimator of the integrated volatility can be constructed as
Σ t P A = Δ θ ψ 2 s = 0 n - k n + 1 δ ¯ s ( p ˜ ) 2 - ψ 1 Δ 2 θ 2 ψ 2 s = 1 n δ s ( p ˜ ) 2 ,
where the pre-averaged return process is given by
δ ¯ s ( p ˜ ) : = r = 1 k n g r k n δ s + r ( p ˜ ) = 1 k n j = k n / 2 k n - 1 p ˜ t , s + j - j = 0 k n / 2 - 1 p ˜ t , s + j ,
θ = k n Δ , ψ 1 = 1 and ψ 2 = 1 / 12 , corresponding to the "hat" weight function g ( x ) = x ( 1 - x ) . The Fourier-Fejer estimator [14] is given by
Σ t F = ( 2 π ) 2 N + 1 s = - N N 1 - | s | N c s ( d p ˜ n ) c - s ( d p ˜ n ) ,
where c k ( d p ˜ n ) = 1 2 π i = 1 n exp ( - i k t i - 1 ) δ i ( p ˜ ) .
Finally, the jump-robust Bipower Variation [21] is defined as
Σ t B V = μ 1 - 2 j = 2 n | δ j - 1 ( p ˜ ) | | δ j ( p ˜ ) | ,
with μ 1 0 . 7979 .
Finite sample MSE-based optimal rules for choosing the parameters employed by these estimators are discussed in [15,22,23,24]. Here, we proceed according to the following rules: a simple approximation of the optimal sampling frequency for the Realized Volatility and Bipower Variation estimators is to choose the number of observations approximately equal to n * = ( Q / 4 E [ η 2 ] 2 ) 1 / 3 , where Q is the integrated quarticity estimated by means of low frequency returns. The optimal number of sub-grids S is given by c * n 2 / 3 , where c * = ( Q / 48 E [ η 2 ] 2 ) - 1 / 3 . For the Kernel estimator, we apply the optimal mean square error bandwidth selection suggested by [23] and get H = c * ξ 4 / 5 n 3 / 5 , where c * = ( 144 / 0 . 269 ) 1 / 5 , ξ 2 = E [ η 2 ] / Q . In the case of the Pre-averaging estimator, inspired by [23], we choose k n = c * ξ 4 / 5 n 3 / 5 . Finally, for the Fourier-Fejer estimator, the optimal cutting frequency N can be easily obtained by direct minimization of the estimated MSE given by Theorem 3 in [24].

3. Empirical Evidence

The aim of this section is to provide empirical evidence that alternative non parametric equity volatility estimates constructed from high-frequency equity return data may have different impact on a firm default probability evaluation. The idea is to consider high frequency equity prices of a US company and the corresponding time series of CDS quotes. We choose to use the CDS premium as a direct measure of a firm’s default probability as already done in [9]. The CDS is one of the most popular instrument in the credit derivatives market. Under a CDS contract, the protection seller promises to buy the reference bond at its par value when a predefined default event occurs. In return, the protection buyer makes periodic payments to the seller until the maturity date of the contract or until a credit event occurs. This periodic payment is called the CDS spread and provides a pure measure of the default risk of the reference entity. Indeed, if compared with corporate bond spreads, CDS spreads provide relatively ”pure” pricing of the default risk of the underlying entity and are typically traded on standardized terms. Moreover, bond spreads are more likely to be affected by differences in contractual arrangements, such as differences related to seniority, coupon rates, embedded options, guarantees and liquidity factors, which do not necessarily reflect the ”pure” default risk of the underlying asset. Finally, CDS spreads tend to respond more quickly to changes in credit conditions in the short run, which may be partly due to the absence of funding and short-sale restrictions in the derivatives market.
We adopt equity volatility measures as proxies for the time variation in volatility risk of a firm in order to examine the relationship between equity returns and credit spreads. The empirical results in [9] suggest that the volatility risk alone can predict up to 48% of the variation in CDS spread levels, particularly for investment-grade entities. Their analysis was based on a sample of US corporates over the period January 2001–December 2003.
Here we focus on 5-year CDS contracts for the Dow Jones Alcoa company. The sample covers the period from 10 November 2015 to 23 May 2016, for a total of 134 trading days. Our CDS quotes are daily and are provided by Bloomberg. The corresponding equity price sample contains 52,166 observations at a 1-min frequency for a 6.5 h daily trading period. Table 1 and Table 2 report, respectively, descriptive statistics of high frequency equity prices and CDS premium in our dataset.
High-frequency equity returns are contaminated by transaction costs, bid-ask bounce effects and possibly jumps leading to biases in the variance measures. Figure 1 shows the time plot of the 1-min log-returns and the autocorrelation function. Row data exhibit a strongly significant negative first order autocorrelation and higher order autocorrelations remain significant up to lag 12. Sparse sampling at 5 min frequency eliminates most of the microstructure effects. Figure 2 shows the histogram of the 1-min equity log-returns. We can see that the distribution of returns is far from being normal. Skewness is equal to –1.2078, while kurtosis is 73.1947.
Jumps have been identified and measured using the Threshold Bipower Variation method (TBV) of [25], which is based on the joint use of Bipower variation and threshold estimation. This method provides a powerful test for jump detection, which is employed at the significance level of 99.99%. This procedure allows to identify a total of 25 days with jumps which are removed from the sample. Since the TBV estimator is not robust to microstructure noise, we compute this measure of the integrated volatility in the presence of jumps using sparse sampling.
A first empirical evidence of the correlation between CDS spreads and equity volatility measures can be obtained by visually inspecting Figure 3, that shows the time plot of the daily 5-years CDS spreads in basis points (bottom panel) and the corresponding equity Realized Volatility (top panel) for our sample of 109 days after jump removal. The final sample covers the period 10 November 2015–23 May 2016. As we can observe from the figure, CDS spreads increased substantially in the middle of the considered period and gradually declined at the end of the sample period. The (annualized) daily return volatility varies between 0.2016 and 1.0838 following a similar pattern. From a statistical point view, the CDS spread and Realized Volatility series are correlated with positive correlation of 61 . 7 % that exactly reflects their common qualitative behavior.
In order to examine the explanatory power of the different equity volatility measures introduced in Section 2, for each estimated series of equity integrated volatility, we project the CDS spread on day t on a constant and the corresponding estimated daily volatility measure Σ t s . The Mincer-Zarnowitz evaluation regression takes the form
C D S t = ϕ 0 + ϕ 1 Σ t s + ϵ t ,
where t = 1 , , 109 . The R 2 from these regressions provides a direct assessment of the variability in the CDS premium that is explained by the particular estimates in the regressions. The R 2 can therefore be interpreted as a simple gauge of the explanatory power of the volatility estimate and hence of its potential economic significance.
Table 3 shows the OLS estimates from regressions (1). All the estimators are constructed using 1-min returns. The only exceptions are the Realized Volatility and the Bipower Variation, obtained from sparse sampled 5-min returns. All the estimated coefficients are significant at the 0.1% significance level. The signs of coefficients are correct, since a higher equity volatility estimate raises credit spread, and the magnitudes are sensible: for instance, a volatility shock of 1% in the Realized Volatility estimate raises the credit spread by 4%.
As we can observe from the results reported in Table 3, a significant portion of the CDS premium can be explained by the equity volatility estimates. This empirical evidence confirms the key role played by high-frequency equity prices in the evaluation of a firm default probability. The general outcome of our empirical study confirms the results already stated by [9] highlighting the key role of the information content of equity volatility and extends them from an econometric point of view, covering a bigger set of equity volatility estimators besides the 5-min Realized Volatility. Regardless of the specific method implemented to estimate the volatility, we show that the portion of CDS premium variability explained by equity volatility is between 36% and 42%. This smaller percentage w.r.t. the 48% found in [9] should be linked to the differences in market conditions between the two samples considered in the empirical case study.
The commonly used 5-min Realized Volatility yields a R 2 of 38%, while estimators robust to microstructure effects can yield higher R 2 using all available data. In particular, the Fourier estimator and the Two-scale estimator provide the highest R 2 results at 42% and 41%, respectively, with a gain of 10% over the Realized Volatility.

4. Firm Value Model in the Presence of Market Microstructure Noise

This Section describes our structural modeling framework and the model we use to describe trading noise on high-frequency equity data. We consider a stylized Merton-type [1] structural approach that incorporates stochastic volatility and jumps in the presence of market microstructure noise. Firm’s asset value process is defined through a jump-diffusion model, including Bates [18] and Heston [17] models as described in Section 4.1. Market microstructure is considered in our setting through the assumption of equity prices being contaminated by trading noise under two different scenarios: a) independent log-Gaussian distributed noise process, b) trading noise process correlated with intra-day equity log-returns; Section 4.2 gives a full description of market microstructure.

4.1. Assets Value Dynamics and Default Assumptions

We assume the structural credit risk modeling approach defined by Merton [1] and follow the intuition of [3] introducing stochastic volatility and jumps in the underlying firm’s asset value dynamics as
d A t A t = ( μ - δ - λ μ J ) d t + V t d W t 1 + J t d q t d V t = k ( θ - V t ) d t + σ V t d W t 2 ,
where A t is firm value, μ is the instantaneous asset return, δ the asset payout ratio. The asset jump component has a Poisson mixing Gaussian distribution, where q t is a Poisson process with intensity λ and log ( 1 + J t ) are Normally distributed jump sizes in the logarithmic firm value, with expectation log ( 1 + μ J ) - σ J 2 / 2 and variance σ J 2 . The Poisson process q t is independent of the Brownian motions and the jump sizes.The asset return variance, V t , follows a square-root process with long-run mean θ, mean reversion speed k and volatility-of-volatility parameter σ. Finally, ρ is the correlation between Brownian motions W t 1 , W t 2 . This dynamics has been studied in the equity-option price literature ([18,26]), while the general credit risk modeling framework has been analyzed by [10]. The model (2) includes the Heston model [17], Merton setting with a jump component [19] and also the original Merton framework [1] by assuming V t constant through time (i.e., k , θ , σ = 0 ). The firm has two classes of outstanding claims: equity and a zero-coupon debt with promised payment B at maturity T. To be suitable for pricing corporate debt, we adopt the following bankruptcy assumptions from the classical setting [1]: (i) there is only one outstanding bond, i.e., a zero-coupon bond with a promised payment B at maturity; (ii) default can occur only at bond maturity with debt face value as default boundary; and (iii) when default occurs, the absolute priority rule prevails. The payoffs to debt holders and equity holders at time T become, respectively
D T = min ( A T , B ) , S T = max ( A T - B , 0 ) .
From now on, the focus of our study is concentrated on equity value and default probabilities in order to develop our computational econometric analysis. Using no-arbitrage arguments, equity price, namely S t , can be computed as the price of a European call option written on the underlying asset A t with maturity T
S t = A t F 1 * - B e - r ( T - t ) F 2 * ,
where F i * = F i * ( A t , V t ; T , B ) i = 1 , 2 are the so called risk-neutral probabilities, computed using the corresponding characteristic function (see [17,18]).
We apply Itô’s lemma and get the following specification for equity price dynamics
d S t S t = 1 S t μ t ( · ) d t + A t S t S t A t V t d W t 1 + 1 S t S t V t σ V t d W t 2 + 1 S t [ S t ( A t ( 1 + J t ) , V t ; Ω ) - S t ( A t , V t ; Ω ) ] d q t ,
where μ t ( · ) is the instantaneous equity return under the historical probability measure, Ω is the parameter vector. Let Σ t s indicate the instantaneous volatility of log-equity price given by
Σ t s = A t S t 2 S t A t 2 V t + σ S t 2 S t V t 2 V t + 2 ρ σ A t S t 2 S t A t S t V t V t .
Notice that the evolution of equity volatility is driven by two time-varying factors, A t and V t , while the underlying asset volatility is only driven by V t . As particular case, if asset volatility is constant ( v : = V ), then Equation (4) reduces to the standard Merton (see [1]) formula Σ t s = A t S t S t A t v . The jump size of log-equity price is given by
J t s = log [ S t ( A t ( 1 + J t ) , V t ; Ω ) ] - log [ S t ( A t , V t ; Ω ) ] .
The firm’s probability of default at maturity T is the probability of A T being below the constant barrier represented by the face value of debt B. Under the physical probability measure this can be computed through
P ( A T < B | A t ) = 1 - P ( A T B | A t ) = 1 - F 2 ( A t , V t ; T , B ) ,
where F 2 ( A t , V t ; T , B ) will be evaluated with the historical drift for asset returns.

4.2. Market Microstructure Noise on High-Frequency Equity Prices

Let us consider an individual firm: we can obtain a time series of equity prices { S t j , j = 1 , , n } , with a given constant sampling interval h = t j - t j - 1 . If one could observe the efficient (i.e., not contaminated by noise) equity price, than equity volatility could be estimated by the well known Realized Volatility estimator [27] at any desired accuracy level using high frequency data. However, as stressed by [7], the relationship between the unobserved asset and the observed equity value predicted by the pricing formula (3) may be masked by trading noise in reality. Econometric literature suggests that observed equity prices can diverge from their equilibrium value due to illiquidity, asymmetric information, price discreteness and other measurement errors. Thus, the approach used to estimate equity volatility must take into account the presence of trading noise affecting the sample of observations.
In this paper we assume an additive error structure for the trading noise on the logarithmic equity value as follows
log S ˜ t j = log S t j + η t j ,
where the random shocks η t j , for 0 j n are i.i.d. random variables with zero mean and bounded fourth moment and independent from the efficient log-return process. The assumption of independence can be relaxed by considering a particular form of dependent noise, given by [12], with market microstructure noise that is time-dependent in tick time and correlated with efficient returns
η ˜ j : = α [ log S t j - log S t j - 1 ] + η j ,
where α is a real constant and η ˜ j and η j are the shorten notation for η ˜ t j and η t j . The case α = 0 corresponds to the case of independent noise assumption. The basic idea of our paper is that, using suitable volatility estimators, we can infer the volatility process Σ t s of equity returns from noisy high-frequency data. Then, equity volatility estimate can be used to back out the (unobservable) asset volatility v t = V t in order to match exactly the 5-years actual default probability given in Equation (5). The resulting calibrated asset volatility finally allows to compute default probabilities for different maturities.

5. Equity Volatility Estimation and Default Probability Computation

This section describes in detail our Monte Carlo simulation analysis and gives numerical results about the performance of different equity volatility estimators in their capability of: (i) filtering the microstructure noise and extracting equity volatility correctly; (ii) backing out asset volatility and (iii) predicting default probabilities.
We perform Monte Carlo simulations by generating the underlying asset dynamic according to Heston and Bates models, respectively, taking calibrated parameters from [9] for the base case. The full set of parameters values used in our study is contained in Table 4.
In order to avoid arbitrage opportunities, all securities written on the underlying firm value A t must have the same Sharpe ratio, see [1] Equation (6). The same underlying assumption on the Sharpe ratio is made for the asset volatility calibration procedure described in [9]. Therefore, the following equation must be satisfied:
π a V t = π e Σ t s ,
where π a : = μ - r + δ is the total asset risk premium, π e is the historical equity risk premium and we replace Σ t s by the value Σ ^ t , that is any of the equity volatility estimates described in Section 2. Equation (8) is the technical tool we use in order to back out the asset volatility v t : = V t that exactly fits the 5-years probability of default given in Equation (5) according to the given equity volatility estimate. More precisely, we plug relation (8) into (5) and solve for v ^ t = V t . This step enables us to tie down the risk-neutral distribution with the objective distribution without resorting to a certain arbitrary assumption or noisy estimate of the asset return process. Once noisy equity prices are generated, we compare the performance of different volatility estimators in their ability of inferring the ”true” equity volatility Σ t s . Then, we analyze the influence of this results on the underlying asset volatility calibration and default probability computation.
In order to understand the influence of alternative equity volatility estimators on both asset volatility calibration and default probability computation, we preliminarily show in Figure 4 the theoretical relationship between these quantities. The idea is to take as base case the pair of parameters Σ t s = 37 . 21 % , v t = 25 . 01 % and evaluate the impact of different equity volatility values Σ t E [ 0 . 3 , 0 . 45 ] through our calibration device on both calibrated asset volatility and default probability, as if Σ t E was the equity volatility estimate obtained through the non-parametric estimator E. The theoretical dependence of the calibrated asset volatility v t E and the default probability is shown in Figure 4 for the Bates model, but the qualitative behaviour is the same also for the Heston model: as we can see, the equity volatility estimate is a key variable determining both the calibration of asset volatility and the default probability. The impact on default probability is measured in terms of the absolute error defined as D P A b s E r r : = D P ( Σ t E ) - D P ( Σ t s ) . The calibrated asset volatility depicted in Subplot (i) is obtained by matching the 5-years default probability P D 5 = 1 . 37 % , then this value is used to compute the absolute error on the 3-years default probability reported in Subplot (ii). To generalize the outcome of the analysis, the relationship between calibrated asset volatility and equity volatility estimation is nonlinear (convex) monotonous decreasing: overestimating equity volatility generates an underestimation of the calibrated asset volatility. The same monotonous decreasing relation holds between Σ t E and the absolute error D P A b s E r r , but the absolute magnitude of the error changes may depend on whether we are under-estimating or over-estimating Σ t s and which time horizon we are considering. As a general feature, the influence of a non-correct equity volatility estimation can be not negligible. Thus, from a credit risk perspective, a robust equity volatility estimator is crucial to properly evaluate default probability at each maturity.

5.1. Equity Volatility Estimation with High-Frequency Data

We perform a Monte Carlo simulation by generating a sample of 1000 trading days and equispaced intra-day data at the frequency h = 1 min (i.e., n = 1440 ). Once the underlying asset value dynamics is generated by model (2), high-frequency equity prices are obtained through the no arbitrage method given by Equation (3). Market microstructure noise is considered, alternatively, for both cases described by Equations (6) and (7). The random shocks η j are i.i.d Gaussian random variables with zero mean and standard deviation equal to 1 . 4 times the log-equity return standard deviation. We set α = 0 . 5 in the dependent noise case (7).
In the case of Bates model, the generated dataset is cleaned for jumps according to the jump removal procedure sketched in Section 3. This procedure allows to eliminate around 7.5%–10% of days from our sample. Nevertheless, some days in the cleaned sample may still exhibit small jumps that can affect volatility estimation. In particular, days in our sample can present different statistical autocorrelation structures. Figure 5 shows the intraday log-returns, normalized to have mean zero and standard deviation one, for the model with Gaussian uncorrelated noise. Panel (A) refers to a day with no jumps and first order autocorrelation structure; panel (B) shows a big spike corresponding to a jump in the equity log-price process, still retaining a first order autocorrelation structure; (C) refers to a day with no jumps and autocorrelation at lags 3, 5, 9 and 15.
Next we plot the histograms of the intraday returns of equity log-prices. Figure 6 displays the histogram along with the standard normal density function, which is essentially confined within ( - 3 , 3 ) . Clearly, the histograms display a slightly high peak and asymmetric heavy tails. The distribution in case (B) is leptokurtic, i.e., more outlier-prone than the normal distribution. Moreover, skewness is positive, i.e., the data are spread out more to the right of the mean than to the left and the returns have heavier right tail.
Table 5 and Table 6 present numerical evidence for each equity volatility estimator introduced in Section 2 and used in our comparison when (a) trading noise is independent of intra-day equity log-returns and (b) trading noise is correlated with intra-day equity log-returns, respectively; both Heston [17] and Bates [18] dynamics are considered. All the estimators are applied on the equity log-price series after jump removal. The only exception is for the jump-robust Bipower Variation which is applied to the original dataset including jumps after optimal sparse sampling.
Table 5 and Table 6 list the mean squared error (MSE) and bias (BIAS) achieved by the different non-parametric equity volatility estimators. Σ t R V represents the sparse sampled Realized Volatility estimator based on 5-min returns, while Σ t R V S S refers to the sparse sampled Realized Volatility estimator where the sampling frequency is optimized in order to filter the microstructure effects, as explained in Section 2. Our results strongly confirm well known stylized facts documented by the econometric literature and highlighted by [24]: Σ t R V estimates are strongly affected by noise and sparse sampling can only moderately provide efficient estimates. The 5-min Realized Volatility estimator has the worst performance in terms of both MSE and BIAS among all estimators for Heston and Bates settings, with a MSE of order e - 02 and BIAS of order e - 01 . This performance can be even worst when considering the 1-min Realized Volatility estimator: results are not reported in the table, but equity volatility estimates are completely swamped by trading noise in such a case. The performance of the Bipower Variation estimator Σ t B V for equity volatility is comparable to the one obtained through the Realized Volatility estimator based on sparse sampling Σ t R V S S . This is particularly relevant in the Bates model simulation of Table 6. In fact, jumps have been removed from our samples using the procedure sketched in Section 2. Although some of the jumps may still be present after jump removal, this procedure seems to be effective in eliminating the most relevant jumps.
As an alternative to sparse sampling, the first order correction of Σ t H L can reduce the bias due to the spurious first order autocorrelation in equity returns introduced by the trading noise: overall, this estimator performs much better than Σ t R V S S in both models and for both assumptions on the trading noise. However, the BIAS reduction is more effective in the Heston model.
The best results in terms of MSE are obtained when considering the set of estimators specifically designed to handle microstructure effects, i.e., Σ t T S , Σ t K , Σ t P A and Σ t F .
If we want to try to rank the estimators according to the noise provision, we can state that both the Two-Scale Σ t T S and the Kernel estimators Σ t K have the best performance in Heston setting, with Σ t K having a smaller BIAS. The performance of the pre-averaging estimator Σ t P A and the Fourier estimator Σ t F are comparable in both models.

5.2. Asset Volatility Calibration

This section analyzes how the choice of different volatility estimators affects the calibration of the underlying asset volatility. In order to study the influence of different equity volatility estimators on the default probability predicted by Merton model (see [1]), we proceed by developing the following calibration exercise for the underlying asset volatility. For each equity volatility estimator, generically denoted by E, and each day t in our sample, we find the corresponding daily asset volatility estimate v t ^ E by matching the 5-years default probability coming from our equity volatility estimate Σ t E with the one evaluated through the model using the parameter values of Table 4. This is done for each day in the sample. In this calibration exercise we act as if we did not know the asset values to mimic the real-life estimation situation and we conduct inference only based on observable quantities such as measurable equity volatility and the 5-years default probabilities. Default probabilities are computed, for any maturity, according to Equation (5). In order to avoid arbitrage opportunities, following [9], we consider as key assumption that all securities written on the underlying firm value A t must have the same Sharpe ratio, see (8) (cfr. Equation (6) in [1]). This consideration enables us to express the instantaneous asset return μ as a function of the unknown asset volatility, given each equity volatility estimate Σ t E , and then to solve Equation (5) for the 5-years default probability with respect to the asset volatility, to obtain the corresponding asset volatility estimate v t ^ E . Once asset volatility is known, we can compute default probabilities for any other maturity according to Equation (5).
Table 7 shows descriptive statistics of calibrated asset volatility for the Heston model, Table 8 contains statistics for the Bates case. We report results obtained by matching 5-years default probabilities for each equity volatility estimator E. Panel (a) of both tables refers to a trading noise of the form (6); panel (b) refers to results obtained for a trading noise of the form (7). The set of equity volatility estimators represented by Two-scale, Kernel, Fourier and Pre-Average estimators provides the most accurate estimation of the resulting calibrated daily asset volatility for Heston case, showing the smallest standard deviation. The same behavior is confirmed by results under Bates setting in Table 8. These results are confirmed by statistics for the ratio v t ^ E / v t , where the variable v t is the underlying asset volatility calibrated by considering the model daily equity volatility Σ t s on each day t in the Monte Carlo sample. In particular, it is evident how the high-frequency Realized Volatility procedure largely underestimates the underlying asset volatility: v t ^ R V / v t is strongly biased due to microstructure effects, while the optimized v t ^ R V S S is less biased, even if a slightly larger variance appears.

5.3. Default Probability Computation

Once asset volatility is calibrated, we can turn to default probability computation. The average (over the daily Monte Carlo replications) default probabilities for different maturities obtained from alternative equity volatility estimation procedures are plotted in Figure 7 for Heston case, Gaussian noise. For each day in the sample and for each volatility estimator E , we use the calibrated asset volatility v t ^ E in order to compute the default probabilities for any maturity (from 1 to 5 years) through Equation (5). The figure highlights how the Realized Volatility approach based on high frequency noisy data drastically underestimates default probabilities, thus sparse sampling becomes mandatory when equity data are affected by microstructure effects. On the whole, all the other procedures seem to provide sensible results and only a deeper analysis reveals differences among different specific estimators, as shown in Table 9 and Table 10 where relative errors are reported.
Comparing the performance of alternative equity volatility estimators on default probabilities is done for different maturities (from 1 to 5 years) through the mean relative error between the estimated default probability and the theoretical one. For each maturity, we consider the following measure>
D P E r r E : = E D P ( Σ t E ) - D P ( Σ t s ) D P ( Σ t s ) · 100 ,
where D P ( Σ t E ) is the default probability evaluated through Equation (5) when equity and asset volatility are estimated through a generic couple of estimators Σ t E , v t ^ E ; D P ( Σ t s ) is the corresponding theoretical default probability when equity and asset volatility are, respectively, Σ t s and v t (model daily volatilities). Table 9 and Table 10 summarize the corresponding numerical results for both Heston and Bates cases, suggesting that the choice of the volatility estimator largely affects default probabilities.
Table 9 highlights that for the Heston model the smallest overall error on default probabilities is provided by the Fourier estimator, followed by the Two-Scales and the Kernel ones. For maturities greater than 2.5 years, even the Σ t H L and Pre-Average estimators give good results. The optimally sampled Bipower Variation Σ t B V provides results that are comparable with Σ t R V S S , underestimating the default risk. The classical 5-min Realized Volatility approach can severely underestimate risk, leading to an error of more than 90 % in absolute value for maturities before 1 y and around 60%−70% in absolute value for maturities up to 2 years.
A similar ranking among estimators holds in Bates case, as reported in Table 10. While for maturities greater than 3 years the Kernel estimator has the best performance, i.e., the smallest average error on default probabilities, overall the Fourier estimator reveals to be the more stable one, with relative errors across all maturities being below 5 . 6 % in absolute value. Differently from the Fourier case, all other estimators show a peak in the default probability relative error for the 1 y time horizon.
For a better visual insight, Figure 8 and Figure 9 show the mean relative error comparing different calibration procedures for Bates model, with, alternatively, Gaussian (a) and Correlated noise (b). A negative (positive) error reveals that the calibration procedure underestimates (overestimates) default probabilities. What emerges from the reported plots of default probability relative errors is that the Fourier estimator (red line in the Figures) is the most efficient, although it slightly underestimates default probabilities. This is due to the convexity of the non linear relationship between default probability error and asset volatility, as depicted in Figure 4 and discussed in Section 6. All the other estimators produce relative errors on default probabilities that are more inaccurate (and more volatile) as function of the maturity, reaching high absolute values for very short time horizons. The Realized Volatility estimator has a poor performance from 1y to 5y, always providing a significant underestimation of the default probability. This qualitative behavior reflects also results obtained for Heston setting, for both cases of Gaussian and Correlated trading noise and we omit them for brevity.

6. Sensitivity Analysis

In this Section we study the sensitivity of the asset volatility calibration and default probability estimation w.r.t. different values of the key parameters considered in our stochastic volatility setting: the correlation between volatility and returns (leverage effect) in Heston model, captured by parameter ρ, and the average percentage jump size in Bates model, captured by parameter μ J .
We consider as base case the setting described in the above analysis and define alternative sets of parameters. This is done by keeping unchanged the 5-years default probability P D 5 y = 1 . 37 % . The leverage effect is analyzed within Heston model in Section 6.1, while the sensitivity w.r.t. the jump parameter in Bates model is studied in Section 6.2. In both cases we consider a sample of 1000 days for the Monte Carlo simulation and a market microstructure effect described by a Correlated noise given in Equation (7). The simulations are based on equispaced intra-day data generated on a 1-min frequency. The analysis is conducted by exploring how the performance of different equity volatility estimators changes for different levels of the leverage effect and jump amplitude; then, the impact of equity volatility estimates on calibrated asset volatility and default probability is studied. From a credit risk perspective, this analysis is useful since it can suggest which equity volatility estimator is more appropriate depending on the specificity of the real dataset. The base cases presented below are the same already discussed in the previous section.

6.1. Leverage Effect

In this subsection we develop a sensitivity analysis on Heston model by considering alternative values of the correlation parameter ρ which captures the leverage effect. The idea is to quantify the impact of this key variable on both asset volatility calibration and default probability estimation. Table 11 summarizes the results of this analysis: starting from the base case described in the previous section, we compare the results associated to a correlation ρ = - 24 . 02 % with those obtained (i) for a higher (in absolute value) negative correlation ρ = - 30 % and (ii) for the uncorrelated case ρ = 0 . As well documented in the econometric literature, the leverage effect is known to be negative in financial markets, that is why we consider ρ = 0 as extreme scenario. The initial asset volatility considered is the one corresponding to a 5-years default probability P D 5 y = 1 . 37 % (see Table 11).
Table 12 reports MSE and BIAS of each equity volatility estimator under different levels of ρ. The outcome of our analysis highlights that the performance of equity volatility estimators is not much affected by the leverage parameter. Changing parameter ρ does not seem to generate a strong and significant directional impact on MSE and BIAS of the considered equity volatility estimators. Nevertheless, Table 12 reveals that the volatility estimators improve slightly their performance in terms of MSE as the negative leverage effect increases in absolute value; the Two-Scales, Kernel, Pre-Average and Fourier estimators give the best performance in terms of MSE, with an order of magnitude of e - 04 in all cases.
As highlighted by Table 13, the leverage effect impacts the underlying asset volatility calibration: increasing the correlation level enhances the variance of the calibrated asset volatility, as the last column of the table shows. The 5-min Realized volatility estimator worsens its performance when ρ ( < 0 ) increases in absolute value, reaching an average of asset volatility ratio of 71 % . The same happens for the Bipower variation estimator and the Realized volatility estimator with sparse sampling: B V worsen its asset volatility mean ratio from 95 . 2 % to 94 . 1 % going from ρ = 0 to ρ = - 30 % , while the R V S S estimator shows results from 94 % to 92 . 8 % in the corresponding case studies. Results are even worst if we consider the 1-min Realized Volatility estimator, completely swamped by trading noise, and reaching less than 40 % in terms of mean asset volatility ratio. We recall that results for the 1-min Realized Volatility estimator are not reported to have a more consistent benchmark with all other estimators. In terms of standard deviation, the performance of Σ t T S , Σ t H L , Σ t K , Σ t P A , Σ t F is comparable, going from 4 . 9 % to 6 . 8 % depending on the level of the leverage ratio.
Table 14 reports the results for the default probability errors under the three alternative settings associated to different leverage effects: all estimators (except the Fourier one) perform better in the uncorrelated case and get worse as the negative leverage effect increases in absolute value. From a credit risk perspective, this sensitivity analysis suggests that when dealing with a sample characterized by a correlation parameter ρ close to zero, the Kernel estimator should be preferred in order to evaluate the default probability, by reaching PD errors lower than 1 . 2 % for 3 y and lower than 0 . 16 % for 4 . 5 y maturity. The Σ t T S , Σ t H L , Σ t P A , Σ t F also provide good estimated for maturities greater than 3 . 5 y ; for shorter maturities, the smallest errors are obtained by the Fourier estimator with an error of - 1 . 75 % for 1.5y and a maximum of 2 . 8 % for the 1y time horizon, significantly lower than errors obtained for the same short term maturity with all other estimators. At the same time, the Two-Scale, Kernel and Fourier estimators represent the best choice for a sample of data characterized by a more realistic high negative leverage effect: here the extreme case considered is ρ = - 30 % , and Σ t T S , Σ t K , Σ t F show a good performance. For ρ = - 30 % , the Fourier estimator still has the smallest volatility of default probability errors if compared with the performance of all other estimators.
By looking at asset volatility and equity volatility estimation, we can observe what follows: the performance of the three noise-robust estimators Σ t T S , Σ t K , Σ t F is comparable, and in some cases the Two-Scale and Kernel estimators can even do better than the Fourier one. Nevertheless, in terms of default probability relative error, the Fourier estimator provides better results, in particular for short maturities. The intuition behind the observed behavior is linked to the convexity of the non-linear relationship between default probability error and underlying asset volatility, as depicted in Figure 4. As we can observe, the relationship between default probability error and equity volatility estimates is non-linear and flattening as the equity volatility estimate increases. The same applies for the relationship between asset volatility and equity volatility. While the Two-Scale and Kernel estimators bring to an average overestimation of the underlying asset volatility, the Fourier estimates reveals an overall underestimation of the asset volatility. As a consequence, Fourier results are associated to a slight underestimation of the default probability while the Two-Scale and Kernel estimators can bring to higher default probability overestimation. This effect is amplified for short maturities, since the target variable to calibrate is the 5-years default probability.

6.2. Average Jumps Size

This section compares the performance of alternative non-parametric estimators for different values of the average jump size, measuring their sensitivity to parameter μ J . The aim of this analysis is to study how alternative scenarios with a different jump amplitude can affect asset volatility calibration and default probability estimation inside a setting described by the Bates model.
Table 15 summarizes the cases we consider in this analysis: starting from the base parameters associated to μ J = 2 . 35 % (as in the previous Section, Table 4), we consider two alternative cases: one with higher average jump size μ J = 3 . 35 % , one with lower average jump amplitude μ J = 1 . 35 % . The initial asset volatility considered in each case is the one associated to a 5-years default probability P D 5 y = 1 . 37 % (see Table 15). Here we consider a sample of 1000 days for our Monte Carlo simulation and a Correlated trading noise as defined in Equation (7). The simulation is again based on equispaced intra-day data generated on a 1-min frequency.
Results for the equity volatility estimators are reported in Table 16: the Two-Scale estimator is the one having the smallest BIAS in the three cases analyzed, with a performance comparable to the Kernel estimator for the case μ J = 3 . 35 . In case μ J = 2 . 35 % or μ J = 1 . 35 % , the Two-Scale estimator still exhibits an order of magnitude on the BIAS of e - 04 , while the BIAS for the Kernel goes to e - 03 . A subset of estimators, namely Σ t B V , Σ t R V S S , Σ t K , Σ t P A , Σ t F , shows an improved performance in terms of MSE when going from μ J = 3 . 35 % to μ J = 1 . 35 % , thus for a decreasing average jumps size.
The jumps amplitude effect on asset volatility calibration is shown in Table 17: the general outcome is that the average jumps size has a quite negligible effect on asset volatility ratio for all estimators except the Bipower variation, showing the highest standard deviation of more than 15 % . The higher jumps amplitude μ J = 3 . 35 % is associated to an average asset volatility ratio for the Bipower variation of 93 % with a 15 . 44 % standard deviation. The Σ t T S , Σ t H L , Σ t K , Σ t P A , Σ t F estimators show the best performance: from our results, we cannot say that the jumps amplitude has a strong directional effect on asset volatility calibration. The 5-min Realized volatility estimator has the worst performance among the whole set of estimators, for each value of the jumps amplitude, showing an average asset volatility ratio of around 70 % . The performance is even worst if we consider the 1-min Realized volatility, providing estimates below 50 % in terms of mean asset volatility ratio for the Bates case. A possible explanation of the observed behavior for all estimators under Bates setting can be related to the fact that we are filtering out days with jumps: if the procedure to remove jumps gives good results for each value of the parameter μ J , this can explain why we do not see a huge impact on calibrated asset volatility. On the contrary, the Bipower variation estimates are computed on the full sample before jump removal and therefore are more sensitive to the jump size.
Finally, from Table 18 we can see the sensitivity of errors on default probability for different levels of the parameter μ J : the 5-min Realized volatility estimator still has a poor performance for each value of μ j , since it always underestimates default probability with errors reaching more than - 91 % in all cases for short maturities. The Realized volatility estimator with sparse sampling on the contrary improves its performance when μ J increases: this effect is greater for maturities up to 2.5 years.
The Fourier estimator gives the best results for the high average jumps size level μ J = 3 . 35 % . Under this scenario, also the Kernel and Two-Scale estimators perform quite well from 2.5 years onwards, with errors below 2 . 7 % .
From a credit risk perspective, the choice of the volatility estimator is affected by the level of the jumps amplitude: for intermediate and low levels of the average jump amplitude, the best performance in terms of relative errors is given by the Kernel estimator for maturities greater than 2.5 years, while for short maturities the Fourier one has to be preferred. When the dataset is instead characterized by a high average jump size, i.e., μ j = 3 . 35 % in our test case, the Fourier estimator provides the smallest standard deviation of errors on the default probability estimation, thus it seems to be the best choice from a risk management perspective under this scenario. If we consider single time horizons estimates, the Two-Scale and Kernel can have estimates which are better than Fourier, but the overall performance is dominated by the Fourier estimator with errors always below 3 % for all maturities.

7. Conclusions

In this paper we address the problem of calibrating the underlying asset volatility and estimating the corresponding default probability inside a structural credit risk model. This is done by means of a non-parametric estimation of equity volatility coming from high-frequency intra-day equity prices.
We consider a stylized Merton-type [1] structural model for defaultable bonds that incorporates stochastic volatility and jumps in the presence of market microstructure noise. The model framework for the unlevered firm’s assets value process is defined through, alternatively, Bates [18] and Heston [17] models. Two market microstructure assumptions are considered, namely independent Gaussian and Correlated (to equity returns) trading noise.
We show that by exploiting the information content of intra-day high-frequency prices and by filtering out microstructure effects, it is possible to efficiently retrieve the underlying asset volatility and thus evaluate the corresponding default probabilities under the historical measure. The simulation analysis is conducted for bonds with rating A and the calibrated underlying assets volatility is obtained as result of a specific calibration procedure defined by matching the 5 y default probability for this rating class and considering equity volatility estimation as input. This technique allows us to overcome the problem of the non-observability of the underlying asset volatility even if working in a structural model.
A regression analysis conducted between CDS daily quotes of a US company and the corresponding series of equity volatility measure extracted from high frequency equity prices highlights that alternative non parametric equity volatility estimators can have a different impact in terms of explanatory power of CDS premium variability and thus on default probability evaluation. These empirical evidences are confirmed by the simulation study under both Heston and Bates settings.
The outcome of our Monte Carlo simulation analysis highlights that asset volatility and default probabilities are deeply affected by the choice of the non-parametric equity volatility estimator. The commonly used Realized Volatility estimator is unable to provide reliable estimates for equity volatility in the presence of market microstructure noise, leading to a significant underestimation of both asset volatility and default probabilities. This confirms the results obtained by [16] and extend their validity to the case of stochastic volatility models with jump components for the firm’s assets.
A sensitivity analysis is also provided to evaluate to what extent the choice of the non-parametric estimators for equity volatility affects the calibrated asset volatility and default probabilities when data are characterized by alternative leverage parameters (Heston model) or alternative average jumps size (Bates model). The analysis reveals interesting results from a credit risk point of view, suggesting that (i) the characteristics of the dataset are crucial to determine which is the proper estimator to consider and (ii) the convexity of the non-linear relationship between default probability error and underlying asset volatility plays a key role. As an example, when the dataset is characterized by a high negative leverage effect, the Fourier estimator has a performance which is comparable with other estimators in retrieving the underlying asset volatility but has the best performance in terms of default probability, especially for short maturities, due to the non-linear relationship between default probability error and underlying asset volatility.

Acknowledgments

The authors would like to thank the two anonymous referees for their helpful suggestions and comments.

Author Contributions

The two authors contributed equally to the work.

Conflicts of Interest

The authors declare no conflict of interest. The corresponding author, Flavia Barsotti, declares that the views expressed in this paper are those of the author and should not be attributed to UniCredit Group or to the author as representative or employee of UniCredit Group.

References

  1. R.C. Merton. “On the pricing of corporate debt: the risk structure of interest rates.” J. Finance 29 (1974): 449–470. [Google Scholar] [CrossRef]
  2. P.E. Jones, P.M. Scott, and E. Rosenfeld. “Contingent claims analysis of corporate capital structures: An empirical investigation.” J. Finance 39 (1984): 611–625. [Google Scholar] [CrossRef]
  3. J. Huang, and M. Huang. “How much of the corporate-treasury yield spread is due to credit risk? ” Rev. Asset Pricing Stud. 2 (2012): 153–202. [Google Scholar] [CrossRef]
  4. Y.H. Eom, J. Helwege, and J. Huang. “Structural models of corporate bond pricing: An empirical analysis.” Rev. Financ. Stud. 17 (2008): 499–544. [Google Scholar] [CrossRef]
  5. F.M. Bandi, and J.R. Russel. “Separating market microstructure noise from volatility.” J. Financ. Econ. 79 (2006): 655–692. [Google Scholar] [CrossRef]
  6. O.E. Barndorff-Nielsen, P.R. Hansen, A. Lunde, and N. Shephard. “Designing realised kernels to measure the ex-post variation of equity prices in the presence of noise.” Econometrica 76 (2008): 1481–1536. [Google Scholar] [CrossRef]
  7. J.C. Duan, and A. Fulop. “Estimating the structural credit risk model when equity prices are contaminated by trading noises.” J. Econom. 150 (2009): 288–296. [Google Scholar] [CrossRef]
  8. J.C. Duan. “Maximum likelihood estimation using price data of the derivative contract.” Math. Finance 4 (1994): 155–167. [Google Scholar] [CrossRef]
  9. B.Y. Zhang, H. Zhou, and H. Zhu. “Explaining Credit Default Swap Spreads with the Equity Volatility and Jump Risks of Individual Firms.” Rev. Financ. Stud. 22 (2009): 5099–5131. [Google Scholar] [CrossRef]
  10. J. Huang. Affine Structural Models of Corporate Bond Pricing. Working Paper; State College, PA, USA: Penn State University, 2005. [Google Scholar]
  11. Y. Aït-Sahalia, P. Mykland, and L. Zhang. “How often to sample a continuous-time process in the presence of market microstructure noise.” Rev. Financ. Stud. 18 (2005): 351–416. [Google Scholar] [CrossRef]
  12. P.R. Hansen, and A. Lunde. “Realized variance and market microstructure noise (with discussions).” J. Bus. Econ. Stat. 24 (2006): 127–218. [Google Scholar] [CrossRef]
  13. J. Jacod, Y. Li, P. Mykland, M. Podolskij, and M. Vetter. “Microstructure noise in the continuous case: The pre-averaging approach.” Stoch. Process. Their Appl. 119 (2009): 2249–2276. [Google Scholar] [CrossRef]
  14. P. Malliavin, and M.E. Mancino. “A Fourier transform method for nonparametric estimation of multivariate volatility.” Ann. Stat. 37 (2009): 1983–2010. [Google Scholar] [CrossRef]
  15. L. Zhang, P. Mykland, and Y. Aït-Sahalia. “A tale of two time scales: Determining integrated volatility with noisy high frequency data.” J. Am. Stat. Assoc. 100 (2005): 1394–1411. [Google Scholar] [CrossRef]
  16. F. Barsotti, and S. Sanfelici. “Firm’s volatility risk under microstructure noise.” In Mathematical and Statistical Methods for Actuarial Sciences and Finance. Edited by M. Corazza and C. Pizzi. Milan, Italy: Springer-Verlag, 2013, pp. 55–68. [Google Scholar]
  17. S. Heston. “A closed-form solution for options with stochastic volatility with applications to bond and currency options.” Rev. Financ. Stud. 6 (2003): 327–343. [Google Scholar] [CrossRef]
  18. D.S. Bates. “Jumps and Stochastic Volatility: Exchange Rate process Implicit in Deutsche Mark Options.” Rev. Financ. Stud. 9 (1996): 69–107. [Google Scholar] [CrossRef]
  19. R.C. Merton. “Option Pricing when Underlying Stock Returns are Discontinuous.” J. Financ. Econ. 3 (1976): 125–144. [Google Scholar] [CrossRef]
  20. O.E. Barndorff-Nielsen, P.R. Hansen, A. Lunde, and N. Shephard. “Subsampling realised kernels.” J. Econom. 160 (2010): 204–219. [Google Scholar] [CrossRef]
  21. O.E. Barndorff-Nielsen, and N. Shephard. “Power and bipower variation with stochastic volatility and jumps.” J. Financ. Econom. 2 (2004): 1–37. [Google Scholar] [CrossRef]
  22. F.M. Bandi, and J.R. Russel. “Market microstructure noise, integrated variance estimators, and the accuracy of asymptotic approximations.” J. Econom. 160 (2006): 145–159. [Google Scholar] [CrossRef]
  23. O.E. Barndorff-Nielsen, P.R. Hansen, A. Lunde, and N. Shephard. “Multivariate realised kernels: Consistent positive semi-definite estimators of the covariation of equity prices with noise and non-synchronous trading.” J. Econom. 162 (2008): 149–169. [Google Scholar] [CrossRef]
  24. M.E. Mancino, and S. Sanfelici. “Robustness of Fourier estimator of integrated volatility in the presence of microstructure noise.” Comput. Stat. Data Anal. 52 (2008): 2966–2989. [Google Scholar] [CrossRef]
  25. F. Corsi, D. Pirino, and R. Renó. “Threshold bipower variation and the impact of jumps on volatility forecasting.” J. Econom. 159 (2010): 276–288. [Google Scholar] [CrossRef]
  26. G. Bakshi, C. Cao, and Z. Chen. “Empirical performance of alternative option pricing models.” J. Finance 52 (1997): 2003–2049. [Google Scholar] [CrossRef]
  27. O.E. Barndorff-Nielsen, and N. Shephard. “Econometric analysis of realized volatility and its use in estimating stochastic volatility models.” J. R. Stat. Soc. Ser. B 64 (2002): 253–280. [Google Scholar] [CrossRef]
Figure 1. Time plot of the 1-min equity log-returns for Alcoa company over the period 10 November 2015 to 23 May 2016.
Figure 1. Time plot of the 1-min equity log-returns for Alcoa company over the period 10 November 2015 to 23 May 2016.
Econometrics 04 00031 g001
Figure 2. Histogram of the 1-min equity log-returns for Alcoa company over the period 10 November 2015 to 23 May 2016.
Figure 2. Histogram of the 1-min equity log-returns for Alcoa company over the period 10 November 2015 to 23 May 2016.
Econometrics 04 00031 g002
Figure 3. Time plot of the 5-years CDS spreads in basis points (lower panel) and corresponding equity volatility (upper panel) for our sample of 109 days after jump removal.
Figure 3. Time plot of the 5-years CDS spreads in basis points (lower panel) and corresponding equity volatility (upper panel) for our sample of 109 days after jump removal.
Econometrics 04 00031 g003
Figure 4. Asset volatility calibration and default probability computation. The plot shows how the estimated equity volatility Σ t E can affect the resulting asset volatility calibration and default probability computation. Subplot (i) reports the underlying asset volatility calibrated for different values of Σ t E ; subplot (ii) shows the absolute error on default probability when a 3y horizon is considered. The base case is Σ t s = 37 . 21 % , v t = 25 . 01 % (red dashed line in Subplot (i)) and the underlying asset volatility is calibrated by matching the 5y default probability to highlight the behaviour of our calibration algorithm. The absolute error on the 3y-default probability is defined through the variable D P A b s E r r : = D P ( Σ t E ) - D P ( Σ t s ) and reported in percentage (%).
Figure 4. Asset volatility calibration and default probability computation. The plot shows how the estimated equity volatility Σ t E can affect the resulting asset volatility calibration and default probability computation. Subplot (i) reports the underlying asset volatility calibrated for different values of Σ t E ; subplot (ii) shows the absolute error on default probability when a 3y horizon is considered. The base case is Σ t s = 37 . 21 % , v t = 25 . 01 % (red dashed line in Subplot (i)) and the underlying asset volatility is calibrated by matching the 5y default probability to highlight the behaviour of our calibration algorithm. The absolute error on the 3y-default probability is defined through the variable D P A b s E r r : = D P ( Σ t E ) - D P ( Σ t s ) and reported in percentage (%).
Econometrics 04 00031 g004
Figure 5. Equity log-returns. Time plot of the tick-by-tick normalized equity log-returns and the ACF on three sample days. (A) day with no jumps and first order autocorrelation structure; (B) day with jumps and first order autocorrelation structure; (C) day with no jumps and autocorrelation up to lag 15.
Figure 5. Equity log-returns. Time plot of the tick-by-tick normalized equity log-returns and the ACF on three sample days. (A) day with no jumps and first order autocorrelation structure; (B) day with jumps and first order autocorrelation structure; (C) day with no jumps and autocorrelation up to lag 15.
Econometrics 04 00031 g005
Figure 6. Intraday returns of equity log-prices. Comparison of the histogram of the normalized intraday returns and the density of N ( 0 , 1 ) . (A) day with no jumps and first order autocorrelation structure; (B) day with jumps and first order autocorrelation structure; (C) day with no jumps and autocorrelation up to lag 15.
Figure 6. Intraday returns of equity log-prices. Comparison of the histogram of the normalized intraday returns and the density of N ( 0 , 1 ) . (A) day with no jumps and first order autocorrelation structure; (B) day with jumps and first order autocorrelation structure; (C) day with no jumps and autocorrelation up to lag 15.
Econometrics 04 00031 g006
Figure 7. Heston model: default probability. The plot reports the average results obtained by computing default probabilities according to Equation (5). For each day in the sample and for each volatility estimator E we use the asset volatility v t ^ E , calibrated as explained in the text, in order to compute the default probabilities for any maturity (from 1 to 5 years). Results are based on 1000 daily Monte Carlo simulations; trading noise is defined by Equation (6).
Figure 7. Heston model: default probability. The plot reports the average results obtained by computing default probabilities according to Equation (5). For each day in the sample and for each volatility estimator E we use the asset volatility v t ^ E , calibrated as explained in the text, in order to compute the default probabilities for any maturity (from 1 to 5 years). Results are based on 1000 daily Monte Carlo simulations; trading noise is defined by Equation (6).
Econometrics 04 00031 g007
Figure 8. Bates model: default probability relative error. The plot shows default probability mean relative error given in Equation (9) for maturities from 1 to 5 years. Results are based on a sample of 934 days for our Monte Carlo simulations (due to jumps removal) and refer to trading noise (6).
Figure 8. Bates model: default probability relative error. The plot shows default probability mean relative error given in Equation (9) for maturities from 1 to 5 years. Results are based on a sample of 934 days for our Monte Carlo simulations (due to jumps removal) and refer to trading noise (6).
Econometrics 04 00031 g008
Figure 9. Bates model: default probability relative error. The plot shows default probability mean relative error given in Equation (9) for maturities from 1 to 5 years. Results are based on a sample of 901 days for our Monte Carlo simulations (due to jumps removal) and refer to trading noise (7).
Figure 9. Bates model: default probability relative error. The plot shows default probability mean relative error given in Equation (9) for maturities from 1 to 5 years. Results are based on a sample of 901 days for our Monte Carlo simulations (due to jumps removal) and refer to trading noise (7).
Econometrics 04 00031 g009
Table 1. Summary statistics of the 1-min transaction data on equity prices for the Alcoa company in the period from 10 November 2015 to 23 May 2016 (52,166 trades). ”Std. Dev.” denotes the sample standard deviation of the variable.
Table 1. Summary statistics of the 1-min transaction data on equity prices for the Alcoa company in the period from 10 November 2015 to 23 May 2016 (52,166 trades). ”Std. Dev.” denotes the sample standard deviation of the variable.
VariableMeanStd. Dev.MinMax
AA US Equity price9.021.056.1511.49
log-return (%)1.68 × 10 - 4 1.78 × 10 - 1 −4.66 × 10 0 3.73 × 10 0
Table 2. Summary statistics of the daily CDS premium data (basis points) for Alcoa company in the period from 10 November 2015 to 23 May 2016. ”Std. Dev.” denotes the sample standard deviation of the variable.
Table 2. Summary statistics of the daily CDS premium data (basis points) for Alcoa company in the period from 10 November 2015 to 23 May 2016. ”Std. Dev.” denotes the sample standard deviation of the variable.
VariableMeanStd. Dev.MinMax
CDS premium433.8791.09329.24661.63
Table 3. OLS estimates from regressions (1) of the 5-year CDS spread (%) on day t over a constant and each corresponding equity volatility estimate obtained by estimators defined in Section 2. Standard deviation and T-statistics are reported in brackets.
Table 3. OLS estimates from regressions (1) of the 5-year CDS spread (%) on day t over a constant and each corresponding equity volatility estimate obtained by estimators defined in Section 2. Standard deviation and T-statistics are reported in brackets.
Estimator ϕ 0 ^ ϕ 1 ^ R 2
Σ t R V 2.54764.09150.3810
Std. Dev.(0.2312)(0.5042)
T-statistics(11.0173)(8.1155)
Σ t T S 2.06154.96570.4118
Std. Dev.(0.2715)(0.5737)
T-statistics(7.5916)(8.6556)
Σ t H L 2.42684.41080.3884
Std. Dev.(0.2418)(0.5350)
T-statistics(10.0355)(8.2447)
Σ t K 2.52934.23300.3797
Std. Dev.(0.2340)(0.5229)
T-statistics(10.8107)(8.0944)
Σ t P A 2.61404.32760.3635
Std. Dev.(0.2314)(0.5535)
T-statistics(11.2956)(7.8184)
Σ t F 1.97355.13980.4236
Std. Dev.(0.2749)(0.5796)
T-statistics(7.1788)(8.8677)
Σ t B V 2.58914.12630.3801
Std. Dev.(0.2267)(0.5094)
T-statistics(11.4188)(8.1010)
Table 4. Base case parameters values.
Table 4. Base case parameters values.
ParameterBatesHeston
A 0 100100
r 5 % 5 %
Σ e 37 . 21 % 37 . 21 %
P D 5 y 1 . 37 % 1 . 37 %
B 43 . 13 43 . 13
δ 2 % 2 %
π a 4 . 04 % 3 . 49 %
v 0 25 . 01 % 21 . 65 %
k 0 . 87 0 . 74
θ 3 . 8 % 4 . 24 %
σ 4 . 42 % 4 . 01 %
ρ - 24 . 04 % - 24 . 02 %
λ 8 . 88 %
μ J 2 . 35 %
σ J 14 . 79 %
Note: The Table describes base case parameters values considered in our Monte Carlo simulation analysis. Σ e denotes the equity volatility, P D 5 y the 5-years default probability and π a the historical total asset risk premium. The values reported for the initial underlying asset v 0 are consistent with an equity risk premium π E = 6 % .
Table 5. Heston model: performance of equity volatility estimators.
Table 5. Heston model: performance of equity volatility estimators.
Heston Σ t BV Σ t RV Σ t RVSS Σ t TS Σ t HL Σ t K Σ t PA Σ t F
(a)MSE1.79 × 10 - 3 1.31 × 10 - 2 1.65 × 10 - 3 3.43 × 10 - 4 1.05 × 10 - 3 4.37 × 10 - 4 5.10 × 10 - 4 5.09 × 10 - 4
BIAS2.35 × 10 - 2 1.13 × 10 - 1 2.62 × 10 - 2 –1.47 × 10 - 3 –9.48 × 10 - 4 −7.69 × 10 - 4 –4.77 × 10 - 3 4.57 × 10 - 3
(b)MSE2.42 × 10 - 3 2.28 × 10 - 2 2.30 × 10 - 3 3.58 × 10 - 4 1.66 × 10 - 3 4.64 × 10 - 4 5.05 × 10 - 4 5.71 × 10 - 4
BIAS2.84 × 10 - 2 1.48 × 10 - 1 3.21 × 10 - 2 −7.66 × 10 - 4 5.10 × 10 - 4 –3.72 × 10 - 4 −4.20 × 10 - 3 6.06 × 10 - 3
Note: The table shows mean squared error (MSE) and BIAS for different equity volatility estimators. Panel (a) refers to the trading noise given in Equation (6), panel (b) to Equation (7), with α = 0 . 5 . The underlying dynamics follows Heston [17] model.
Table 6. Bates model: performance of equity volatility estimators.
Table 6. Bates model: performance of equity volatility estimators.
Bates Σ t BV Σ t RV Σ t RVSS Σ t TS Σ t HL Σ t K Σ t PA Σ t F
(a)MSE2.77 × 10 - 3 1.63 × 10 - 2 2.36 × 10 - 3 4.69 × 10 - 4 1.23 × 10 - 3 5.44 × 10 - 4 5.50 × 10 - 4 6.75 × 10 - 4
BIAS3.16 × 10 - 2 1.24 × 10 - 1 3.19 × 10 - 2 7.18 × 10 - 4 2.08 × 10 - 3 2.10 × 10 - 3 -2.45 × 10 - 3 8.40 × 10 - 3
(b)MSE3.36 × 10 - 3 2.78 × 10 - 2 3.12 × 10 - 3 6.37 × 10 - 4 1.88 × 10 - 3 7.18 × 10 - 4 7.51 × 10 - 4 8.91 × 10 - 4
BIAS3.40 × 10 - 2 1.63 × 10 - 1 3.71 × 10 - 2 –6.38 × 10 - 4 1.54 × 10 - 3 2.43 × 10 - 3 –1.65 × 10 - 3 1.03 × 10 - 2
Note: The table shows mean squared error (MSE) and BIAS for different equity volatility estimators. Panel (a) refers to the trading noise given in Equation (6), panel (b) to Equation (7), with α = 0 . 5 . The underlying dynamics follows [18] model.
Table 7. Heston model: calibrated asset volatility.
Table 7. Heston model: calibrated asset volatility.
Heston
MeanMedian10 perc.90 perc.MinMaxStd Dev.
(a) v t ^ B V 0.20820.20560.18330.23500.16140.29600.0205
v t ^ R V 0.17040.16990.15990.18150.14530.20270.0085
v t ^ R V S S 0.20630.20450.18380.23080.15830.29220.0182
v t ^ T S 0.22140.22020.20350.24100.18310.27180.0151
v t ^ H L 0.22290.21950.19500.25540.17640.33850.0243
v t ^ K 0.22100.22020.20210.24060.18020.27360.0155
v t ^ P A 0.22380.22280.20290.24530.18020.27740.0166
v t ^ F 0.21780.21690.19840.23750.17760.26840.0156
v t ^ B V / v t 0.96170.94980.84641.08530.74541.36730.0948
v t ^ R V / v t 0.78700.78450.73870.83830.67090.93640.0392
v t ^ R V S S / v t 0.95280.94480.84891.06590.73101.34950.0840
v t ^ T S / v t 1.02261.01720.94001.11340.84561.25540.0696
v t ^ H L / v t 1.02971.01400.90051.17960.81491.56370.1123
v t ^ K / v t 1.02101.01690.93351.11130.83261.26370.0714
v t ^ P A / v t 1.03371.02900.93711.13290.83261.28150.0766
v t ^ F / v t 1.00611.00180.91641.09700.82031.23960.0722
(b) v t ^ B V 0.20620.20400.17910.23680.14850.29020.0227
v t ^ R V 0.16020.16000.15030.17050.13610.18460.0079
v t ^ R V S S 0.20360.20230.18020.23000.15710.29420.0196
v t ^ T S 0.22090.21930.20210.24130.17810.27770.0153
v t ^ H L 0.22340.21800.19090.26370.16580.37790.0309
v t ^ K 0.22090.21930.20070.24290.17820.29420.0164
v t ^ P A 0.22350.22260.20200.24620.17910.29590.0172
v t ^ F 0.21700.21560.19670.23970.17600.28490.0165
v t ^ B V / v t 0.93910.92820.82001.07570.70791.27510.0984
v t ^ R V / v t 0.72980.72960.69470.76350.63080.82850.0273
v t ^ R V S S / v t 0.92760.92160.82271.03570.73961.30230.0838
v t ^ T S / v t 1.00611.00230.93481.07890.86731.20230.0553
v t ^ H L / v t 1.01760.99110.87571.19510.78261.80070.1353
v t ^ K / v t 1.00601.00060.92831.09140.84231.24110.0627
v t ^ P A / v t 1.01761.01210.93591.10620.85701.26060.0662
v t ^ F / v t 0.98830.98340.91031.07340.83881.21480.0642
Note: The table shows descriptive statistics of calibrated asset volatility for different equity volatility estimators. Results are based on 1000 daily Monte Carlo simulations. Panel (a) refers to the trading noise given in Equation (6), panel (b) to Equation (7), with α = 0 . 5 .
Table 8. Bates model: calibrated asset volatility.
Table 8. Bates model: calibrated asset volatility.
Bates
MeanMedian10 perc.90 perc.MinMaxStd Dev.
(a) v t ^ B V 0.25260.24970.20760.30120.15580.37920.0381
v t ^ R V 0.20420.20380.18930.22110.15970.25020.0125
v t ^ R V S S 0.25020.24770.22020.28390.18260.35030.0254
v t ^ T S 0.27030.27010.24590.29620.20180.33600.0203
v t ^ H L 0.27140.26730.23610.31400.18570.37390.0304
v t ^ K 0.26940.26730.24440.29790.20070.34110.0215
v t ^ P A 0.27300.27100.24700.30180.20190.34980.0223
v t ^ F 0.26480.26240.23960.29280.19800.33400.0215
v t ^ B V / v t 0.93910.92410.76051.12600.53921.48020.1501
v t ^ R V / v t 0.75790.75730.71690.80340.58010.88810.0355
v t ^ R V S S / v t 0.92830.92080.82621.04730.64911.21190.0860
v t ^ T S / v t 1.00271.00200.93121.08330.78401.20040.0600
v t ^ H L / v t 1.00680.99310.88841.14670.74371.41640.1036
v t ^ K / v t 0.99940.99310.92221.08910.77941.23730.0659
v t ^ P A / v t 1.01271.00710.93151.10730.78521.24100.0684
v t ^ F / v t 0.98250.97710.90271.07450.77271.21130.0668
(b) v t ^ B V 0.25210.24850.20670.30550.16160.37940.0387
v t ^ R V 0.19090.19090.17870.20430.15110.22520.0106
v t ^ R V S S 0.24740.24500.21600.28250.17270.36270.0270
v t ^ T S 0.27120.26940.24510.29900.20190.34800.0212
v t ^ H L 0.27260.26820.23180.32280.18470.37970.0347
v t ^ K 0.26920.26860.24240.29780.19180.34560.0225
v t ^ P A 0.27260.27210.24460.30290.19210.35680.0238
v t ^ F 0.26360.26330.23680.29240.18800.34810.0224
v t ^ B V / v t 0.93780.92220.75701.14820.54351.48660.1514
v t ^ R V / v t 0.70890.70960.67100.74940.56720.82060.0322
v t ^ R V S S / v t 0.91870.91010.81121.03890.62861.40600.0940
v t ^ T S / v t 1.00701.00240.92461.09720.73061.28130.0700
v t ^ H L / v t 1.01290.99810.87111.18640.68011.47730.1275
v t ^ K / v t 0.99970.99710.91151.09180.69011.25400.0724
v t ^ P A / v t 1.01201.00860.92221.10730.69111.29460.0766
v t ^ F / v t 0.97860.97670.89521.06810.67661.26320.0719
Note: The table shows descriptive statistics of calibrated asset volatility for different equity volatility estimators. Results are based on a sample of 934 days in case (a) and 901 days in case (b) for our Monte Carlo simulations (due to jumps removal). Panel (a) refers to the trading noise given in Equation (6), panel (b) to Equation (7), with α = 0 . 5 .
Table 9. Heston model: default probability relative error.
Table 9. Heston model: default probability relative error.
Heston DP Err BV DP Err RV DP Err RVSS DP Err TS DP Err HL DP Err K DP Err PA DP Err F
(a)1 y−13.3336−92.8243−28.511221.993878.971323.447742.20607.7918
1.5 y−17.5285−78.4652−24.56419.059426.26139.140118.37080.2651
2 y−14.7317−62.4385−18.88724.789512.37084.602010.1900−1.1713
2.5 y−11.3451−47.6616−14.00022.80636.70212.58226.2117−1.3314
3 y−8.2927−34.8476−10.02111.70673.84121.51183.9022−1.1385
3.5 y−5.6774−23.9395−6.77151.03302.21240.88722.4191−0.8426
4 y−3.4884−14.7071−4.12150.57001.16960.47371.3652−0.5485
4.5 y−1.6097−6.8180−1.88940.24380.48320.19780.5945−0.2621
(b)1 y−13.4191−96.3065−31.870820.9369168.709025.542842.40586.8763
1.5 y−19.3928−85.0507−28.22788.395144.90369.722718.0795−0.8700
2 y−16.6626−69.9185−21.95034.322519.17464.81879.9133−2.1040
2.5 y−12.9684−54.5511−16.37672.47939.74222.67836.0040−2.0250
3 y−9.5387−40.4839−11.76901.48195.31111.55813.7556−1.6324
3.5 y−6.5555−28.1046−7.97720.88452.92800.91072.3219−1.1722
4 y−4.0377−17.3914−4.86190.48141.49190.48431.3077−0.7493
4.5 y−1.8674−8.1086−2.23290.20370.59520.20150.5683−0.3538
Note: The table shows default probability mean relative error given in Equation (9) for maturities from 1 to 4.5 years. Results are based on 1000 daily Monte Carlo simulations using different equity volatility estimators. Default probability relative errors are in percentage (%). The simulation is based on equi-spaced intra-day data generated on a 1-min frequency. Panel (a) refers to a trading noise defined by Equation (6), panel (b) to Equation (7).
Table 10. Bates model: default probability relative error.
Table 10. Bates model: default probability relative error.
Bates DP Err BV DP Err RV DP Err RVSS DP Err TS DP Err HL DP Err K DP Err PA DP Err F
(a)1 y14.4309−87.1895−29.929311.213336.544710.647422.8589−2.1273
1.5 y−6.3833−75.8094−25.65935.056815.01334.119811.6211−4.3575
2 y−10.3506−62.5783−20.50192.66517.55501.81406.9034−4.1742
2.5 y−10.1010−49.6830−15.84341.50974.14380.81914.3846−3.4849
3 y−8.4835−37.7336−11.78090.88102.34540.35682.8397−2.6973
3.5 y−6.4324−26.8797−8.26760.49381.28130.11941.7782−1.9495
4 y−4.2495−17.0514−5.17430.25560.64470.01891.0198−1.2420
4.5 y−2.0866−8.1234−2.44130.10090.2468−0.01160.4469−0.5943
(b)1 y13.4257−91.9420−31.616918.498759.624013.096425.4460−3.3672
1.5 y−6.7622−82.7271−28.05548.786724.57315.157112.4494−5.5782
2 y−10.6074−70.1849−22.71474.889612.50412.32037.1783−5.2046
2.5 y−10.3035−56.8420−17.68592.92216.95221.07764.4463−4.3048
3 y−8.6408−43.8048−13.21621.79023.98220.48532.8177−3.3196
3.5 y−6.5510−31.5429−9.30841.06242.21630.18031.7301−2.3910
4 y−4.3266−20.1699−5.84250.57981.13330.04360.9752−1.5208
4.5 y−2.1247−9.6830−2.76210.24220.4423−0.00500.4206−0.7267
Note: The table shows default probability mean relative error given in Equation (9) for maturities from 1 to 4.5 years. Due to jumps removal, results are based on a sample of 934 days in case (a) and 901 days in case (b) for our Monte Carlo simulations using different equity volatility estimators. Default probability relative errors are in percentage (%). The simulation is based on equispaced intra-day data generated on a 1-min frequency. Panel (a) refers to a trading noise being defined by Equation (6).
Table 11. Test cases.
Table 11. Test cases.
Heston
ρ v 0
ρ = - 30 . 00 % v 0 = 20 . 92 %
ρ = - 24 . 02 % v 0 = 21 . 65 %
ρ = 0 . 0 v 0 = 24 . 20 %
Note: The table summarizes the values considered in the sensitivity analysis for three alternative test cases associated with a 5-years default probability P D 5 y = 1 . 37 % .
Table 12. Heston model: leverage effect on equity volatility estimators.
Table 12. Heston model: leverage effect on equity volatility estimators.
Heston Σ t BV Σ t RV Σ t RVSS Σ t TS Σ t HL Σ t K Σ t PA Σ t F
ρ = - 30 %
MSE2.25 × 10 - 3 2.18 × 10 - 2 2.04 × 10 - 3 3.53 × 10 - 4 1.45 × 10 - 3 4.23 × 10 - 4 4.65 × 10 - 4 5.28 × 10 - 4
BIAS2.57 × 10 - 2 1.46 × 10 - 1 2.92 × 10 - 2 −5.06 × 10 - 4 −7.05 × 10 - 4 −1.96 × 10 - 5 −4.01 × 10 - 3 5.84 × 10 - 3
ρ = - 24 . 02 %
MSE2.42 × 10 - 3 2.28 × 10 - 2 2.30 × 10 - 3 3.58 × 10 - 4 1.66 × 10 - 3 4.64 × 10 - 4 5.05 × 10 - 4 5.71 × 10 - 4
BIAS2.84 × 10 - 2 1.48 × 10 - 1 3.21 × 10 - 2 −7.66 × 10 - 4 5.10 × 10 - 4 −3.72 × 10 - 4 −4.20 × 10 - 3 6.06 × 10 - 3
ρ = 0
MSE2.86 × 10 - 3 2.70 × 10 - 2 2.64 × 10 - 3 4.95 × 10 - 4 2.11 × 10 - 3 4.93 × 10 - 4 5.52 × 10 - 4 6.25 × 10 - 4
BIAS3.03 × 10 - 2 1.62 × 10 - 1 3.50 × 10 - 2 −1.62 × 10 - 3 −1.24 × 10 - 3 −1.64 × 10 - 4 −4.61 × 10 - 3 6.55 × 10 - 3
Note: The table shows mean squared error (MSE) and BIAS for different equity volatility estimators. The panel refers to the trading noise given in Equation (7), with α = 0 . 5 . Alternative values for the correlation ρ are considered, i.e., ρ { - 30 % , - 24 . 02 % , 0 } .
Table 13. Heston model: leverage effect on calibrated asset volatility.
Table 13. Heston model: leverage effect on calibrated asset volatility.
Heston
MeanMedian10 perc.90 perc.MinMaxStd Dev.
ρ = - 30 %
v t ^ B V / v t 0.94140.92870.80991.08200.71701.68300.1100
v t ^ R V / v t 0.71120.71050.67530.74810.63330.81610.0285
v t ^ R V S S / v t 0.92860.91860.81801.04010.72801.56500.0909
v t ^ T S / v t 1.00551.00290.93221.08150.85021.28040.0598
v t ^ H L / v t 1.02090.99770.87451.18640.75051.76500.1370
v t ^ K / v t 1.00480.99910.92501.08810.82961.24800.0651
v t ^ P A / v t 1.01801.01230.93221.10790.82981.26400.0689
v t ^ F / v t 0.98750.98220.90611.07580.79971.22540.0672
ρ = - 24 . 02 %
v t ^ B V / v t 0.93910.92820.82001.07570.70791.27510.0984
v t ^ R V / v t 0.72980.72960.69470.76350.63080.82850.0273
v t ^ R V S S / v t 0.92760.92160.82271.03570.73961.30230.0838
v t ^ T S / v t 1.00611.00230.93481.07890.86731.20230.0553
v t ^ H L / v t 1.01760.99110.87571.19510.78261.80070.1353
v t ^ K / v t 1.00601.00060.92831.09140.84231.24110.0627
v t ^ P A / v t 1.01761.01210.93591.10620.85701.26060.0662
v t ^ F / v t 0.98830.98340.91031.07340.83881.21480.0642
ρ = 0
v t ^ B V / v t 0.95180.94200.85451.06140.74861.26000.0809
v t ^ R V / v t 0.77750.77680.74960.80610.70930.86360.0218
v t ^ R V S S / v t 0.94030.93710.85981.03160.78101.16550.0653
v t ^ T S / v t 1.00751.00310.94811.06870.87981.23560.0497
v t ^ H L / v t 1.02010.99720.89811.17180.80001.58790.1145
v t ^ K / v t 1.00440.99920.94401.06890.88411.18710.0494
v t ^ P A / v t 1.01471.01210.94971.08120.88801.21850.0529
v t ^ F / v t 0.99050.98790.92841.05350.86281.18340.0509
Note: The table shows descriptive statistics of calibrated asset volatility for different equity volatility estimators. The panel refers to the trading noise given in Equation (7), with α = 0 . 5 . Alternative values for the correlation ρ are considered, i.e., ρ { - 30 % , - 24 . 02 % , 0 } .
Table 14. Heston model: leverage effect on default probability estimation.
Table 14. Heston model: leverage effect on default probability estimation.
Heston DP Err BV DP Err RV DP Err RVSS DP Err TS DP Err HL DP Err K DP Err PA DP Err F
ρ = - 30 %
1 y4.7996−96.9042−23.740922.3603164.522924.327542.80466.8151
1.5 y−14.3485−86.5006−25.64298.553145.94919.060418.3006−0.9973
2 y−14.2628−71.7115−20.58324.243420.16674.370210.0212−2.2244
2.5 y−11.6035−56.2816−15.53882.358710.47612.36516.0606−2.1285
3 y−8.7107−41.9399−11.23251.37175.82611.33783.7860−1.7152
3.5 y−6.0728−29.2146−7.65460.79143.25560.75702.3270−1.2429
4 y−3.7607−18.1143−4.66810.42491.69580.39791.3155−0.7859
4.5 y−1.7491−8.4585−2.14840.17610.68650.16170.5710−0.3716
ρ = - 24 . 02 %
1 y−13.4191−96.3065−31.870820.9369168.709025.542842.40586.8763
1.5 y−19.3928−85.0507-28.22788.395144.90369.722718.0795−0.8700
2 y−16.6626-69.9185-21.95034.322519.17464.81879.9133−2.1040
2.5 y−12.9684−54.5511−16.37672.47939.74222.67836.0040−2.0250
3 y−9.5387−40.4839−11.76901.48195.31111.55813.7556−1.6324
3.5 y−6.5555−28.1046−7.97720.88452.92800.91072.3219−1.1722
4 y−4.0377−17.3914−4.86190.48141.49190.48431.3077−0.7493
4.5 y−1.8674−8.1086−2.23290.20370.59520.20150.5683−0.3538
ρ = 0
1y−15.5280−94.5066−34.758422.5433132.970018.009733.33032.8602
1.5y−18.0677−81.4410−27.51599.416040.91347.159615.0375−1.7585
2y−14.9654−65.8099−20.67775.034118.84183.62908.5020−2.2432
2.5y−11.4706−50.7851−15.17412.983310.12422.05825.2617−1.9516
3y−8.3639−37.4239−10.79861.83395.78291.21933.3458−1.5084
3.5y−5.7269−25.8687−7.28101.11173.30500.71652.0835−1.0672
4y−3.5034−15.9461−4.40690.61891.75620.38801.1869−0.6675
4.5y−1.6112−7.3986−2.01280.26920.72880.16650.5234−0.3087
Note: The table shows default probability mean relative error given in Equation (9) for maturities from 1 to 4.5 years. Results are based on a sample of 1000 days for our Monte Carlo simulations. Default probability relative errors are in percentage (%). The panel refers to the trading noise given in Equation (7), with α = 0 . 5 . Alternative values for the correlation ρ are considered, i.e., ρ { - 30 % , - 24 . 02 % , 0 } .
Table 15. Test cases.
Table 15. Test cases.
Bates
μ J v 0
μ J = 3 . 35 % v 0 = 25 . 04 %
μ J = 2 . 35 % v 0 = 25 . 01 %
μ J = 1 . 35 % v 0 = 24 . 94 %
Note: The table summarizes the values considered in this sensitivity analysis for three alternative test cases associated with a 5-years default probability P D 5 y = 1 . 37 % .
Table 16. Bates model: effect of average jumps size on equity volatility estimators.
Table 16. Bates model: effect of average jumps size on equity volatility estimators.
Bates Σ t BV Σ t RV Σ t RVSS Σ t TS Σ t HL Σ t K Σ t PA Σ t F
μ J = 3 . 35 %
MSE3.55 × 10 - 3 2.80 × 10 - 2 2.94 × 10 - 3 6.29 × 10 - 4 1.84 × 10 - 3 6.23 × 10 - 4 6.70 × 10 - 4 7.61 × 10 - 4
BIAS3.36 × 10 - 2 1.64 × 10 - 1 3.50 × 10 - 2 2.58 × 10 - 4 2.12 × 10 - 3 5.23 × 10 - 4 −3.99 × 10 - 3 7.70 × 10 - 3
μ J = 2 . 35 %
MSE3.36 × 10 - 3 2.78 × 10 - 2 3.12 × 10 - 3 6.37 × 10 - 4 1.88 × 10 - 3 7.18 × 10 - 4 7.51 × 10 - 4 8.91 × 10 - 4
BIAS3.40 × 10 - 2 1.63 × 10 - 1 3.71 × 10 - 2 −6.38 × 10 - 4 1.54 × 10 - 3 2.43 × 10 - 3 −1.65 × 10 - 3 1.03 × 10 - 2
μ J = 1 . 35 %
MSE3.32 × 10 - 3 2.81 × 10 - 2 2.74 × 10 - 3 5.31 × 10 - 4 1.81 × 10 - 3 5.42 × 10 - 4 5.85 × 10 - 4 6.83 × 10 - 4
BIAS3.51 × 10 - 2 1.65 × 10 - 1 3.58 × 10 - 2 4.81 × 10 - 4 2.37 × 10 - 3 1.15 × 10 - 3 −3.10 × 10 - 3 8.59 × 10 - 3
Note: The table shows mean squared error (MSE) and BIAS for different equity volatility estimators. The panel refers to the trading noise given in Equation (7), with α = 0 . 5 . Alternative values for the jump parameter μ J are considered: μ J { 1 . 35 % , 2 . 35 % , 3 . 35 % } .
Table 17. Bates model: effect of average jumps size on calibrated asset volatility.
Table 17. Bates model: effect of average jumps size on calibrated asset volatility.
Bates
MeanMedian10 perc.90 perc.MinMaxStd Dev.
μ J = 3 . 35 %
v t ^ B V / v t 0.93920.92670.76021.14960.41261.44780.1544
v t ^ R V / v t 0.71020.70980.67470.75010.54550.81800.0311
v t ^ R V S S / v t 0.92470.91620.81291.04810.65221.29650.0943
v t ^ T S / v t 1.00500.99980.92261.09320.77001.28980.0705
v t ^ H L / v t 1.01130.99140.87031.18600.68871.48650.1253
v t ^ K / v t 1.00471.00020.92481.09210.74241.28300.0695
v t ^ P A / v t 1.01821.01580.92881.10870.74041.27690.0733
v t ^ F / v t 0.98540.98400.90331.07200.71901.22040.0698
μ J = 2 . 35 %
v t ^ B V / v t 0.93780.92220.75701.14820.54351.48660.1514
v t ^ R V / v t 0.70890.70960.67100.74940.56720.82060.0322
v t ^ R V S S / v t 0.91870.91010.81121.03890.62861.40600.0940
v t ^ T S / v t 1.00701.00240.92461.09720.73061.28130.0700
v t ^ H L / v t 1.01290.99810.87111.18640.68011.47730.1275
v t ^ K / v t 0.99970.99710.91151.09180.69011.25400.0724
v t ^ P A / v t 1.01201.00860.92221.10730.69111.29460.0766
v t ^ F / v t 0.97860.97670.89521.06810.67661.26320.0719
μ J = 1 . 35 %
v t ^ B V / v t 0.93170.91960.75121.12450.55981.48280.1508
v t ^ R V / v t 0.70570.70520.66880.74430.57630.81890.0303
v t ^ R V S S / v t 0.92000.91690.81271.03600.68701.25060.0871
v t ^ T S / v t 1.00390.99920.92281.08750.76331.32450.0657
v t ^ H L / v t 1.01180.98920.86601.19910.70201.45250.1273
v t ^ K / v t 1.00210.99940.91731.08500.76101.23450.0658
v t ^ P A / v t 1.01491.01400.92621.09920.76051.27290.0699
v t ^ F / v t 0.98190.98220.89871.06690.74531.23600.0661
Note: The table shows descriptive statistics of calibrated asset volatility for different equity volatility estimators. The panel refers to the trading noise given in Equation (7), with α = 0 . 5 . Alternative value for the jump parameter μ J are considered: μ J { 1 . 35 % , 2 . 35 % , 3 . 35 % } . Results are based on samples of 905 , 901 , 904 days respectively (due to jumps removal over a sample of 1000 days).
Table 18. Bates model: effect of average jumps size on default probability estimation.
Table 18. Bates model: effect of average jumps size on default probability estimation.
Bates DP Err BV DP Err RV DP Err RVSS DP Err TS DP Err HL DP Err K DP Err PA DP Err F
μ J = 3 . 35 %
1 y15.4331-92.3935−28.662817.479157.903816.434529.88760.8040
1.5 y−5.5499−82.9953−25.83037.856623.32467.423515.4241−2.6013
2 y−9.8511−70.3362−21.00374.197611.72653.98509.3216−2.9963
2.5 y−9.8283−56.9149−16.37582.42246.46312.31126.0103−2.6650
3 y−8.3473−43.8289−12.24201.43843.67641.37963.9382−2.1293
3.5 y−6.3752−31.5402−8.61990.82532.03040.79682.4951−1.5712
4 y−4.2349−20.1568−5.40920.43681.03080.42461.4449−1.0138
4.5 y−2.0879−9.6728−2.55610.17680.39940.17310.6388−0.4898
μ J = 2 . 35 %
1 y13.4257−91.9420−31.616918.498759.624013.096425.4460−3.3672
1.5 y−6.7622−82.7271−28.05548.786724.57315.157112.4494−5.5782
2 y−10.6074−70.1849−22.71474.889612.50412.32037.1783−5.2046
2.5 y−10.3035−56.8420−17.68592.92216.95221.07764.4463−4.3048
3 y−8.6408−43.8048−13.21621.79023.98220.48532.8177−3.3196
3.5 y−6.5510−31.5429−9.30841.06242.21630.18031.7301−2.3910
4 y−4.3266−20.1699−5.84250.57981.13330.04360.9752−1.5208
4.5 y−2.1247−9.6830−2.76210.24220.4423−0.00500.4206−0.7267
μ J = 1 . 35 %
1 y8.1961−91.5704−33.420314.175857.091512.316124.6874−3.0626
1.5 y−9.5962−82.6645−28.52016.464323.55545.401812.8792−4.7706
2 y−12.5535−70.2623−22.75043.451411.95532.74557.7550−4.4219
2.5 y−11.7049−56.9596−17.56331.98386.63311.49144.9708−3.6450
3 y−9.6430−43.9142−13.04881.17343.79600.82993.2388−2.8032
3.5 y−7.2316−31.6256−9.15010.67072.11230.43982.0399−2.0157
4 y−4.7442−20.2222−5.72380.35331.08060.21371.1752−1.2808
4.5 y−2.3185−9.7054−2.69870.14230.42220.07810.5170−0.6114
Note: The table shows default probability mean relative error given in Equation (9) for maturities from 1 to 4.5 years. Results are based on a sample of 905, 901, 904 days (due to jumps removal over a sample of 1000 days), respectively, in the three cases depicted below for our Monte Carlo simulations. Default probability relative errors are in percentage (%). The panel refers to a trading noise defined by Equation (7). Alternative values for the jump parameter μ J are considered: μ J { 1.35 % , 2.35 % , 3.35 % } .

Share and Cite

MDPI and ACS Style

Barsotti, F.; Sanfelici, S. Market Microstructure Effects on Firm Default Risk Evaluation. Econometrics 2016, 4, 31. https://doi.org/10.3390/econometrics4030031

AMA Style

Barsotti F, Sanfelici S. Market Microstructure Effects on Firm Default Risk Evaluation. Econometrics. 2016; 4(3):31. https://doi.org/10.3390/econometrics4030031

Chicago/Turabian Style

Barsotti, Flavia, and Simona Sanfelici. 2016. "Market Microstructure Effects on Firm Default Risk Evaluation" Econometrics 4, no. 3: 31. https://doi.org/10.3390/econometrics4030031

APA Style

Barsotti, F., & Sanfelici, S. (2016). Market Microstructure Effects on Firm Default Risk Evaluation. Econometrics, 4(3), 31. https://doi.org/10.3390/econometrics4030031

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop