Source author record

Michel Mandjes

Michel Mandjes 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

32works
8topics
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

32 published item(s)

preprint2023arXiv

Bankruptcy probabilities under non-Poisson inspection

This paper concerns an insurance firm's surplus process observed at renewal inspection times, with a focus on assessing the probability of the surplus level dropping below zero. For various types of inter-inspection time distributions, an explicit expression for the corresponding transform is given. In addition, Cramér-Lundberg type asymptotics are established. Also, an importance sampling based Monte Carlo algorithm is proposed, and is shown to be logarithmically efficient.

preprint2022arXiv

A relative approach to opinion formation

Formal models of opinion formation commonly represent an individual's opinion by a value on a fixed opinion interval. We propose an alternative modeling method wherein interpretation is only provided to the relative positions of opinions vis-à-vis each other. This method is then considered in a similar setting as the discrete-time Altafini model (an extension of the well-known DeGroot model), but with more general influence weights. Even in a linear framework, the model can describe, in the long run, polarization, dynamics with a periodic pattern, and (modulus) consensus formation. In addition, in our alternative approach key characteristics of the opinion dynamic can be derived from real-valued square matrices of influence weights, which immediately allows one to transfer matrix theory insights to the field of opinion formation dynamics under more relaxed conditions than in the DeGroot or discrete-time Altafini models. A few specific themes are covered: (i) We demonstrate how stable patterns in relative opinion dynamics are identified which are hidden when opinions are considered in an absolute opinion framework. (ii) For the two-agent case, we provide an exhaustive closed-form description of the relative opinion model's dynamic in the long run. (iii) We explore group dynamics analytically, in particular providing a non-trivial condition under which a subgroup's asymptotic behavior carries over to the entire population.

preprint2022arXiv

Estimating customer impatience in a service system with unobserved balking

This paper studies a service system in which arriving customers are provided with information about the delay they will experience. Based on this information they decide to wait for service or to leave the system. Specifically, every customer has a patience threshold and they balk if the observed delay is above the threshold. The main objective is to estimate the parameters of the customers' patience-level distribution and the corresponding potential arrival rate, using knowledge of the actual queue-length process only. The main complication, and distinguishing feature of our setup, lies in the fact that customers who decide not to join are not observed, remarkably, we manage to devise a procedure to estimate the underlying patience and arrival rate parameters. The model is a multi-server queue with a Poisson stream of customers, enabling evaluation of the corresponding likelihood function of the state-dependent effective arrival process. We establish strong consistency of the MLE and derive the asymptotic distribution of the estimation error. Several applications and extensions of the method are discussed. The performance is further assessed through a series of numerical experiments. By fitting parameters of hyperexponential and generalized-hyperexponential distributions our method provides a robust estimation framework for any continuous patience-level distribution.

preprint2022arXiv

Optimal departure time advice in road networks with stochastic disruptions

Due to recurrent (e.g. daily or weekly) patterns and non-recurrent disruptions (e.g. caused by incidents), travel times in road networks are time-dependent and inherently random. This is challenging for travelers planning a future trip, aiming to ensure on-time arrival at the destination, while also trying to limit the total travel-time budget spent. The focus of this paper lies on determining their optimal departure time: the latest time of departure for which a chosen on-time arrival probability can be guaranteed. To model the uncertainties in the network, a Markovian background process is used, tracking events affecting the driveable vehicle speeds on the links, thus enabling us to incorporate both recurrent and non-recurrent effects. It allows the evaluation of the travel-time distribution, given the state of this process at departure, on each single link. Then, a computationally efficient algorithm is devised that uses these individual link travel-time distributions to obtain the optimal departure time for a given path or origin-destination pair. Since the conditions in the road network, and thus the state of the background process, may change between the time of request and the advised time of departure, we consider an online version of this procedure as well, in which the traveler receives departure time updates while still at the origin. Finally, numerical experiments are conducted to exemplify a selection of properties of the optimal departure time and, moreover, quantify the performance of the presented algorithms in an existing road network -- the Dutch highway network.

