Well-known inequalities for the spectral gap of a discrete-time Markov chain, such as Poincaré's and Cheeger's inequalities, do not perform well if the transition graph of the Markov chain is strongly connected. For example in the case of nearest-neighbour random walk on the n-dimensional cube Poincaré's and Cheeger's inequalities are off by a factor n. Using a coupling technique and a contraction principle lower bounds on the spectral gap can be derived. Finally, we show that using the contraction principle yields a sharp estimate for nearest-neighbour random walk on the n-dimensional cube.

We extend large exceedence results for i.i.d. -valued random variables to a class of uniformly recurrent Markov-additive processes and stationary strong-mixing processes. As in the i.i.d. case, the results are proved via large deviations estimates.

An optimal coupling is a bivariate distribution with specified marginals achieving maximal correlation. We show that optimal couplings are totally positive and, in fact, satisfy a strictly stronger condition we call the nonintersection property. For discrete distributions we illustrate the equivalence between optimal coupling and a certain transportation problem. Specifically, the optimal solutions of greedily-solvable transportation problems are totally positive, and even nonintersecting, through a rearrangement of matrix entries that results in a Monge sequence. In coupling continuous random variables or random vectors, we exploit a characterization of optimal couplings in terms of subgradients of a closed convex function to establish a generalization of the nonintersection property. We argue that nonintersection is not only stronger than total positivity, it is the more natural concept for the singular distributions that arise in coupling continuous random variables.

We consider the problem of stochastic comparison of general GARCH-like processes for different parameters and different distributions of the innovations. We identify several stochastic orders that are propagated from the innovations to the GARCH process itself, and we discuss their interpretations. We focus on the convex order and show that in the case of symmetric innovations it is also propagated to the cumulated sums of the GARCH process. More generally, we discuss multivariate comparison results related to the multivariate convex and supermodular orders. Finally, we discuss ordering with respect to the parameters in the GARCH(1, 1) case.

Consider a number of events in a probability space. Let X be a random variable that is the number of events that occur. Given some of the moments of the distribution of X, it is possible to obtain bounds on the probability that at least one event occurs. The best possible bounds are given here. If there are many equiprobable events that are d- wise independent, and d is even, then the probability that at least one event happens is at least 1 — O(µ–d/2), where μ = E(X).

It is shown that methods used by the author to approximate limit probabilities for Markov processes from their Q-matrices extend to semi-Markov processes. The limit probabilities for semi-Markov processes can be approximated using only truncations of the embedded Markov chain transition matrix and the vector of mean holding times.

Given a stationary Markov process with s transient states and r absorbing states, a matrix infinite series solution is presented for the variance of duration of stay in state j within the interval [0, t), given initial state i. Closed forms are derived for absorbing states, and for transient states if eigenvalues are real and distinct. Several relationships among Markov matrices are presented.

In a superfair red and black gambling house where the player must bet at least 2 units at each stage, a gambler wishes to maximize the probability of reaching a goal integer N before reaching zero. For win probability p > 1/2, when N is even, an optimal strategy is to bet 3 units when you have 3 units or N − 3 units and to bet 2 units otherwise. When N is odd, there are two strategies which are optimal depending on the value of the win probability p. When p is smaller than a certain value, p*, the above strategy is optimal, and when p is larger than this value, the timid strategy is optimal.

A summation formula is given for computing the spectrum of intervals of a geometric branching Poisson process from the cumulant generating function of the associated counting process.

This paper presents a degree-two probability lower bound for the union of an arbitrary set of events in an arbitrary probability space. The bound is designed in terms of the first-degree Bonferroni summation and pairwise joint probabilities of events, which are represented as weights of edges in a Hamilton-type circuit in a connected graph. The proposed lower bound strengthens the Dawson–Sankoff lower bound in the same way that Hunter and Worsley's degree-two upper bound improves the degree-two Bonferroni-type optimal upper bound. It can be applied to statistical inference in time series and outlier diagnoses as well as the study of dose response curves.

