Source author record

José E. Figueroa-López

José E. Figueroa-López appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2022arXiv

Efficient Volatility Estimation for Lévy Processes with Jumps of Unbounded Variation

Statistical inference for stochastic processes based on high-frequency observations has been an active research area for more than a decade. One of the most well-known and widely studied problems is that of estimation of the quadratic variation of the continuous component of an Itô semimartingale with jumps. Several rate- and variance-efficient estimators have been proposed in the literature when the jump component is of bounded variation. However, to date, very few methods can deal with jumps of unbounded variation. By developing new high-order expansions of the truncated moments of a Lévy process, we construct a new rate- and variance-efficient estimator for a class of Lévy processes of unbounded variation, whose small jumps behave like those of a stable Lévy process with Blumenthal-Getoor index less than $8/5$. The proposed method is based on a two-step debiasing procedure for the truncated realized quadratic variation of the process. Our Monte Carlo experiments indicate that the method outperforms other efficient alternatives in the literature in the setting covered by our theoretical framework.

preprint2022arXiv

Estimation of Tempered Stable Lévy Models of Infinite Variation

We propose a new method for the estimation of a semiparametric tempered stable Lévy model. The estimation procedure combines iteratively an approximate semiparametric method of moment estimator, Truncated Realized Quadratic Variations (TRQV), and a newly found small-time high-order approximation for the optimal threshold of the TRQV of tempered stable processes. The method is tested via simulations to estimate the volatility and the Blumenthal-Getoor index of the generalized CGMY model as well as the integrated volatility of a Heston-type model with CGMY jumps. The method outperforms other efficient alternatives proposed in the literature when working with a Lévy process (i.e., the volatility is constant), or when the index of jump intensity $Y$ is larger than $3/2$ in the presence of stochastic volatility.

preprint2022arXiv

Kernel Estimation of Spot Volatility with Microstructure Noise Using Pre-Averaging

We first revisit the problem of estimating the spot volatility of an Itô semimartingale using a kernel estimator. We prove a Central Limit Theorem with optimal convergence rate for a general two-sided kernel. Next, we introduce a new pre-averaging/kernel estimator for spot volatility to handle the microstructure noise of ultra high-frequency observations. We prove a Central Limit Theorem for the estimation error with an optimal rate and study the optimal selection of the bandwidth and kernel functions. We show that the pre-averaging/kernel estimator's asymptotic variance is minimal for exponential kernels, hence, justifying the need of working with kernels of unbounded support as proposed in this work. We also develop a feasible implementation of the proposed estimators with optimal bandwidth. Monte Carlo experiments confirm the superior performance of the devised method.

preprint2021arXiv

Market Making with Stochastic Liquidity Demand: Simultaneous Order Arrival and Price Change Forecasts

We provide an explicit characterization of the optimal market making strategy in a discrete-time Limit Order Book (LOB). In our model, the number of filled orders during each period depends linearly on the distance between the fundamental price and the market maker's limit order quotes, with random slope and intercept coefficients. The high-frequency market maker (HFM) incurs an end-of-the-day liquidation cost resulting from linear price impact. The optimal placement strategy incorporates in a novel and parsimonious way forecasts about future changes in the asset's fundamental price. We show that the randomness in the demand slope reduces the inventory management motive, and that a positive correlation between demand slope and investors' reservation prices leads to wider spreads. Our analysis reveals that the simultaneous arrival of buy and sell market orders (i) reduces the shadow cost of inventory, (ii) leads the HFM to reduce price pressures to execute larger flows, and (iii) introduces patterns of nonlinearity in the intraday dynamics of bid and ask spreads. Our empirical study shows that the market making strategy outperforms those which ignores randomness in demand, simultaneous arrival of buy and sell market orders, and local drift in the fundamental price.

preprint2020arXiv

Optimal Iterative Threshold-Kernel Estimation of Jump Diffusion Processes

