In this paper we propose a stochastic process, which is a Cox-Ingersoll-Ross process with Hawkes jumps. It can be seen as a generalization of the classical Cox-Ingersoll-Ross process and the classical Hawkes process with exponential exciting function. Our model is a special case of the affine point processes. We obtain Laplace transforms and limit theorems, including the law of large numbers, central limit theorems, and large deviations.

We derive in this paper closed formulae for the joint probability generating function of the number of customers in the two FIFO queues of a generalized processor-sharing (GPS) system with two classes of customers arriving according to Poisson processes and requiring exponential service times. In contrast to previous studies published on the GPS system, we show that it is possible to establish explicit expressions for the generating functions of the number of customers in each queue without calling for the formulation of a Riemann–Hilbert problem. We specifically prove that the problem of determining the unknown functions due to the reflecting conditions on the boundaries of the positive quarter plane can be reduced to a Poisson equation. The explicit formulae are then used to derive some characteristics of the GPS system (in particular the tails of the probability distributions of the numbers of customers in each queue).

The ergodic properties of birth and death processes are studied. We obtain some explicit estimates for the rate of convergence by the methods of theory of differential equations.

Consider the random uniform placement of a finite number of arcs on the circle, where the arc lengths are sampled from a distribution on (0, 1). We provide exact formulae for the probability that the circle is completely covered and for the distribution of the number of uncovered gaps, extending Stevens's (1939) formulae for the case of fixed equal arc lengths. A special class of arc length distributions is considered, and exact probabilities of coverage are tabulated for the uniform distribution on (0, 1). Some asymptotic results for the number of gaps are also given.

We present a new inspection policy useful when testing is needed to detect failures of a single-unit system. It is supposed that tests may fail and give an erroneous result. The inspection policy minimizing cost per unit of time for an infinite time span is also discussed. In addition, we study the behaviour of the optimum policy for some time to failure distributions often assumed in reliability: exponential and Pareto.

The partially observed control problem is considered for stochastic processes with control entering into the diffusion and the observation. The maximum principle is proved for the partially observable optimal control. A pure probabilistic approach is used, and the adjoint processes are characterized as solutions of related backward stochastic differential equations in finite-dimensional spaces. Most of the derivation is identified with that of the completely observable case.

An autoregressive process ARP(1) with Pareto-distributed inputs, analogous to those of Lawrance and Lewis (1977), (1980), is defined and its properties developed. It is shown that the stationary distributions are Pareto. Further, the maximum and minimum processes are asymptotically Weibull, and the ARP(1) process is shown to be closed under maximization or minimization when the number of terms is geometrically distributed. The ARP(1) process leads naturally to an extremal process in the sense of Lamperti (1964). Statistical inference for the ARP(1) process is developed. An absolutely continuous variant of the Pareto process is described.

We consider initially two parallel routes, each of two queues in tandem, with arriving customers choosing the route giving them the shortest expected time in the system, given the queue lengths at the customer's time of arrival. All interarrival and service times are exponential.

We then augment this network to obtain a Wheatstone bridge, in which customers may cross from one route to the other between queues, again choosing the route giving the shortest expected time in the system, given the queue lengths ahead of them.

We find that Braess's paradox can occur: namely in equilibrium the expected transit time in the augmented network, for some service rates, can be greater than in the initial network.

Let X1, X2, ·· ·be stationary normal random variables with ρn = cov(X0, Xn). The asymptotic joint distribution of and is derived under the condition ρn log n → γ [0,∞). It is seen that the two statistics are asymptotically independent only if γ = 0.

This paper studies inspected systems with non-self-announcing failures where the rate of deterioration is governed by a Markov chain. We compute the lifetime distribution and availability when the system is inspected according to a periodic inspection policy. In doing so, we expose the role of certain transient distributions of the environment.

Discrete renewal theory is generalized to study the occurrence of a collection of patterns in random sequences, where a renewal is defined to be the occurrence of one of the patterns in the collection which does not overlap an earlier renewal. The action of restriction enzymes on DNA sequences provided motivation for this work. Related results of Guibas and Odlyzko are discussed.