We extend the work of Antunović et al. (2016) on competing types in preferential attachment models to include cases where the types have different fitnesses, which may be either multiplicative or additive. We show that, depending on the values of the parameters of the models, there are different possible limiting behaviours depending on the zeros of a certain function. In particular, we show the existence of choices of the parameters where one type is favoured both by having higher fitness and by the type of attachment mechanism, but the other type has a positive probability of dominating the network in the limit.

The cumulative distribution of the finite sum of the binary sequence of order k is studied and some of its applications discussed. Certain properties of this sequence are investigated and uniformly superior bounds for the cumulative distribution under minimal information on the ‘success' probabilities are derived.

An alternative proof of a result of Blomqvist on the “heavy-traffic” problem for the two-barrier random walks using weak convergence methods is indicated.

We study a simple random process in which vertices of a connected graph reach consensus through pairwise interactions. We compute outcome probabilities, which do not depend on the graph structure, and consider the expected time until a consensus is reached. In some cases we are able to show that this is minimised by Kn. We prove an upper bound for the p=0 case and give a family of graphs which asymptotically achieve this bound. In order to obtain the mean of the waiting time we also study a gambler's ruin process with delays. We give the mean absorption time and prove that it monotonically increases with p∈[0,1∕2] for symmetric delays.

Let (Xt) and (Yt) be continuous-time Markov chains with countable state spaces E and F and let K be an arbitrary subset of E x F. We give necessary and sufficient conditions on the transition rates of (Xt) and (Yt) for the existence of a coupling which stays in K. We also show that when such a coupling exists, it can be chosen to be Markovian and give a way to construct it. In the case E=F and K ⊆ E x E, we see how the problem of construction of the coupling can be simplified. We give some examples of use and application of our results, including a new concept of lumpability in Markov chains.

Birth–death processes form a natural class where ideas and results on large deviations can be tested. We derive a large-deviation principle under an assumption that the rate of jump down (death) grows asymptotically linearly with the population size, while the rate of jump up (birth) grows sublinearly. We establish a large-deviation principle under various forms of scaling of the underlying process and the corresponding normalization of the logarithm of the large-deviation probabilities. The results show interesting features of dependence of the rate functional upon the parameters of the process and the forms of scaling and normalization.

We show in this paper that if a stationary traffic source is regulated by a leaky bucket with leak rate ρ and bucket size σ, then the amount of information generated in successive time intervals is dominated, in the increasing convex ordering sense, by that of a Poisson arrival process with rate ρ/σ, with each arrival bringing an amount of information equal to σ. By exploiting this property, we then show that the mean value in the stationary regime of the content of a buffer drained at constant rate and fed with the superposition of regulated flows is less than the mean value of the same buffer fed with an adequate Poisson process, whose characteristics depend upon the regulated input flows.

We study a Pólya-type urn model defined as follows. Start at time 0 with a single ball of some colour. Then, at each time n≥1, choose a ball from the urn uniformly at random. With probability ½<p<1, return the ball to the urn along with another ball of the same colour. With probability 1−p, recolour the ball to a new colour and then return it to the urn. This is equivalent to the supercritical case of a random graph model studied by Backhausz and Móri (2015), (2016) and Thörnblad (2015). We prove that, with probability 1, there is a dominating colour, in the sense that, after some random but finite time, there is a colour that always has the most number of balls. A crucial part of the proof is the analysis of an urn model with two colours, in which the observed ball is returned to the urn along with another ball of the same colour with probability p, and removed with probability 1−p. Our results here generalise a classical result about the Pólya urn model (which corresponds to p=1).

We construct a process with inverse gamma increments and an asymptotically self-similar limit. This construction supports the use of long-range-dependent t subordinator models for actual financial data as advocated in Heyde and Leonenko (2005), in that it allows for noninteger-valued model parameters, as is found empirically in model estimation from data.

We consider an r component system having an arbitrary binary monotone structure function. We suppose that shocks occur according to a point process and that, independent of what has already occurred, each new shock is one of r different types, with respective probabilities p1, …, pr. We further suppose that there are given integers n1, …, nr such that component i fails (and remains failed) when there have been a total of ni type-i shocks. Letting L be the time at which the system fails, we are interested in using simulation to estimate E[L], E[L2], and P(L > t). We show how to efficiently accomplish this when the point process is (i) a Poisson, (ii) a renewal, and (iii) a Hawkes process.