In this paper, we propose a new threshold-kernel jump-detection method for jump-diffusion processes, which iteratively applies thresholding and kernel methods in an approximately optimal way to achieve improved finite-sample performance. We use the expected number of jump misclassifications as the objective function to optimally select the threshold parameter of the jump detection scheme. We prove that the objective function is quasi-convex and obtain a new second-order infill approximation of the optimal threshold in closed form. The approximate optimal threshold depends not only on the spot volatility, but also the jump intensity and the value of the jump density at the origin. Estimation methods for these quantities are then developed, where the spot volatility is estimated by a kernel estimator with thresholding and the value of the jump density at the origin is estimated by a density kernel estimator applied to those increments deemed to contain jumps by the chosen thresholding criterion. Due to the interdependency between the model parameters and the approximate optimal estimators built to estimate them, a type of iterative fixed-point algorithm is developed to implement them. Simulation studies for a prototypical stochastic volatility model show that it is not only feasible to implement the higher-order local optimal threshold scheme but also that this is superior to those based only on the first order approximation and/or on average values of the parameters over the estimation time period.

preprint2016arXiv

One-level limit order book models with memory and variable spread

We propose a new model for the level I of a Limit Order Book (LOB), which incorporates the information about the standing orders at the opposite side of the book after each price change and the arrivals of new orders within the spread. Our main result gives a diffusion approximation for the mid-price process. To illustrate the applicability of the considered framework, we also propose a feasible method to compute several quantities of interest}, such as the distribution of the time span between price changes and the probability of consecutive price increments conditioned on the current state of the book. The proposed method is used to develop an efficient simulation scheme for the price dynamics, which is then applied to assess numerically the accuracy of the diffusion approximation.

preprint2016arXiv

Optimal Kernel Estimation of Spot Volatility of Stochastic Differential Equations

Kernel Estimation is one of the most widely used estimation methods in non-parametric Statistics, having a wide-range of applications, including spot volatility estimation of stochastic processes. The selection of bandwidth and kernel function is of great importance, especially for the finite sample settings commonly encountered in econometric applications. In the context of spot volatility estimation, most of the proposed selection methods are either largely heuristic or just formally stated without any feasible implementation. In this work, an objective method of bandwidth and kernel selection is proposed, under some mild conditions on the volatility, which not only cover classical Brownian motion driven dynamics but also some processes driven by long-memory fractional Brownian motions or other Gaussian processes. We characterize the leading order terms of the Mean Squared Error, which are also ratified by central limit theorems for the estimation error. As a byproduct, an approximated optimal bandwidth is then obtained in closed form. This result allows us to develop a feasible plug-in type bandwidth selection procedure, for which, as a sub-problem, we propose a new estimator of the volatility of volatility. The optimal selection of kernel function is also discussed. For Brownian Motion type volatilities, the optimal kernel function is proved to be the exponential kernel. For fractional Brownian motion type volatilities, numerical results to compute the optimal kernel are devised and, for the deterministic volatility case, explicit optimal kernel functions of different orders are derived. Simulation studies further confirm the good performance of the proposed methods.

preprint2015arXiv

Short-time asymptotics for the implied volatility skew under a stochastic volatility model with Lévy jumps

The implied volatility skew has received relatively little attention in the literature on short-term asymptotics for financial models with jumps, despite its importance in model selection and calibration. We rectify this by providing high-order asymptotic expansions for the at-the-money implied volatility skew, under a rich class of stochastic volatility models with independent stable-like jumps of infinite variation. The case of a pure-jump stable-like Lévy model is also considered under the minimal possible conditions for the resulting expansion to be well defined. Unlike recent results for "near-the-money" option prices and implied volatility, the results herein aid in understanding how the implied volatility smile near expiry is affected by important features of the continuous component, such as the leverage and vol-of-vol parameters. As intermediary results we obtain high-order expansions for at-the-money digital call option prices, which furthermore allow us to infer analogous results for the delta of at-the-money options. Simulation results indicate that our asymptotic expansions give good fits for options with maturities up to one month, underpinning their relevance in practical applications, and an analysis of the implied volatility skew in recent S&P500 options data shows it to be consistent with the infinite variation jump component of our models.