preprint2021arXiv

A diffusion-based analysis of a multi-class road traffic network

This paper studies a stochastic model that describes the evolution of vehicle densities in a road network. It is consistent with the class of (deterministic) kinematic wave models, which describe traffic flows on the basis of conservation laws that incorporate the macroscopic fundamental diagram (a functional relationship between vehicle density and flow). Our setup is capable of handling multiple types of vehicle densities, with general macroscopic fundamental diagrams, on a network with arbitrary topology. Interpreting our system as a spatial population process, we derive, under a natural scaling, fluid and diffusion limits. More specifically, the vehicle density process can be approximated with a suitable Gaussian process, which yield accurate normal approximations to the joint (in the spatial and temporal sense) vehicle density process. The corresponding means and variances can be computed efficiently. Along the same lines, we develop an approximation to the vehicles' travel-time distribution between any given origin and destination pair. Finally, we present a series of numerical experiments that demonstrate the accuracy of the approximations and illustrate the usefulness of the results.

preprint2020arXiv

A Gaussian segment-based traffic flow model for the design and control of transport networks

In the setting of a recently developed cellular stochastic traffic flow model, it has shown that the joint per-cell vehicle densities, as a function of time, can be accurately approximated by a Gaussian process, which has the attractive feature that its means and (spatial and temporal) covariances can be efficiently evaluated. The present article demonstrates the rich potential of this methodology in the context of road traffic control and transportation network design. To solidly provide empirical backing for the use of a multivariate Gaussian approximation, we rely on a detailed historical dataset that contains traffic flow data. Then, in the remainder of the paper, we provide a sequence of design and control related example questions that can be analyzed using the Gaussian methodology. These cover the following topics: (i) evaluation of stationary performance measures, (ii) route selection, (iii) control of traffic flows, and (iv) performance of traffic networks with arbitrary topology. In discussing the setup, results, and applications of these examples, we stress the appropriateness of our {stochastic} traffic model over a deterministic counterpart.

preprint2020arXiv

A Multiplicative Version of the Lindley Recursion

This paper presents an analysis of the stochastic recursion $W_{i+1} = [V_iW_i+Y_i]^+$ that can be interpreted as an autoregressive process of order 1, reflected at 0. We start our exposition by a discussion of the model's stability condition. Writing $Y_i=B_i-A_i$, for independent sequences of non-negative i.i.d.\ random variables $\{A_i\}_{i\in N_0}$ and $\{B_i\}_{i\in N_0}$, and assuming $\{V_i\}_{i\in N_0}$ is an i.i.d. sequence as well (independent of $\{A_i\}_{i\in N_0}$ and $\{B_i\}_{i\in N_0}$), we then consider three special cases: (i) $V_i$ attains negative values only and $B_i$ has a rational LST, (ii) $V_i$ equals a positive value $a$ with certain probability $p\in (0,1)$ and is negative otherwise, and both $A_i$ and $B_i$ have a rational LST, (iii) $V_i$ is uniformly distributed on $[0,1]$, and $A_i$ is exponentially distributed. In all three cases we derive transient and stationary results, where the transient results are in terms of the transform at a geometrically distributed epoch.

preprint2020arXiv

A roundabout model with on-ramp queues: exact results and scaling approximations

This paper introduces a general model of a single-lane roundabout, represented as a circular lattice that consists of $L$ cells, with Markovian traffic dynamics. Vehicles enter the roundabout via on-ramp queues that have stochastic arrival processes, remain on the roundabout a random number of cells, and depart via off-ramps. Importantly, the model does not oversimplify the dynamics of traffic on roundabouts, while various performance-related quantities (such as delay and queue length) allow an analytical characterization. In particular, we present an explicit expression for the marginal stationary distribution of each cell on the lattice. Moreover, we derive results that give insight on the dependencies between parts of the roundabout, and on the queue distribution. Finally, we find scaling limits that allow, for every partition of the roundabout in segments, to approximate 1) the joint distribution of the occupation of these segments by a multivariate Gaussian distribution; and 2) the joint distribution of their total queue lengths by a collection of independent Poisson random variables. To verify the scaling limit statements, we develop a novel way to empirically assess convergence in distribution of random variables.

