Recently, the properties of ESSs for random payoff matrices, and stable polymorphisms for random fitness matrices, have been investigated. These problems are very closely related, and some similarities and differences are examined here. In contrast to some earlier work, non-robustness of the conclusions, to changes in the distribution of the random matrices, are found in two areas: in the asymptotic number of polymorphisms (or ESSs) whose support is a given size, and in the location of a polymorphism (or ESS) whose size of support is 2 or 3, in a large matrix.

We consider an M/GI/1 queue with two types of customers, positive and negative, which cancel each other out. The server provides service to either a positive customer or a negative customer. In such a system, the queue length can be either positive or negative and an arrival either joins the queue, if it is of the same sign, or instantaneously removes a customer of the opposite sign at the end of the queue or in service. This study is a generalization of Gelenbe's original concept of a queue with negative customers, where only positive customers need services and negative customers arriving at an empty system are lost or need no service. In this paper, we derive the transient and the stationary probability distributions for the major performance measures in terms of generating functions and Laplace transforms. It has been shown that the previous results for the system with negative arrivals of zero service time are special cases of our model. In addition, we obtain the stationary waiting time distribution of this system in terms of a Laplace transform.

A generating function technique is used to determine the probability that the deviation between two empirical distributions drawn from the same population lies within a given band a specified number of times. We also treat the asymptotic problem of very large sample size, and obtain explicit expressions when the relative number of failures is very small or very large.

Using a representation in terms of a two-type branching particle system, we prove that positive solutions of the system remain bounded for suitable bounded initial conditions, provided A and B generate processes with independent increments and one of the processes is transient with a uniform power decay of its semigroup. For the case of symmetric stable processes on R1,this answers a question raised in [4].

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.

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.

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.

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.

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 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 lumping of a Markov chain is a coordinatewise projection of the chain. We characterise the entropy rate preservation of a lumping of an aperiodic and irreducible Markov chain on a finite state space by the random growth rate of the cardinality of the realisable preimage of a finite-length trajectory of the lumped chain and by the information needed to reconstruct original trajectories from their lumped images. Both are purely combinatorial criteria, depending only on the transition graph of the Markov chain and the lumping function. A lumping is strongly k-lumpable, if and only if the lumped process is a kth-order Markov chain for each starting distribution of the original Markov chain. We characterise strong k-lumpability via tightness of stationary entropic bounds. In the sparse setting, we give sufficient conditions on the lumping to both preserve the entropy rate and be strongly k-lumpable.

In this paper we are interested in a nonlinear parabolic evolution equation occurring in rheology. We give a probabilistic interpretation to this equation by associating a nonlinear martingale problem with it. We prove the existence of a unique solution, P, to this martingale problem. For any t, the time marginal of P at time t admits a density ρ(t,x) with respect to the Lebesgue measure, where the function ρ is the unique weak solution to the evolution equation in a well-chosen energy space. Next we introduce a simulable system of n interacting particles and prove that the empirical measure of this system converges to P as n tends to ∞. This propagation-of-chaos result ensures that the solution to the equation of interest can be approximated using a Monte Carlo method. Finally, we illustrate the convergence in some numerical experiments.

The autopsy data of a coherent system is the pair (T, I), where T is the time to the failure of the system, and I is the set of components failed by that time. A basic question is under which conditions the distribution of the autopsy data determines the joint life distribution of the components. This problem has been studied under various distributional and structural assumptions (always under the assumption of independent component life distributions) by Meilijson (1981), Nowik (1990) and Antoine, Doss and Hollander (1993). In this paper we study the same problem in the situation when the autopsy data is supplemented with data from lifemonitoring and conditional lifemonitoring of some components, as introduced by Meilijson (1994).

We consider a space-time random field on ${{\mathbb{R}^d} \times {\mathbb{R}}}$ given as an integral of a kernel function with respect to a Lévy basis with a convolution equivalent Lévy measure. The field obeys causality in time and is thereby not continuous along the time axis. For a large class of such random fields we study the tail behaviour of certain functionals of the field. It turns out that the tail is asymptotically equivalent to the right tail of the underlying Lévy measure. Particular examples are the asymptotic probability that there is a time point and a rotation of a spatial object with fixed radius, in which the field exceeds the level x, and that there is a time interval and a rotation of a spatial object with fixed radius, in which the average of the field exceeds the level x.

We study the genealogy of so-called immortal branching processes, i.e. branching processes where each individual upon death is replaced by at least one new individual, and conclude that their marginal distributions are compound geometric. The result also implies that the limiting distributions of properly scaled supercritical branching processes are compound geometric. We exemplify our results with an expression for the marginal distribution for a class of branching processes that have recently appeared in the theory of coalescent processes and continuous stable random trees. The limiting distribution can be expressed in terms of the Fox H-function, and in special cases by the Meijer G-function.

Let X, B, and Y be the Dirichlet, Bernoulli, and beta-independent random variables such that X ~ D(a0, …, ad), Pr(B = (0, …, 0, 1, 0, …, 0)) = ai / a with a = ∑i=0dai, and Y ~ β(1, a). Then, as proved by Sethuraman (1994), X ~ X(1 - Y) + BY. This gives the stationary distribution of a simple Markov chain on a tetrahedron. In this paper we introduce a new distribution on the tetrahedron called a quasi-Bernoulli distribution Bk(a0, …, ad) with k an integer such that the above result holds when B follows Bk(a0, …, ad) and when Y ~ β(k, a). We extend it even more generally to the case where X and B are random probabilities such that X is Dirichlet and B is quasi-Bernoulli. Finally, the case where the integer k is replaced by a positive number c is considered when a0 = · · · = ad = 1.

We apply the known formulae of the RESTART problem to Markov models of software (and many other) systems, and derive new equations. We show how checkpoints might be included, with their resultant performance under RESTART. The result is a complete procedure for finding the mean, variance, and tail behavior of the job completion time as a function of the failure rate. We also provide a detailed example.

We propose a method to generate a max-stable process in C[0, 1] from a max-stable random vector in Rd by generalizing the max-linear model established by Wang and Stoev (2011). For this purpose, an interpolation technique that preserves max-stability is proposed. It turns out that if the random vector follows some finite-dimensional distribution of some initial max-stable process, the approximating processes converge uniformly to the original process and the pointwise mean-squared error can be represented in a closed form. The obtained results carry over to the case of generalized Pareto processes. The introduced method enables the reconstruction of the initial process only from a finite set of observation points and, thus, a reasonable prediction of max-stable processes in space becomes possible. A possible extension to arbitrary dimensions is outlined.

We present an elegant uniqueness criterion for the weighted Markov branching process in the potentially explosive case.