BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

John Schoenmakers

John Schoenmakers appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.PRq-fin.CPq-fin.PRmath.NAmath.OCmath.STStatistics TheoryNumerical AnalysisMachine LearningMethodologyq-fin.MF

Catalog footprint

What is connected

14works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

14 published item(s)

preprint2022arXiv

Reinforced optimal control

Least squares Monte Carlo methods are a popular numerical approximation method for solving stochastic control problems. Based on dynamic programming, their key feature is the approximation of the conditional expectation of future rewards by linear least squares regression. Hence, the choice of basis functions is crucial for the accuracy of the method. Earlier work by some of us [Belomestny, Schoenmakers, Spokoiny, Zharkynbay. Commun.~Math.~Sci., 18(1):109-121, 2020](arXiv:1808.02341) proposes to reinforce the basis functions in the case of optimal stopping problems by already computed value functions for later times, thereby considerably improving the accuracy with limited additional computational cost. We extend the reinforced regression method to a general class of stochastic control problems, while considerably improving the method's efficiency, as demonstrated by substantial numerical examples as well as theoretical analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

From optimal martingales to randomized dual optimal stopping

In this article we study and classify optimal martingales in the dual formulation of optimal stopping problems. In this respect we distinguish between weakly optimal and surely optimal martingales. It is shown that the family of weakly optimal and surely optimal martingales may be quite large. On the other hand it is shown that the Doob-martingale, that is, the martingale part of the Snell envelope, is in a certain sense the most robust surely optimal martingale under random perturbations. This new insight leads to a novel randomized dual martingale minimization algorithm that doesn't require nested simulation. As a main feature, in a possibly large family of optimal martingales the algorithm efficiently selects a martingale that is as close as possible to the Doob martingale. As a result, one obtains the dual upper bound for the optimal stopping problem with low variance.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

From rough path estimates to multilevel Monte Carlo

New classes of stochastic differential equations can now be studied using rough path theory (e.g. Lyons et al. [LCL07] or Friz--Hairer [FH14]). In this paper we investigate, from a numerical analysis point of view, stochastic differential equations driven by Gaussian noise in the aforementioned sense. Our focus lies on numerical implementations, and more specifically on the saving possible via multilevel methods. Our analysis relies on a subtle combination of pathwise estimates, Gaussian concentration, and multilevel ideas. Numerical examples are given which both illustrate and confirm our findings.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Affine LIBOR models with multiple curves: theory, examples and calibration

We introduce a multiple curve framework that combines tractable dynamics and semi-analytic pricing formulas with positive interest rates and basis spreads. Negatives rates and positive spreads can also be accommodated in this framework. The dynamics of OIS and LIBOR rates are specified following the methodology of the affine LIBOR models and are driven by the wide and flexible class of affine processes. The affine property is preserved under forward measures, which allows us to derive Fourier pricing formulas for caps, swaptions and basis swaptions. A model specification with dependent LIBOR rates is developed, that allows for an efficient and accurate calibration to a system of caplet prices.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Forward-reverse EM algorithm for Markov chains: convergence and numerical analysis

We develop a forward-reverse EM (FREM) algorithm for estimating parameters that determine the dynamics of a discrete time Markov chain evolving through a certain measurable state space. As a key tool for the construction of the FREM method we develop forward-reverse representations for Markov chains conditioned on a certain terminal state. These representations may be considered as an extension of the earlier work Bayer and Schoenmakers [2013] on conditional diffusions. We proof almost sure convergence of our algorithm for a Markov chain model with curved exponential family structure. On the numerical side we give a complexity analysis of the forward-reverse algorithm by deriving its expected cost. Two application examples are discuss to demonstrate the scope of possible applications ranging from models based on continuous time processes to discrete time Markov chain models.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Generalized Post-Widder inversion formula with application to statistics

In this work we derive an inversion formula for the Laplace transform of a density observed on a curve in the complex domain, which generalizes the well known Post-Widder formula. We establish convergence of our inversion method and derive the corresponding convergence rates for the case of a Laplace transform of a smooth density. As an application we consider the problem of statistical inference for variance-mean mixture models. We construct a nonparametric estimator for the mixing density based on the generalized Post-Widder formula, derive bounds for its root mean square error and give a brief numerical example.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Option pricing in affine generalized Merton models

In this article we consider affine generalizations of the Merton jump diffusion model [Merton, J. Fin. Econ., 1976] and the respective pricing of European options. On the one hand, the Brownian motion part in the Merton model may be generalized to a log-Heston model, and on the other hand, the jump part may be generalized to an affine process with possibly state dependent jumps. While the characteristic function of the log-Heston component is known in closed form, the characteristic function of the second component may be unknown explicitly. For the latter component we propose an approximation procedure based on the method introduced in [Belomestny et al., J. Func. Anal., 2009]. We conclude with some numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulation of forward-reverse stochastic representations for conditional diffusions