preprint2020arXiv

Large deviations for acyclic networks of queues with correlated Gaussian inputs

We consider an acyclic network of single-server queues with heterogeneous processing rates. It is assumed that each queue is fed by the superposition of a large number of i.i.d. Gaussian processes with stationary increments and positive drifts, which can be correlated across different queues. The flow of work departing from each server is split deterministically and routed to its neighbors according to a fixed routing matrix, with a fraction of it leaving the network altogether. We study the exponential decay rate of the probability that the steady-state queue length at any given node in the network is above any fixed threshold, also referred to as the "overflow probability". In particular, we first leverage Schilder's sample-path large deviations theorem to obtain a general lower bound for the limit of this exponential decay rate, as the number of Gaussian processes goes to infinity. Then, we show that this lower bound is tight under additional technical conditions. Finally, we show that if the input processes to the different queues are non-negatively correlated, non short-range dependent fractional Brownian motions, and if the processing rates are large enough, then the asymptotic exponential decay rates of the queues coincide with the ones of isolated queues with appropriate Gaussian inputs.

preprint2020arXiv

Performance of the smallest-variance-first rule in appointment sequencing

A classical problem in appointment scheduling, with applications in health care, concerns the determination of the patients' arrival times that minimize a cost function that is a weighted sum of mean waiting times and mean idle times. One aspect of this problem is the sequencing problem, which focuses on ordering the patients. We assess the performance of the smallest-variance-first (SVF) rule, which sequences patients in order of increasing variance of their service durations. While it was known that SVF is not always optimal, it has been widely observed that it performs well in practice and simulation. We provide a theoretical justification for this observation by proving, in various settings, quantitative worst-case bounds on the ratio between the cost incurred by the SVF rule and the minimum attainable cost. We also show that, in great generality, SVF is asymptotically optimal, i.e., the ratio approaches 1 as the number of patients grows large. While evaluating policies by considering an approximation ratio is a standard approach in many algorithmic settings, our results appear to be the first of this type in the appointment scheduling literature.

preprint2020arXiv

Single-server queues under overdispersion in the heavy-traffic regime

This paper addresses the analysis of the queue-length process of single-server queues under overdispersion, i.e., queues fed by an arrival process for which the variance of the number of arrivals in a given time window exceeds the corresponding mean. Several variants are considered, using concepts as mixing and Markov modulation, resulting in different models with either endogenously triggered or exogenously triggered random environments. Only in special cases explicit expressions can be obtained, e.g. when the random arrival and/or service rate can attain just finitely many values. While for more general model variants exact analysis is challenging, one ${\it can}$ derive limit theorems in the heavy-traffic regime. In some of our derivations we rely on evaluating the relevant Laplace transform in the heavy-traffic scaling using Taylor expansions, whereas other results are obtained by applying the continuous mapping theorem.

preprint2019arXiv

Queues on a dynamically evolving graph

This paper considers a population process on a dynamically evolving graph, which can be alternatively interpreted as a queueing network. The queues are of infinite-server type, entailing that at each node all customers present are served in parallel. The links that connect the queues have the special feature that they are unreliable, in the sense that their status alternates between 'up' and 'down'. If a link between two nodes is down, with a fixed probability each of the clients attempting to use that link is lost; otherwise the client remains at the origin node and reattempts using the link (and jumps to the destination node when it finds the link restored). For these networks we present the following results: (a) a system of coupled partial differential equations that describes the joint probability generating function corresponding to the queues' time-dependent behavior (and a system of ordinary differential equations for its stationary counterpart), (b) an algorithm to evaluate the (time-dependent and stationary) moments, and procedures to compute user-perceived performance measures which facilitate the quantification of the impact of the links' outages, (c) a diffusion limit for the joint queue length process.We include explicit results for a series relevant special cases, such as tandem networks and symmetric fully connected networks.

