BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Yuliya Mishura 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
Asymptotic expansion is presented for an estimator of the Hurst coefficient of a fractional Brownian motion. For this, a recently developed theory of asymptotic expansion of the distribution of Wiener functionals is applied. The effects of the asymptotic expansion are demonstrated by numerical studies.
In this paper, we analyze the drift-implicit (or backward) Euler numerical scheme for a class of stochastic differential equations with unbounded drift driven by an arbitrary $λ$-Hölder continuous process, $λ\in(0,1)$. We prove that, under some mild moment assumptions on the Hölder constant of the noise, the $L^r(Ω;L^\infty([0,T]))$-rate of convergence is equal to $λ$. To exemplify, we consider numerical schemes for the generalized Cox--Ingersoll-Ross and Tsallis--Stariolo--Borland models. The results are illustrated by simulations.
We study convexity properties of Rényi entropy as function of $α>0$ on finite alphabets. We also describe robustness of the Rényi entropy on finite alphabets, and it turns out that the rate of respective convergence depends on initial alphabet. We establish convergence of the disturbed entropy when the initial distribution is uniform but the number of events increases to $\infty$ and prove that limit of Rényi entropy of binomial distribution is equal to Rényi entropy of Poisson distribution.
The paper contains sufficient conditions on the function $f$ and the stochastic process $X$ that supply the rate of divergence of the integral functional $\int_0^Tf(X_t)^2dt$ at the rate $T^{1-ε}$ as $T\to\infty$ for every $ε>0$. These conditions include so called small ball estimates which are discussed in detail. Statistical applications are provided.
The paper is devoted to the approximate solutions of the Fredholm integral equations of the second kind with the weak singular kernel that can have additional singularity in the numerator. We describe two problems that lead to such equations. They are the problem of minimization of small deviation and the entropy minimization problem. Both of them appear when considering dynamical system involving mixed fractional Brownian motion. In order to deal with the kernel with additional singularity applying well-known methods for weakly singular kernels, we prove the theorem on the approximation of solution of integral equation with the kernel containing additional singularity by the solutions of the integral equations whose kernels are weakly singular but the numerator is continuous. We demonstrate numerically how our methods work being applied to our specific integral equations.
We study the asymptotic behaviour of a properly normalized time changed Wiener processes. The time change reflects the fact that we consider the Laplace operator (which generates a Wiener process) multiplied by a possibly degenerate state-space dependent intensity $λ(x)$. Applying a functional limit theorem for the superposition of stochastic processes, we prove functional limit theorems for the normalized time changed Wiener process. The normalization depends on the asymptotic behaviour of the intensity function $λ$. One of the possible limits is a skew Brownian motion.
We study boundary non-crossing probabilities $$ P_{f,u} := \mathrm P\big(\forall t\in \mathbb T\ X_t + f(t)\le u(t)\big) $$ for continuous centered Gaussian process $X$ indexed by some arbitrary compact separable metric space $\mathbb T$. We obtain both upper and lower bounds for $P_{f,u}$. The bounds are matching in the sense that they lead to precise logarithmic asymptotics for the large-drift case $P_{y f,u}$, $y \to+\infty$, which are two-term approximations (up to $o(y)$). The asymptotics are formulated in terms of the solution $\tilde f$ to the constrained optimization problem $$ \|h\|_{\mathbb H_X}\to \min, \quad h\in \mathbb H_X, h\ge f $$ in the reproducing kernel Hilbert space $\mathbb H_X$ of $X$. Several applications of the results are further presented.
We find the best approximation of the fractional Brownian motion with the Hurst index $H\in (0,1/2)$ by Gaussian martingales of the form $\int _0^ts^γdW_s$, where $W$ is a Wiener process, $γ>0$.
In this paper the fractional Cox-Ingersoll-Ross process on $\mathbb{R}_+$ for $H<1/2$ is defined as a square of a pointwise limit of the processes $Y_{\varepsilon}$, satisfying the SDE of the form $d Y_{\varepsilon}(t)=( \frac{k}{ Y_{\varepsilon}(t)\mathbb{1}_{\{ Y_{\varepsilon}(t)>0\}}+\varepsilon}-a Y_{\varepsilon}(t))dt+σdB^H(t)$, as $\varepsilon\downarrow0$. Properties of such limit process are considered. SDE for both the limit process and the fractional Cox-Ingersoll-Ross process are obtained.
In this paper we find fractional Riemann-Liouville derivatives for the Takagi-Landsberg functions. Moreover, we introduce their generalizations called weighted Takagi-Landsberg functions which have arbitrary bounded coefficients in the expansion under Schauder basis. The class of the weighted Takagi- Landsberg functions of order $H>0$ on $[0,1]$ coincides with the $H$-Hölder continuous functions on $[0,1]$. Based on computed fractional integrals and derivatives of the Haar and Schauder functions, we get a new series representation of the fractional derivatives of a Hölder continuous function. This result allows to get the new formula of a Riemann-Stieltjes integral. The application of such series representation is the new method of numerical solution of the Volterra and linear integral equations driven by a Hölder continuous function.
We study Volterra processes $X_t = \int_0^t K(t,s) dW_s$, where $W$ is a standard Wiener process, and the kernel has the form $K(t,s) = a(s) \int_s^t b(u) c(u-s) du$. This form generalizes the Volterra kernel for fractional Brownian motion (fBm) with Hurst index $H>1/2$. We establish smoothness properties of $X$, including continuity and Holder property. It happens that its Holder smoothness is close to well-known Holder smoothness of fBm but is a bit worse. We give a comparison with fBm for any smoothness theorem. Then we investigate the problem of inverse representation of $W$ via $X$ in the case where $c\in L^1[0,T]$ creates a Sonine pair, i.e. there exists $h\in L^1[0,T]$ such that $c * h = 1$. It is a natural extension of the respective property of fBm that generates the same filtration with the underlying Wiener process. Since the inverse representation of the Gaussian processes under consideration are based on the properties of Sonine pairs, we provide several examples of Sonine pairs, both well-known and new. Key words: Gaussian process, Volterra process, Sonine pair, continuity, Holder property, inverse representation.
The present paper investigates the effects of tempering the power law kernel of moving average representation of a fractional Brownian motion (fBm) on some local and global properties of this Gaussian stochastic process. Tempered fractional Brownian motion (TFBM) and tempered fractional Brownian motion of the second kind (TFBMII) are the processes that are considered in order to investigate the role of tempering. Tempering does not change the local properties of fBm including the sample paths and p-variation, but it has a strong impact on the Breuer-Major theorem, asymptotic behavior of the 3rd and 4th cumulants of fBm and the optimal fourth moment theorem.
We study the existence and uniqueness of solutions to stochastic differential equations with Volterra processes driven by Lévy noise. For this purpose, we study in detail smoothness properties of these processes. Special attention is given to two kinds of Volterra-Gaussian processes that generalize the compact interval representation of fractional Brownian motion and to stochastic equations with such processes.
We define a time-changed fractional Ornstein-Uhlenbeck process by composing a fractional Ornstein-Uhlenbeck process with the inverse of a subordinator. Properties of the moments of such process are investigated and the existence of the density is shown. We also provide a generalized Fokker-Planck equation for the density of the process.
We consider a stochastic differential equation of the form $dr_t = (a - b r_t) dt + σ\sqrt{r_t}dW_t$, where $a$, $b$ and $σ$ are positive constants. The solution corresponds to the Cox-Ingersoll-Ross process. We study the estimation of an unknown drift parameter $(a,b)$ by continuous observations of a sample path $\{r_t,t\in[0,T]\}$. First, we prove the strong consistency of the maximum likelihood estimator. Since this estimator is well-defined only in the case $2a>σ^2$, we propose another estimator that is defined and strongly consistent for all positive $a$, $b$, $σ$. The quality of the estimators is illustrated by simulation results.
We construct the least-square estimator for the unknown drift parameter in the multifractional Ornstein-Uhlenbeck model and establish its strong consistency in the non-ergodic case. The proofs are based on the asymptotic bounds with probability 1 for the rate of the growth of the trajectories of multifractional Brownian motion (mBm) and of some other functionals of mBm, including increments and fractional derivatives. As the auxiliary results having independent interest, we produce the asymptotic bounds with probability 1 for the rate of the growth of the trajectories of the general Gaussian process and some functionals of it, in terms of the covariance function of its increments.
We introduce a pathwise integration for Volterra processes driven by Lévy noise or martingale noise. These processes are widely used in applications to turbulence, signal processes, biology, and in environmental finance. Indeed they constitute a very flexible class of models, which include fractional Brownian and Lévy motions and it is part of the so-called ambit fields. A pathwise integration with respect of such Volterra processes aims at producing a framework where modelling is easily understandable from an information perspective. The techniques used are based on fractional calculus and in this there is a bridging of the stochastic and deterministic techniques. The present paper aims at setting the basis for a framework in which further computational rules can be devised. Our results are general in the choice of driving noise. Additionally we propose some further details in the relevant context subordinated Wiener processes.
We consider the fractional Ornstein-Uhlenbeck process with an unknown drift parameter and known Hurst parameter $H$. We propose a new method to test the hypothesis of the sign of the parameter and prove the consistency of the test. Contrary to the previous works, our approach is applicable for all $H\in(0,1)$. We also study the estimators for drift parameter for continuous and discrete observations and prove their strong consistency for all $H\in(0,1)$.
We consider a discrete-time approximation of paths of an Ornstein--Uhlenbeck process as a mean for estimation of a price of European call option in the model of financial market with stochastic volatility. The Euler--Maruyama approximation scheme is implemented. We determine the estimates for the option price for predetermined sets of parameters. The rate of convergence of the price and an average volatility when discretization intervals tighten are determined. Discretization precision is analyzed for the case where the exact value of the price can be derived.
We consider the pricing problem related to payoffs that can have discontinuities of polynomial growth. The asset price dynamic is modeled within the Black and Scholes framework characterized by a stochastic volatility term driven by a fractional Ornstein-Uhlenbeck process. In order to solve the aforementioned problem, we consider three approaches. The first one consists in a suitable transformation of the initial value of the asset price, in order to eliminate possible discontinuities. Then we discretize both the Wiener process and the fractional Brownian motion and estimate the rate of convergence of the related discretized price to its real value, the latter one being impossible to be evaluated analytically. The second approach consists in considering the conditional expectation with respect to the entire trajectory of the fractional Brownian motion (fBm). Then we derive a closed formula which involves only integral functional depending on the fBm trajectory, to evaluate the price; finally we discretize the fBm and estimate the rate of convergence of the associated numerical scheme to the option price. In both cases the rate of convergence is the same and equals $n^{-rH}$, where $n$ is a number of the points of discretization, $H$ is the Hurst index of fBm, and $r$ is the Hölder exponent of volatility function. The third method consists in calculating the density of the integral functional depending on the trajectory of the fBm via Malliavin calculus also providing the option price in terms of the associated probability density.
The paper deals with the expected maxima of continuous Gaussian processes $X = (X_t)_{t\ge 0}$ that are Hölder continuous in $L_2$-norm and/or satisfy the opposite inequality for the $L_2$-norms of their increments. Examples of such processes include the fractional Brownian motion and some of its "relatives" (of which several examples are given in the paper). We establish upper and lower bounds for $E \max_{0\le t\le 1}X_t$ and investigate the rate of convergence to that quantity of its discrete approximation $E \max_{0\le i\le n}X_{i/n}$. Some further properties of these two maxima are established in the special case of the fractional Brownian motion.
We consider Langevin equation involving fractional Brownian motion with Hurst index $H\in(0,\frac12)$. Its solution is the fractional Ornstein-Uhlenbeck process and with unknown drift parameter $θ$. We construct the estimator that is similar in form to maximum likelihood estimator for Langevin equation with standard Brownian motion. Observations are discrete in time. It is assumed that the interval between observations is $n^{-1}$, i.e. tends to zero (high frequency data) and the number of observations increases to infinity as $n^m$ with $m>1$. It is proved that for positive $θ$ the estimator is strongly consistent for any $m>1$ and for negative $θ$ it is consistent when $m>\frac{1}{2H}$.
We construct an estimator of the unknown drift parameter $θ\in {\mathbb{R}}$ in the linear model \[X_t=θt+σ_1B^{H_1}(t)+σ_2B^{H_2}(t),\;t\in[0,T],\] where $B^{H_1}$ and $B^{H_2}$ are two independent fractional Brownian motions with Hurst indices $H_1$ and $H_2$ satisfying the condition $\frac{1}{2}\leq H_1<H_2<1.$ Actually, we reduce the problem to the solution of the integral Fredholm equation of the 2nd kind with a specific weakly singular kernel depending on two power exponents. It is proved that the kernel can be presented as the product of a bounded continuous multiplier and weak singular one, and this representation allows us to prove the compactness of the corresponding integral operator. This, in turn, allows us to establish an existence--uniqueness result for the sequence of the equations on the increasing intervals, to construct accordingly a sequence of statistical estimators, and to establish asymptotic consistency.
We construct the maximum likelihood estimator (MLE) of the unknown drift parameter $θ\in \mathbb{R}$ in the linear model $X_t=θt+σB^{H_1}(t)+B^{H_2}(t),\;t\in[0,T],$ where $B^{H_1}$ and $B^{H_2}$ are two independent fractional Brownian motions with Hurst indices $\frac12<H_1<H_2<1.$ The formula for MLE is based on the solution of the integral equation with weak polar kernel.
Explicit solution of an infinite horizon optimal stopping problem for a Levy processes with a polynomial reward function is given, in terms of the overall supremum of the process, when the solution of the problem is one-sided. The results are obtained via the generalization of known results about the averaging function associated with the problem. This averaging function can be directly computed in case of polynomial rewards. To illustrate this result, examples for general quadratic and cubic polynomials are discussed in case the process is Brownian motion, and the optimal stopping problem for a quartic polynomial and a Kou's process is solved.
We deal with a generalization of the classical risk model when an insurance company gets additional funds whenever a claim arrives and consider some practical approaches to the estimation of the ruin probability. In particular, we get an upper exponential bound and construct an analogue to the De Vylder approximation for the ruin probability. We compare results of these approaches with statistical estimates obtained by the Monte Carlo method for selected distributions of claim sizes and additional funds.
We consider the Black--Scholes model of financial market modified to capture the stochastic nature of volatility observed at real financial markets. For volatility driven by the Ornstein--Uhlenbeck process, we establish the existence of equivalent martingale measure in the market model. The option is priced with respect to the minimal martingale measure for the case of uncorrelated processes of volatility and asset price, and an analytic expression for the price of European call option is derived. We use the inverse Fourier transform of a characteristic function and the Gaussian property of the Ornstein--Uhlenbeck process.
We show that small ball estimates together with Holder continuity assumption allow to obtain new representation results in models with long memory. In order to apply these results, we establish small ball probability estimates for Gaussian processes whose incremental variance admits two-sided estimates and the incremental covariance preserves sign. As a result, we obtain small ball estimates for integral transforms of Wiener processes and of fractional Brownian motion with Volterra kernels.
We establish the rate of convergence of distributions of sums of independent identically distributed random variables to the Gaussian distribution in terms of truncated pseudomoments by implementing the idea of Yu. Studnyev for getting estimates of the rate of convergence of the order higher than $n^{-1/2}$.
We consider a utility maximization problem in a broad class of markets. Apart from traditional semimartingale markets, our class of markets includes processes with long memory, fractional Brownian motion and related processes, and, in general, Gaussian processes satisfying certain regularity conditions on their covariance functions. Our choice of markets is motivated by the well-known phenomena of the so-called "constant" and "variable depth" memory observed in real world price processes, for which fractional and multifractional models are the most adequate descriptions. We introduce the notion of a Wiener-transformable Gaussian process, and give examples of such processes, and their representations. The representation for the solution of the utility maximization problem in our specific setting is presented for various utility functions.
Let $W_i=\{W_i(t), t\in \mathbb{R}_+\}, i=1,2$ be two Wiener processes and $W_3=\{W_3(\mathbf{t}), \mathbf{t}\in \mathbb{R}_+^2\}$ be a two-parameter Brownian sheet, all three processes being mutually independent. We derive upper and lower bounds for the boundary non-crossing probability $$P_f=P\{W_1(t_1)+W_2(t_2)+W_3(\mathbf{t})+h(\mathbf{t})\leq u(\mathbf{t}), \mathbf{t}\in\mathbb{R}_+^2\},$$ where $h, u: \mathbb{R}_+^2\rightarrow \mathbb{R}_+$ are two measurable functions. We show further that for large trend functions $γf>0$ asymptotically when $γ\to \infty$ we have that $\ln P_{γf}$ is the same as $\ln P_{γ\underline{f}}$ where $\underline{f}$ is the projection of $f$ on some closed convex set of the reproducing kernel Hilbert Space of $W$. It turns out that our approach is applicable also for the additive Brownian pillow.
The paper is concerned with a mixed stochastic delay differential equation involving both a Wiener process and a $γ$-Hölder continuous process with $γ>1/2$ (e.g. a fractional Brownian motion with Hurst parameter greater than $1/2$). It is shown that its solution depends continuously on the coefficients and the initial data. Two applications of this result are given: the convergence of solutions to equations with vanishing delay to the solution of equation without delay and the convergence of Euler approximations for mixed stochastic differential equations. As a side result of independent interest, the integrability of solution to mixed stochastic delay differential equations is established.
We consider anisotropic self-similar random fields, in particular, the fractional Brownian sheet. This Gaussian field is an extension of fractional Brownian motion. We prove some properties of covariance function for self-similar fields with rectangular increments. Using Lamperti transformation we obtain properties of covariance function for the corresponding stationary fields. We present an example of a Gaussian self-similar field with stationary rectangular increments that is not a fractional Brownian sheet.
We consider a generalization of the classical risk model when the premium intensity depends on the current surplus of an insurance company. All surplus is invested in the risky asset, the price of which follows a geometric Brownian motion. We get an exponential bound for the infinite-horizon ruin probability. To this end, we allow the surplus process to explode and investigate the question concerning the probability of explosion of the surplus process between claim arrivals.
We obtain results on both weak and almost sure asymptotic behaviour of power variations of a linear combination of independent Wiener process and fractional Brownian motion. These results are used to construct strongly consistent parameter estimators in mixed models.
We consider a problem of statistical estimation of an unknown drift parameter for a stochastic differential equation driven by fractional Brownian motion. Two estimators based on discrete observations of solution to the stochastic differential equations are constructed. It is proved that the estimators converge almost surely to the parameter value as the observation interval expands and the distance between observations vanishes. A bound for the rate of convergence is given and numerical simulations are presented. As an auxilliary result of independent interest we establish global estimates for fractional derivative of fractional Brownian motion.
In this paper we investigate the boundary non-crossing probabilities of a fractional Brownian motion considering some general deterministic trend function. We derive bounds for non-crossing probabilities and discuss the case of a large trend function. As a by-product we solve a minimization problem related to the norm of the trend function.
In this article we study a homogeneous transient diffusion process $X$. We combine the theories of differential equations and of stochastic processes to obtain new results for homogeneous diffusion processes, generalizing the results of Salminen and Yor. The distribution of local time of $X$ is found in a closed form. To this end, a second order differential equation corresponding to the generator of $X$ is considered, and properties of its monotone solutions as functions of a parameter are established using their probabilistic representations. We also provide expressions and upper bounds for moments, exponential moments, and potentials of integral functionals of $X$.
We consider the homogeneous stochastic differential equation with unknown parameter to be estimated. We prove that the standard maximum likelihood estimate is strongly consistent under very mild conditions. There are also established the conditions for strong consistency of the discretized estimator.
We study the problem of optimal approximation of a fractional Brownian motion by martingales. We prove that there exist a unique martingale closest to fractional Brownian motion in a specific sense. It shown that this martingale has a specific form. Numerical results concerning the approximation problem are given.
For a mixed stochastic differential equation involving standard Brownian motion and an almost surely Hölder continuous process $Z$ with Hölder exponent $γ>1/2$, we establish a new result on its unique solvability. We also establish an estimate for difference of solutions to such equations with different processes $Z$ and deduce a corresponding limit theorem. As a by-product, we obtain a result on existence of moments of a solution to a mixed equation under an assumption that $Z$ has certain exponential moments.
We show that a pathwise stochastic integral with respect to fractional Brownian motion with an adapted integrand $g$ can have any prescribed distribution, moreover, we give both necessary and sufficient conditions when random variables can be represented in this form. We also prove that any random variable is a value of such integral in some improper sense. We discuss some applications of these results, in particular, to fractional Black--Scholes model of financial market.
Assume that $X$ is a continuous square integrable process with zero mean, defined on some probability space $(Ω,\mathrm {F},\mathrm {P})$. The classical characterization due to P. Lévy says that $X$ is a Brownian motion if and only if $X$ and $X_t^2-t$, $t\ge0,$ are martingales with respect to the intrinsic filtration $\mathrm {F}^X$. We extend this result to fractional Brownian motion.
We consider a stochastic differential equation involving standard and fractional Brownian motion with unknown drift parameter to be estimated. We investigate the standard maximum likelihood estimate of the drift parameter, two non-standard estimates and three estimates for the sequential estimation. Model strong consistency and some other properties are proved. The linear model and Ornstein-Uhlenbeck model are studied in detail. As an auxiliary result, an asymptotic behavior of the fractional derivative of the fractional Brownian motion is established.
We consider a mixed stochastic differential equation driven by possibly dependent fractional Brownian motion and Brownian motion. Under mild regularity assumptions on the coefficients, it is proved that the equation has a unique solution.
Fractional Brownian motion can be represented as an integral of a deterministic kernel w.r.t. an ordinary Brownian motion either on infinite or compact interval. In previous literature fractional Lévy processes are defined by integrating the infinite interval kernel w.r.t. a general Lévy process. In this article we define fractional Lévy processes using the compact interval representation. We prove that the fractional Lévy processes presented via different integral transformations have the same finite dimensional distributions if and only if they are fractional Brownian motions. Also, we present relations between different fractional Lévy processes and analyze the properties of such processes. A financial example is introduced as well.
The paper focuses on discrete-type approximations of solutions to non-homogeneous stochastic differential equations (SDEs) involving fractional Brownian motion (fBm). We prove that the rate of convergence for Euler approximations of solutions of pathwise SDEs driven by fBm with Hurst index $H>1/2$ can be estimated by $O(δ^{2H-1})$ ($δ$ is the diameter of partition). For discrete-time approximations of Skorohod-type quasilinear equation driven by fBm we prove that the rate of convergence is $O(δ^H)$.