In this paper we analyze a tollbooth tandem queueing problem with an infinite number of servers. A customer starts service immediately upon arrival but cannot leave the system before all customers who arrived before him/her have left, i.e. customers depart the system in the same order as they arrive. Distributions of the total number of customers in the system, the number of departure-delayed customers in the system, and the number of customers in service at time t are obtained in closed form. Distributions of the sojourn times and departure delays of customers are also obtained explicitly. Both transient and steady state solutions are derived first for Poisson arrivals, and then extended to cases with batch Poisson and nonstationary Poisson arrival processes. Finally, we report several stochastic ordering results on how system performance measures are affected by arrival and service processes.

A gambler starts with fortune f < 1 and plays in a Vardi casino with infinitely many tables indexed by their odds, r ≥ 0. In addition, all tables return the same expected winnings per dollar, c < 0, and a discount factor is applied after each round. We determine the optimal probability of reaching fortune 1, as well as an optimal strategy that is different from bold play for fortunes larger than a critical value depending exclusively on c and 1 + a, the discount factor. The general result is computed explicitly for some relevant special cases. The question of whether bold play is an optimal strategy is discussed for various choices of the parameters.

Consider two sequences of bounded random variables, a value and a timing process, that satisfy the large deviation principle (LDP) with rate function J(⋅,⋅) and whose cumulative process satisfies the LDP with rate function I(⋅). Under mixing conditions, an LDP for estimates of I constructed by transforming an estimate of J is proved. For the case of a cumulative renewal process it is demonstrated that this approach is favourable to a more direct method, as it ensures that the laws of the estimates converge weakly to a Dirac measure at I.

We consider a class of optimal stopping problems involving both the running maximum as well as the prevailing state of a linear diffusion. Instead of tackling the problem directly via the standard free boundary approach, we take an alternative route and present a parameterized family of standard stopping problems of the underlying diffusion. We apply this family to delineate circumstances under which the original problem admits a unique, well-defined solution. We then develop a discretized approach resulting in a numerical algorithm for solving the considered class of stopping problems. We illustrate the use of the algorithm in both a geometric Brownian motion and a mean reverting diffusion setting.

We consider a discrete-time population growth system called the Bienaymé–Galton–Watson stochastic branching system. We deal with a noncritical case, in which the per capita offspring mean $m\neq1$. The famous Kolmogorov theorem asserts that the expectation of the population size in the subcritical case $m<1$ on positive trajectories of the system asymptotically stabilizes and approaches ${1}/\mathcal{K}$, where $\mathcal{K}$ is called the Kolmogorov constant. The paper is devoted to the search for an explicit expression of this constant depending on the structural parameters of the system. Our argumentation is essentially based on the basic lemma describing the asymptotic expansion of the probability-generating function of the number of individuals. We state this lemma for the noncritical case. Subsequently, we find an extended analogue of the Kolmogorov constant in the noncritical case. An important role in our discussion is also played by the asymptotic properties of transition probabilities of the Q-process and their convergence to invariant measures. Obtaining the explicit form of the extended Kolmogorov constant, we refine several limit theorems of the theory of noncritical branching systems, showing explicit leading terms in the asymptotic expansions.

A set of data with positive values follows a Pareto distribution if the log–log plot of value versus rank is approximately a straight line. A Pareto distribution satisfies Zipf’s law if the log–log plot has a slope of $-1$. Since many types of ranked data follow Zipf’s law, it is considered a form of universality. We propose a mathematical explanation for this phenomenon based on Atlas models and first-order models, systems of strictly positive continuous semimartingales with parameters that depend only on rank. We show that the stationary distribution of an Atlas model will follow Zipf’s law if and only if two natural conditions, conservation and completeness, are satisfied. Since Atlas models and first-order models can be constructed to approximate systems of time-dependent rank-based data, our results can explain the universality of Zipf’s law for such systems. However, ranked data generated by other means may follow non-Zipfian Pareto distributions. Hence, our results explain why Zipf’s law holds for word frequency, firm size, household wealth, and city size, while it does not hold for earthquake magnitude, cumulative book sales, and the intensity of wars, all of which follow non-Zipfian Pareto distributions.

We propose a stochastic model for the failure times of items subject to two external random shocks occurring as events in an underlying bivariate counting process. This is a special formulation of the competing risks model, which is of interest in reliability theory and survival analysis. Specifically, we assume that a system, or an item, fails when the sum of the two types of shock reaches a critical random threshold. In detail, the two kinds of shock occur according to a bivariate space-fractional Poisson process, which is a two-dimensional vector of independent homogeneous Poisson processes time-changed by an independent stable subordinator. Various results are given, such as analytic hazard rates, failure densities, the probability that the failure occurs due to a specific type of shock, and the survival function. Some special cases and ageing notions related to the NBU characterization are also considered. In this way we generalize certain results in the literature, which can be recovered when the underlying process reduces to the homogeneous Poisson process.

We consider two random sequential packing processes in which spheres of unit radius are randomly attached to the surface of a fixed unit sphere. Independent random spheres are generated and added successively, provided there is no overlap with previous spheres. In model 1, the process stops when a trial sphere intersects one of the previously-accepted spheres. In model 2, random sequential packing, any such overlapping trial sphere is discarded and the next random sphere is tried, until it is impossible to add any further spheres.

Previous workers have conjectured convincingly that no exact analytical solution is possible for this type of problem. We use Monte Carlo simulation methods to estimate transition probabilities for the two models. Because some probabilities are extremely small, a simulation using independent repetitions of the model would be inefficient. We designed a branching process of conditionally binomial trials, and performed over 108 trials on a supercomputer.

We introduce the concept of the absorbing process for analysing a state process. Our aim is to show the existence of the absorbing process with probability one. This process is shown to be stationary, asymptotically stationary, periodic or a.m.s., if the input distribution has such properties. The real process is absorbed into this process so that its stability and some other properties are easily derived.

