A Markov-modulated Poisson process (MMPP) is a Poisson process whose rate is a finite Markov chain. The Poisson process is a simple MMPP. An MMPP/M/1 queue is a queue with MMPP arrivals, an infinite capacity, and a single exponential server. We prove that the output of an MMPP/M/1 queue is not an MMPP process unless the input is Poisson. We derive this result by analyzing the structure of the non-linear filter of the state given the departure process of the queue. The practical relevance of the result is that it rules out the existence of simple finite descriptions of queueing networks with MMPP inputs.

In this paper we prove a quenched functional central limit theorem for a biased random walk on a supercritical Galton–Watson tree with leaves. This extends a result of Peres and Zeitouni (2008) where the case without leaves was considered. A conjecture of Ben Arous and Fribergh (2016) suggests an upper bound on the bias which we observe to be sharp.

The aim of this paper is to study genealogical processes in a geographically structured population with weak migration. The coalescence time for sampled genes from different colonies diverges to infinity as the migration rates among colonies are close to zero. We investigate the moment generating functions of the coalescence time, the number of segregating sites and the number of allele types in sampled genes when there is low migration. Employing a perturbation method, we obtain a system of recurrence relations for the approximate solutions of these moment generating functions and solve them in some cases.

Christ and Avi-Itzhak (2002) analyzed a queueing system with two competing servers who determine their service rates so as to optimize their individual utilities. The system is formulated as a two-person game; Christ and Avi-Itzhak proved the existence of a unique Nash equilibrium which is symmetric. In this paper, we explore globally optimal solutions. We prove that the unique Nash equilibrium is generally strictly inferior to a globally optimal solution and that optimal solutions are symmetric and require the servers to adopt service rates that are smaller than those occurring in equilibrium. Furthermore, given a symmetric globally optimal solution, we show how to impose linear penalties on the service rates so that the given optimal solution becomes a unique Nash equilibrium. When service rates are not observable, we show how the same effect is achieved by imposing linear penalties on a corresponding signal.

A Markov chain X with finite state space {0,…,N} and tridiagonal transition matrix is considered, where transitions from i to i-1 occur with probability (i/N)(1-p(i/N)) and transitions from i to i+1 occur with probability (1-i/N)p(i/N). Here p:[0,1]→[0,1] is a given function. It is shown that if p is continuous with p(x)≤p(1) for all x∈[0,1] then, for each N, a dual process Y to X (with respect to a specific duality function) exists if and only if 1-p is completely monotone with p(0)=0. A probabilistic interpretation of Y in terms of an ancestral process of a mixed multitype Moran model with a random number of types is presented. It is shown that, under weak conditions on p, the process Y, properly time and space scaled, converges to an Ornstein–Uhlenbeck process as N tends to ∞. The asymptotics of the stationary distribution of Y is studied as N tends to ∞. Examples are presented involving selection mechanisms. results.

For independent exponentially distributed random variables $X_i$ , $i\in {\mathcal{N}}$ , with distinct rates ${\lambda}_i$ we consider sums $\sum_{i\in\mathcal{A}} X_i$ for $\mathcal{A}\subseteq {\mathcal{N}}$ which follow generalized exponential mixture distributions. We provide novel explicit results on the conditional distribution of the total sum $\sum_{i\in {\mathcal{N}}}X_i$ given that a subset sum $\sum_{j\in \mathcal{A}}X_j$ exceeds a certain threshold value $t>0$ , and vice versa. Moreover, we investigate the characteristic tail behavior of these conditional distributions for $t\to\infty$ . Finally, we illustrate how our probabilistic results can be applied in practice by providing examples from both reliability theory and risk management.

In an incomplete financial market in which the dynamics of the asset prices is driven by a d-dimensional continuous semimartingale X, we consider the problem of pricing European contingent claims embedded in a power utility framework. This problem reduces to identifying the p-optimal martingale measure, which can be given in terms of the solution to a semimartingale backward equation. We use this characterization to examine two extreme cases. In particular, we find a necessary and sufficient condition, written in terms of the mean-variance trade-off, for the p-optimal martingale measure to coincide with the minimal martingale measure. Moreover, if and only if an exponential function of the mean-variance trade-off is a martingale strongly orthogonal to the asset price process, the p-optimal martingale measure can be simply expressed in terms of a Doléans-Dade exponential involving X.

Block, Savits and Singh previously presented a preservation result for bathtub-shaped functions. In this paper, we strengthen and extend their results to a class of functions which allow for more than one change of monotonicity. Various applications are also considered.

