The intensity of a stationary process of k-dimensional affine subspaces (k-flats) of ℝd with directional distribution from a given family R is estimated by observing the process in a compact window. To this end we introduce a type of unbiased estimator (the R-estimator) using the available information about the directional distribution.

Special cases are estimators for the intensity of stationary k-flat processes (1) with known directional distribution, (2) with directional distribution invariant with respect to a subgroup of the group of rotations in ℝd and (3) with unknown directional distribution.

We give sufficient conditions for the R-estimator to be the uniformly best unbiased estimator for the intensity of stationary Poisson k-flat processes with directional distribution in R. Equivalent statements for certain types of stationary Cox flat processes can be deduced directly from the results in the Poisson case.

Moreover, we consider stationary ergodic flat processes with directional distribution in R and general stationary flat processes with unknown directional distribution, all with a non-degeneracy property. In both cases our estimator turns out to be the uniformly best unbiased estimator from a restricted set of estimators. The result for general stationary flat processes is proved with the help of a factorization result for the second factorial moment measure.

We define nearest-neighbour point processes on graphs with Euclidean edges and linear networks. They can be seen as analogues of renewal processes on the real line. We show that the Delaunay neighbourhood relation on a tree satisfies the Baddeley‒Møller consistency conditions and provide a characterisation of Markov functions with respect to this relation. We show that a modified relation defined in terms of the local geometry of the graph satisfies the consistency conditions for all graphs with Euclidean edges that do not contain triangles.

We consider a threshold epidemic model on a clustered random graph model obtained from local transformations in an alternating branching process that approximates a bipartite graph. In other words, our epidemic model is such that an individual becomes infected as soon as the proportion of his/her infected neighbors exceeds the threshold q of the epidemic. In our random graph model, each individual can belong to several communities. The distributions for the community sizes and the number of communities an individual belongs to are arbitrary. We consider the case where the epidemic starts from a single individual, and we prove a phase transition (when the parameter q of the model varies) for the appearance of a cascade, i.e. when the epidemic can be propagated to an infinite part of the population. More precisely, we show that our epidemic is entirely described by a multi-type (and alternating) branching process, and then we apply Sevastyanov's theorem about the phase transition of multi-type Galton-Watson branching processes. In addition, we compute the entries of the mean progeny matrix corresponding to the epidemic. The phase transition for the contagion is given in terms of the largest eigenvalue of this matrix.

It is proved that a distribution F is strongly unimodal iff any two quantiles of the convolution of F with any other distribution are further apart than the corresponding quantiles of F itself. Some related characterizations of strong unimodality are also given.

In a branching process with variable lifetime, introduced by Fildes (1972) define Yjk(t) as the number of particles alive in generation k at time t when the initial particle is born in generation j. A limit law similar to that derived in the Bellman-Harris process is proved where it is shown that Yjk(t) suitably normalised converges in mean square to a random variable which is the limit random variable of Znm–n in the Galton-Watson process (m is the mean number of particles born).

We study the problem of maximizing the long-run average growth of total wealth for a logarithmic utility function under the existence of fixed and proportional transaction costs. The market model consists of one riskless asset and d risky assets. Impulsive control theory is applied to this problem. We derive a quasivariational inequality (QVI) of ‘ergodic’ type and obtain a weak solution for the inequality. Using this solution, we obtain an optimal investment strategy to achieve the optimal growth.

