This paper introduces a convenient class of spatiotemporal random field models that can be interpreted as the mean-square solutions of stochastic fractional evolution equations.

Random vectors in the positive orthant whose distributions possess hidden regular variation are a subclass of those whose distributions are multivariate regularly varying with asymptotic independence. The concept is an elaboration of the coefficient of tail dependence of Ledford and Tawn. We show that the rank transform that brings unequal marginals to the standard case also preserves the hidden regular variation. We discuss applications of the results to two examples, one involving flood risk and the other Internet data.

We consider a multi-type branching random walk on d-dimensional Euclidian space. The~uniform convergence, as n goes to infinity, of a scaled version of the Laplace transform of the point process given by the nth generation particles of each type is obtained. Similar results in the one-type case, where the transform gives a martingale, were obtained in Biggins (1992) and Barral (2001). This uniform convergence of transforms is then used to obtain limit results for numbers in the underlying point processes. Supporting results, which are of interest in their own right, are obtained on (i) ‘Perron-Frobenius theory’ for matrices that are smooth functions of a variable λ∈L and are nonnegative when λ∈L−⊂L, where L is an open set in ℂd, and (ii) saddlepoint approximations of multivariate distributions. The saddlepoint approximations developed are strong enough to give a refined large deviation theorem of Chaganty and Sethuraman (1993) as a by-product.

We consider risk processes that locally behave like Brownian motion with some drift and variance, these both depending on an underlying Markov chain that is also used to generate the claims arrival process. Thus, claims arrive according to a renewal process with waiting times of phase type. Positive claims (downward jumps) are always possible but negative claims (upward jumps) are also allowed. The claims are assumed to form an independent, identically distributed sequence, independent of everything else. As main results, the joint Laplace transform of the time to ruin and the undershoot at ruin, as well as the probability of ruin, are explicitly determined under the assumption that the Laplace transform of the positive claims is a rational function. Both the joint Laplace transform and the ruin probability are decomposed according to the type of ruin: ruin by jump or ruin by continuity. The methods used involve finding certain martingales by first finding partial eigenfunctions for the generator of the Markov process composed of the risk process and the underlying Markov chain. We also use certain results from complex function theory as important tools.

In this paper, we study the extremal behavior of stationary mixed moving average processes of the form Y(t)=∫ℝ+×ℝf(r,t-s) dΛ(r,s), t∈ℝ, where f is a deterministic function and Λ is an infinitely divisible, independently scattered random measure whose underlying driving Lévy process is regularly varying. We give sufficient conditions for the stationarity of Y and compute the tail behavior of certain functionals of Y. The extremal behavior is modeled by marked point processes on a discrete-time skeleton chosen properly by the jump times of the underlying driving Lévy process and the extremes of the kernel function. The sequences of marked point processes converge weakly to a cluster Poisson random measure and reflect extremes of Y at a high level. We also show convergence of the partial maxima to the Fréchet distribution. Our models and results cover short- and long-range dependence regimes.

Variograms and covariance functions are the fundamental tools for modeling dependent data observed over time, space, or space-time. This paper aims at constructing nonseparable spatio-temporal variograms and covariance models. Special attention is paid to an intrinsically stationary spatio-temporal random field whose covariance function is of Schoenberg-Lévy type. The correlation structure is studied for its increment process and for its partial derivative with respect to the time lag, as well as for the superposition over time of a stationary spatio-temporal random field. As another approach, we investigate the permissibility of the linear combination of certain separable spatio-temporal covariance functions to be a valid covariance, and obtain a subclass of stationary spatio-temporal models isotropic in space.

Assuming that the forward rates ftu are semimartingales, we give conditions on their components under which the discounted bond prices are martingales. To achieve this, we give sufficient conditions for the integrated processes ftu=∫0uftvdv to be semimartingales, and identify their various components. We recover the no-arbitrage conditions in models well known in the literature and, finally, we formulate a new random field model for interest rates and give its equivalent martingale measure (no-arbitrage) condition.

We prove the existence and uniqueness of the solution to certain reflected backward stochastic differential equations (RBSDEs) with one continuous barrier and deterministic terminal time, under monotonicity, and general increasing growth conditions on the associated coefficient. As an application, we obtain, in some constraint cases, the price of an American contingent claim as the unique solution of such an RBSDE.

Diffusion models in economics, finance, queueing, mathematical biology, and electrical engineering often involve reflecting barriers. In this paper, we study the analytical representation of transition densities for reflected one-dimensional diffusions in terms of their associated Sturm-Liouville spectral expansions. In particular, we provide explicit analytical expressions for transition densities of Brownian motion with drift, the Ornstein-Uhlenbeck process, and affine (square-root) diffusion with one or two reflecting barriers. The results are easily implementable on a personal computer and should prove useful in applications.

Let (Xt)t≥0 be a continuous-time irreducible Markov chain on a finite state space E, let v be a map v: E→ℝ\{0}, and let (φt)t≥0 be an additive functional defined by φt=∫0tv(Xs)d s. We consider the case in which the process (φt)t≥0 is oscillating and that in which (φt)t≥0 has a negative drift. In each of these cases, we condition the process (Xt,φt)t≥0 on the event that (φt)t≥0 is nonnegative until time T and prove weak convergence of the conditioned process as T→∞.

In this paper, we analyze the diffusion limit of a discrete-time queueing system with constant service rate and connections that randomly enter and depart from the system. Each connection generates periodic traffic while it is active, and a connection's lifetime has finite mean. This can model a time division multiple access system with constant bit-rate connections. The diffusion scaling retains semiperiodic behavior in the limit, allowing for both short-time analysis (within one frame) and long-time analysis (over multiple frames). Weak convergence of the cumulative arrival process and the stationary buffer-length distribution is proved. It is shown that the limit of the cumulative arrival process can be viewed as a discrete-time stationary-increment Gaussian process interpolated by Brownian bridges. We present bounds on the overflow probability of the limit queueing process as functions of the arrival rate and the connection lifetime distribution. Also, numerical and simulation results are presented for geometrically distributed connection lifetimes.