preprint2015arXiv

Small-time expansions for state-dependent local jump-diffusion models with infinite jump activity

In this article, we consider a Markov process X, starting from x and solving a stochastic differential equation, which is driven by a Brownian motion and an independent pure jump component exhibiting state-dependent jump intensity and infinite jump activity. A second order expansion is derived for the tail probability P[X(t)>x+y] in small time t, for y>0. As an application of this expansion and a suitable change of the underlying probability measure, a second order expansion, near expiration, for out-of-the-money European call option prices is obtained when the underlying stock price is modeled as the exponential of the jump-diffusion process X under the risk-neutral probability measure.

preprint2014arXiv

High-order short-time expansions for ATM option prices of exponential Lévy models

In the present work, a novel second-order approximation for ATM option prices is derived for a large class of exponential Lévy models with or without Brownian component. The results hereafter shed new light on the connection between both the volatility of the continuous component and the jump parameters and the behavior of ATM option prices near expiration. In the presence of a Brownian component, the second-order term, in time-$t$, is of the form $d_{2}\,t^{(3-Y)/2}$, with $d_{2}$ only depending on $Y$, the degree of jump activity, on $σ$, the volatility of the continuous component, and on an additional parameter controlling the intensity of the "small" jumps (regardless of their signs). This extends the well known result that the leading first-order term is $σt^{1/2}/\sqrt{2π}$. In contrast, under a pure-jump model, the dependence on $Y$ and on the separate intensities of negative and positive small jumps are already reflected in the leading term, which is of the form $d_{1}t^{1/Y}$. The second-order term is shown to be of the form $\tilde{d}_{2} t$ and, therefore, its order of decay turns out to be independent of $Y$. The asymptotic behavior of the corresponding Black-Scholes implied volatilities is also addressed. Our approach is sufficiently general to cover a wide class of Lévy processes which satisfy the latter property and whose Lévy densitiy can be closely approximated by a stable density near the origin. Our numerical results show that the first-order term typically exhibits rather poor performance and that the second-order term can significantly improve the approximation's accuracy, particularly in the absence of a Brownian component.

preprint2014arXiv

Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility

In Figueroa-López et al. (2013), a second order approximation for at-the-money (ATM) option prices is derived for a large class of exponential Lévy models, with or without a Brownian component. The purpose of this article is twofold. First, we relax the regularity conditions imposed in Figueroa-López et al. (2013) on the Lévy density to the weakest possible conditions for such an expansion to be well defined. Second, we show that the formulas extend both to the case of "close-to-the-money" strikes and to the case where the continuous Brownian component is replaced by an independent stochastic volatility process with leverage.

preprint2014arXiv

Small-time asymptotics of stopped Lévy bridges and simulation schemes with controlled bias

We characterize the small-time asymptotic behavior of the exit probability of a Lévy process out of a two-sided interval and of the law of its overshoot, conditionally on the terminal value of the process. The asymptotic expansions are given in the form of a first-order term and a precise computable error bound. As an important application of these formulas, we develop a novel adaptive discretization scheme for the Monte Carlo computation of functionals of killed Lévy processes with controlled bias. The considered functionals appear in several domains of mathematical finance (e.g., structural credit risk models, pricing of barrier options, and contingent convertible bonds) as well as in natural sciences. The proposed algorithm works by adding discretization points sampled from the Lévy bridge density to the skeleton of the process until the overall error for a given trajectory becomes smaller than the maximum tolerance given by the user.

preprint2014arXiv

Small-time expansions for local jump-diffusion models with infinite jump activity

