Modern Stochastics: Theory and Applications

A nonlinear stochastic differential equation with the order of nonlinearity higher than one, with several discrete and distributed delays and time varying coefficients is considered. It is shown that the sufficient conditions for exponential mean square stability of the linear part of the considered nonlinear equation also are sufficient conditions for stability in probability of the initial nonlinear equation. Some new sufficient condition of stability in probability for the zero solution of the considered nonlinear non-autonomous stochastic differential equation is obtained which can be considered as a multi-condition of stability because it allows to get for one considered equation at once several different complementary of each other sufficient stability conditions. The obtained results are illustrated with numerical simulations and figures.

In this paper we present some new limit theorems for power variations of stationary increment Lévy driven moving average processes. Recently, such asymptotic results have been investigated in [Ann. Probab. 45(6B) (2017), 4477–4528, Festschrift for Bernt Øksendal, Stochastics 81(1) (2017), 360–383] under the assumption that the kernel function potentially exhibits a singular behaviour at 0. The aim of this work is to demonstrate how some of the results change when the kernel function has multiple singularity points. Our paper is also related to the article [Stoch. Process. Appl. 125(2) (2014), 653–677] that studied the same mathematical question for the class of Brownian semi-stationary models.

We consider the infinite divisibility of distributions of some well-known inverse subordinators. Using a tail probability bound, we establish that distributions of many of the inverse subordinators used in the literature are not infinitely divisible. We further show that the distribution of a renewal process time-changed by an inverse stable subordinator is not infinitely divisible, which in particular implies that the distribution of the fractional Poisson process is not infinitely divisible.

We investigate the pricing of cliquet options in a jump-diffusion model. The considered option is of monthly sum cap style while the underlying stock price model is driven by a drifted Lévy process entailing a Brownian diffusion component as well as compound Poisson jumps. We also derive representations for the density and distribution function of the emerging Lévy process. In this setting, we infer semi-analytic expressions for the cliquet option price by two different approaches. The first one involves the probability distribution function of the driving Lévy process whereas the second draws upon Fourier transform techniques. With view on sensitivity analysis and hedging purposes, we eventually deduce representations for several Greeks while putting emphasis on the Vega.

A fractional advection-dispersion equation (fADE) has been advocated for heavy-tailed flows where the usual Brownian diffusion models fail. A stochastic differential equation (SDE) driven by a stable Lévy process gives a forward equation that matches the space-fractional advection-dispersion equation and thus gives the stochastic framework of particle tracking for heavy-tailed flows. For constant advection and dispersion coefficient functions, the solution to such SDE itself is a stable process and can be derived easily by least square parameter fitting from the observed flow concentration data. However, in a more generalized scenario, a closed form for the solution to a stable SDE may not exist. We propose a numerical method for solving/generating a stable SDE in a general set-up. The method incorporates a discretized finite volume scheme with the characteristic line to solve the fADE or the forward equation for the Markov process that solves the stable SDE. Then we use a numerical scheme to generate the solution to the governing SDE using the fADE solution. Also, often the functional form of the advection or dispersion coefficients are not known for a given plume concentration data to start with. We use a Levenberg–Marquardt (L-M) regularization method to estimate advection and dispersion coefficient function from the observed data (we present the case for a linear advection) and proceed with the SDE solution construction described above.

Confidence ellipsoids for linear regression coefficients are constructed by observations from a mixture with varying concentrations. Two approaches are discussed. The first one is the nonparametric approach based on the weighted least squares technique. The second one is an approximate maximum likelihood estimation with application of the EM-algorithm for the estimates calculation.

This paper deals with a homoskedastic errors-in-variables linear regression model and properties of the total least squares (TLS) estimator. We partly revise the consistency results for the TLS estimator previously obtained by the author [18]. We present complete and comprehensive proofs of consistency theorems. A theoretical foundation for construction of the TLS estimator and its relation to the generalized eigenvalue problem is explained. Particularly, the uniqueness of the estimate is proved. The Frobenius norm in the definition of the estimator can be substituted by the spectral norm, or by any other unitarily invariant norm; then the consistency results are still valid.

The effect that weighted summands have on each other in approximations of $S={w_{1}}{S_{1}}+{w_{2}}{S_{2}}+\cdots +{w_{N}}{S_{N}}$ is investigated. Here, ${S_{i}}$’s are sums of integer-valued random variables, and ${w_{i}}$ denote weights, $i=1,\dots ,N$. Two cases are considered: the general case of independent random variables when their closeness is ensured by the matching of factorial moments and the case when the ${S_{i}}$ has the Markov Binomial distribution. The Kolmogorov metric is used to estimate the accuracy of approximation.

A continuous-time regression model with a jointly strictly sub-Gaussian random noise is considered in the paper. Upper exponential bounds for probabilities of large deviations of the least squares estimator for the regression parameter are obtained.

In the paper we consider time-changed Poisson processes where the time is expressed by compound Poisson-Gamma subordinators $G(N(t))$ and derive the expressions for their hitting times. We also study the time-changed Poisson processes where the role of time is played by the processes of the form $G(N(t)+at)$ and by the iteration of such processes.