Source author record

Jevgenijs Ivanovs

Jevgenijs Ivanovs appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

26 published item(s)

preprint2022arXiv

A series expansion formula of the scale matrix with applications in change-point detection

We introduce a new Levy fluctuation theoretic method to analyze the cumulative sum (CUSUM) procedure in sequential change-point detection. When observations are phase-type distributed and the post-change distribution is given by exponential tilting of its pre-change distribution, the first passage analysis of the CUSUM statistic is reduced to that of a certain Markov additive process. We develop a novel series expansion formula of the scale matrix for Markov additive processes of finite activity, and apply it to derive exact expressions of the average run length, average detection delay, and false alarm probability under the CUSUM procedure.

preprint2021arXiv

Fluctuation theory for one-sided Lévy processes with a matrix-exponential time horizon

There is an abundance of useful fluctuation identities for one-sided Lévy processes observed up to an independent exponentially distributed time horizon. We show that all the fundamental formulas generalize to time horizons having matrix exponential distributions, and the structure is preserved. Essentially, the positive killing rate is replaced by a matrix with eigenvalues in the right half of the complex plane which, in particular, applies to the positive root of the Laplace exponent and the scale function. Various fundamental properties of thus obtained matrices and functions are established, resulting in an easy to use toolkit. An important application concerns deterministic time horizons which can be well approximated by concentrated matrix exponential distributions. Numerical illustrations are also provided.

preprint2021arXiv

On scale functions for Lévy processes with negative phase-type jumps

We provide a novel expression of the scale function for a Lévy processes with negative phase-type jumps. It is in terms of a certain transition rate matrix which is explicit up to a single positive number. A monotone iterative scheme for the calculation of the latter is presented and it is shown that the error decays exponentially fast. Our numerical examples suggest that this algorithm allows to employ phase-type distributions with a hundred of phases, which is problematic when using the known formula for the scale function in terms of roots. Extensions to other distributions, such as matrix-exponential and infinite-dimensional phase-type, can be anticipated.

preprint2020arXiv

Discretization of the Lamperti representation of a positive self-similar Markov process

This paper considers discretization of the Lévy process appearing in the Lamperti representation of a strictly positive self-similar Markov process. Limit theorems for the resulting approximation are established under some regularity assumptions on the given Lévy process. Additionally, the scaling limit of a positive self-similar Markov process at small times is provided. Finally, we present an application to simulation of self-similar Lévy processes conditioned to stay positive.

preprint2020arXiv

Optimal estimation of some random quantities of a Lévy process

In this paper we present new theoretical results on optimal estimation of certain random quantities based on high frequency observations of a Lévy process. More specifically, we investigate the asymptotic theory for the conditional mean and conditional median estimators of the supremum/infimum of a linear Brownian motion and a stable Lévy process. Another contribution of our article is the conditional mean estimation of the local time and the occupation time measure of a linear Brownian motion. We demonstrate that the new estimators are considerably more efficient compared to the classical estimators. Furthermore, we discuss pre-estimation of the parameters of the underlying models, which is required for practical implementation of the proposed statistics.

preprint2020arXiv

Recovering Brownian and jump parts from high-frequency observations of a Lévy process

We introduce two general non-parametric methods for recovering paths of the Brownian and jump components from high-frequency observations of a Lévy process. The first procedure relies on reordering of independently sampled normal increments and thus avoids tuning parameters. The functionality of this method is a consequence of the small time predominance of the Brownian component, the presence of exchangeable structures, and fast convergence of normal empirical quantile functions. The second procedure amounts to filtering the increments and compensating with the final value. It requires a carefully chosen threshold, in which case both methods yield the same rate of convergence. This rate depends on the small-jump activity and is given in terms of the Blumenthal-Getoor index. Finally, we discuss possible extensions, including the multidimensional case, and provide numerical illustrations.

preprint2020arXiv

Sparse Structures for Multivariate Extremes