In this paper we derive stochastic representations for the finite dimensional distributions of a multidimensional diffusion on a fixed time interval, conditioned on the terminal state. The conditioning can be with respect to a fixed point or more generally with respect to some subset. The representations rely on a reverse process connected with the given (forward) diffusion as introduced in Milstein, Schoenmakers and Spokoiny [Bernoulli 10 (2004) 281-312] in the context of a forward-reverse transition density estimator. The corresponding Monte Carlo estimators have essentially root-$N$ accuracy, and hence they do not suffer from the curse of dimensionality. We provide a detailed convergence analysis and give a numerical example involving the realized variance in a stochastic volatility asset model conditioned on a fixed terminal value of the asset.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Statistical Skorohod embedding problem and its generalizations

Given a Lévy process $L$, we consider the so-called statistical Skorohod embedding problem of recovering the distribution of an independent random time $T$ based on i.i.d. sample from $L_{T}.$ Our approach is based on the genuine use of the Mellin and Laplace transforms. We propose a consistent estimator for the density of $T,$ derive its convergence rates and prove their optimality. It turns out that the convergence rates heavily depend on the decay of the Mellin transform of $T.$ We also consider the application of our results to the problem of statistical inference for variance-mean mixture models and for time-changed Lévy processes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Uniform approximation of the Cox-Ingersoll-Ross process

The Doss-Sussmann (DS) approach is used for uniform simulation of the Cox-Ingersoll-Ross (CIR) process. The DS formalism allows to express trajectories of the CIR process through solutions of some ordinary differential equation (ODE) depending on realizations of a Wiener process involved. By simulating the first-passage times of the increments of the Wiener process to the boundary of an interval and solving the ODE, we uniformly approximate the trajectories of the CIR process. In this respect special attention is payed to simulation of trajectories near zero. From a conceptual point of view the proposed method gives a better quality of approximation (from a path-wise point of view) than standard, or even exact simulation of the SDE at some discrete time grid.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Efficient and accurate log-Lévy approximations to Lévy driven LIBOR models

The LIBOR market model is very popular for pricing interest rate derivatives, but is known to have several pitfalls. In addition, if the model is driven by a jump process, then the complexity of the drift term is growing exponentially fast (as a function of the tenor length). In this work, we consider a Lévy-driven LIBOR model and aim at developing accurate and efficient log-Lévy approximations for the dynamics of the rates. The approximations are based on truncation of the drift term and Picard approximation of suitable processes. Numerical experiments for FRAs, caps, swaptions and sticky ratchet caps show that the approximations perform very well. In addition, we also consider the log-Lévy approximation of annuities, which offers good approximations for high volatility regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Optimal dual martingales, their analysis and application to new algorithms for Bermudan products

In this paper we introduce and study the concept of optimal and surely optimal dual martingales in the context of dual valuation of Bermudan options, and outline the development of new algorithms in this context. We provide a characterization theorem, a theorem which gives conditions for a martingale to be surely optimal, and a stability theorem concerning martingales which are near to be surely optimal in a sense. Guided by these results we develop a framework of backward algorithms for constructing such a martingale. In turn this martingale may then be utilized for computing an upper bound of the Bermudan product. The methodology is pure dual in the sense that it doesn't require certain input approximations to the Snell envelope. In an Itô-Lévy environment we outline a particular regression based backward algorithm which allows for computing dual upper bounds without nested Monte Carlo simulation. Moreover, as a by-product this algorithm also provides approximations to the continuation values of the product, which in turn determine a stopping policy. Hence, we may obtain lower bounds at the same time. In a first numerical study we demonstrate the backward dual regression algorithm in a Wiener environment at well known benchmark examples. It turns out that the method is at least comparable to the one in Belomestny et. al. (2009) regarding accuracy, but regarding computational robustness there are even several advantages.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dual representations for general multiple stopping problems

In this paper, we study the dual representation for generalized multiple stopping problems, hence the pricing problem of general multiple exercise options. We derive a dual representation which allows for cashflows which are subject to volume constraints modeled by integer valued adapted processes and refraction periods modeled by stopping times. As such, this extends the works by Schoenmakers (2010), Bender (2011a), Bender (2011b), Aleksandrov and Hambly (2010), and Meinshausen and Hambly (2004) on multiple exercise options, which either take into consideration a refraction period or volume constraints, but not both simultaneously. We also allow more flexible cashflow structures than the additive structure in the above references. For example some exponential utility problems are covered by our setting. We supplement the theoretical results with an explicit Monte Carlo algorithm for constructing confidence intervals for the price of multiple exercise options and exemplify it by a numerical study on the pricing of a swing option in an electricity market.

0 citations0 reviewsNo code linkedOpen access