Researcher profile

Lauri Viitasaari

Lauri Viitasaari contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 published item(s)

preprint2023arXiv

On Lamperti transformation and characterisations of discrete random fields

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.

preprint2022arXiv

Geometric Characterization of the Eyring-Kramers Formula

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.

preprint2022arXiv

On sharp rate of convergence for discretisation of integrals driven by fractional Brownian motions and related processes with discontinuous integrands

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.

preprint2021arXiv

Sobolev regularity of occupation measures and paths, variability and compositions

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.

preprint2020arXiv

Continuous time Gaussian process dynamical models in gene regulatory network inference

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.

preprint2020arXiv

Latent Model Extreme Value Index Estimation

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.

preprint2020arXiv

Local times and sample path properties of the Rosenblatt process

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.

preprint2020arXiv

Modeling temporally uncorrelated components for complex-valued stationary processes

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.

preprint2020arXiv

Quantitative normal approximations for the stochastic fractional heat equation

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.

preprint2020arXiv

Stochastic differential equations with noise perturbations and Wong-Zakai approximation of fractional Brownian motion

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)$.