preprint2016arXiv

Convergence rates of Laplace-transform based estimators

This paper considers the problem of estimating probabilities of the form $\mathbb{P}(Y \leq w)$, for a given value of $w$, in the situation that a sample of i.i.d.\ observations $X_1, \ldots, X_n$ of $X$ is available, and where we explicitly know a functional relation between the Laplace transforms of the non-negative random variables $X$ and $Y$. A plug-in estimator is constructed by calculating the Laplace transform of the empirical distribution of the sample $X_1, \ldots, X_n$, applying the functional relation to it, and then (if possible) inverting the resulting Laplace transform and evaluating it in $w$. We show, under mild regularity conditions, that the resulting estimator is weakly consistent and has expected absolute estimation error $O(n^{-1/2} \log(n+1))$. We illustrate our results by two examples: in the first we estimate the distribution of the workload in an M/G/1 queue from observations of the input in fixed time intervals, and in the second we identify the distribution of the increments when observing a compound Poisson process at equidistant points in time (usually referred to as `decompounding').

preprint2016arXiv

Detecting Markov Chain Instability: A Monte Carlo Approach

We devise a Monte Carlo based method for detecting whether a non-negative Markov chain is stable for a given set of parameter values. More precisely, for a given subset of the parameter space, we develop an algorithm that is capable of deciding whether the set has a subset of positive Lebesgue measure for which the Markov chain is unstable. The approach is based on a variant of simulated annealing, and consequently only mild assumptions are needed to obtain performance guarantees. The theoretical underpinnings of our algorithm are based on a result stating that the stability of a set of parameters can be phrased in terms of the stability of a single Markov chain that searches the set for unstable parameters. Our framework leads to a procedure that is capable of performing statistically rigorous tests for instability, which has been extensively tested using several examples of standard and non-standard queueing networks.

preprint2016arXiv

Functional central limit theorems for Markov-modulated infinite-server systems

In this paper we study the Markov-modulated M/M/$\infty$ queue, with a focus on the correlation structure of the number of jobs in the system. The main results describe the system's asymptotic behavior under a particular scaling of the model parameters in terms of a functional central limit theorem. More specifically, relying on the martingale central limit theorem, this result is established, covering the situation in which the arrival rates are sped up by a factor $N$ and the transition rates of the background process by $N^α$, for some $α>0$. The results reveal an interesting dichotomy, with crucially different behavior for $α>1$ and $α<1$, respectively. The limiting Gaussian process, which is of the Ornstein-Uhlenbeck type, is explicitly identified, and it is shown to be in accordance with explicit results on the mean, variances and covariances of the number of jobs in the system.

preprint2016arXiv

Lévy-driven GPS queues with heavy-tailed input

In this paper we derive exact large-buffer asymptotics for a two-class Generalized Processor Sharing (GPS) model, under the assumption that the input traffic streams generated by both classes correspond to heavy-tailed Lévy processes. Four scenarios need to be distinguished, which differ in terms of (i)~the level of heavy-tailedness of the driving Lévy processes as well as (ii)~the values of the corresponding mean rates relative to the GPS weights. The derived results are illustrated by two important special cases, in which the queues' inputs are modeled by heavy-tailed compound Poisson processes and by $α$-stable Lévy motions.

preprint2016arXiv

Refined large deviations asymptotics for Markov-modulated infinite-server systems

Many networking-related settings can be modeled by Markov-modulated infinite-server systems. In such models, the customers' arrival rates and service rates are modulated by a Markovian background process, additionally, there are infinitely many servers (and consequently the resulting model is often used as a proxy for the corresponding many-server model). The Markov-modulated infinite-server model hardly allows any explicit analysis, apart from results in terms of systems of (ordinary or partial) differential equations for the underlying probability generating functions, and recursions to obtain all moments. As a consequence, recent research efforts have pursued an asymptotic analysis in various limiting regimes, notably the central-limit regime (describing fluctuations around the average behavior) and the large-deviations regime (focusing on rare events). Many of these results use the property that the number of customers in the system obeys a Poisson distribution with a random parameter. The objective of this paper is to develop techniques to accurately approximate tail probabilities in the large-deviations regime. We consider the scaling in which the arrival rates are inflated by a factor $N$, and we are interested in the probability that the number of customers exceeds a given level $Na$. Where earlier contributions focused on so-called $logarithmic$ $asymptotics$ of this exceedance probability (which are inherently imprecise), the present paper improves upon those results in that $exact$ $asymptotics$ are established. These are found in two steps: first the distribution of the random parameter of the Poisson distribution is characterized, and then this knowledge is used to identify the exact asymptotics. The paper is concluded by a set of numerical experiments, in which the accuracy of the asymptotic results is assessed.

preprint2016arXiv

Scaling limits for infinite-server systems in a random environment

This paper studies the effect of an overdispersed arrival process on the performance of an infinite-server system. In our setup, a random environment is modeled by drawing an arrival rate $Λ$ from a given distribution every $Δ$ time units, yielding an i.i.d. sequence of arrival rates $Λ_1,Λ_2, \ldots$. Applying a martingale central limit theorem, we obtain a functional central limit theorem for the scaled queue length process. We proceed to large deviations and derive the logarithmic asymptotics of the queue length's tail probabilities. As it turns out, in a rapidly changing environment (i.e., $Δ$ is small relative to $Λ$) the overdispersion of the arrival process hardly affects system behavior, whereas in a slowly changing random environment it is fundamentally different; this general finding applies to both the central limit and the large deviations regime. We extend our results to the setting where each arrival creates a job in multiple infinite-server queues.

preprint2015arXiv

Analysis of Markov-modulated infinite-server queues in the central-limit regime

This paper focuses on an infinite-server queue modulated by an independently evolving finite-state Markovian background process, with transition rate matrix $Q\equiv(q_{ij})_{i,j=1}^d$. Both arrival rates and service rates are depending on the state of the background process. The main contribution concerns the derivation of central limit theorems for the number of customers in the system at time $t\ge 0$, in the asymptotic regime in which the arrival rates $λ_i$ are scaled by a factor $N$, and the transition rates $q_{ij}$ by a factor $N^α$, with $α\in \mathbb R^+$. The specific value of $α$ has a crucial impact on the result: (i) for $α>1$ the system essentially behaves as an M/M/$\infty$ queue, and in the central limit theorem the centered process has to be normalized by $\sqrt{N}$; (ii) for $α<1$, the centered process has to be normalized by $N^{{1-}α/2}$, with the deviation matrix appearing in the expression for the variance.

preprint2015arXiv

Extremes of multidimensional Gaussian processes

This paper considers extreme values attained by a centered, multidimensional Gaussian process $X(t)= (X_1(t),\ldots,X_n(t))$ minus drift $d(t)=(d_1(t),\ldots,d_n(t))$, on an arbitrary set $T$. Under mild regularity conditions, we establish the asymptotics of \[\log\mathbb P\left(\exists{t\in T}:\bigcap_{i=1}^n\left\{X_i(t)-d_i(t)>q_iu\right\}\right),\] for positive thresholds $q_i>0$, $i=1,\ldots,n$, and $u\to\infty$. Our findings generalize and extend previously known results for the single-dimensional and two-dimensional cases. A number of examples illustrate the theory.

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

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.

preprint2015arXiv

Transient analysis of one-sided Lévy-driven queues

In this paper we analyze the transient behavior of the workload process in a Lévy input queue. We are interested in the value of the workload process at a random epoch; this epoch is distributed as the sum of independent exponential random variables. We consider both cases of spectrally one-sided Lévy input processes, for which we succeed in deriving explicit results. As an application we approximate the mean and the Laplace transform of the workload process after a deterministic time.

preprint2014arXiv

Logarithmic asymptotics for multidimensional extremes under non-linear scalings

Let $\boldsymbol W=\{\boldsymbol W_n:n\in\mathbb N\}$ be a sequence of random vectors in $\mathbb R^d$, $d\ge 1$. This paper considers the logarithmic asymptotics of the extremes of $\boldsymbol W$, that is, for any vector $\boldsymbol q>\boldsymbol 0$ in $\mathbb R^d$, we find $$\log\mathbb P\left(\exists{n\in\mathbb N}:\boldsymbol W_n> u \boldsymbol q\right), \quad\text{as} u\to\infty.$$ We follow the approach of the restricted large deviation principle introduced in Duffy et al. \textit{Logarithmic asymptotics for the supremum of a stochastic process} (Ann. Appl. Probab., 13:430--445, 2003). That is, we assume that, for every $\boldsymbol q\ge\boldsymbol 0$, and some scalings $\{a_n\},\{v_n\}$, $\frac{1}{v_n}\log\mathbb P\left(\boldsymbol W_n/a_n\ge u \boldsymbol q\right)$ has a, continuous in $\boldsymbol q$, limit $J_{\boldsymbol W}(\boldsymbol q)$. We allow the scalings $\{a_n\}$ and $\{v_n\}$ to be regularly varying with a positive index. This approach is general enough to incorporate sequences $\boldsymbol W$, such that the probability law of $\boldsymbol W_n/a_n$ satisfies the large deviation principle with continuous, not necessarily convex, rate functions. The formula for these asymptotics agrees with the seminal papers on this topic.

preprint2014arXiv

On the tail asymptotics of the area swept under the Brownian storage graph

In this paper, the area swept under the workload graph is analyzed: with $\{Q(t) : t\ge0\}$ denoting the stationary workload process, the asymptotic behavior of \[π_{T(u)}(u):={\mathbb{P}}\biggl(\int_0^ {T(u)}Q(r)\,\mathrm{d}r>u\biggr)\] is analyzed. Focusing on regulated Brownian motion, first the exact asymptotics of $π_{T(u)}(u)$ are given for the case that $T(u)$ grows slower than $\sqrt{u}$, and then logarithmic asymptotics for (i) $T(u)=T\sqrt{u}$ (relying on sample-path large deviations), and (ii) $\sqrt{u}=\mathrm{o}(T(u))$ but $T(u)=\mathrm{o}(u)$. Finally, the Laplace transform of the residual busy period are given in terms of the Airy function.

preprint2013arXiv

Limit theorems for reflected Ornstein-Uhlenbeck processes

This paper studies one-dimensional Ornstein-Uhlenbeck processes, with the distinguishing feature that they are reflected on a single boundary (put at level 0) or two boundaries (put at levels 0 and d>0). In the literature they are referred to as reflected OU (ROU) and doubly-reflected OU (DROU) respectively. For both cases, we explicitly determine the decay rates of the (transient) probability to reach a given extreme level. The methodology relies on sample-path large deviations, so that we also identify the associated most likely paths. For DROU, we also consider the `idleness process' $L_t$ and the `loss process' $U_t$, which are the minimal nondecreasing processes which make the OU process remain $\geqslant 0$ and $\leqslant d$, respectively. We derive central limit theorems for $U_t$ and $L_t$, using techniques from stochastic integration and the martingale central limit theorem.

preprint2011arXiv

Gaussian queues in light and heavy traffic

In this paper we investigate Gaussian queues in the light-traffic and in the heavy-traffic regime. The setting considered is that of a centered Gaussian process $X\equiv\{X(t):t\in\mathbb R\}$ with stationary increments and variance function $σ^2_X(\cdot)$, equipped with a deterministic drift $c>0$, reflected at 0: \[Q_X^{(c)}(t)=\sup_{-\infty<s\le t}(X(t)-X(s)-c(t-s)).\] We study the resulting stationary workload process $Q^{(c)}_X\equiv\{Q_X^{(c)}(t):t\ge0\}$ in the limiting regimes $c\to 0$ (heavy traffic) and $c\to\infty$ (light traffic). The primary contribution is that we show for both limiting regimes that, under mild regularity conditions on the variance function, there exists a normalizing function $δ(c)$ such that $Q^{(c)}_X(δ(c)\cdot)/σ_X(δ(c))$ converges to a non-trivial limit in $C[0,\infty)$.

preprint2011arXiv

On the infimum attained by a reflected Lévy process

This paper considers a Lévy-driven queue (i.e., a Lévy process reflected at 0), and focuses on the distribution of $M(t)$, that is, the minimal value attained in an interval of length $t$ (where it is assumed that the queue is in stationarity at the beginning of the interval). The first contribution is an explicit characterization of this distribution, in terms of Laplace transforms, for spectrally one-sided Lévy processes (i.e., either only positive jumps or only negative jumps). The second contribution concerns the asymptotics of $\prob{M(T_u)> u}$ (for different classes of functions $T_u$ and $u$ large); here we have to distinguish between heavy-tailed and light-tailed scenarios.

preprint2011arXiv

Quasi-stationary workload in a Lévy-driven storage system

In this paper we analyze the quasi-stationary workload of a Lévy-driven storage system. More precisely, assuming the system is in stationarity, we study its behavior conditional on the event that the busy period $T$ in which time 0 is contained has not ended before time $t$, as $t\to\infty$. We do so by first identifying the double Laplace transform associated with the workloads at time 0 and time $t$, on the event $\{T>t\}.$ This transform can be explicitly computed for the case of spectrally one-sided jumps. Then asymptotic techniques for Laplace inversion are relied upon to find the corresponding behavior in the limiting regime that $t\to\infty.$ Several examples are treated; for instance in the case of Brownian input, we conclude that the workload distribution at time 0 and $t$ are both Erlang(2).

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.

preprint2009arXiv

Sample-path Large Deviations in Credit Risk

The event of large losses plays an important role in credit risk. As these large losses are typically rare, and portfolios usually consist of a large number of positions, large deviation theory is the natural tool to analyze the tail asymptotics of the probabilities involved. We first derive a sample-path large deviation principle (LDP) for the portfolio's loss process, which enables the computation of the logarithmic decay rate of the probabilities of interest. In addition, we derive exact asymptotic results for a number of specific rare-event probabilities, such as the probability of the loss process exceeding some given function.

preprint2004arXiv

Exact asymptotics for fluid queues fed by multiple heavy-tailed on-off flows

We consider a fluid queue fed by multiple On-Off flows with heavy-tailed (regularly varying) On periods. Under fairly mild assumptions, we prove that the workload distribution is asymptotically equivalent to that in a reduced system. The reduced system consists of a ``dominant'' subset of the flows, with the original service rate subtracted by the mean rate of the other flows. We describe how a dominant set may be determined from a simple knapsack formulation. The dominant set consists of a ``minimally critical'' set of On-Off flows with regularly varying On periods. In case the dominant set contains just a single On-Off flow, the exact asymptotics for the reduced system follow from known results. For the case of several On-Off flows, we exploit a powerful intuitive argument to obtain the exact asymptotics. Combined with the reduced-load equivalence, the results for the reduced system provide a characterization of the tail of the workload distribution for a wide range of traffic scenarios.