Extreme value statistics provides accurate estimates for the small occurrence probabilities of rare events. While theory and statistical tools for univariate extremes are well-developed, methods for high-dimensional and complex data sets are still scarce. Appropriate notions of sparsity and connections to other fields such as machine learning, graphical models and high-dimensional statistics have only recently been established. This article reviews the new domain of research concerned with the detection and modeling of sparse patterns in rare events. We first describe the different forms of extremal dependence that can arise between the largest observations of a multivariate random vector. We then discuss the current research topics including clustering, principal component analysis and graphical modeling for extremes. Identification of groups of variables which can be concomitantly extreme is also addressed. The methods are illustrated with an application to flood risk assessment.

preprint2020arXiv

Zooming-in on a Lévy process: Failure to observe threshold exceedance over a dense grid

For a Lévy process $X$ on a finite time interval consider the probability that it exceeds some fixed threshold $x>0$ while staying below $x$ at the points of a regular grid. We establish exact asymptotic behavior of this probability as the number of grid points tends to infinity. We assume that $X$ has a zooming-in limit, which necessarily is $1/α$-self-similar Lévy process with $α\in(0,2]$, and restrict to $α>1$. Moreover, the moments of the difference of the supremum and the maximum over the grid points are analyzed and their asymptotic behavior is derived. It is also shown that the zooming-in assumption implies certain regularity properties of the ladder process, and the decay rate of the left tail of the supremum distribution is determined.

preprint2016arXiv

Exit identities for Lévy processes observed at Poisson arrival times

For a spectrally one-sided Lévy process, we extend various two-sided exit identities to the situation when the process is only observed at arrival epochs of an independent Poisson process. In addition, we consider exit problems of this type for processes reflected either from above or from below. The resulting Laplace transforms of the main quantities of interest are in terms of scale functions and turn out to be simple analogues of the classical formulas.

preprint2016arXiv

Time inhomogeneity in longest gap and longest run problems

Consider an inhomogeneous Poisson process and let $D$ be the first of its epochs which is followed by a gap of size $\ell>0$. We establish a criterion for $D<\infty$ a.s., as well as for $D$ being long-tailed and short-tailed, and obtain logarithmic tail asymptotics in various cases. These results are translated into the discrete time framework of independent non-stationary Bernoulli trials where the analogue of $D$ is the waiting time for the first run of ones of length $\ell$. A main motivation comes from computer reliability, where $D+\ell$ represents the actual execution time of a program or transfer of a file of size $\ell$ in presence of failures (epochs of the process) which necessitate restart.

preprint2015arXiv

A bivariate risk model with mutual deficit coverage

We consider a bivariate Cramer-Lundberg-type risk reserve process with the special feature that each insurance company agrees to cover the deficit of the other. It is assumed that the capital transfers between the companies are instantaneous and incur a certain proportional cost, and that ruin occurs when neither company can cover the deficit of the other. We study the survival probability as a function of initial capitals and express its bivariate transform through two univariate boundary transforms, where one of the initial capitals is fixed at 0. We identify these boundary transforms in the case when claims arriving at each company form two independent processes. The expressions are in terms of Wiener-Hopf factors associated to two auxiliary compound Poisson processes. The case of non-mutual (reinsurance) agreement is also considered.

preprint2015arXiv

Lévy-driven polling systems and continuous-state branching processes

In this paper we consider a ring of $N\ge 1$ queues served by a single server in a cyclic order. After having served a queue (according to a service discipline that may vary from queue to queue), there is a switch-over period and then the server serves the next queue and so forth. This model is known in the literature as a \textit{polling model}. Each of the queues is fed by a non-decreasing Lévy process, which can be different during each of the consecutive periods within the server's cycle. The $N$-dimensional Lévy processes obtained in this fashion are described by their (joint) Laplace exponent, thus allowing for non-independent input streams. For such a system we derive the steady-state distribution of the joint workload at embedded epochs, i.e. polling and switching instants. Using the Kella-Whitt martingale, we also derive the steady-state distribution at an arbitrary epoch. Our analysis heavily relies on establishing a link between fluid (Lévy input) polling systems and multi-type Jiřina processes (continuous-state discrete-time branching processes). This is done by properly defining the notion of the \textit{branching property} for a discipline, which can be traced back to Fuhrmann and Resing. This definition is broad enough to contain the most important service disciplines, like exhaustive and gated.

