We consider the distributional equation Z =^D ∑_k=1^NA_kZ(k), where N is a random variable taking value in N₀ = {0, 1, . . .}, A₁, A₂, . . . are p x p nonnegative random matrices, and Z, Z(1), Z(2), . . ., are independent and identically distributed random vectors in R₊^p with R₊ = [0, ∞), which are independent of (N, A₁, A₂, . . .). Let {Y_n} be the multidimensional Mandelbrot martingale defined as sums of products of random matrices indexed by nodes of a Galton-Watson tree plus an appropriate vector. Its limit Y is a solution of the equation above. For α > 1, we show a sufficient condition for E|Y|^α ∈ (0, ∞). Then for a nondegenerate solution Z of the distributional equation above, we show the decay rates of Ee^-t∙Z as |t| → ∞ and those of the tail probability P(y ∙ Z ≤ x) as x → 0 for given y = (y¹, . . ., y^p) ∈ R₊^p, and the existence of the harmonic moments of y ∙ Z. As an application, these results concerning the moments (of positive and negative orders) of Y are applied to a special multitype branching random walk. Moreover, for the case where all the vectors and matrices of the equation above are complex, a sufficient condition for the L^α convergence and the αth-moment of the Mandelbrot martingale {Y_n} are also established.

A reflected Ornstein-Uhlenbeck process is a process that returns continuously and immediately to the interior of the state space when it attains a certain boundary. It is an extended model of the traditional Ornstein-Uhlenbeck process being extensively used in finance as a one-factor short-term interest rate model. Under some mild conditions, this paper is devoted to the study of the analogue of the Cramer-Rao lower bound of a general class of parameter estimation of the unknown parameter in reflected Ornstein-Uhlenbeck processes.

The purpose of this paper is twofold. On the one hand, we provide sufficient conditions for the excess wealth order. These conditions are based on properties of the quantile functions which are useful when the dispersive order does not hold. On the other hand, we study sufficient conditions for the comparison in the increasing convex order of spacings of generalized order statistics. These results will be combined to show how we can provide comparisons of quantities of interest in reliability and insurance.

In this paper we build the increasing convex (concave) order for the scalar product of random vectors with an upper (lower) tail permutation decreasing joint density. As applications, we revisit allocations of portfolio risks in financial engineering and of coverage limits and deductibles in insurance. Some related results in the literature are substantially updated.

Consider a random directed graph on n vertices with independent and identically distributed outdegrees with distribution F having mean μ, and destinations of arcs selected uniformly at random. We show that if μ > 1 then for large n there is very likely to be a unique giant strong component with proportionate size given as the product of two branching process survival probabilities, one with offspring distribution F and the other with Poisson offspring distribution with mean μ. If μ ≤ 1 there is very likely to be no giant strong component. We also extend this to allow for F varying with n.

We investigate default probabilities and default correlations of Merton-type credit portfolio models in stress scenarios where a common risk factor is truncated. For elliptically distributed asset variables, the asymptotic limits of default probabilities and default correlations depend on the max-domain of attraction of the asset variables. In the regularly varying case, we derive an integral representation for multivariate default probabilities, which turn out to be strictly smaller than 1. Default correlations are in (0, 1). In the rapidly varying case, asymptotic multivariate default probabilities are 1 and asymptotic default correlations are 0.

We analyze the dynamics of a stochastic process driven by binomial random variables, where the probability of success depends on the past realization. We study the limit behavior when the group size is fixed but the number of iterations increases. It will become apparent that the so-called policy function and its fixed point play an outstanding role. Some applications to a statistical analysis of gender bias are also briefly discussed.

We provide an algorithm to find the value and an optimal strategy of the Ten Thousand dice game solitaire variant in the framework of Markov control processes. Once an optimal critical threshold is found, the set of nonstopping states of the game becomes finite and the solution is found by a backwards algorithm that gives the values for each one of these states of the game. The algorithm is finite and exact. The strategy to find the critical threshold comes from the continuous pasting condition used in optimal stopping problems for continuous-time processes with jumps.

