Interest has been shown in Markovian sequences of geometric shapes. Mostly the equations for invariant probability measures over shape space are extremely complicated and multidimensional. This paper deals with rectangles which have a simple one-dimensional shape descriptor. We explore the invariant distributions of shape under a variety of randomised rules for splitting the rectangle into two sub-rectangles, with numerous methods for selecting the next shape in sequence. Many explicit results emerge. These help to fill a vacant niche in shape theory, whilst contributing at the same time, new distributions on [0,1] and interesting examples of Markov processes or, in the language of another discipline, of stochastic dynamical systems.

We consider d-dimensional stochastic processes which take values in (R+)d These processes generalize Galton-Watson branching processes, but the main assumption of branching processes, independence between particles, is dropped. Instead, we assume for some

Here τ:(R+)d→R +, |x| = σ1d |x(i)|, A {x ∈(R+)d: |x| 1} and T: A→A. Under various assumptions on the maps τ and T it is shown that with probability one there exists a ρ > 1, a fixed point p ∈ A of T and a random variable w such that limn→∞Zn|ρnwp. This result is a generalization of the main limit theorem for supercritical branching processes; note, however, that in the present situation both ρ and ρ are random as well. The results are applied to a population genetical model for zygotic selection without mutation at one locus.

Various types of time and event sampling of a stationary and orderly bivariate point process are considered. Fundamental relations between inter-event intervals and the event counting process are derived. Relations between first forward recurrence times and their moments for different types of sampling are obtained.

Generalized perturbation analysis (PA) estimates to study sensitivity of performance measures of discrete events dynamic systems for discontinuous sample functions are introduced. Their convergence conditions and rate of convergence are given. It is shown that the PA estimates based on a single sample path always converge faster to the unknown sensitivity parameter (vector of parameters) than their counterpart—crude Monte Carlo ones.

The purpose of this note is to point out the connection between the invariance property of M/M/1 and GI/G/1 queues (which has been reported in several papers) and the interchangeability and reversibility properties of tandem queues. This enables us to gain new insights for both problems and obtain stronger invariance results for M/M/1, GI/G/1, as well as loss systems M/M/1/N, GI/G/1/N and tandem systems.

We consider the problem of numerical integration when the sampling nodes form a stationary point process on the real line. In previous papers it was argued that a naïve Riemann sum approach can cause a severe variance inflation when the sampling points are not equidistant. We show that this inflation can be avoided using a higher-order Newton–Cotes quadrature rule which exploits smoothness properties of the integrand. Under mild assumptions, the resulting estimator is unbiased and its variance asymptotically obeys a power law as a function of the mean point distance. If the Newton–Cotes rule is of sufficiently high order, the exponent of this law turns out to only depend on the point process through its mean point distance. We illustrate our findings with the stereological estimation of the volume of a compact object, suggesting alternatives to the well-established Cavalieri estimator.

We consider the statistical estimation problem of potential functions of Gibbs models on the plane lattice. We assume that the area on which a random point pattern is observed is sufficiently large and take an asymptotic point of view. The main result is to show the locally asymptotic normality of the Gibbs model under certain conditions. From this result we can show the optimality of the maximum likelihood estimator employing known results about locally asymptotic normal families, though a practical computation of the maximum likelihood estimator presents difficulties. An estimation procedure based on the method of moments is also proposed.

We consider a stochastic control model driven by a fractional Brownian motion. This model is a formal approximation to a queueing network with an ON-OFF input process. We study stochastic control problems associated with the long-run average cost, the infinite-horizon discounted cost, and the finite-horizon cost. In addition, we find a solution to a constrained minimization problem as an application of our solution to the long-run average cost problem. We also establish Abelian limit relationships among the value functions of the above control problems.

We consider a stochastic control model for a queueing system driven by a two-dimensional fractional Brownian motion with Hurst parameter 0 < H < 1. In particular, when H > ½, this model serves to approximate a controlled two-station tandem queueing model with heavy-tailed ON/OFF sources in heavy traffic. We establish the weak convergence results for the distribution of the state process and construct an explicit stationary state process associated with given controls. Based on suitable coupling arguments, we show that each state process couples with its stationary counterpart and we use it to represent the long-run average cost functional in terms of the stationary process. Finally, we establish the existence result of an optimal control, which turns out to be independent of the initial data.