preprint2015arXiv

Sparre-Andersen identity and the last passage time

It is shown that the celebrated result of Sparre Andersen for random walks and Lévy processes has intriguing consequences when the last time of the process in $(-\infty,0]$, say $σ$, is added to the picture. In the case of no positive jumps this leads to six random times, all of which have the same distribution - the uniform distribution on $[0,σ]$. Surprisingly, this result does not appear in the literature, even though it is based on some classical observations concerning exchangeable increments.

preprint2015arXiv

Splitting and time reversal for Markov additive processes

We consider a Markov additive process with a finite phase space and study its path decompositions at the times of extrema, first passage and last exit. For these three families of times we establish splitting conditional on the phase, and provide various relations between the laws of post- and pre-splitting processes using time reversal. These results offer valuable insight into behavior of the process, and while being structurally similar to the Lévy process case, they demonstrate various new features. As an application we formulate the Wiener-Hopf factorization, where time is counted in each phase separately and killing of the process is phase dependent. Assuming no positive jumps, we find concise formulas for these factors, and also characterize the time of last exit from the negative half-line using three different approaches, which further demonstrates applicability of path decomposition theory.

preprint2015arXiv

Strikingly simple identities relating exit problems for Lévy processes under continuous and Poisson observations

We consider exit problems for general Lévy processes, where the first passage over a threshold is detected either immediately or at an epoch of an independent homogeneous Poisson process. It is shown that the two corresponding one-sided problems are related through a surprisingly simple identity. Moreover, we identify a simple link between two-sided exit problems with one continuous and one Poisson exit. Finally, Poisson exit of a reflected process is connected to the continuous exit of a process reflected at Poisson epochs, and a link between some Parisian type exit problems is established. With the appropriate perspective, the proofs of all these relations turn out to be quite elementary. For spectrally one-sided Lévy processes this approach enables alternative proofs for a number of previously established identities, providing additional insight.

preprint2015arXiv

Transient analysis of a stationary Lévy-driven queue

In this paper we study a queue with Lévy input, without imposing any a priori assumption on the jumps being one-sided. The focus is on computing the transforms of all sorts of quantities related to the transient workload, assuming the workload is in stationarity at time 0. The results are simple expressions that are in terms of the bivariate Laplace exponents of ladder processes. In particular, we derive the transform of the minimum workload attained over an exponentially distributed interval.

preprint2014arXiv

A Lévy process on the real line seen from its supremum and max-stable processes

We consider a process $Z$ on the real line composed from a Lévy process and its exponentially tilted version killed with arbitrary rates and give an expression for the joint law of $Z$ seen from its supremum, the supremum $\overline Z$ and the time $T$ at which the supremum occurs. In fact, it is closely related to the laws of the original and the tilted Lévy processes conditioned to stay negative and positive. The result is used to derive a new representation of stationary particle systems driven by Lévy processes. In particular, this implies that a max-stable process arising from Lévy processes admits a mixed moving maxima representation with spectral functions given by the conditioned Lévy processes.

preprint2014arXiv

Power identities for Lévy risk models under taxation and capital injections

In this paper we study a spectrally negative Lévy process which is refracted at its running maximum and at the same time reflected from below at a certain level. Such a process can for instance be used to model an insurance surplus process subject to tax payments according to a loss-carry-forward scheme together with the flow of minimal capital injections required to keep the surplus process non-negative. We characterize the first passage time over an arbitrary level and the cumulative amount of injected capital up to this time by their joint Laplace transform, and show that it satisfies a simple power relation to the case without refraction. It turns out that this identity can also be extended to a certain type of refraction from below. The net present value of tax collected before the cumulative injected capital exceeds a certain amount is determined, and a numerical illustration is provided.

preprint2013arXiv