This paper provides a comparison of the successive lumping (SL) methodology developed in Katehakis et al. (2015) with the popular lattice path counting (Mohanty (1979)) in obtaining rate matrices for queueing models, satisfying the specific quasi birth and death structure as in Van Leeuwaarden et al. (2009) and Van Leeuwaarden and Winands (2006). The two methodologies are compared both in terms of applicability requirements and numerical complexity by analyzing their performance for the same classical queueing models considered in Van Leeuwaarden et al. (2009). The main findings are threefold. First, when both methods are applicable, the SL-based algorithms outperform the lattice path counting algorithm (LPCA). Second, there are important classes of problems (for example, models with (level) nonhomogenous rates or with finite state spaces) for which the SL methodology is applicable and for which the LPCA cannot be used. Third, another main advantage of SL algorithms over lattice path counting is that the former includes a method to compute the steady state distribution using this rate matrix.

One often observed empirical regularity is a power-law behavior of the tails of some distribution of interest. We propose a limit law for normalized random means that exhibits such heavy tails irrespective of the distribution of the underlying sampling units: the limit is a t-distribution if the random variables have finite variances. The generative scheme is then extended to encompass classic limit theorems for random sums. The resulting unifying framework has wide empirical applicability which we illustrate by considering two empirical regularities in two different fields. First, we turn to urban geography and explain why city-size growth rates are approximately t-distributed, using a model of random sector growth based on the central place theory. Second, turning to an issue in finance, we show that high-frequency stock index returns can be modeled as a generalized asymmetric Laplace process. These empirical illustrations elucidate the situations in which heavy tails can arise.

By considering the extreme behavior of bivariate random vectors with a polar representation R(u(T), v(T)), it is commonly assumed that the radial component R and the angular component T are stochastically independent. We investigate how to relax this rigid independence assumption such that conditional limit theorems can still be deduced. For this purpose, we introduce a novel measure for the dependence structure and present convenient criteria for validity of limit theorems possessing a geometrical meaning. Thus, our results verify a stability of the available limit results, which is essential in applications where the independence of the polar components is not necessarily present or exactly fulfilled.

For the directed edge preferential attachment network growth model studied by Bollobás et al. (2003) and Krapivsky and Redner (2001), we prove that the joint distribution of in-degree and out-degree has jointly regularly varying tails. Typically, the marginal tails of the in-degree distribution and the out-degree distribution have different regular variation indices and so the joint regular variation is nonstandard. Only marginal regular variation has been previously established for this distribution in the cases where the marginal tail indices are different.

We consider the drawing with and without replacement models from numerical populations. For order statistics based on the above drawing schemes, we provide sharp upper bounds for their variances, expressed in the single observation variance units. We also characterize the populations for which the bounds are attained.

A cascading failure is a failure in a system of interconnected parts in which the failure of a part can trigger the failure of successive parts. Although an initial and introductory approach for probabilistic modeling and analysis of the cascading failures was suggested in the literature, any general framework and fundamental results have yet to be reported. In this paper, applying the point process approach, we suggest a general framework for modeling and analysis of the cascading failures. Furthermore, a new concept of 'information-based residual lifetime' will be defined and discussed.

In this paper we derive sufficient conditions for the permanence and ergodicity of a stochastic predator-prey model with a Beddington-DeAngelis functional response. The conditions obtained are in fact very close to the necessary conditions. Both nondegenerate and degenerate diffusions are considered. One of the distinctive features of our results is that they enable the characterization of the support of a unique invariant probability measure. It proves the convergence in total variation norm of the transition probability to the invariant measure. Comparisons to the existing literature and matters related to other stochastic predator-prey models are also given.

