We present and investigate a general model for inhomogeneous random digraphs with labeled vertices, where the arcs are generated independently, and the probability of inserting an arc depends on the labels of its endpoints and on its orientation. For this model, the critical point for the emergence of a giant component is determined via a branching process approach.

We prove that a stochastic process of pure coagulation has at any time t ≥ 0 a time-dependent Gibbs distribution if and only if the rates ψ(i, j) of single coagulations are of the form ψ(i; j) = if(j) + jf(i), where f is an arbitrary nonnegative function on the set of positive integers. We also obtain a recurrence relation for weights of these Gibbs distributions that allow us to derive the general form of the solution and the explicit solutions in three particular cases of the function f. For the three corresponding models, we study the probability of coagulation into one giant cluster by time t > 0.

In a Galton-Watson branching process that is not extinct by the nth generation and has at least two individuals, pick two individuals at random by simple random sampling without replacement. Trace their lines of descent back in time till they meet. Call that generation Xn a pairwise coalescence time. Similarly, let Yn denote the coalescence time for the whole population of the nth generation conditioned on the event that it is not extinct. In this paper the distributions of Xn and Yn, and their limit behaviors as n → ∞ are discussed for both the critical and subcritical cases.

We focus on supercritical decomposable (reducible) multitype branching processes. Types are partitioned into irreducible equivalence classes. In this context, extinction of some classes is possible without the whole process becoming extinct. We derive criteria for the almost-sure extinction of the whole process, as well as of a specific class, conditionally given the class of the initial particle. We give sufficient conditions under which the extinction of a class implies the extinction of another class or of the whole process. Finally, we show that the extinction probability of a specific class is the minimal nonnegative solution of the usual extinction equation but with added constraints.

We study generalized branching random walks on the real line R that allow time dependence and local dependence between siblings. Specifically, starting from one particle at time 0, the system evolves such that each particle lives for one unit amount of time, gives birth independently to a random number of offspring according to some branching law, and dies. The offspring from a single particle are assumed to move to new locations on R according to some joint displacement distribution; the branching laws and displacement distributions depend on time. At time n, Fn(·) is used to denote the distribution function of the position of the rightmost particle in generation n. Under appropriate tail assumptions on the branching laws and offspring displacement distributions, we prove that Fn(· - Med(Fn)) is tight in n, where Med(Fn) is the median of Fn. The main part of the argument is to demonstrate the exponential decay of the right tail 1 - Fn(· - Med(Fn)).

We analyse the behaviour of supercritical super-Brownian motion with a barrier through the pathwise backbone embedding of Berestycki, Kyprianou and Murillo-Salas (2011). In particular, by considering existing results for branching Brownian motion due to Harris and Kyprianou (2006) and Maillard (2011), we obtain, with relative ease, conclusions regarding the growth in the right-most point in the support, analytical properties of the associated one-sided Fisher-Kolmogorov-Petrovskii-Piscounov wave equation, as well as the distribution of mass on the exit measure associated with the barrier.

Nearly fifty years after the introduction of skew Brownian motion by Itô and McKean (1963), the first passage time distribution remains unknown. In this paper we first generalize results of Pitman and Yor (2011) and Csáki and Hu (2004) to derive formulae for the distribution of ranked excursion heights of skew Brownian motion, and then use these results to derive the first passage time distribution.

Let (X, J) denote a Markov-modulated Brownian motion (MMBM) and denote its supremum process by S. For some a > 0, let σ(a) denote the time when the reflected process Y := S - X first surpasses the level a. Furthermore, let σ−(a) denote the last time before σ(a) when X attains its current supremum. In this paper we shall derive the joint distribution of Sσ(a), σ−(a), and σ(a), where the latter two will be given in terms of their Laplace transforms. We also provide some remarks on scale matrices for MMBMs with strictly positive variation parameters. This extends recent results for spectrally negative Lévy processes to MMBMs. Due to well-known fluid embedding and state-dependent killing techniques, the analysis applies to Markov additive processes with phase-type jumps as well. The result is of interest to applications such as the dividend problem in insurance mathematics and the buffer overflow problem in queueing theory. Examples will be given for the former.

In this paper we show that fractional Brownian motion with H < ½ can arise as a limit of a simple class of traffic processes that we call ‘scheduled traffic models’. To our knowledge, this paper provides the first simple traffic model leading to fractional Brownnian motion with H < ½. We also discuss some immediate implications of this result for queues fed by scheduled traffic, including a heavy-traffic limit theorem.

We study quasistationary distributions on a drifted Brownian motion killed at 0, when +∞ is an entrance boundary and 0 is an exit boundary. We prove the existence of a unique quasistationary distribution and of the Yaglom limit.

Suppose that both you and your friend toss an unfair coin n times, for which the probability of heads is equal to α. What is the probability that you obtain at least d more heads than your friend if you make r additional tosses? We obtain asymptotic and monotonicity/convexity properties for this competing probability as a function of n, and demonstrate surprising phase transition phenomenon as the parameters d, r, and α vary. Our main tools are integral representations based on Fourier analysis.

Using Stein's method, we derive explicit upper bounds on the total variation distance between a Poisson-binomial law (the distribution of a sum of independent but not necessarily identically distributed Bernoulli random variables) and a Pólya distribution with the same support, mean, and variance; a nonuniform bound on the pointwise distance between the probability mass functions is also given. A numerical comparison of alternative distributional approximations on a somewhat representative collection of case studies is also exhibited. The evidence proves that no single one is uniformly most accurate, though it suggests that the Pólya approximation might be preferred in several parameter domains encountered in practice.