In this paper, in work strongly related with that of Coffman et al. [5], Bruss and Robertson [2], and Rhee and Talagrand [15], we focus our interest on an asymptotic distributional comparison between numbers of ‘smallest’ i.i.d. random variables selected by either on-line or off-line policies. Let X1,X2,… be a sequence of i.i.d. random variables with distribution function F(x), and let X1,n,…,Xn,n be the sequence of order statistics of X1,…,Xn. For a sequence (cn)n≥1 of positive constants, the smallest fit off-line counting random variable is defined by Ne(cn) := max {j ≤ n : X1,n + … + Xj,n ≤ cn}. The asymptotic joint distributional comparison is given between the off-line count Ne(cn) and on-line counts Nnτ for ‘good’ sequential (on-line) policies τ satisfying the sum constraint ∑j≥1XτjI(τj≤n) ≤ cn. Specifically, for such policies τ, under appropriate conditions on the distribution function F(x) and the constants (cn)n≥1, we find sequences of positive constants (Bn)n≥1, (Δn)n≥1 and (Δ'n)n≥1 such that

for some non-degenerate random variables W and W'. The major tools used in the paper are convergence of point processes to Poisson random measure and continuous mapping theorems, strong approximation results of the normalized empirical process by Brownian bridges, and some renewal theory.

We consider the problem of estimation in a partially observed discrete-time Galton-Watson branching process, focusing on the first two moments of the offspring distribution. Our study is motivated by modelling the counts of new cases at the onset of a stochastic epidemic, allowing for the facts that only a part of the cases is detected, and that the detection mechanism may affect the evolution of the epidemic. In this setting, the offspring mean is closely related to the spreading potential of the disease, while the second moment is connected to the variability of the mean estimators. Inference for branching processes is known for its nonstandard characteristics, as compared with classical inference. When, in addition, the true process cannot be directly observed, the problem of inference suffers significant further perturbations. We propose nonparametric estimators related to those used when the underlying process is fully observed, but suitably modified to take into account the intricate dependence structure induced by the partial observation and the interaction scheme. We show consistency, derive the limiting laws of the estimators, and construct asymptotic confidence intervals, all valid conditionally on the explosion set.

We consider a branching random walk on a multitype (with Q types of particles), supercritical Galton–Watson tree which satisfies the Kesten–Stigum condition. We assume that the displacements associated with the particles of type Q have regularly varying tails of index $\alpha$, while the other types of particles have lighter tails than the particles of type Q. In this paper we derive the weak limit of the sequence of point processes associated with the positions of the particles in the nth generation. We verify that the limiting point process is a randomly scaled scale-decorated Poisson point process using the tools developed by Bhattacharya, Hazra, and Roy (2018). As a consequence, we obtain the asymptotic distribution of the position of the rightmost particle in the nth generation.

In the literature on active redundancy allocation, the redundancy lifetimes are usually postulated to be independent of the component lifetimes for the sake of technical convenience. However, this unrealistic assumption leads to a risk of inaccurately evaluating system reliability, because it overlooks the statistical dependence of lifetimes due to common stresses. In this study, for k-out-of-n:F systems with component and redundancy lifetimes linked by the Archimedean copula, we show that (i) allocating more homogeneous redundancies to the less reliable components tends to produce a redundant system with stochastically larger lifetime, (ii) the reliability of the redundant system can be uniformly maximized through balancing the allocation of homogeneous redundancies in the context of homogeneous components, and (iii) allocating a more reliable matched redundancy to a less reliable component produces a more reliable system. These novel results on k-out-of-n:F systems in which component and redundancy lifetimes are statistically dependent are more applicable to the complicated engineering systems that arise in real practice. Some numerical examples are also presented to illustrate these findings.

A necessary condition for a rational Laplace–Stieltjes transform to correspond to a matrix exponential distribution is that the pole of maximal real part is real and negative. Given a rational Laplace–Stieltjes transform with such a pole, we present a method to determine whether or not the numerator polynomial admits a transform that corresponds to a matrix exponential distribution. The method relies on the minimization of a continuous function of one variable over the nonnegative real numbers. Using this approach, we give an alternative characterization for all matrix exponential distributions of order three.

We consider a model for a time series of spatial locations, in which points are placed sequentially at random into an initially empty region of ℝd, and given the current configuration of points, the likelihood at location x for the next particle is proportional to a specified function βk of the current number (k) of points within a specified distance of x. We show that the maximum likelihood estimator of the parameters βk (assumed to be zero for k exceeding some fixed threshold) is consistent in the thermodynamic limit where the number of points grows in proportion to the size of the region.

We consider the operation of an insurer with a large initial surplus x>0. The insurer's surplus process S(t) (with S(0)=x) evolves in the range S(t)≥ 0 as a generalized renewal process with positive mean drift and with jumps at time epochs T1,T2,…. At the time Tη(x) when the process S(t) first becomes negative, the insurer's ruin (in the ‘classical’ sense) occurs, but the insurer can borrow money via a line of credit. After this moment the process S(t) behaves as a solution to a certain stochastic differential equation which, in general, depends on the indebtedness, -S(t). This behavior of S(t) lasts until the time θ(x,y) at which the indebtedness reaches some ‘critical’ level y>0. At this moment the line of credit will be closed and the insurer's absolute ruin occurs with deficit -S(θ(x,y)). We find the asymptotics of the absolute ruin probability and the limiting distributions of η(x), θ(x,y), and -S(θ(x,y)) as x → ∞, assuming that the claim distribution is regularly varying. The second-order approximation to the absolute ruin probability is also obtained. The abovementioned results are obtained by using limiting theorems for the joint distribution of η(x) and -S(Tη(x)).

Conditions for a correlated cumulative shock model under which the system failure time is HNBUE are given. It is shown that the proof of a theorem given by Sumita and Shanthikumar (1985) relative to this property is not correct and a correct proof of the theorem is given.

We study the asymptotics of a Markovian system of N ≥ 3 particles in [0, 1]d in which, at each step in discrete time, the particle farthest from the current centre of mass is removed and replaced by an independent U[0, 1]d random particle. We show that the limiting configuration contains N − 1 coincident particles at a random location ξN ∈ [0, 1]d. A key tool in the analysis is a Lyapunov function based on the squared radius of gyration (sum of squared distances) of the points. For d = 1, we give additional results on the distribution of the limit ξN, showing, among other things, that it gives positive probability to any nonempty interval subset of [0, 1], and giving a reasonably explicit description in the smallest nontrivial case, N = 3.

We illustrate a technique for computing certain integrals that arise in probability theory by giving a new derivation of a formula of Kingman. This formula contains the joint distribution of the processes F(t) = inf {s: X(t + s) = b} and B(t) = inf{s: X(t - s) = b} where X is a time homogeneous, continuous parameter, Markov process and b is a fixed point in its state space. We then extend this formula to the situation in which b is replaced by a finite set {b1, …, bn}.

The steepest increase property of phase-type (PH) distributions was first proposed in O’Cinneide (1999) and proved in O’Cinneide (1999) and Yao (2002), but since then has received little attention in the research community. In this work we demonstrate that the steepest increase property can be applied for proving previously unknown moment bounds of PH distributions with infinite or finite support. Of special interest are moment bounds free of specific PH representations except the size of the representation. For PH distributions with infinite support, it is shown that such a PH distribution is stochastically smaller than or equal to an Erlang distribution of the same size. For PH distributions with finite support, a class of distributions which was introduced and investigated in Ramaswami and Viswanath (2014), it is shown that the squared coefficient of variation of a PH distribution with finite support is greater than or equal to 1/(m(m + 2)), where m is the size of its PH representation.

In the high-level operations of computer vision it is taken for granted that image features have been reliably detected. This paper addresses the problem of feature extraction by scale-space methods. There has been a strong development in scale-space theory and its applications to low-level vision in the last couple of years. Scale-space theory for continuous signals is on a firm theoretical basis. However, discrete scale-space theory is known to be quite tricky, particularly for low levels of scale-space smoothing. The paper is based on two key ideas: to investigate the stochastic properties of scale-space representations and to investigate the interplay between discrete and continuous images. These investigations are then used to predict the stochastic properties of sub-pixel feature detectors.

The modeling of image acquisition, image interpolation and scale-space smoothing is discussed, with particular emphasis on the influence of random errors and the interplay between the discrete and continuous representations. In doing so, new results are given on the stochastic properties of discrete and continuous random fields. A new discrete scale-space theory is also developed. In practice this approach differs little from the traditional approach at coarser scales, but the new formulation is better suited for the stochastic analysis of sub-pixel feature detectors.

The interpolated images can then be analysed independently of the position and spacing of the underlying discretisation grid. This leads to simpler analysis of sub-pixel feature detectors. The analysis is illustrated for edge detection and correlation. The stochastic model is validated both by simulations and by the analysis of real images.

For several common parent laws, the number of vertices of a sample convex hull follows a kind of law of large numbers. We exhibit an example of a parent law which contradicts a general conjecture about this matter.