We study continuous-time birth-death type processes, where individuals have independent and identically distributed lifetimes, according to a random variable Q, with E[Q] = 1, and where the birth rate if the population is currently in state (has size) n is α(n). We focus on two important examples, namely α(n) = λ n being a branching process, and α(n) = λn(N - n) / N which corresponds to an SIS (susceptible → infective → susceptible) epidemic model in a homogeneously mixing community of fixed size N. The processes are assumed to start with a single individual, i.e. in state 1. Let T, A_n, C, and S denote the (random) time to extinction, the total time spent in state n, the total number of individuals ever alive, and the sum of the lifetimes of all individuals in the birth-death process, respectively. We give expressions for the expectation of all these quantities and show that these expectations are insensitive to the distribution of Q. We also derive an asymptotic expression for the expected time to extinction of the SIS epidemic, but now starting at the endemic state, which is not independent of the distribution of Q. The results are also applied to the household SIS epidemic, showing that, in contrast to the household SIR (susceptible → infective → recovered) epidemic, its threshold parameter R_* is insensitive to the distribution of Q.

Starting from a Markov chain with a finite or a countable infinite alphabet, we consider the chain obtained when all but one symbol are indistinguishable for the practitioner. We study conditions on the transition matrix of the Markov chain ensuring that the image chain has continuous or discontinuous transition probabilities with respect to the past.

We consider a class of discrete-time Markov chains with state space [0, 1] and the following dynamics. At each time step, first the direction of the next transition is chosen at random with probability depending on the current location. Then the length of the jump is chosen independently as a random proportion of the distance to the respective end point of the unit interval, the distributions of the proportions being fixed for each of the two directions. Chains of that kind were the subjects of a number of studies and are of interest for some applications. Under simple broad conditions, we establish the ergodicity of such Markov chains and then derive closed-form expressions for the stationary densities of the chains when the proportions are beta distributed with the first parameter equal to 1. Examples demonstrating the range of stationary distributions for processes described by this model are given, and an application to a robot coverage algorithm is discussed.

We consider a continuous, infinitely divisible random field in R^d given as an integral of a kernel function with respect to a Lévy basis with convolution equivalent Lévy measure. For a large class of such random fields we compute the asymptotic probability that the supremum of the field exceeds the level x as x → ∞. Our main result is that the asymptotic probability is equivalent to the right tail of the underlying Lévy measure.

We describe an Euler scheme to approximate solutions of Lévy driven stochastic differential equations (SDEs) where the grid points are given by the arrival times of a Poisson process and thus are random. This result extends the previous work of Ferreiro-Castilla et al. (2014). We provide a complete numerical analysis of the algorithm to approximate the terminal value of the SDE and prove that the mean-square error converges with rate O(n^-1/2). The only requirement of the methodology is to have exact samples from the resolvent of the Lévy process driving the SDE. Classical examples, such as stable processes, subclasses of spectrally one-sided Lévy processes, and new families, such as meromorphic Lévy processes (Kuznetsov et al. (2012), are examples for which our algorithm provides an interesting alternative to existing methods, due to its straightforward implementation and its robustness with respect to the jump structure of the driving Lévy process.

We establish the moderate deviation principle for the solutions of a class of stochastic partial differential equations with non-Lipschitz continuous coefficients. As an application, we derive the moderate deviation principle for two important population models: super-Brownian motion and the Fleming-Viot process.

We present a very short derivation of the integral representation of the two-sided Skorokhod reflection Z of a continuous function X of bounded variation, which is a generalization of the integral representation of the one-sided map featured in Anantharam and Konstantopoulos (2011) and Konstantopoulos et al. (1996). We also show that Z satisfies a simpler integral representation when additional conditions are imposed on X.

A computable necessary and sufficient condition of separation cutoff is obtained for a sequence of continuous-time upward skip-free chains with the stochastically monotone time-reversals.

Let {Z_n, n = 0, 1, 2, . . .} be a supercritical branching process, {N_t, t ≥ 0} be a Poisson process independent of {Z_n, n = 0, 1, 2, . . .}, then {Z_Nt, t ≥ 0} is a supercritical Poisson random indexed branching process. We show a law of large numbers, central limit theorem, and large and moderate deviation principles for log Z_Nt.