BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Lauri Viitasaari appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
In this article we characterise discrete time stationary fields by difference equations involving stationary increment fields and self-similar fields. This gives connections between stationary fields, stationary increment fields and, through Lamperti transformation, self-similar fields. Our contribution is a natural generalisation of recently proved results covering the case of stationary processes.
In this paper we consider the mean transition time of an over-damped Brownian particle between local minima of a smooth potential. When the minima and saddles are non-degenerate this is in the low noise regime exactly characterized by the so called Eyring-Kramers law and gives the mean transition time as a quantity depending on the curvature of the minima and the saddle. In this paper we find an extension of the Eyring-Kramers law giving an upper bound on the mean transition time when both the minima/saddles are degenerate (flat) while at the same time covering multiple saddles at the same height. Our main contribution is a new sharp characterization of the capacity of two local minimas as a ratio of two geometric quantities, i.e., the smallest separating surface and the geodesic distance.
We consider equidistant approximations of stochastic integrals driven by Hölder continuous Gaussian processes of order $H>\frac12$ with discontinuous integrands involving bounded variation functions. We give exact rate of convergence in the $L^1$-distance and provide examples with different drivers. It turns out that the exact rate of convergence is proportional to $n^{1-2H}$ that is twice better compared to the best known results in the case of discontinuous integrands, and corresponds to the known rate in the case of smooth integrands. The novelty of our approach is that, instead of using multiplicative estimates for the integrals involved, we apply change of variables formula together with some facts on convex functions allowing us to compute expectations explicitly.
We prove a result on the fractional Sobolev regularity of composition of paths of low fractional Sobolev regularity with functions of bounded variation. The result relies on the notion of variability, proposed by us in the previous article [43, arXiv:2003.11698]. Here we work under relaxed hypotheses, formulated in terms of Sobolev norms, and we can allow discontinuous paths, which is new. The result applies to typical realizations of certain Gaussian or Lévy processes, and we use it to show the existence of Stieltjes type integrals involving compositions.
One of the focus areas of modern scientific research is to reveal mysteries related to genes and their interactions. The dynamic interactions between genes can be encoded into a gene regulatory network (GRN), which can be used to gain understanding on the genetic mechanisms behind observable phenotypes. GRN inference from time series data has recently been a focus area of systems biology. Due to low sampling frequency of the data, this is a notoriously difficult problem. We tackle the challenge by introducing the so-called continuous-time Gaussian process dynamical model, based on Gaussian process framework that has gained popularity in nonlinear regression problems arising in machine learning. The model dynamics are governed by a stochastic differential equation, where the dynamics function is modelled as a Gaussian process. We prove the existence and uniqueness of solutions of the stochastic differential equation. We derive the probability distribution for the Euler discretised trajectories and establish the convergence of the discretisation. We develop a GRN inference method called BINGO, based on the developed framework. BINGO is based on MCMC sampling of trajectories of the GPDM and estimating the hyperparameters of the covariance function of the Gaussian process. Using benchmark data examples, we show that BINGO is superior in dealing with poor time resolution and it is computationally feasible.
We propose a novel strategy for multivariate extreme value index estimation. In applications such as finance, volatility and risk present in the components of a multivariate time series are often driven by the same underlying factors, such as the subprime crisis in the US. To estimate the latent risk, we apply a two-stage procedure. First, a set of independent latent series is estimated using a method of latent variable analysis. Then, univariate risk measures are estimated individually for the latent series to assess their contribution to the overall risk. As our main theoretical contribution, we derive conditions under which the effect of the first step to the asymptotic behavior of the risk estimators is negligible. Simulations demonstrate the theory under both i.i.d. and dependent data, and an application into financial data illustrates the usefulness of the method in extracting joint sources of risk in practice.
Let $Z = (Z_t)_{t \geq 0}$ be the Rosenblatt process with Hurst index $H \in (1/2, 1)$. We prove joint continuity for the local time of $Z$, and establish Hölder conditions for the local time. These results are then used to study the irregularity of the sample paths of $Z$. Based on analogy with similar known results in the case of fractional Brownian motion, we believe our results are sharp. A main ingredient of our proof is a rather delicate spectral analysis of arbitrary linear combinations of integral operators, which arise from the representation of the Rosenblatt process as an element in the second chaos.
We consider a complex-valued linear mixture model, under discrete weakly stationary processes. We recover latent components of interest, which have undergone a linear mixing. We study asymptotic properties of a classical unmixing estimator, that is based on simultaneous diagonalization of the covariance matrix and an autocovariance matrix with lag $τ$. Our main contribution is that our asymptotic results can be applied to a large class of processes. In related literature, the processes are typically assumed to have weak correlations. We extend this class and consider the unmixing estimator under stronger dependency structures. In particular, we analyze the asymptotic behavior of the unmixing estimator under both, long- and short-range dependent complex-valued processes. Consequently, our theory covers unmixing estimators that converge slower than the usual $\sqrt{T}$ and unmixing estimators that produce non-Gaussian asymptotic distributions. The presented methodology is a powerful prepossessing tool and highly applicable in several fields of statistics. Complex-valued processes are frequently encountered in, for example, biomedical applications and signal processing. In addition, our approach can be applied to model real-valued problems that involve temporally uncorrelated pairs. These are encountered in, for example, applications in finance.
In this article we present a {\it quantitative} central limit theorem for the stochastic fractional heat equation driven by a a general Gaussian multiplicative noise, including the cases of space-time white noise and the white-colored noise with spatial covariance given by the Riesz kernel or a bounded integrable function. We show that the spatial average over a ball of radius $R$ converges, as $R$ tends to infinity, after suitable renormalization, towards a Gaussian limit in the total variation distance. We also provide a functional central limit theorem. As such, we extend recently proved similar results for stochastic heat equation to the case of the fractional Laplacian and to the case of general noise.
In this article we study effects that small perturbations in the noise have to the solution of differential equations driven by Hölder continuous functions of order $H>\frac12$. As an application, we consider stochastic differential equations driven by a fractional Brownian motion. We introduce a Wong--Zakai type stationary approximation to the fractional Brownian motions and prove that it converges in a suitable space. Moreover, we provide sharp results on the rate of convergence in the $p$-norm. Our stationary approximation is suitable for all values of $H\in (0,1)$.
In this paper, we study linear backward stochastic differential equations driven by a class of centered Gaussian non-martingales, including fractional Brownian motion with Hurst parameter $H\in (0,1)\setminus \{\frac12\}$. We show that, for every choice of deterministic coefficient functions, there is a square integrable terminal condition such that the equation has no solution.
This paper deals with stochastic integrals of form $\int_0^T f(X_u)d Y_u$ in a case where the function $f$ has discontinuities, and hence the process $f(X)$ is usually of unbounded $p$-variation for every $p\geq 1$. Consequently, integration theory introduced by Young or rough path theory introduced by Lyons cannot be applied directly. In this paper we prove the existence of such integrals in a pathwise sense provided that $X$ and $Y$ have suitably regular paths together with some minor additional assumptions. In many cases of interest, our results extend the celebrated results by Young.
It was shown in Mishura et al. (Stochastic Process. Appl. 123 (2013) 2353-2369), that any random variable can be represented as improper pathwise integral with respect to fractional Brownian motion. In this paper, we extend this result to cover a wide class of Gaussian processes. In particular, we consider a wide class of processes that are Hölder continuous of order $α>1/2$ and show that only local properties of the covariance function play role for such results.
We first study the drift parameter estimation of the fractional Ornstein-Uhlenbeck process (fOU) with periodic mean for every $\frac{1}{2}<H<1$. More precisely, we extend the consistency proved in \cite{DFW} for $\frac{1}{2}<H<\frac{3}{4}$ to the strong consistency for any $\frac{1}{2}<H<1$ on the one hand, and on the other, we also discuss the asymptotic normality given in \cite{DFW}. In the second main part of the paper, we study the strong consistency and the asymptotic normality of the fOU of the second kind with periodic mean for any $\frac{1}{2}<H<1$.
We study the Langevin equation with stationary-increment Gaussian noise. We show the strong consistency and the asymptotic normality with Berry--Esseen bound of the so-called alternative estimator of the mean reversion parameter. The conditions and results are stated in terms of the variance function of the noise. We consider both the case of continuous and discrete observations. As examples we consider fractional and bifractional Ornstein--Uhlenbeck processes. Finally, we discuss the maximum likelihood and the least squares estimators.
We show that every separable Gaussian process with integrable variance function admits a Fredholm representation with respect to a Brownian motion. We extend the Fredholm representation to a transfer principle and develop stochastic analysis by using it. We show the convenience of the Fredholm representation by giving applications to equivalence in law, bridges, series expansions, stochastic differential equations and maximum likelihood estimations.
We study asymptotic normality of the randomized periodogram estimator of quadratic variation in the mixed Brownian--fractional Brownian model. In the semimartingale case, that is, where the Hurst parameter $H$ of the fractional part satisfies $H\in(3/4,1)$, the central limit theorem holds. In the nonsemimartingale case, that is, where $H\in(1/2,3/4]$, the convergence toward the normal distribution with a nonzero mean still holds if $H=3/4$, whereas for the other values, that is, $H\in(1/2,3/4)$, the central convergence does not take place. We also provide Berry--Esseen estimates for the estimator.
We show that every multiparameter Gaussian process with integrable variance function admits a Wiener integral representation of Fredholm type with respect to the Brownian sheet. The Fredholm kernel in the representation can be constructed as the unique symmetric square root of the covariance. We analyze the equivalence of multiparameter Gaussian processes by using the Fredholm representation and show how to construct series expansions for multiparameter Gaussian processes by using the Fredholm kernel.
Let $W_t$ be a standard Brownian motion. It is well-known that the Langevin equation $d U_t = -θU_td t + d W_t$ defines a stationary process called Ornstein-Uhlenbeck process. Furthermore, Langevin equation can be used to construct other stationary processes by replacing Brownian motion $W_t$ with some other process $G$ with stationary increments. In this article we prove that the converse also holds and all continuous stationary processes arise from a Langevin equation with certain noise $G=G_θ$. Discrete analogies of our results are given and applications are discussed.
Quadratic variations of Gaussian processes play important role in both stochastic analysis and in applications such as estimation of model parameters, and for this reason the topic has been extensively studied in the literature. In this article we study the problem for general Gaussian processes and we provide sufficient and necessary conditions for different types of convergence which include convergence in probability, almost sure convergence, $L^p$-convergence as well as convergence in law. Furthermore, we study general Gaussian vectors from which different interesting cases including first or second order quadratic variations can be studied by appropriate choice of the underlying vector. Finally, we provide a practical and simple approach to attack the problem which simplifies the existing methodology considerably.
We study integral representations of random variables with respect to general Hölder continuous processes and with respect to two particular cases; fractional Brownian motion and mixed fractional Brownian motion. We prove that arbitrary random variable can be represented as an improper integral, and that the stochastic integral can have any distribution. If in addition the random variable is a final value of an adapted Hölder continuous process, then it can be represented as a proper integral. It is also shown that in the particular case of mixed fractional Brownian motion, any adapted random variable can be represented as a proper integral.
We show that if a random variable is a final value of an adapted Holder continuous process, then it can be represented as a stochastic integral with respect to fractional Brownian motion, and the integrand is an adapted process, continuous up to the final point.
In this article, we consider European options of type $h(X^1_T, X^2_T,\ldots, X^n_T)$ depending on several underlying assets. We study how such options can be valued in terms of simple vanilla options in non-specified market models. We consider different approaches related to static hedging and derive several pricing formulas for a wide class of payoff functions $h:\R_+^n\rightarrow \R$. We also give new relations between prices of different options both in one dimensional and multidimensional case.
The continuity of Gaussian processes is extensively studied topic and it culminates in the Talagrand's notion of majorizing measures that gives complicated necessary and sufficient conditions. In this note we study the Hölder continuity of Gaussian processes. It turns out that necessary and sufficient conditions can be stated in a simple form that is a variant of the celebrated Kolmogorov-Čentsov condition.
In this note, we consider European options of type $h(X^1_T, X^2_T,\ldots, X^n_T)$ depending on several underlying assets. We give a multidimensional version of the result of Breeden and Litzenberger \cite{Breeden} on the relation between derivatives of the call price and the risk-neutral density of the underlying asset. The pricing measure is assumed to be absolutely continuous with respect to the Lebesgue measure on the state space.
Fractional Ornstein-Uhlenbeck process of the second kind $(\text{fOU}_{2})$ is solution of the Langevin equation $\mathrm{d}X_t = -θX_t\,\mathrm{d}t+\mathrm{d}Y_t^{(1)}, \ θ>0$ with Gaussian driving noise $ Y_t^{(1)} := \int^t_0 e^{-s} \,\mathrm{d}B_{a_s}$, where $ a_t= H e^{\frac{t}{H}}$ and $B$ is a fractional Brownian motion with Hurst parameter $H \in (0,1)$. In this article, we consider the case $H>\frac{1}{2}$. Then using the ergodicity of $\text{fOU}_{2}$ process, we construct consistent estimators of drift parameter $θ$ based on discrete observations in two possible cases: $(i)$ the Hurst parameter $H$ is known and $(ii)$ the Hurst parameter $H$ is unknown. Moreover, using Malliavin calculus technique, we prove central limit theorems for our estimators which is valid for the whole range $H \in (\frac{1}{2},1)$.
We prove the Ito-Tanaka formula and the existence of pathwise stochastic integrals for a wide class of Gaussian processes. Motivated by financial applications, we define the stochastic integrals as forward-type pathwise integrals introduced by Föllmer and as pathwise generalized Lebesgue-Stieltjes integrals introduced by Zähle. As an application, we illustrate the importance of Ito-Tanaka formula for pricing and hedging of financial derivatives.
In this article we study existence of pathwise stochastic integrals with respect to a general class of $n$-dimensional Gaussian processes and a wide class of adapted integrands. More precisely, we study integrands which are functions that are of locally bounded variation with respect to all variables. Moreover, multidimensional Itô formula is derived.
In this article, we study the rate of convergence of prices when a model is approximated by some simplified model. We also provide a method how explicit error formula for more general options can be obtained if such formula is available for digital option prices. We illustrate our results by considering convergence of binomial prices to Black-Scholes prices. We also consider smooth convergence in which the approximation does not oscillate for general class of payoff functions.
There exist several methods how more general options can be priced with call prices. In this article, we extend these results to cover a wider class of options and market models. In particular, we introduce a new pricing formula which can be used to price more general options if prices for call options and digital options are known for every strike price. Moreover, we derive similar results for barrier type options. As a consequence, we obtain a static hedging for general options in the general class of models. Our result can be utilised in several significant applications. As a simple example, we derive an upper bound for the value of a general American option with convex payoff and characterise conditions under which the value of this option equals to the value of the corresponding European option.
In this article, an uniform discretization of stochastic integrals $\int_{0}^{1} f'_-(B_t)\ud B_t$, with respect to fractional Brownian motion with Hurst parameter $H \in (1/2,1)$, for a large class of convex functions $f$ is considered. In Statistics & Decisions, 27, 129-143, for any convex function $f$, the almost sure convergence of uniform discretization to such stochastic integral is proved. Here we prove $L^r$- convergence of uniform discretization to stochastic integral. In addition, we obtain a rate of convergence. It turns out that the rate of convergence can be brought as closely as possible to $H - 1/2$.