arXiv research

New methods for $Z$-transform inversion and Wiener-Hopf factorization.

This paper considers the valuation of exotic path-dependent options in Lévy models, in particular options on the supremum and the infimum of the asset price process. Using the Wiener--Hopf factorization, we derive expressions for the analytically extended characteristic function of the supremum and the infimum of a Lév…

New pricing methods for $α$-quantile and early-exercise options using Spitzer identities.

Fast method developed for pricing barrier options and joint Lévy process distributions.

This paper describes a new method of bond portfolio optimization based on stochastic string models of correlation structure in bond returns. The paper shows how to approximate correlation function of bond returns, compute the optimal portfolio allocation using Wiener-Hopf factorization, and check whether a collection o…

In this note we apply the recently established Wiener-Hopf Monte Carlo (WHMC) simulation technique for Levy processes from Kuznetsov et al. [17] to path functionals, in particular first passage times, overshoots, undershoots and the last maximum before the passage time. Such functionals have many applications, for inst…

We show that the jumps correlation matrix of a multivariate Hawkes process is related to the Hawkes kernel matrix through a system of Wiener-Hopf integral equations. A Wiener-Hopf argument allows one to prove that this system (in which the kernel matrix is the unknown) possesses a unique causal solution and consequentl…

Efficient methods for Lévy models using SINH-regular processes.

In this paper we consider the problem of pricing a perpetual American put option in an exponential regime-switching Lévy model. For the case of the (dense) class of phase-type jumps and finitely many regimes we derive an explicit expression for the value function. The solution of the corresponding first passage problem…

In this paper we develop an algorithm to calculate the prices and Greeks of barrier options in a hyper-exponential additive model with piecewise constant parameters. We obtain an explicit semi-analytical expression for the first-passage probability. The solution rests on a randomization and an explicit matrix Wiener-Ho…

The Wiener-Hopf factorization is obtained in closed form for a phase type approximation to the CGMY Lévy process. This allows, for the approximation, exact computation of first passage times to barrier levels via Laplace transform inversion. Calibration of the CGMY model to market option prices defines the risk neutral…

The paper derives formulas for option pricing and random walk expectations.

Lewis and Mordecki have computed the Wiener-Hopf factorization of a Lévy process whose restriction on $]0,+\infty[$ of their Lévy measure has a rational Laplace transform. That allows to compute the distribution of $(X_t,\inf_{0\leq s\leq t}X_s)$. For the same class of Lévy processes, we compute the distribution of $ (…

We present a numerical scheme to calculate fluctuation identities for exponential Lévy processes in the continuous monitoring case. This includes the Spitzer identities for touching a single upper or lower barrier, and the more difficult case of the two-barriers exit problem. These identities are given in the Fourier-L…

In this paper, we argue that, once the costs of maintaining the hedging portfolio are properly taken into account, semi-static portfolios should more properly be thought of as separate classes of derivatives, with non-trivial, model-dependent payoff structures. We derive new integral representations for payoffs of exot…

In equity and foreign exchange markets the risk-neutral dynamics of the underlying asset are commonly represented by stochastic volatility models with jumps. In this paper we consider a dense subclass of such models and develop analytically tractable formulae for the prices of a range of first-generation exotic derivat…

We derive a new equation for the optimal investment boundary of a general irreversible investment problem under exponential Lévy uncertainty. The problem is set as an infinite time-horizon, two-dimensional degenerate singular stochastic control problem. In line with the results recently obtained in a diffusive setting,…

Improved iterative methods for risk parity portfolio weights.

We describe a novel optimization method for finite sums (such as empirical risk minimization problems) building on the recently introduced SAGA method. Our method achieves an accelerated convergence rate on strongly convex smooth problems. Our method has only one parameter (a step size), and is radically simpler than o…

A new method combines Laplace and Variational Bayes for scalable inference.

Unified framework for model explanation methods based on feature removal.

This work reviews and evaluates methods for predicting prediction intervals in regression problems.

Derives kernel PCA with Nyström method for scalability.

In this paper, the author considers the numerical computation of CVA for large systems by Mote Carlo methods. He introduces two types of stochastic mesh methods for the computations of CVA. In the first method, stochastic mesh method is used to obtain the future value of the derivative contracts. In the second method, …

New method combines spectral and sparse methods for Gaussian processes.

A comprehensive benchmark of 15 scRNA-seq imputation methods across various datasets and analyses.

New methods using natural gradient for structured optimization.

Recently, {\it stochastic momentum} methods have been widely adopted in training deep neural networks. However, their convergence analysis is still underexplored at the moment, in particular for non-convex optimization. This paper fills the gap between practice and theory by developing a basic convergence analysis of t…

We investigate methods for pricing American options under the variance gamma model. The variance gamma process is a pure jump process which is constructed by replacing the calendar time by the gamma time in a Brownian motion with drift, which makes it a time-changed Brownian motion. In general, the finite difference me…

A new method speeds up deep neural network training.

We propose a new stochastic dual coordinate ascent technique that can be applied to a wide range of regularized learning problems. Our method is based on Alternating Direction Multiplier Method (ADMM) to deal with complex regularization functions such as structured regularizations. Although the original ADMM is a batch…

NCG methods improve shape optimization efficiency.

We propose two localized Radial Basis Function (RBF) methods, the Radial Basis Function Partition of Unity method (RBF-PUM) and the Radial Basis Function generated Finite Differences method (RBF-FD), for solving financial derivative pricing problems arising from market models with multiple stochastic factors. We demons…

Proposes UTC method for stock price prediction with uncertainty quantification.

Survey of spectral, probabilistic, and deep metric learning methods.

Various approaches to gene selection for cancer classification based on microarray data can be found in the literature and they may be grouped into two categories: univariate methods and multivariate methods. Univariate methods look at each gene in the data in isolation from others. They measure the contribution of a p…

A novel weighted feature selection method using fuzzy sets improves classification accuracy and stability.

New method improves accuracy in computing implied volatility.

The paper examines Wiener process for LID estimation methods.

We evaluate how modern outlier detection methods perform in identifying outliers in e-commerce conversion rate data. Based on the limitations identified, we then present a novel method to detect outliers in e-commerce conversion rate. This unsupervised method is made more business relevant by letting it automatically a…

Spectral methods of moments provide a powerful tool for learning the parameters of latent variable models. Despite their theoretical appeal, the applicability of these methods to real data is still limited due to a lack of robustness to model misspecification. In this paper we present a hierarchical approach to methods…

We propose an optimization method for minimizing the finite sums of smooth convex functions. Our method incorporates an accelerated gradient descent (AGD) and a stochastic variance reduction gradient (SVRG) in a mini-batch setting. Unlike SVRG, our method can be directly applied to non-strongly and strongly convex prob…

We discuss the relevance of the recent Machine Learning (ML) literature for economics and econometrics. First we discuss the differences in goals, methods and settings between the ML literature and the traditional econometrics and statistics literatures. Then we discuss some specific methods from the machine learning l…

R package for counterfactual explanation methods.

We generalize Newton-type methods for minimizing smooth functions to handle a sum of two convex functions: a smooth function and a nonsmooth function with a simple proximal mapping. We show that the resulting proximal Newton-type methods inherit the desirable convergence behavior of Newton-type methods for minimizing s…

The field of medical image reconstruction has seen roughly four types of methods. The first type tended to be analytical methods, such as filtered back-projection (FBP) for X-ray computed tomography (CT) and the inverse Fourier transform for magnetic resonance imaging (MRI), based on simple mathematical models for the …

Derives new optimization methods using variational integrators.

A new method for faster optimization in high dimensions.