A risk model with an observer in a Markov environment

We consider a spectrally-negative Markov additive process as a model of a risk process in random environment. Following recent interest in alternative ruin concepts, we assume that ruin occurs when an independent Poissonian observer sees the process negative, where the observation rate may depend on the state of the environment. Using an approximation argument and spectral theory we establish an explicit formula for the resulting survival probabilities in this general setting. We also discuss an efficient evaluation of the involved quantities and provide a numerical illustration.

preprint2013arXiv

Exact boundaries in sequential testing for phase-type distributions

We consider Wald's sequential probability ratio test for deciding whether a sequence of independent and identically distributed observations comes from a specified phase-type distribution or from an exponentially tilted alternative distribution. In this setting, we derive exact decision boundaries for given Type I and Type II errors by establishing a link with ruin theory. Information on the mean sample size of the test can be retrieved as well. The approach relies on the use of matrix-valued scale functions associated to a certain one-sided Markov additive process. By suitable transformations the results also apply to other types of distributions including some distributions with regularly varying tail.

preprint2013arXiv

Potential measures of one-sided Markov additive processes with reflecting and terminating barriers

Consider a one-sided Markov additive process with an upper and a lower barrier, where each can be either reflecting or terminating. For both defective and non-defective processes and all possible scenarios we identify the corresponding potential measures, which generalizes a number of results for one-sided Levy processes. The resulting rather neat formulas have various applications, and in particular they lead to quasi-stationary distributions of the corresponding processes.

preprint2013arXiv

Two coupled Levy queues with independent input

We consider a pair of coupled queues driven by independent spectrally-positive Levy processes. With respect to the bi-variate workload process this framework includes both the coupled processor model and the two-server fluid network with independent Levy inputs. We identify the joint transform of the stationary workload distribution in terms of Wiener-Hopf factors corresponding to two auxiliary Levy processes with explicit Laplace exponents. We reinterpret and extend the ideas of Cohen and Boxma (1983) to provide a general and uniform result with a neat transform expression.

preprint2011arXiv

Occupation densities in solving exit problems for Markov additive processes and their reflections

This paper solves exit problems for spectrally negative Markov additive processes and their reflections. A so-called scale matrix, which is a generalization of the scale function of a spectrally negative \levy process, plays a central role in the study of exit problems. Existence of the scale matrix was shown in Thm. 3 of Kyprianou and Palmowski (2008). We provide a probabilistic construction of the scale matrix, and identify the transform. In addition, we generalize to the MAP setting the relation between the scale function and the excursion (height) measure. The main technique is based on the occupation density formula and even in the context of fluctuations of spectrally negative Lévy processes this idea seems to be new. Our representation of the scale matrix $W(x)=e^{-Λx}\eL(x)$ in terms of nice probabilistic objects opens up possibilities for further investigation of its properties.

preprint2010arXiv

First passage process of a Markov additive process, with applications to reflection problems

In this paper we consider the first passage process of a spectrally negative Markov additive process (MAP). The law of this process is uniquely characterized by a certain matrix function, which plays a crucial role in fluctuation theory. We show how to identify this matrix using the theory of Jordan chains associated with analytic matrix functions. Importantly, our result also provides us with a technique, which can be used to derive various further identities. We then proceed to show how to compute the stationary distribution associated with a one-sided reflected (at zero) MAP for both the spectrally positive and spectrally negative cases as well as for the two sided reflected Markov-modulated Brownian motion; these results can be interpreted in terms of queues with MAP input.

preprint2010arXiv

Markov-modulated Brownian motion with two reflecting barriers

We consider a Markov-modulated Brownian motion reflected to stay in a strip [0,B]. The stationary distribution of this process is known to have a simple form under some assumptions. We provide a short probabilistic argument leading to this result and explaining its simplicity. Moreover, this argument allows for generalizations including the distribution of the reflected process at an independent exponentially distributed epoch. Our second contribution concerns transient behavior of the reflected system. We identify the joint law of the processes t,X(t),J(t) at inverse local times.