Source author record

Salvador Ortiz-Latorre

Salvador Ortiz-Latorre appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.PR q-fin.PR math.NA Numerical Analysis Artificial Intelligence Machine Learning q-fin.CP q-fin.RM

Catalog footprint

What is connected

8works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2023arXiv

SPDE bridges with observation noise and their spatial approximation

This paper introduces SPDE bridges with observation noise and contains an analysis of their spatially semidiscrete approximations. The SPDEs are considered in the form of mild solutions in an abstract Hilbert space framework suitable for parabolic equations. They are assumed to be linear with additive noise in the form of a cylindrical Wiener process. The observational noise is also cylindrical and SPDE bridges are formulated via conditional distributions of Gaussian random variables in Hilbert spaces. A general framework for the spatial discretization of these bridge processes is introduced. Explicit convergence rates are derived for a spectral and a finite element based method. It is shown that for sufficiently rough observation noise, the rates are essentially the same as those of the corresponding discretization of the original SPDE.

preprint2022arXiv

An application of the splitting-up method for the computation of a neural network representation for the solution for the filtering equations

The filtering equations govern the evolution of the conditional distribution of a signal process given partial, and possibly noisy, observations arriving sequentially in time. Their numerical approximation plays a central role in many real-life applications, including numerical weather prediction, finance and engineering. One of the classical approaches to approximate the solution of the filtering equations is to use a PDE inspired method, called the splitting-up method, initiated by Gyongy, Krylov, LeGland, among other contributors. This method, and other PDE based approaches, have particular applicability for solving low-dimensional problems. In this work we combine this method with a neural network representation. The new methodology is used to produce an approximation of the unnormalised conditional distribution of the signal process. We further develop a recursive normalisation procedure to recover the normalised conditional distribution of the signal process. The new scheme can be iterated over multiple time steps whilst keeping its asymptotic unbiasedness property intact. We test the neural network approximations with numerical approximation results for the Kalman and Benes filter.

preprint2021arXiv

Pathwise approximations for the solution of the non-linear filtering problem

We consider high order approximations of the solution of the stochastic filtering problem, derive their pathwise representation in the spirit of the earlier work of Clark and Davis and prove their robustness property. In particular, we show that the high order discretised filtering functionals can be represented by Lipschitz continuous functions defined on the observation path space. This property is important from the practical point of view as it is in fact the pathwise version of the filtering functional that is sought in numerical applications. Moreover, the pathwise viewpoint will be a stepping stone into the rigorous development of machine learning methods for the filtering problem. This work is a continuation of a recent work by two of the authors where a discretisation of the solution of the filtering problem of arbitrary order has been established. We expand the previous work by showing that robust approximations can be derived from the discretisations therein.

preprint2020arXiv

A decomposition formula for fractional Heston jump diffusion models

We present an option pricing formula for European options in a stochastic volatility model. In particular, the volatility process is defined using a fractional integral of a diffusion process and both the stock price and the volatility processes have jumps in order to capture the market effect known as leverage effect. We show how to compute a martingale representation for the volatility process. Finally, using Itô calculus for processes with discontinuous trajectories, we develop a first order approximation formula for option prices. There are two main advantages in the usage of such approximating formulas to traditional pricing methods. First, to improve computational effciency, and second, to have a deeper understanding of the option price changes in terms of changes in the model parameters.

preprint2020arXiv

Variance and interest rate risk in unit-linked insurance policies

One of the risks derived from selling long term policies that any insurance company has, arises from interest rates. In this paper we consider a general class of stochastic volatility models written in forward variance form. We also deal with stochastic interest rates to obtain the risk-free price for unit-linked life insurance contracts, as well as providing a perfect hedging strategy by completing the market. We conclude with a simulation experiment, where we price unit-linked policies using Norwegian mortality rates. In addition we compare prices for the classical Black-Scholes model against the Heston stochastic volatility model with a Vasicek interest rate model.

preprint2014arXiv

A change of measure preserving the affine structure in the BNS model for commodity markets

For a commodity spot price dynamics given by an Ornstein-Uhlenbeck process with Barndorff-Nielsen and Shephard stochastic volatility, we price forwards using a class of pricing measures that simultaneously allow for change of level and speed in the mean reversion of both the price and the volatility. The risk premium is derived in the case of arithmetic and geometric spot price processes, and it is demonstrated that we can provide flexible shapes that is typically observed in energy markets. In particular, our pricing measure preserves the affine model structure and decomposes into a price and volatility risk premium, and in the geometric spot price model we need to resort to a detailed analysis of a system of Riccati equations, for which we show existence and uniqueness of solution and asymptotic properties that explains the possible risk premium profiles. Among the typical shapes, the risk premium allows for a stochastic change of sign, and can attain positive values in the short end of the forward market and negative in the long end.

preprint2013arXiv

A pricing measure to explain the risk premium in power markets

In electricity markets, it is sensible to use a two-factor model with mean reversion for spot prices. One of the factors is an Ornstein-Uhlenbeck (OU) process driven by a Brownian motion and accounts for the small variations. The other factor is an OU process driven by a pure jump Lévy process and models the characteristic spikes observed in such markets. When it comes to pricing, a popular choice of pricing measure is given by the Esscher transform that preserves the probabilistic structure of the driving Lévy processes, while changing the levels of mean reversion. Using this choice one can generate stochastic risk premiums (in geometric spot models) but with (deterministically) changing sign. In this paper we introduce a pricing change of measure, which is an extension of the Esscher transform. With this new change of measure we also can slow down the speed of mean reversion and generate stochastic risk premiums with stochastic non constant sign, even in arithmetic spot models. In particular, we can generate risk profiles with positive values in the short end of the forward curve and negative values in the long end. Finally, our pricing measure allows us to have a stationary spot dynamics while still having randomly fluctuating forward prices for contracts far from maturity.

preprint2012arXiv

Optimal simulation schemes for Lévy driven stochastic differential equations

We consider a general class of high order weak approximation schemes for stochastic differential equations driven by Lévy processes with infinite activity. These schemes combine a compound Poisson approximation for the jump part of the Lévy process with a high order scheme for the Brownian driven component, applied between the jump times. The overall approximation is analyzed using a stochastic splitting argument. The resulting error bound involves separate contributions of the compound Poisson approximation and of the discretization scheme for the Brownian part, and allows, on one hand, to balance the two contributions in order to minimize the computational time, and on the other hand, to study the optimal design of the approximating compound Poisson process. For driving processes whose Lévy measure explodes near zero in a regularly varying way, this procedure allows to construct discretization schemes with arbitrary order of convergence.