Determinantal point processes (DPPs) enable the modeling of repulsion: they provide diverse sets of points. The repulsion is encoded in a kernel K that can be seen, in a discrete setting, as a matrix storing the similarity between points. The main exact algorithm to sample DPPs uses the spectral decomposition of K, a computation that becomes costly when dealing with a high number of points. Here we present an alternative exact algorithm to sample in discrete spaces that avoids the eigenvalues and the eigenvectors computation. The method used here is innovative, and numerical experiments show competitive results with respect to the initial algorithm.

This paper proposes a class of complete financial market models, the benchmark models, with security price processes that exhibit intensity-based jumps. The benchmark or reference unit is chosen to be the growth-optimal portfolio. Primary security account prices, when expressed in units of the benchmark, turn out to be local martingales. In the proposed framework an equivalent risk-neutral measure need not exist. Benchmarked fair derivative prices are obtained as conditional expectations of future benchmarked prices under the real-world probability measure. This concept of fair pricing generalizes the classical risk-neutral approach and the actuarial present-value pricing methodology.

We consider a finite-capacity queueing system where arriving customers offer rewards which are paid upon acceptance into the system. The gatekeeper, whose objective is to ‘maximize’ rewards, decides if the reward offered is sufficient to accept or reject the arriving customer. Suppose the arrival rates, service rates, and system capacity are changing over time in a known manner. We show that all bias optimal (a refinement of long-run average reward optimal) policies are of threshold form. Furthermore, we give sufficient conditions for the bias optimal policy to be monotonic in time. We show, via a counterexample, that if these conditions are violated, the optimal policy may not be monotonic in time or of threshold form.

For stationary sequences X = {Xn}n≥1 we relate τ, the limiting mean number of exceedances of high levels un by X1,…,Xn, and ν, the limiting mean number of upcrossings of the same level, through the expression θ = (ν/τ)η, where θ is the extremal index of X and η is a new parameter here called the upcrossings index. The upcrossings index is a measure of the clustering of upcrossings of u by variables in X, and the above relation extends the known relation θ = ν/τ, which holds under the mild-oscillation local restriction D″(u) on X. We present a new family of local mixing conditions D̃(k)(u) under which we prove that (a) the intensity of the limiting point process of upcrossings and η can both be computed from the k-variate distributions of X; and (b) the cluster size distributions for the exceedances are asymptotically equivalent to those for the lengths of one run of exceedances or the lengths of several consecutive runs which are separated by at most k − 2 nonexceedances and, except for the last one, each contain at most k − 2 exceedances.

It can be seen that a mixture of an exponential distribution and a gamma distribution with increasing failure rate for the right choice of parameters can yield a distribution with a bathtub-shaped failure rate. In this paper we consider a continuous mixture of exponentials and a continuous mixture of gammas with increasing failure rates and show that the resulting mixture has a bathtub-shaped failure rate.

We present bounds on the decay parameter for absorbing birth–death processes adapted from results of Chen (2000), (2001). We address numerical issues associated with computing these bounds, and assess their accuracy for several models, including the stochastic logistic model, for which estimates of the decay parameter have been obtained previously by Nåsell (2001).

In this paper we investigate the number and maximum severity of the ruin excursion of the insurance portfolio reserve process in the Cramér–Lundberg model with and without tax payments. We also provide a relation of the Cramér–Lundberg risk model with the G/G/∞ queue and use it to derive some explicit ruin probability formulae. Finally, the renewal risk model with tax is considered, and an asymptotic identity is derived that in some sense extends the tax identity of the Cramér– Lundberg risk model.

Existing results on the asymptotic behaviour of solutions to two dynamic programming equations involving non-negative matrices are reviewed and strengthened in certain directions. The results are then applied to strategies for harvesting of a small population so as to optimise its survival potential in a limited environment.

Inspired by a PDE–ODE system of aggregation developed in the biomathematical literature, we investigate an interacting particle system representing aggregation at the level of individuals. We prove that the empirical density of the individual converges to the solution of the PDE–ODE system.

Let {W(t), 0 ≦ t < ∞} be the standard Wiener process. The probabilities of the type P[sup0≦t ≦ TW(t) − f(t) ≧ 0] have been extensively studied when f(t) is a deterministic function. This paper discusses the probabilities of the type P{sup0≦t ≦ TW(t) − [f(t) + X(t)] ≧ 0} when X(t) is a stochastic process. By taking compound Poisson processes as X(t), the paper gives procedures for finding such probabilities.

Motivated by a problem in neural encoding, we introduce an adaptive (or real-time) parameter estimation algorithm driven by a counting process. Despite the long history of adaptive algorithms, this kind of algorithm is relatively new. We develop a finite-time averaging analysis which is nonstandard partly because of the point process setting and partly because we have sought to avoid requiring mixing conditions. This is significant since mixing conditions often place restrictive history-dependent requirements on algorithm convergence.