In this paper we build the increasing convex (concave) order for the scalar product of random vectors with an upper (lower) tail permutation decreasing joint density. As applications, we revisit allocations of portfolio risks in financial engineering and of coverage limits and deductibles in insurance. Some related results in the literature are substantially updated.

A multiclass c-server retrial queueing system in which customers arrive according to a class-dependent Markovian arrival process (MAP) is considered. Service and retrial times follow class-dependent phase-type (PH) distributions with the further assumption that PH distributions of retrial times are acyclic. A necessary and sufficient condition for ergodicity is obtained from criteria based on drifts. The infinite state space of the model is truncated with an appropriately chosen Lyapunov function. The truncated model is described as a multidimensional Markov chain, and a Kronecker representation of its generator matrix is numerically analyzed.

We propose a class of models of random walks in a random environment where an exact solution can be given for a stationary distribution. The environment is cast in terms of a Jackson/Gordon–Newell network although alternative interpretations are possible. The main tool is the detailed balance equations. The difference compared to earlier works is that the position of the random walk influences the transition intensities of the network environment and vice versa, creating strong correlations. The form of the stationary distribution is closely related to the well-known product formula.

Consider a d-type (d<∞) Galton–Watson branching process, conditioned on the event that there are at least k≥2 individuals in the nth generation, pick k individuals at random from the nth generation and trace their lines of descent backward in time till they meet. In this paper, the limit behaviors of the distributions of the generation number of the most recent common ancestor of any k chosen individuals and of the whole population are studied for both critical and subcritical cases. Also, we investigate the limit distribution of the joint distribution of the generation number and their types.

In this paper we study a reflected AR(1) process, i.e. a process (Z n )n obeying the recursion Z n +1= max{aZ n +X n ,0}, with (X n )n a sequence of independent and identically distributed (i.i.d.) random variables. We find explicit results for the distribution of Z n (in terms of transforms) in case X n can be written as Y n −B n , with (B n )n being a sequence of independent random variables which are all Exp(λ) distributed, and (Y n )n i.i.d.; when |a|<1 we can also perform the corresponding stationary analysis. Extensions are possible to the case that (B n )n are of phase-type. Under a heavy-traffic scaling, it is shown that the process converges to a reflected Ornstein–Uhlenbeck process; the corresponding steady-state distribution converges to the distribution of a normal random variable conditioned on being positive.

We consider a directed graph on the integers with a directed edge from vertex i to j present with probability n-1, whenever i < j, independently of all other edges. Moreover, to each edge (i, j) we assign weight n-1(j - i). We show that the closure of vertex 0 in such a weighted random graph converges in distribution to the Poisson-weighted infinite tree as n → ∞. In addition, we derive limit theorems for the length of the longest path in the subgraph of the Poisson-weighted infinite tree which has all vertices at weighted distance of at most ρ from the root.

A computable necessary and sufficient condition of separation cutoff is obtained for a sequence of continuous-time upward skip-free chains with the stochastically monotone time-reversals.

We consider a system of N parallel queues with identical exponential service rates and a single dispatcher where tasks arrive as a Poisson process. When a task arrives, the dispatcher always assigns it to an idle server, if there is any, and to a server with the shortest queue among d randomly selected servers otherwise (1≤d≤N). This load balancing scheme subsumes the so-called join-the-idle queue policy (d=1) and the celebrated join-the-shortest queue policy (d=N) as two crucial special cases. We develop a stochastic coupling construction to obtain the diffusion limit of the queue process in the Halfin‒Whitt heavy-traffic regime, and establish that it does not depend on the value of d, implying that assigning tasks to idle servers is sufficient for diffusion level optimality.

Classical Jackson networks are a well-established tool for the analysis of complex systems. In this paper we analyze Jackson networks with the additional features that (i) nodes may have an infinite supply of low priority work and (ii) nodes may be unstable in the sense that the queue length at these nodes grows beyond any bound. We provide the limiting distribution of the queue length distribution at stable nodes, which turns out to be of product form. A key step in establishing this result is the development of a new algorithm based on adjusted traffic equations for detecting unstable nodes. Our results complement the results known in the literature for the subcases of Jackson networks with either infinite supply nodes or unstable nodes by providing an analysis of the significantly more challenging case of networks with both types of nonstandard node present. Building on our product-form results, we provide closed-form solutions for common customer and system oriented performance measures.

We present a very short derivation of the integral representation of the two-sided Skorokhod reflection Z of a continuous function X of bounded variation, which is a generalization of the integral representation of the one-sided map featured in Anantharam and Konstantopoulos (2011) and Konstantopoulos et al. (1996). We also show that Z satisfies a simpler integral representation when additional conditions are imposed on X.

We consider three different schemes for signal routeing on a tree. The vertices of the tree represent transceivers that can transmit and receive signals, and are equipped with independent and identically distributed weights representing the strength of the transceivers. The edges of the tree are also equipped with independent and identically distributed weights, representing the costs for passing the edges. For each one of our schemes, we derive sharp conditions on the distributions of the vertex weights and the edge weights that determine when the root can transmit a signal over arbitrarily large distances.