Let X1, X2,… be real-valued random variables. For u>0, define the time of ruin T = T(u) by T = inf{n: X1+⋯+Xn>u} or T=∞ if X1+⋯+Xn≤u for every n = 1,2,…. We are interested in the ruin probabilities of general processes {Xn} for large u. In the presence of heavy tails, one often finds power estimates. Our objective is to specify the associated powers and provide the crude estimate P(T≤xu)≈u−R(x) for large u, for a given x∈ℝ. The rate R(x) will be described by means of tails of partial sums and maxima of {Xn}. We also extend our results to the case of the infinite time horizon.

Let (Xt)t≥0 be a continuous-time irreducible Markov chain on a finite state space E, let v: E→ℝ\{0}, and let (φt)t≥0 be defined by φt=∫0tv(Xs)d s. We consider the case in which the process (φt)t≥0 is oscillating and that in which (φt)t≥0 has a negative drift. In each of these cases, we condition the process (Xt,φt)t≥0 on the event that (φt)t≥0 hits level y before hitting 0 and prove weak convergence of the conditioned process as y→∞. In addition, we show the relationship between the conditioning of the process (φt)t≥0 with a negative drift to oscillate and the conditioning of it to stay nonnegative for a long time, and the relationship between the conditioning of (φt)t≥0 with a negative drift to drift to ∞ and the conditioning of it to hit large levels before hitting 0.

Consider a continuous-time Markov chain evolving in a random environment. We study certain forms of interaction between the process of interest and the environmental process, under which the stationary joint distribution is tractable. Moreover, we obtain necessary and sufficient conditions for a product-form stationary distribution. A number of examples that illustrate the applicability of our results in queueing and population growth models are also included.

We consider a service system (QS) that operates according to the first-come-first-served (FCFS) discipline, and in which the service rate is an increasing function of the queue length. Customers arrive sequentially at the system, and decide whether or not to join using decision rules based upon the queue length on arrival. Each customer is interested in selecting a rule that meets a certain optimality criterion with regard to their expected sojourn time in the system; as a consequence, the decision rules of other customers must be taken into account. Within a particular class of decision rules for an associated infinite-player game, the structure of the Nash equilibrium routeing policies is characterized. We prove that, within this class, there exist a finite number of Nash equilibria, and that at least one of these is nonrandomized. Finally, with the aid of simulation experiments, we explore the extent to which the Nash equilibria are characteristic of customer joining behaviour under a learning rule based on system-wide data.

We study stationary processes given as solutions to stochastic differential equations driven by fractional Brownian motion. This rich class includes the fractional Ornstein-Uhlenbeck process and those processes that can be obtained from it by state space transformations. An explicit formula in terms of Euler's Γ-function describes the asymptotic behaviour of the covariance function of the fractional Ornstein-Uhlenbeck process near zero, which, by an application of Berman's condition, guarantees that this process is in the maximum domain of attraction of the Gumbel distribution. Necessary and sufficient conditions on the state space transforms are stated to classify the maximum domain of attraction of solutions to stochastic differential equations driven by fractional Brownian motion.

Let X1,X2,… be a sequence of independent and identically distributed random variables with some continuous distribution function F. Let L(n) and X(n) denote the nth record time and the nth record value, respectively. We refer to the variables Xi as near-nth-record observations if Xi∈(X(n)-a,X(n)], with a>0, and L(n)<i<L(n+1). In this work we study asymptotic properties of the number of near-record observations. We also discuss sums of near-record observations.

In this paper, we study the stationary distributions for reflected diffusions with jumps in the positive orthant. Under the assumption that the stationary distribution possesses a density in R+n that satisfies certain finiteness conditions, we characterize the Fokker-Planck equation. We then provide necessary and sufficient conditions for the existence of a product-form distribution for diffusions with oblique boundary reflections and jumps. To do so, we exploit a recent characterization of the boundary properties of such reflected processes. In particular, we show that the conditions generalize those for semimartingale reflecting Brownian motions and reflected Lévy processes. We provide explicit results for some models of interest.

We derive necessary and sufficient conditions for the existence of bounded or summable solutions to systems of linear equations associated with Markov chains. This substantially extends a famous result of G. E. H. Reuter, which provides a convenient means of checking various uniqueness criteria for birth-death processes. Our result allows chains with much more general transition structures to be accommodated. One application is to give a new proof of an important result of M. F. Chen concerning upwardly skip-free processes. We then use our generalization of Reuter's lemma to prove new results for downwardly skip-free chains, such as the Markov branching process and several of its many generalizations. This permits us to establish uniqueness criteria for several models, including the general birth, death, and catastrophe process, extended branching processes, and asymptotic birth-death processes, the latter being neither upwardly skip-free nor downwardly skip-free.

In this paper, we provide a novel approach to studying the heavy-tailed asymptotics of the stationary probability vector of a Markov chain of GI/G/1 type, whose transition matrix is constructed from two matrix sequences referred to as a boundary matrix sequence and a repeating matrix sequence, respectively. We first provide a necessary and sufficient condition under which the stationary probability vector is heavy tailed. Then we derive the long-tailed asymptotics of the R-measure in terms of the RG-factorization of the repeating matrix sequence, and a Wiener-Hopf equation for the boundary matrix sequence. Based on this, we are able to provide a detailed analysis of the subexponential asymptotics of the stationary probability vector.