An infinite sequence (Y1, Y2,…) of independent Bernoulli random variables with P(Yi = 1) = a / (a + b + i - 1), i = 1, 2,…, where a > 0 and b ≥ 0, will be called a Bern(a, b) sequence. Consider the counts Z1, Z2, Z3,… of occurrences of patterns or strings of the form {11}, {101}, {1001},…, respectively, in this sequence. The joint distribution of the counts Z1, Z2,… in the infinite Bern(a, b) sequence has been studied extensively. The counts from the initial finite sequence (Y1, Y2,…, Yn) have been studied by Holst (2007), (2008b), who obtained the joint factorial moments for Bern(a, 0) and the factorial moments of Z1, the count of the string {1, 1}, for a general Bern(a, b). We consider stopping the Bernoulli sequence at a random time and describe the joint distribution of counts, which extends Holst's results. We show that the joint distribution of counts from a class of randomly stopped Bernoulli sequences possesses the mixture of independent Poissons property: there is a random vector conditioned on which the counts are independent Poissons. To obtain these results, we extend the conditional marked Poisson process technique introduced in Huffer, Sethuraman and Sethuraman (2009). Our results avoid previous combinatorial and induction methods which generally only yield factorial moments.

Hamilton's method is a natural and common method to distribute seats proportionally between states (or parties) in a parliament. In the USA it has been abandoned due to some drawbacks, in particular the possibility of the Alabama paradox, but it is still in use in many other countries. In this paper we give, under certain assumptions, a closed formula for the asymptotic probability, as the number of seats tends to infinity, that the Alabama paradox occurs given the vector p1,…, pm of relative sizes of the states. From the formula we deduce a number of consequences. For example, the expected number of states that will suffer from the Alabama paradox is asymptotically bounded above by 1 / e and on average approximately 0.123.

Estimator algorithms in learning automata are useful tools for adaptive, real-time optimization in computer science and engineering applications. In this paper we investigate theoretical convergence properties for a special case of estimator algorithms - the pursuit learning algorithm. We identify and fill a gap in existing proofs of probabilistic convergence for pursuit learning. It is tradition to take the pursuit learning tuning parameter to be fixed in practical applications, but our proof sheds light on the importance of a vanishing sequence of tuning parameters in a theoretical convergence analysis.

Let (Xt)0 ≤ t ≤ T be a one-dimensional stochastic process with independent and stationary increments, either in discrete or continuous time. In this paper we consider the problem of stopping the process (Xt) ‘as close as possible’ to its eventual supremum MT := sup0 ≤ t ≤ TXt, when the reward for stopping at time τ ≤ T is a nonincreasing convex function of MT - Xτ. Under fairly general conditions on the process (Xt), it is shown that the optimal stopping time τ takes a trivial form: it is either optimal to stop at time 0 or at time T. For the case of a random walk, the rule τ ≡ T is optimal if the steps of the walk stochastically dominate their opposites, and the rule τ ≡ 0 is optimal if the reverse relationship holds. An analogous result is proved for Lévy processes with finite Lévy measure. The result is then extended to some processes with nonfinite Lévy measure, including stable processes, CGMY processes, and processes whose jump component is of finite variation.

The classical secretary problem for selecting the best item is studied when the actual values of the items are observed with noise. One of the main appeals of the secretary problem is that the optimal strategy is able to find the best observation with a nontrivial probability of about 0.37, even when the number of observations is arbitrarily large. The results are strikingly different when the qualities of the secretaries are observed with noise. If there is no noise then the only information that is needed is whether an observation is the best among those already observed. Since the observations are assumed to be independent and identically distributed, the solution to this problem is distribution free. In the case of noisy data, the results are no longer distribution free. Furthermore, we need to know the rank of the noisy observation among those already observed. Finally, the probability of finding the best secretary often goes to 0 as the number of observations, n, goes to ∞. The results heavily depend on the behavior of pn, the probability that the observation that is best among the noisy observations is also best among the noiseless observations. Results involving optimal strategies if all that is available is noisy data are described and examples are given to elucidate the results.

In this paper we propose a class of financial market models which are based on telegraph processes with alternating tendencies and jumps. It is assumed that the jumps have random sizes and that they occur when the tendencies are switching. These models are typically incomplete, but the set of equivalent martingale measures can be described in detail. We provide additional suggestions which permit arbitrage-free option prices as well as hedging strategies to be obtained.

In this paper we study the distribution of the location, at time t, of a particle moving U time units upwards, V time units downwards, and W time units of no movement (idle). These are repeated cyclically, according to independent alternating renewals. The distributions of U, V, and W are absolutely continuous. The velocities are v = +1 upwards, v = -1 downwards, and v = 0 during idle periods. Let Y+(t), Y−(t), and Y0(t) denote the total time in (0, t) of movements upwards, downwards, and no movements, respectively. The exact distribution of Y+(t) is derived. We also obtain the probability law of X(t) = Y+(t) - Y−(t), which describes the particle's location at time t. Explicit formulae are derived for the cases of exponential distributions with equal rates, with different rates, and with linear rates (leading to damped processes).

Tankov (2011) improved the Fréchet bounds for a bivariate copula when its values on a compact subset of [0, 1]2 are given. He showed that the best possible bounds are quasi-copulas and gave a sufficient condition for these bounds to be copulas. In this note we give weaker sufficient conditions to ensure that the bounds are copulas. We also show how this can be useful in portfolio selection. It turns out that finding a copula as a lower bound plays a key role in determining optimal investment strategies explicitly for investors with some type of state-dependent constraints.