The aim of this paper is to extend the existing theory of second-order self-similar processes as defined by Cox (1984) from the univariate case to higher dimensions. Multivariate self-similar processes defined in terms of second-order theory for stationary time series can be used as models for long-range dependent observations when the marginal observations are long-range dependent. An interesting question concerns the correlation structure within the processes when the marginal processes are correlated. We show that the self-similarity requirement, as defined in this article, implies a cross-correlation structure similar to that for the marginal processes. This occurs both in the time domain and in the frequency domain. This fact can be used to obtain generalized least squares estimates for the long-range dependence parameters. We discuss some difficulties concerning estimation based on simulations.

We derive two kinds of rate conservation laws for describing the time-dependent behavior of a process defined with a stationary marked point process and starting at time 0. These formulas are called TRCLs (time-dependent rate conservation laws). It is shown that TRCLs are useful to study the transient behaviors of risk and storage processes with stationary claim and supply processes and with a general premium and release rates, respectively. Detailed discussions are given for the severity for the risk process, and for the workload process of a single-server queue.

Through the study of a simple embedded martingale we obtain an extension of the Kesten–Stigum theorem and prove a central limit theorem for controlled Galton-Watson processes.

We study the problem of maximizing the probability of stopping at an object which is best in at least one of a given set of criteria, using only stopping rules based on the knowledge of whether the current object is relatively best in each of the criteria. The asymptotic results for the case of independent criteria are shown to hold in certain cases where the componentwise maxima are, pairwise, either asymptotically independent or asymptotically full dependent.

An example of the former is a random sample from a bivariate correlated normal distribution; thus our results settle a question posed recently by T. S. Ferguson.

Let X, X1, X2, … be i.i.d. Sn =Σ1 n Xj , E|X| > 0, E(X) = 0 and τ = inf {n ≥ 1: Sn ≥ 0}. By Wald's equation, E(τ) =∞. If E(X 2) <∞, then by a theorem of Burkholder and Gundy (1970), E(τ1/2) =∞. In this paper, we prove that if E((X– )2 ) <∞, then E(τ1/2) =∞. When X is integer-valued and X ≥ −1 a.s., a necessary and sufficient condition for E(τ1–1/p ) <∞, p > 1, is Σn–1–1p E|Sn| <∞.

In this article, we generalize results by Dimitrov and Khalil (1990), Khalil et al. (1991), and van Harn and Steutel (1991) and obtain some characterizations of the exponential and geometric laws.

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.

This paper considers the absorption of a non-decreasing compound Poisson process of finite order in a general upper boundary. The problem is relevant in fields such as risk theory, Kolmogorov–Smirnov statistics and sequential analysis. The probability of absorption and first-passage times are given in terms of a generating function which depends on the boundary only and can be computed readily. Absorption is certain or not as the asymptotic slope of the boundary is greater or less than the expected increase of the process in unit time. The case of the linear boundary is considered in detail.

If X is a Brownian motion with drift and γ = inf{t > 0: Mt = t} we derive the joint density of the triple {U, γ, Δ}, where and Δ= γ —Xγ . In the case δ = 0 it follows easily from this that Δ has an Exp(2) distribution and this in turn implies the rather surprising result that if τ= inf{t > 0: Xt = Mt = t}, then Pr{τ = 0} = 0 and . We also derive various other distributional results involving the pair (X, M), including for example the distribution of ; in particular we show that, in case δ. = 1, when Pr{0 < τ < ∞} = 1, the ratio τ+/τ has the arc-sine distribution.

Exceedances of a non-stationary sequence above a boundary define certain point processes, which converge in distribution under mild mixing conditions to Poisson processes. We investigate necessary and sufficient conditions for the convergence of the point process of exceedances, the point process of upcrossings and the point process of clusters of exceedances. Smooth regularity conditions, as smooth oscillation of the non-stationary sequence, imply that these point processes converge to the same Poisson process. Since exceedances are asymptotically rare, the results are extended to triangular arrays of rare events.

In this paper, optimal stopping problems for semi-Markov processes are studied in a fairly general setting. In such a process transitions are made from state to state in accordance with a Markov chain, but the amount of time spent in each state is random. The times spent in each state follow a general renewal process. They may depend on the present state as well as on the state into which the next transition is made.

Our goal is to maximize the expected net return, which is given as a function of the state at time t minus some cost function. Discounting may or may not be considered. The main theorems (Theorems 3.5 and 3.11) are expressions for the optimal stopping time in the undiscounted and discounted case. These theorems generalize results of Zuckerman [16] and Boshuizen and Gouweleeuw [3]. Applications are given in various special cases.

