We consider two independent homogeneous Poisson processes Π0 and Π1 in the plane with intensities λ0 and λ1, respectively. We study additive functionals of the set of Π0-particles within a typical Voronoi Π1-cell. We find the first and the second moments of these variables as well as upper and lower bounds on their distribution functions, implying an exponential asymptotic behavior of their tails. Explicit formulae are given for the number and the sum of distances from Π0-particles to the nucleus within a typical Voronoi Π1-cell.

In design-based stereology, fixed parameters (such as volume, surface area, curve length, feature number, connectivity) of a non-random geometrical object are estimated by intersecting the object with randomly located and oriented geometrical probes (e.g. test slabs, planes, lines, points). Estimation accuracy may in principle be increased by increasing the number of probes, which are usually laid in a systematic pattern. An important prerequisite to increase accuracy, however, is that the relevant estimators are unbiased and consistent. The purpose of this paper is therefore to give sufficient conditions for the unbiasedness and strong consistency of design-based stereological estimators obtained by systematic sampling. Relevant mechanisms to increase sample size, compatible with stereological practice, are considered.

A rectangular tessellation is a covering of the plane by non-overlapping rectangles. A basic theory for general homogeneous random rectangular tessellations is developed, and it is shown that many first-order mean values may be expressed in terms of just three basic quantities. Corresponding values for independent superpositions of two or more such tessellations are derived. The most interesting homogeneous rectangular tessellations are those with only T-vertices (i.e. no X-vertices). Gilbert's (1967) isotropic model adapted to this two-orthogonal-orientations case, although simply specified, appears theoretically intractable, due to a complex ‘blocking' effect. However, the approximating penetration model, also introduced by Gilbert, is found to be both tractable and informative about the true model. A multi-stage method for simulating the model is developed, and the distributions of important characteristics estimated.

Let θ (a) be the first time when the range (Rn; n ≧ 0) is equal to a, Rn being equal to the difference of the maximum and the minimum, taken at time n, of a simple random walk on ℤ. We compute the g.f. of θ (a); this allows us to compute the distributions of θ (a) and Rn. We also investigate the asymptotic behaviour of θ (n), n going to infinity.

We investigate the ‘clumping versus local finiteness' behavior in the infinite backward tree for a class of branching particle systems in ℝd with symmetric stable migration and critical ‘genuine multitype' branching. Under mild assumptions on the branching we establish, by analysing certain ergodic properties of the individual ancestral process, a critical dimension dc such that the (measure-valued) tree-top is almost surely locally finite if and only if d > dc. This result is used to obtain L1-norm asymptotics of a corresponding class of systems of non-linear partial differential equations.

In phylogenetic analysis it is useful to study the distribution of the parsimony length of a tree under the null model, by which the leaves are independently assigned letters according to prescribed probabilities. Except in one special case, this distribution is difficult to describe exactly. Here we analyze this distribution by providing a recursive and readily computable description, establishing large deviation bounds for the parsimony length of a fixed tree on a single site and for the minimum length (maximum parsimony) tree over several sites. We also show that, under very general conditions, the former distribution converges asymptotically to the normal, thereby settling a recent conjecture. Furthermore, we show how the mean and variance of this distribution can be efficiently calculated. The proof of normality requires a number of new and recent results, as the parsimony length is not directly expressible as a sum of independent random variables, and so normality does not follow immediately from a standard central limit theorem.

In the first part of this paper a global Kushner–Clark theorem about the convergence of stochastic algorithms is proved: we show that, under some natural assumptions, one can ‘read' from the trajectories of its ODE whether or not an algorithm converges. The classical stochastic optimization results are included in this theorem. In the second part, the above smoothness assumption on the mean vector field of the algorithm is relaxed using a new approach based on a path-dependent Lyapounov functional. Several applications, for non-smooth mean vector fields and/or bounded Lyapounov function settings, are derived. Examples and simulations are provided that illustrate and enlighten the field of application of the theoretical results.