1. In 1923 Eggenberger and Pólya introduced the following ‘urn scheme’ as a model for the development of a contagious phenomenon. A box contains b black and r red balls, and a ball is drawn from it at random with ‘double replacement’ (i.e. whatever ball is drawn, it is returned to the box together with a fresh ball of the same colour); the procedure is then continued indefinitely. A slightly more complicated version with m-fold replacement is sometimes discussed, but it will be sufficient for our purposes to keep m = 2 and it will be convenient further to simplify the scheme by taking b = r = 1 as the initial condition. We shall however generalise the scheme in another direction by allowing an arbitrary number k(≧2) of colours. Thus initially the box will contain k differently coloured balls and successive random drawings will be followed by double replacement as before. We write sn (a k-vector with jth component ) for the numerical composition of the box immediately after the nth replacement, so that and we observe that is a Markov process for which the state-space consists of all ordered k-ads of positive integers, the (constant) transition-probability matrix having elements determined by where Sn is the sum of the components of sn and (e(i))j = δij. We shall calculate the Martin boundary for this Markov process, and point out some applications to stochastic models for population growth.

The similarity between order statistics and record values motivated us, in this paper, to investigate some relationships between the order statistics and the record values. These relationships are employed to characterize a continuous distribution by some moment properties of the spacings of the record values, and hence obtain characterizations of the exponential distribution. Some of the well-known results follow trivially.

A general multivariate non-linear regression model is considered, including as special cases linear regression when the regression matrix is of less than full rank, simultaneous equations systems and regression on an unobservable predetermined variable. Given a time-series of observations at unit intervals we consider the estimation of the parameters, subject to non-linear constraints, by minimizing a criterion based on the Fourier-transformed model. We allow the residuals to be generated by a stationary, linear, process and establish asymptotic properties of our estimates.

In this paper we present a robust numerical method to compute the scale function W(q)(x) of a general spectrally negative Lévy process (X, P). The method is based on the Esscher transform of measure Pν under which X is taken and the scale function is determined. This change of measure makes it possible for the scale function to be bounded and, hence, makes numerical computation easy, fast, and stable. Working under the new measure Pν and using the method of Abate and Whitt (1992) and Choudhury, Lucantoni, and Whitt (1994), we give a fast stable numerical algorithm for the computation of W(q)(x).

Block replacement policy, under which items are replaced at constant intervals of time and at failure, is wasteful because sometimes almost new items are also replaced. In order to eliminate this kind of waste, the following replacement policy is suggested.

(1) Items are exchanged for new items at time kT (k ··· 1, 2, …).

(2) If items fail in [(k – 1) T, kT – v), they are replaced by new items, and if in [kT – v, kT), they are replaced by used items.

The results of this ‘extended block replacement policy with used items' are studied for the Erlang lifetime distribution.

In this paper conditions for the convergence of a class of simulated annealing algorithms for continuous global optimization are given. The previous literature about the subject gives results for the convergence of algorithms in which the next candidate point is generated according to a probability distribution whose support is the whole feasible set. A class of possible cooling schedules has been introduced in order to remove this restriction.

We prove that, under rather general conditions, the law of a continuous Gaussian process represented by a stochastic integral of a deterministic kernel, with respect to a standard Wiener process, can be weakly approximated by the law of some processes constructed from a standard Poisson process. An example of a Gaussian process to which this result applies is the fractional Brownian motion with any Hurst parameter.

In the last few years there has been renewed interest in the classical control problem of de Finetti (1957) for the case where the underlying source of randomness is a spectrally negative Lévy process. In particular, a significant step forward was made by Loeffen (2008), who showed that a natural and very general condition on the underlying Lévy process which allows one to proceed with the analysis of the associated Hamilton-Jacobi-Bellman equation is that its Lévy measure is absolutely continuous, having completely monotone density. In this paper we consider de Finetti's control problem, but with the restriction that control strategies are absolutely continuous with respect to the Lebesgue measure. This problem has been considered by Asmussen and Taksar (1997), Jeanblanc-Picqué and Shiryaev (1995), and Boguslavskaya (2006) in the diffusive case, and Gerber and Shiu (2006) for the case of a Cramér-Lundberg process with exponentially distributed jumps. We show the robustness of the condition that the underlying Lévy measure has a completely monotone density and establish an explicit optimal strategy for this case that envelopes the aforementioned existing results. The explicit optimal strategy in question is the so-called refraction strategy.

It is possible to interpret the classical central limit theorem for sums of independent random variables as a convergence rate result for the law of large numbers. For example, if Xi, i = 1, 2, 3, ··· are independent and identically distributed random variables with EXi = μ, var Xi = σ2 < ∞ and then the central limit theorem can be written in the form This provides information on the rate of convergence in the strong law as . (“a.s.” denotes almost sure convergence.) It is our object in this paper to discuss analogues for the super-critical Galton-Watson process.