We consider a Markov process $X$, which is the solution of a stochastic differential equation driven by a Lévy process $Z$ and an independent Wiener process $W$. Under some regularity conditions, including non-degeneracy of the diffusive and jump components of the process as well as smoothness of the Lévy density of $Z$ outside any neighborhood of the origin, we obtain a small-time second-order polynomial expansion for the tail distribution and the transition density of the process $X$. Our method of proof combines a recent regularizing technique for deriving the analog small-time expansions for a Lévy process with some new tail and density estimates for jump-diffusion processes with small jumps based on the theory of Malliavin calculus, flow of diffeomorphisms for SDEs, and time-reversibility. As an application, the leading term for out-of-the-money option prices in short maturity under a local jump-diffusion model is also derived.

preprint2012arXiv

High-order short-time expansions for ATM option prices under the CGMY model

The short-time asymptotic behavior of option prices for a variety of models with jumps has received much attention in recent years. In the present work, a novel second-order approximation for ATM option prices under the CGMY Lévy model is derived, and then extended to a model with an additional independent Brownian component. Our results shed light on the connection between both the volatility of the continuous component and the jump parameters and the behavior of ATM option prices near expiration. In case of an additional Brownian component, the second-order term, in time-t, is of the form $ d_{2} t^{(3-Y)/2}$, with the coefficient $d_{2}$ depending only on the overall jump intensity parameter C and the tail-heaviness parameter Y. This extends the known result that the leading term is $(σ/\sqrt{2π})t^{1/2}$, where $σ$ is the volatility of the continuous component. In contrast, under a pure-jump CGMY model, the dependence on the two parameters C and Y is already reflected in the leading term, which is of the form $d_{1} t^{1/Y}$. Information on the relative frequency of negative and positive jumps appears only in the second-order term, which is shown to be of the form $d_{2} t$ and whose order of decay turns out to be independent of Y. The third-order asymptotic behavior of the option prices as well as the asymptotic behavior of the corresponding Black-Scholes implied volatilities are also addressed. Our numerical results show that in most cases the second-order term significantly outperform the first-order approximation.

preprint2011arXiv

Sieve-based confidence intervals and bands for Lévy densities

The estimation of the Lévy density, the infinite-dimensional parameter controlling the jump dynamics of a Lévy process, is considered here under a discrete-sampling scheme. In this setting, the jumps are latent variables, the statistical properties of which can be assessed when the frequency and time horizon of observations increase to infinity at suitable rates. Nonparametric estimators for the Lévy density based on Grenander's method of sieves was proposed in Figueroa-López [IMS Lecture Notes 57 (2009) 117--146]. In this paper, central limit theorems for these sieve estimators, both pointwise and uniform on an interval away from the origin, are obtained, leading to pointwise confidence intervals and bands for the Lévy density. In the pointwise case, our estimators converge to the Lévy density at a rate that is arbitrarily close to the rate of the minimax risk of estimation on smooth Lévy densities. In the case of uniform bands and discrete regular sampling, our results are consistent with the case of density estimation, achieving a rate of order arbitrarily close to $\log^{-1/2}(n)\cdot n^{-1/3}$, where $n$ is the number of observations. The convergence rates are valid, provided that $s$ is smooth enough and that the time horizon $T_n$ and the dimension of the sieve are appropriately chosen in terms of $n$.

preprint2006arXiv

Risk bounds for the non-parametric estimation of Lévy processes

Estimation methods for the Lévy density of a Lévy process are developed under mild qualitative assumptions. A classical model selection approach made up of two steps is studied. The first step consists in the selection of a good estimator, from an approximating (finite-dimensional) linear model ${\mathcal{S}}$ for the true Lévy density. The second is a data-driven selection of a linear model ${\mathcal{S}}$, among a given collection $\{\mathcal{S}_m\}_{m\in {\mathcal{M}}}$, that approximately realizes the best trade-off between the error of estimation within ${\mathcal{S}}$ and the error incurred when approximating the true Lévy density by the linear model ${\mathcal{S}}$. Using recent concentration inequalities for functionals of Poisson integrals, a bound for the risk of estimation is obtained. As a byproduct, oracle inequalities and long-run asymptotics for spline estimators are derived. Even though the resulting underlying statistics are based on continuous time observations of the process, approximations based on high-frequency discrete-data can be easily devised.