We study a classical stochastic control problem arising in financial economics: to maximize expected logarithmic utility from terminal wealth and/or consumption. The novel feature of our work is that the portfolio is allowed to anticipate the future, i.e. the terminal values of the prices, or of the driving Brownian motion, are known to the investor, either exactly or with some uncertainty. Results on the finiteness of the value of the control problem are obtained in various setups, using techniques from the so-called enlargement of filtrations. When the value of the problem is finite, we compute it explicitly and exhibit an optimal portfolio in closed form.

A single machine is available to process a collection of jobs J, each of which evolves stochastically under processing. Jobs incur costs while awaiting the machine at a rate which is state dependent and processing must respect a set of precedence constraints Γ. Index policies are optimal in a variety of scenarios. The indices concerned are characterised as values of restart problems with the average reward criterion. This characterisation yields a range of efficient approaches to their computation. Index-based suboptimality bounds are derived for general processing policies. These bounds enable us to develop sensitivity analyses and to evaluate scheduling heuristics.

We derive optimal gambling and investment policies for cases in which the underlying stochastic process has parameter values that are unobserved random variables. For the objective of maximizing logarithmic utility when the underlying stochastic process is a simple random walk in a random environment, we show that a state-dependent control is optimal, which is a generalization of the celebrated Kelly strategy: the optimal strategy is to bet a fraction of current wealth equal to a linear function of the posterior mean increment. To approximate more general stochastic processes, we consider a continuous-time analog involving Brownian motion. To analyze the continuous-time problem, we study the diffusion limit of random walks in a random environment. We prove that they converge weakly to a Kiefer process, or tied-down Brownian sheet. We then find conditions under which the discrete-time process converges to a diffusion, and analyze the resulting process. We analyze in detail the case of the natural conjugate prior, where the success probability has a beta distribution, and show that the resulting limit diffusion can be viewed as a rescaled Brownian motion. These results allow explicit computation of the optimal control policies for the continuous-time gambling and investment problems without resorting to continuous-time stochastic-control procedures. Moreover they also allow an explicit quantitative evaluation of the financial value of randomness, the financial gain of perfect information and the financial cost of learning in the Bayesian problem.

In this paper, we present results for the steady-state system length distributions of the discrete-time GI/G/1 queue. We examine the system at customer arrival epochs (customer departure epochs) and use the residual service time (residual interarrival time) as the supplementary variable. The embedded Markov chain is of GI/M/1 type if the embedding points are arrival epochs and is of M/G/1 type if the embedding points are departure epochs. Using the matrix analytic method, we identify the necessary and sufficient condition for both Markov chains to be positive recurrent. For the GI/M/1 type chain, we derive a matrix-geometric solution for its steady-state distribution and for the M/G/1 type chain, we develop a simple linear transformation that relates it to the GI/M/1 type chain and leads to a simple analytic solution for its steady-state distribution. We also show that the steady-state system length distribution at an arbitrary point in time can be obtained by a simple linear transformation of the matrix-geometric solution for the GI/M/1 type chain. A number of applications of the model to communication systems and numerical examples are also discussed.

This paper studies optimal routing and jockeying policies in a two-station parallel queueing system. It is assumed that jobs arrive to the system in a Poisson stream with rate λand are routed to one of two parallel stations. Each station has a single server and a buffer of infinite capacity. The service times are exponential with server-dependent rates, μ1 and μ2. Jockeying between stations is permitted. The jockeying cost is cij when a job in station i jockeys to station j, i ≠ j. There is no cost when a new job joins either station. The holding cost in station j is hj, h1 ≦ h2, per job per unit time. We characterize the structure of the dynamic routing and jockeying policies that minimize the expected total (holding plus jockeying) cost, for both discounted and long-run average cost criteria. We show that the optimal routing and jockeying controls are described by three monotonically non-decreasing functions. We study the properties of these control functions, their relationships, and their asymptotic behavior. We show that some well-known queueing control models, such as optimal routing to symmetric and asymmetric queues, preemptive or non-preemptive scheduling on homogeneous or heterogeneous servers, are special cases of our system.