Modern Stochastics: Theory and Applications

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 ${\textstyle\int _{0}^{t}}{s^{\gamma }}d{W_{s}}$, where W is a Wiener process, $\gamma >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){1_{\{{Y_{\varepsilon }}(t)>0\}}}+\varepsilon }-a{Y_{\varepsilon }}(t))dt+\sigma d{B^{H}}(t)$, as $\varepsilon \downarrow 0$. 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 define the fractional Cox–Ingersoll–Ross process as $X_{t}:={Y_{t}^{2}}\mathbf{1}_{\{t<\inf \{s>0:Y_{s}=0\}\}}$, where the process $Y=\{Y_{t},t\ge 0\}$ satisfies the SDE of the form $dY_{t}=\frac{1}{2}(\frac{k}{Y_{t}}-aY_{t})dt+\frac{\sigma }{2}d{B_{t}^{H}}$, $\{{B_{t}^{H}},t\ge 0\}$ is a fractional Brownian motion with an arbitrary Hurst parameter $H\in (0,1)$. We prove that $X_{t}$ satisfies the stochastic differential equation of the form $dX_{t}=(k-aX_{t})dt+\sigma \sqrt{X_{t}}\circ d{B_{t}^{H}}$, where the integral with respect to fractional Brownian motion is considered as the pathwise Stratonovich integral. We also show that for $k>0$, $H>1/2$ the process is strictly positive and never hits zero, so that actually $X_{t}={Y_{t}^{2}}$. Finally, we prove that in the case of $H<1/2$ the probability of not hitting zero on any fixed finite interval by the fractional Cox–Ingersoll–Ross process tends to 1 as $k\to \infty $.

We consider a stochastic differential equation of the form with multiplicative stochastic volatility, where Y is some adapted stochastic process. We prove existence–uniqueness results for weak and strong solutions of this equation under various conditions on the process Y and the coefficients a, $\sigma _{1}$, and $\sigma _{2}$. Also, we study the strong consistency of the maximum likelihood estimator for the unknown parameter θ. We suppose that Y is in turn a solution of some diffusion SDE. Several examples of the main equation and of the process Y are provided supplying the strong consistency.

This short note is devoted to the “Fractality and Fractionality” workshop, held on 17–20 May 2016 in Lorentz Center (Leiden, Netherlands).

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 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 construct an estimator of the unknown drift parameter $\theta \in \mathbb{R}$ in the linear model 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}\le 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 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 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.