Modeling of manufacturing lines and data communications makes frequent use of tandem queues with blocking. Here we present and study such a system with a k-stage blocking scheme in which processing in each station requires attendance of the server of that station together with the servers of the next k – 1 stations. This scheme describes the conventional manufacturing and communications blocking schemes but is also representative of a wider range of applications. Explicit expressions for residence times, throughput, queueing times, waiting times and other measures of performance are obtained for the case of communications type flow and just-in-time input. Use of the results in modeling and analyzing a highway traffic flow situation is presented.

We construct a multitype constant-size population model allowing for general selective interactions as well as extreme reproductive events. Our multidimensional model aims for the generality of adaptive dynamics and the tractability of population genetics. It generalises the idea of Krone and Neuhauser [39] and González Casanova and Spanò [29], who represented the selection by allowing individuals to sample several potential parents in the previous generation before choosing the ‘strongest’ one, by allowing individuals to use any rule to choose their parent. The type of the newborn can even not be one of the types of the potential parents, which allows modelling mutations. Via a large population limit, we obtain a generalisation of $\Lambda$ -Fleming–Viot processes, with a diffusion term and a general frequency-dependent selection, which allows for non-transitive interactions between the different types present in the population. We provide some properties of these processes related to extinction and fixation events, and give conditions for them to be realised as unique strong solutions of multidimensional stochastic differential equations with jumps. Finally, we illustrate the generality of our model with applications to some classical biological interactions. This framework provides a natural bridge between two of the most prominent modelling frameworks of biological evolution: population genetics and eco-evolutionary models.

Let Xn(k) be the number of vertices at level k in a random recursive tree with n+1 vertices. We are interested in the asymptotic behavior of Xn(k) for intermediate levels k=kn satisfying kn→∞ and kn=o(logn) as n→∞. In particular, we prove weak convergence of finite-dimensional distributions for the process (Xn ([knu]))u>0, properly normalized and centered, as n→∞. The limit is a centered Gaussian process with covariance (u,v)↦(u+v)−1. One-dimensional distributional convergence of Xn(kn), properly normalized and centered, was obtained with the help of analytic tools by Fuchs et al. (2006). In contrast, our proofs, which are probabilistic in nature, exploit a connection of our model with certain Crump–Mode–Jagers branching processes.

Limit laws and limiting diffusions are obtained for critical branching processes, either Galton-Watson or age-dependent, conditioned on extinction in the interval (T, cT], 1 < c, as T→∞, and also as T→∞ and c ↓ 1.

We consider the convex hull of a sample of n randomly placed points in a unit circle. Obviously not all of the n points are needed to construct the convex hull. We show that asymptotically only the points belonging to a certain ring are used for the convex hull, where the size of the ring decreases at a certain rate such that the error of not completely correct construction of the convex hull by the particular subset of points tends to 0.

The concept of θ-stationarity for a simple second-order point process in R2 is introduced. This concept is closely related to that of isotropy. Some θ-stationary processes are defined. Techniques are given for simulating realisations of these processes. The second-order analysis of these processes which have an obvious point of reference or origin is considered. Methods are suggested for modelling spatial patterns which are realisations of such processes. These methods are illustrated using simulated data. The ideas are extended to multitype point processes.

The Robbins–Monro algorithm with randomly varying truncations for measurements with non-additive noise is considered. Assuming that the function under observation is locally Lipschitz-continuous in its first argument and that the noise is a φ-mixing process, strong consistency of the estimate is shown. Neither growth rate restriction on the function, nor the decreasing rate of the mixing coefficients are required.

In [14] a necessary and sufficient condition was obtained for there to exist uniquely a Q-process with a specified invariant measure, under the assumption that Q is a stable, conservative, single-exit matrix. The purpose of this note is to demonstrate that, for an arbitrary stable and conservative q-matrix, the same condition suffices for the existence of a suitable Q-process, but that this process might not be unique. A range of examples is considered, including pure-birth processes, a birth process with catastrophes, birth-death processes and the Markov branching process with immigration.

Excess life distributions for discrete renewal processes may be computed by using elementary discrete Markov chain concepts involving absorption probabilities. Excess life distributions in general may then be obtained by approximating the renewal process under study by a suitably chosen sequence of discrete renewal processes. The technique is illustrated in the cases of renewal processes with interarrival distributions which are a linear combination of two exponentials and uniform [0, 1]. A related algorithm is described for computer generated approximations of excess life distributions corresponding to continuous interarrival time distributions with an increasing c.d.f. Conditions for convergence of this algorithm are examined.

We consider backward stochastic differential equations in a setting where noise is generated by a countable state, continuous time Markov chain, and the terminal value is prescribed at a stopping time. We show that, given sufficient integrability of the stopping time and a growth bound on the terminal value and BSDE driver, these equations admit unique solutions satisfying the same growth bound (up to multiplication by a constant). This holds without assuming that the driver is monotone in y, that is, our results do not require that the terminal value be discounted at some uniform rate. We show that the conditions are satisfied for hitting times of states of the chain, and hence present some novel applications of the theory of these BSDEs.

A method for the construction of Stein-type covariance identities for a nonnegative continuous random variable is proposed, using a probabilistic analogue of the mean value theorem and weighted distributions. A generalized covariance identity is obtained, and applications focused on actuarial and financial science are provided. Some characterization results for gamma and Pareto distributions are also given. Identities for risk measures which have a covariance representation are obtained; these measures are connected with the Bonferroni, De Vergottini, Gini, and Wang indices. Moreover, under some assumptions, an identity for the variance of a function of a random variable is derived, and its performance is discussed with respect to well-known upper and lower bounds.