Motivated by an application in wireless telecommunication networks, we consider a two-type continuum-percolation problem involving a homogeneous Poisson point process of users and a stationary and ergodic point process of base stations. Starting from a randomly chosen point of the Poisson point process, we investigate the distribution of the minimum number of hops that are needed to reach some point of the base station process. In the supercritical regime of continuum percolation, we use the close relationship between Euclidean and chemical distance to identify the distributional limit of the rescaled minimum number of hops that are needed to connect a typical Poisson point to a point of the base station process as its intensity tends to 0. In particular, we obtain an explicit expression for the asymptotic probability that a typical Poisson point connects to a point of the base station process in a given number of hops.

We consider a class of discrete-time Markov chains with state space [0, 1] and the following dynamics. At each time step, first the direction of the next transition is chosen at random with probability depending on the current location. Then the length of the jump is chosen independently as a random proportion of the distance to the respective end point of the unit interval, the distributions of the proportions being fixed for each of the two directions. Chains of that kind were the subjects of a number of studies and are of interest for some applications. Under simple broad conditions, we establish the ergodicity of such Markov chains and then derive closed-form expressions for the stationary densities of the chains when the proportions are beta distributed with the first parameter equal to 1. Examples demonstrating the range of stationary distributions for processes described by this model are given, and an application to a robot coverage algorithm is discussed.

During the course of a day an individual typically mixes with different groups of individuals. Epidemic models incorporating population structure with individuals being able to infect different groups of individuals have received extensive attention in the literature. However, almost exclusively the models assume that individuals are able to simultaneously infect members of all groups, whereas in reality individuals will typically only be able to infect members of any group they currently reside in. In this paper we develop a model where individuals move between a community and their household during the course of the day, only infecting within their current group. By defining a novel branching process approximation with an explicit expression for the probability generating function of the offspring distribution, we are able to derive the probability of a major epidemic outbreak.

Let r and d be positive integers with r<d. Consider a random d-ary tree constructed as follows. Start with a single vertex, and in each time-step choose a uniformly random leaf and give it d newly created offspring. Let 𝒯d,t be the tree produced after t steps. We show that there exists a fixed δ<1 depending on d and r such that almost surely for all large t, every r-ary subtree of 𝒯d,t has less than t δ vertices. The proof involves analysis that also yields a related result. Consider the following iterative construction of a random planar triangulation. Start with a triangle embedded in the plane. In each step, choose a bounded face uniformly at random, add a vertex inside that face and join it to the vertices of the face. In this way, one face is destroyed and three new faces are created. After t steps, we obtain a random triangulated plane graph with t+3 vertices, which is called a random Apollonian network. We prove that there exists a fixed δ<1, such that eventually every path in this graph has length less than t 𝛿, which verifies a conjecture of Cooper and Frieze (2015).

The shapes of branching trees have been linked to disease transmission patterns. In this paper we use the general Crump‒Mode‒Jagers branching process to model an outbreak of an infectious disease under mild assumptions. Introducing a new class of characteristic functions, we are able to derive a formula for the limit of the frequency of the occurrences of a given shape in a general tree. The computational challenges concerning the evaluation of this formula are in part overcome using the jumping chronological contour process. We apply the formula to derive the limit of the frequency of cherries, pitchforks, and double cherries in the constant-rate birth‒death model, and the frequency of cherries under a nonconstant death rate.

The coalescent was introduced by Kingman (1982a), (1982b) and Tajima (1983) as a continuous-time Markov chain model describing the genealogical relationship among sampled genes from a panmictic population of a species. The random mating in a population is a strict condition and the genealogical structure of the population has a strong influence on the genetic variability and the evolution of the species. In this paper, starting from a discrete-time Markov chain model, we show the weak convergence to a continuous-time Markov chain, called the structured coalescent model, describing the genealogy of the sampled genes from whole population by means of passing the limit of the population size. Herbots (1997) proved the weak convergence to the structured coalescent on the condition of conservative migration and Wright–Fisher-type reproduction. We will give the proof on the condition of general migration rates and exchangeable reproduction.

We define and analyze a coalescent process as a recursive box-filling process whose genealogy is given by an ancestral time-reversed, time-inhomogeneous Bienyamé‒Galton‒Watson process. Special interest is on the expected size of a typical box and its probability of being empty. Special cases leading to exact asymptotic computations are investigated when the coalescing mechanisms are either linear fractional or quadratic.

We obtain a complete description of anisotropic scaling limits of the random grain model on the plane with heavy-tailed grain area distribution. The scaling limits have either independent or completely dependent increments along one or both coordinate axes and include stable, Gaussian, and ‘intermediate’ infinitely divisible random fields. The asymptotic form of the covariance function of the random grain model is obtained. Application to superimposed network traffic is included.

This paper provides a comparison of the successive lumping (SL) methodology developed in Katehakis et al. (2015) with the popular lattice path counting (Mohanty (1979)) in obtaining rate matrices for queueing models, satisfying the specific quasi birth and death structure as in Van Leeuwaarden et al. (2009) and Van Leeuwaarden and Winands (2006). The two methodologies are compared both in terms of applicability requirements and numerical complexity by analyzing their performance for the same classical queueing models considered in Van Leeuwaarden et al. (2009). The main findings are threefold. First, when both methods are applicable, the SL-based algorithms outperform the lattice path counting algorithm (LPCA). Second, there are important classes of problems (for example, models with (level) nonhomogenous rates or with finite state spaces) for which the SL methodology is applicable and for which the LPCA cannot be used. Third, another main advantage of SL algorithms over lattice path counting is that the former includes a method to compute the steady state distribution using this rate matrix.