Let (Y, Z) denote the solution to a forward-backward stochastic differential equation (FBSDE). If one constructs a random walk $B^n$ from the underlying Brownian motion B by Skorokhod embedding, one can show $L_2$-convergence of the corresponding solutions $(Y^n,Z^n)$ to $(Y, Z).$ We estimate the rate of convergence based on smoothness properties, especially for a terminal condition function in $C^{2,\alpha}$. The proof relies on an approximative representation of $Z^n$ and uses the concept of discretized Malliavin calculus. Moreover, we use growth and smoothness properties of the partial differential equation associated to the FBSDE, as well as of the finite difference equations associated to the approximating stochastic equations. We derive these properties by probabilistic methods.

We continue the study of the asymptotic behaviour of a random walk when it exits from a symmetric region of the form {(x, n): |x| ≤ rnb} as r → ∞ which was begun in Part I of this work. In contrast to that paper, we are interested in the case where the probability of exiting at the upper boundary tends to 1. In this scenario we treat the case where the power b lies in the interval [0, 1), and we establish necessary and sufficient conditions for the overshoot to be relatively stable in probability (except for the case ), and for the pth moment of the overshoot to be O(rq) as r → ∞.

We propose and study a number of layer methods for Navier‒Stokes equations (NSEs) with spatial periodic boundary conditions. The methods are constructed using probabilistic representations of solutions to NSEs and exploiting ideas of the weak sense numerical integration of stochastic differential equations. Despite their probabilistic nature, the layer methods are nevertheless deterministic.

We consider de Finetti’s problem for spectrally one-sided Lévy risk models with control strategies that are absolutely continuous with respect to the Lebesgue measure. Furthermore, we consider the version with a constraint on the time of ruin. To characterize the solution to the aforementioned models, we first solve the optimal dividend problem with a terminal value at ruin and show the optimality of threshold strategies. Next, we introduce the dual Lagrangian problem and show that the complementary slackness conditions are satisfied, characterizing the optimal Lagrange multiplier. Finally, we illustrate our findings with a series of numerical examples.

The reliability of many stochastic repairable systems depends on several characteristics that are time dependent. In this paper, we develop general repair models for a repairable system by using auxiliary stochastic processes which describe physical characteristics of the system and derive various properties of resulting models. We also obtain an inference procedure to assess the number of failures and expected number of failures for our proposed models by observing auxiliary processes.

Clearly, our inference procedure is different from the traditional approach, where successive times of occurrence of failures are observed.

A Fourier series is obtained for the variance of the number of points lying in a parallelogram thrown randomly on to a square lattice at a fixed angle.

This paper deals with the thrown string problem posed by Synge [12]. The models studied are those of Willenborg [14] and Kingman [7]. The convergence results obtained are complementary to those of Kingman, although they are derived by quite different methods involving the use of martingales.

Recently, systematic sampling on the circle and the sphere has been studied by Gual-Arnau and Cruz-Orive (2000) from a design-based point of view. In this note, it is shown that their mathematical model for the covariogram is, in a model-based statistical setting, a special case of the p-order shape model suggested by Hobolth, Pedersen and Jensen (2000) and Hobolth, Kent and Dryden (2002) for planar objects without landmarks. Benefits of this observation include an alternative variance estimator, applicable in the original problem of systematic sampling. In a wider perspective, the paper contributes to the discussion concerning design-based versus model-based stereology.

The hyperbolic random geometric graph was introduced by Krioukov et al. (Phys. Rev. E82, 2010). Among many equivalent models for the hyperbolic space, we study the d-dimensional Poincaré ball ( $d\ge 2$ ), with a general connectivity radius. While many phase transitions are known for the expectation asymptotics of certain subgraph counts, very little is known about the second-order results. Two of the distinguishing characteristics of geometric graphs on the hyperbolic space are the presence of tree-like hierarchical structures and the power-law behaviour of the degree distribution. We aim to reveal such characteristics in detail by investigating the behaviour of sub-tree counts. We show multiple phase transitions for expectation and variance in the resulting hyperbolic geometric graph. In particular, the expectation and variance of the sub-tree counts exhibit an intricate dependence on the degree sequence of the tree under consideration. Additionally, unlike the thermodynamic regime of the Euclidean random geometric graph, the expectation and variance may exhibit different growth rates, which is indicative of power-law behaviour. Finally, we also prove a normal approximation for sub-tree counts using the Malliavin–Stein method of Last et al. (Prob. Theory Relat. Fields165, 2016), along with the Palm calculus for Poisson point processes.

The ‘scanning process' Z(t), t ∈ ℝk of the title is a Gaussian random field obtained by associating with Z(t) the value of a set-indexed Brownian motion on the translate t + A0 of some ‘scanning set' A0. We study the basic properties of the random field Z relating, for example, its continuity and other sample path properties to the geometrical properties of A0. We ask if the set A0 determines the scanning process, and investigate when, and how, it is possible to recover the structure of A0 from realisations of the sample paths of the random field Z.