We identify a close relation between stable distributions and the limiting homomorphisms central to the theory of regular variation. In so doing some simplifications are achieved in the direct analysis of these laws in Pitman and Pitman (2016); stable distributions are themselves linked to homomorphy.

In this paper we study a discrete-time optimal switching problem on a finite horizon. The underlying model has a running reward, terminal reward, and signed (positive and negative) switching costs. Using optimal stopping theory for discrete-parameter stochastic processes, we extend a well-known explicit dynamic programming method for computing the value function and the optimal strategy to the case of signed switching costs.

In this paper we introduce the multivariate Brownian semistationary (BSS) process and study the joint asymptotic behaviour of its realised covariation using in-fill asymptotics. First, we present a central limit theorem for general multivariate Gaussian processes with stationary increments, which are not necessarily semimartingales. Then, we show weak laws of large numbers, central limit theorems, and feasible results for BSS processes. An explicit example based on the so-called gamma kernels is also provided.

A general multi-type population model is considered, where individuals live and reproduce according to their age and type, but also under the influence of the size and composition of the entire population. We describe the dynamics of the population as a measure-valued process and obtain its asymptotics as the population grows with the environmental carrying capacity. Thus, a deterministic approximation is given, in the form of a law of large numbers, as well as a central limit theorem. This general framework is then adapted to model sexual reproduction, with a special section on serial monogamic mating systems.

In this paper we study the treewidth of the random geometric graph, obtained by dropping n points onto the square [0,√n]2 and connecting pairs of points by an edge if their distance is at most r=r(n). We prove a conjecture of Mitsche and Perarnau (2014) stating that, with probability going to 1 as n→∞, the treewidth of the random geometric graph is 𝜣(r√n) when lim inf r>rc, where rc is the critical radius for the appearance of the giant component. The proof makes use of a comparison to standard bond percolation and with a little bit of extra work we are also able to show that, with probability tending to 1 as k→∞, the treewidth of the graph we obtain by retaining each edge of the k×k grid with probability p is 𝜣(k) if p>½ and 𝜣(√log k) if p<½.

We derive the large-sample distribution of the number of species in a version of Kingman’s Poisson–Dirichlet model constructed from an $\alpha$ -stable subordinator but with an underlying negative binomial process instead of a Poisson process. Thus it depends on parameters $\alpha\in (0,1)$ from the subordinator and $r>0$ from the negative binomial process. The large-sample distribution of the number of species is derived as sample size $n\to\infty$ . An important component in the derivation is the introduction of a two-parameter version of the Dickman distribution, generalising the existing one-parameter version. Our analysis adds to the range of Poisson–Dirichlet-related distributions available for modeling purposes.

For a class of cell division processes in the Euclidean space ℝd, spatial consistency is investigated. This addresses the problem whether the distribution of the generated structures, restricted to a bounded set V, depends on the choice of a larger region W ⊃ V where the construction of the cell division process is performed. This can also be understood as the problem of boundary effects in the cell division procedure. It is known that the STIT tessellations are spatially consistent. In the present paper it is shown that, within a reasonable wide class of cell division processes, the STIT tessellations are the only ones that are consistent.

We consider regular variation of a Lévy process X := (Xt)t≥0 in with Lévy measure Π, emphasizing the dependence between jumps of its components. By transforming the one-dimensional marginal Lévy measures to those of a standard 1-stable Lévy process, we decouple the marginal Lévy measures from the dependence structure. The dependence between the jumps is modeled by a so-called Pareto Lévy measure, which is a natural standardization in the context of regular variation. We characterize multivariate regularly variation of X by its one-dimensional marginal Lévy measures and the Pareto Lévy measure. Moreover, we define upper and lower tail dependence coefficients for the Lévy measure, which also apply to the multivariate distributions of the process. Finally, we present graphical tools to visualize the dependence structure in terms of the spectral density and the tail integral for homogeneous and nonhomogeneous Pareto Lévy measures.