The results developed in this paper can also be applied to semi-Markov shock models, as considered in Taylor [13], Feldman [6] and Zuckerman [15].

We propose a two-parameter family of conjugate prior distributions for the number of undiscovered objects in a class of Bayesian search models. The family contains the one-parameter Euler and Heine families as special cases. The two parameters may be interpreted respectively as an overall success rate and a rate of depletion of the source of objects. The new family gives enhanced flexibility in modelling.

A reference probability is explicitly constructed under which the signal and observation processes are independent. A simple, explicit recursive form is then obtained for the conditional density of the signal given the observations. Both non-linear and linear filters are considered, as well as two different information patterns.

For two-dimensional spatial data, a spatial unilateral autoregressive moving average (ARMA) model of first order is defined and its properties studied. The spatial correlation properties for these models are explicitly obtained, as well as simple conditions for stationarity and conditional expectation (interpolation) properties of the model. The multiplicative or linear-by-linear first-order spatial models are seen to be a special case which have proved to be of practical use in modeling of two-dimensional spatial lattice data, and hence the more general models should prove to be useful in applications. These unilateral models possess a convenient computational form for the exact likelihood function, which gives proper treatment to the border cell values in the lattice that have a substantial effect in estimation of parameters. Some simulation results to examine properties of the maximum likelihood estimator and a numerical example to illustrate the methods are briefly presented.

The full-information secretary problem in which the objective is to minimize the expected rank is seen to have a value smaller than 7/3 for all n (the number of options). This can be achieved by a simple memoryless threshold rule. The asymptotically optimal value for the class of such rules is about 2.3266. For a large finite number of options, the optimal stopping rule depends on the whole sequence of observations and seems to be intractable. This raises the question whether the influence of the history of all observations may asymptotically fade. We have not solved this problem, but we show that the values for finite n are non-decreasing in n and exhibit a sequence of lower bounds that converges to the asymptotic value which is not smaller than 1.908.

The two-point Markov chain boundary-value problem discussed in this paper is a finite-time version of the quasi-stationary behaviour of Markov chains. Specifically, for a Markov chain {Xt :t = 0, 1, ·· ·}, given the time interval (0, n), the interest is in describing the chain at some intermediate time point r conditional on knowing both the behaviour of the chain at the initial time point 0 and that over the interval (0, n) it has avoided some subset B of the state space. The paper considers both ‘real time' estimates for r = n (i.e. the chain has avoided B since 0), and a posteriori estimates for r < n with at least partial knowledge of the behaviour of Xn. Algorithms to evaluate the distribution of Xr can be as small as O(n 3) (and, for practical purposes, even O(n 2 log n)). The estimates may be stochastically ordered, and the process (and hence, the estimates) may be spatially homogeneous in a certain sense. Maximum likelihood estimates of the sample path are furnished, but by example we note that these ML paths may differ markedly from the path consisting of the expected or average states. The scope for two-point boundary-value problems to have solutions in a Markovian setting is noted.

Several examples are given, together with a discussion and examples of the analogous problem in continuous time. These examples include the basic M/G/k queue and variants that include a finite waiting room, reneging, balking, and Bernoulli feedback, a pure birth process and the Yule process. The queueing examples include Larson's (1990) ‘queue inference engine'.

We propose an AR(1) model that can be used to generate logistic processes. The proposed model has simple probability and correlation structure that can accommodate the full range of attainable correlation. The correlation structure and the joint distribution of the proposed model are given, as well as their conditional mean and variance.

The first-order autoregressive semi-Mittag-Leffler (SMLAR(1)) process is introduced and its properties are studied. As an illustration, we discuss the special case of the first-order autoregressive Mittag-Leffler (MLAR(1)) process.

We define a class of two-dimensional Markov random graphs with I, V, T and Y-shaped nodes (vertices). These are termed polygonal models. The construction extends our earlier work [1]– [5]. Most of the paper is concerned with consistent polygonal models which are both stationary and isotropic and which admit an alternative description in terms of the trajectories in space and time of a one-dimensional particle system with motion, birth, death and branching. Examples of computer simulations based on this description are given.

There are a number of cases in the theories of queues and dams where the limiting distribution of the pertinent processes is geometric with a modified initial term — herein called zero-modified geometric (ZMG). The paper gives a unified treatment of the various cases considered hitherto and some others by using a duality relation between random walks with impenetrable and with absorbing barriers, and deriving the probabilities of absorption by using Waldian identities. Thus the method enables us to distinguish between those cases where the limiting distribution would be ZMG and those where it would not.