Generalizing the classical Banach matchbox problem, we consider the process of removing two types of ‘items’ from a ‘pile’ with selection probabilities for the type of the next item to be removed depending on the current numbers of remaining items, and thus changing sequentially. Under various conditions on the probability pn1,n2 that the next removal will take away an item of type I, given that n1 and n2 are the current numbers of items of the two types, we derive asymptotic formulas (as the initial pile size tends to infinity) for the probability that the items of type I are completely removed first and for the number of items left. In some special cases we also obtain explicit results.

This paper is the third in a series with the general title ‘How to look at the five-dimensional shape space Σ4,3. Parts I and II were concerned with distributions and diffusions respectively. Here we ‘look at’ geodesics.

We consider a model in which players must divide a fixed amount of resource between a number of trials of an underlying contest, with the underlying contest based on the War of Attrition. We are able to find the unique ES set (a simple generalisation of the idea of an evolutionarily stable strategy) in certain circumstances: in particular we find the conditions under which the ES set may contain a pure strategy.

We study the ruin probability where the claim sizes are modeled by a stationary ergodic symmetric α-stable process. We exploit the flow representation of such processes, and we consider the processes generated by conservative flows. We focus on two classes of conservative α-stable processes (one discrete-time and one continuous-time), and give results for the order of magnitude of the ruin probability as the initial capital goes to infinity. We also prove a solidarity property for null-recurrent Markov chains as an auxiliary result, which might be of independent interest.

We consider a new embedded Markov chain for the PH/G/1 queue by recording the queue length, the phase of the arrival process and the number of services completed during the current busy period at the successive departure epochs. Algorithmically tractable matrix formulas are obtained which permit the analysis of the fluctuations of the queue length and waiting times during a typical busy cycle. These are useful in the computation of certain profile curves arising in the statistical analysis of queues. In addition, informative expressions for the mean waiting times in the stable GI/G/1 queue and a simple new algorithm to evaluate the waiting-time distributions for the stationary PH/PH/1 queue are obtained.

Management of a bufferless resource is considered under non-homogeneous demand consisting of one-unit and two-unit requests. Two-unit requests can be served only by a given partition of the resource. Three simple admission policies are evaluated with regard to revenue generation. One policy involves no admission control and two policies involve trunk reservation. A limiting regime in which demand and capacity increase in proportion is considered. It is shown that each policy is asymptotically optimal for a certain range of parameters. Limiting dynamical behavior is obtained via a theory developed by Hunt and Kurtz. The results also point out the remarkable effect of partition constraints.

In this paper we introduce and solve a generalization of the classic average cost Brownian control problem in which a system manager dynamically controls the drift rate of a diffusion process X. At each instant, the system manager chooses the drift rate from a pair {u, v} of available rates and can invoke instantaneous controls either to keep X from falling or to keep it from rising. The objective is to minimize the long-run average cost consisting of holding or delay costs, processing costs, costs for invoking instantaneous controls, and fixed costs for changing the drift rate. We provide necessary and sufficient conditions on the cost parameters to ensure the problem admits a finite optimal solution. When it does, a simple control band policy specifying economic buffer sizes (α, Ω) and up to two switching points is optimal. The controller should invoke instantaneous controls to keep X in the interval (α, Ω). A policy with no switching points relies on a single drift rate exclusively. When there is no cost to change the drift rate, a policy with a single switching point s indicates that the controller should change to the slower drift rate when X exceeds s and use the faster drift rate otherwise. When there is a cost to change the drift rate, a policy with two switching points s < S indicates that the controller should maintain the faster drift rate until X exceeds S and maintain the slower drift rate until X falls below s.

In this paper, we study a large multi-server loss model under the SQ(d) routeing scheme when the service time distributions are general with finite mean. Previous works have addressed the exponential service time case when the number of servers goes to infinity, giving rise to a mean field model. The fixed point of the limiting mean field equations (MFEs) was seen to be insensitive to the service time distribution in simulations, but no proof was available. While insensitivity is well known for loss systems, the models, even with state-dependent inputs, belong to the class of linear Markov models. In the context of SQ(d) routeing, the resulting model belongs to the class of nonlinear Markov processes (processes whose generator itself depends on the distribution) for which traditional arguments do not directly apply. Showing insensitivity to the general service time distributions has thus remained an open problem. Obtaining the MFEs in this case poses a challenge due to the resulting Markov description of the system being in positive orthant as opposed to a finite chain in the exponential case. In this paper, we first obtain the MFEs and then show that the MFEs have a unique fixed point that coincides with the fixed point in the exponential case, thus establishing insensitivity. The approach is via a measure-valued Markov process representation and the martingale problem to establish the mean field limit.