Search

For many practical situations in reliability engineering, components in the system are usually dependent since they generally work in a collaborative environment. In this paper we build sufficient conditions for comparing two coherent systems under different random environments in the sense of the usual stochastic, hazard rate, reversed hazard rate, and likelihood ratio orders. Applications and numerical examples are provided to illustrate all the theoretical results established here.

In this paper we study the Robbins–Monro procedure Xn+1 = Xn - an-1Yn with some fixed number a > 0 and establish the moderate deviation principle of the process {Xn}.

Community detection is considered for a stochastic block model graph of n vertices, with K vertices in the planted community, edge probability p for pairs of vertices both in the community, and edge probability q for other pairs of vertices. The main focus of the paper is on weak recovery of the community based on the graph G, with o(K) misclassified vertices on average, in the sublinear regime n1-o(1) ≤ K ≤ o(n). A critical parameter is the effective signal-to-noise ratio λ = K2(p - q)2 / ((n - K)q), with λ = 1 corresponding to the Kesten–Stigum threshold. We show that a belief propagation (BP) algorithm achieves weak recovery if λ > 1 / e, beyond the Kesten–Stigum threshold by a factor of 1 / e. The BP algorithm only needs to run for log*n + O(1) iterations, with the total time complexity O(|E|log*n), where log*n is the iterated logarithm of n. Conversely, if λ ≤ 1 / e, no local algorithm can asymptotically outperform trivial random guessing. Furthermore, a linear message-passing algorithm that corresponds to applying a power iteration to the nonbacktracking matrix of the graph is shown to attain weak recovery if and only if λ > 1. In addition, the BP algorithm can be combined with a linear-time voting procedure to achieve the information limit of exact recovery (correctly classify all vertices with high probability) for all K ≥ (n / logn) (ρBP + o(1)), where ρBP is a function of p / q.

Mao and Hu (2010) left an open problem about the hazard rate order between the largest order statistics from two samples of n geometric random variables. Du et al. (2012) solved this open problem when n = 2, and Wang (2015) solved for 2 ≤ n ≤ 9. In this paper we completely solve this problem for any value of n.

In this paper we prove a quenched functional central limit theorem for a biased random walk on a supercritical Galton–Watson tree with leaves. This extends a result of Peres and Zeitouni (2008) where the case without leaves was considered. A conjecture of Ben Arous and Fribergh (2016) suggests an upper bound on the bias which we observe to be sharp.

In this paper we propose a new approach to study the Parisian ruin problem for spectrally negative Lévy processes. Since our approach is based on a hybrid observation scheme switching between discrete and continuous observations, we call it a temporal approach as opposed to the spatial approximation approach in the literature. Our approach leads to a unified proof for the underlying processes with bounded or unbounded variation paths, and our result generalizes Loeffen et al. (2013).

In this paper we relate the set of structure-preserving equivalent martingale measures ℳsp for financial models driven by semimartingales with conditionally independent increments to a set of measurable and integrable functions 𝒴. More precisely, we prove that ℳsp ≠ ∅ if and only if 𝒴 ≠ ∅, and connect the sets ℳsp and 𝒴 to the semimartingale characteristics of the driving process. As examples we consider integrated Lévy models with independent stochastic factors and time-changed Lévy models and derive mild conditions for ℳsp ≠ ∅.

In this paper we prove that a parallel system consisting of Weibull components with different scale parameters ages faster than a parallel system comprising Weibull components with equal scale parameters in the convex transform order when the lifetimes of components of both systems have different shape parameters satisfying some restriction. Moreover, while comparing these two systems, we show that the dispersive and the usual stochastic orders, and the right-spread order and the increasing convex order are equivalent. Further, some of the known results in the literature concerning comparisons of k-out-of-n systems in the exponential model are extended to the Weibull model. We also provide solutions to two open problems mentioned by Balakrishnan and Zhao (2013) and Zhao et al. (2016).

While max-stable processes are typically written as pointwise maxima over an infinite number of stochastic processes, in this paper, we consider a family of representations based on ℓp-norms. This family includes both the construction of the Reich–Shaby model and the classical spectral representation by de Haan (1984) as special cases. As the representation of a max-stable process is not unique, we present formulae to switch between different equivalent representations. We further provide a necessary and sufficient condition for the existence of an ℓp-norm-based representation in terms of the stable tail dependence function of a max-stable process. Finally, we discuss several properties of the represented processes such as ergodicity or mixing.

A classical P��lya urn scheme is a Markov process whose evolution is encoded by a replacement matrix $(R_{i,j})_{1\leq i,j\leq d}$. At every discrete time-step, we draw a ball uniformly at random, denote its colour $c$, and replace it in the urn together with $R_{c,j}$ balls of colour $j$ (for all $1\leq j\leq d$). We are interested in multi-drawing P��lya urns, where the replacement rule depends on the random drawing of a set of $m$ balls from the urn (with or without replacement). This generalisation has already been studied in the literature, in particular by Kuba & Mahmoud (ArXiv:1503.09069 and 1509.09053), where second order asymptotic results are proved for $2$-colour urns under the balanced and the affinity assumptions. The main idea of this work is to apply stochastic approximation methods to this problem, which enables us to remove the affinity hypothesis of Kuba & Mahmoud and generalise the result to more-than-two-colour urns. We also give some partial results in the two-colour non-balanced case., This new arxiv version (v6) corrects a mistake that we discovered in the previous versions of this paper (v1-5). The mistake was in Theorem 1$(a)$ and in the last sentence of Theorem 4. In this new version, Theorem 1$(a)$ has been corrected, and Theorem 4 has been deleted

In this paper we consider the problem of simultaneously estimating rare-event probabilities for a class of Gaussian random fields. A conventional rare-event simulation method is usually tailored to a specific rare event and consequently would lose estimation efficiency for different events of interest, which often results in additional computational cost in such simultaneous estimation problems. To overcome this issue, we propose a uniformly efficient estimator for a general family of Hölder continuous Gaussian random fields. We establish the asymptotic and uniform efficiency of the proposed method and also conduct simulation studies to illustrate its effectiveness.

In this paper we obtain a recursive formula for the density of the double-barrier Parisian stopping time. We present a probabilistic proof of the formula for the first few steps of the recursion, and then a formal proof using explicit Laplace inversions. These results provide an efficient computational method for pricing double-barrier Parisian options.

In this paper we present a set of results relating to the occupation time α(t) of a processX(·). The first set of results concerns exact characterizations of α(t), e.g. in terms of its transform up to an exponentially distributed epoch. In addition, we establish a central limit theorem (entailing that a centered and normalized version of α(t)∕tconverges to a zero-mean normal random variable ast→∞) and the tail asymptotics of ℙ(α(t)∕t≥q). We apply our findings to spectrally positive Lévy processes reflected at the infimum and establish various new occupation time results for the corresponding model.

In this paper we investigate the Parisian ruin problem of the general dual Lévy risk model. Unlike the usual concept of ultimate ruin, allowing the surplus level to be negative within a prespecified period indicates that the deficit at Parisian ruin is not necessarily equal to zero. Hence, we consider a Gerber–Shiu type expected discounted penalty function at the Parisian ruin and obtain an explicit expression for this function under the dual Lévy risk model. As particular cases, we calculate the Parisian ruin probability and the expected discountedkth moments of the deficit at the Parisian ruin for the compound Poisson dual risk model and a drift-diffusion model. Numerical examples are given to illustrate the behavior of Parisian ruin and the expected discounted deficit at Parisian ruin.

In this paper we investigate the stochastic properties of the number of failed components of a three-state network. We consider a network made up of n components which is designed for a specific purpose according to the performance of its components. The network starts operating at time t = 0 and it is assumed that, at any time t > 0, it can be in one of states up, partial performance, or down. We further suppose that the state of the network is inspected at two time instants t1 and t2 (t1 < t2). Using the notion of the two-dimensional signature, the probability of the number of failed components of the network is calculated, at t1 and t2, under several scenarios about the states of the network. Stochastic and ageing properties of the proposed failure probabilities are studied under different conditions. We present some optimal age replacement policies to show applications of the proposed criteria. Several illustrative examples are also provided.

In this paper, by using the exact tail asymptotics derived by Debicki, Hashorva and Ji (Ann. Probab. 2014), we proved the Gumbel limit theorem for the maximum of a class of non-homogeneous Gaussian random fields. By using the obtained results, we also derived the Gumbel laws for Shepp statistics of fractional Brownian motion and Gaussian integrated process as well as the Gumbel law for Storage process with fractional Brownian motion as

In this paper we provide a sufficient condition for mean residual life ordering of parallel systems with n ≥ 3 heterogeneous exponential components.

In this paper we study general aggregation of stochastic arrangement increasing random variables, including both the generalized linear combination and the standard aggregation as special cases. In terms of monotonicity, supermodularity, and convexity of the kernel function, we develop several sufficient conditions for the increasing convex order on the generalized aggregations. Some applications in reliability and risks are also presented.

We study the rare-event behavior of the workload process in a transitory queue, where the arrival epochs (or 'points') of a finite number of jobs are assumed to be the ordered statistics of independent and identically distributed (i.i.d.) random variables. The service times (or 'marks') of the jobs are assumed to be i.i.d. random variables with a general distribution, that are jointly independent of the arrival epochs. Under the assumption that the service times are strictly positive, we derive the large deviations principle (LDP) satisfied by the workload process. The analysis leverages the connection between ordered statistics and self-normalized sums of exponential random variables to establish the LDP. In this paper we present the first analysis of rare events in transitory queueing models, supplementing prior work that has focused on fluid and diffusion approximations.

Current literature on stochastic dominance assumes utility/loss functions to be the same across random variables. However, decision models with inconsistent utility functions have been proposed in the literature. The use of inconsistent loss functions when comparing between two random variables can also be appropriate under other problem settings. In this paper we generalize almost stochastic dominance to problems with inconsistent utility/loss functions. In particular, we propose a set of conditions that is necessary and sufficient for clear preferences when the utility/loss functions are allowed to vary across different random variables.

In this paper we study the existence of Lebesgue densities of one-dimensional Lévy processes. Equivalently, we show the absolute continuity of the distributions of one-dimensional Lévy processes. Compared with the previous literature, we consider Lévy processes with Lévy symbols of a logarithmic behavior at ∞.

Batch service has a wide application in manufacturing, communication networks, and cloud computing. In batch service queues with limited resources, one critical issue is to properly schedule the service so as to ensure the quality of service. In this paper we consider an M/G[a,b]/1/N batch service queue with bulking threshold a, max service capacity b, and buffer capacity N, where N can be finite or infinite. Through renewal theory, busy period analysis and decomposition techniques, we demonstrate explicitly how the bulking threshold influences the system performance such as the mean waiting time and time-averaged number of loss customers in batch service queues. We then establish a necessary and sufficient condition on the optimal bulking threshold that minimizes the expected waiting time. Enabled by this condition, we propose a simple algorithm which guarantees to find the optimal threshold in polynomial time. The performance of the algorithm is also demonstrated by numerical examples.

In this paper we consider a diffusive stochastic predator–prey model with a nonlinear functional response and the randomness is assumed to be of Gaussian nature. A large deviation principle is established for solution processes of the considered model by implementing the weak convergence technique.

Drawdown (resp. drawup) of a stochastic process, also referred as the reflected process at its supremum (resp. infimum), has wide applications in many areas including financial risk management, actuarial mathematics and statistics. In this paper, for general time-homogeneous Markov processes, we study the joint law of the first passage time of the drawdown (resp. drawup) process, its overshoot, and the maximum of the underlying process at this first passage time. By using short-time pathwise analysis, under some mild regularity conditions, the joint law of the three drawdown quantities is shown to be the unique solution to an integral equation which is expressed in terms of fundamental two-sided exit quantities of the underlying process. Explicit forms for this joint law are found when the Markov process has only one-sided jumps or is a L\'{e}vy process (possibly with two-sided jumps). The proposed methodology provides a unified approach to study various drawdown quantities for the general class of time-homogeneous Markov processes., Comment: 24 pages, 2 figures

Importance sampling has become an important tool for the computation of extreme quantiles and tail-based risk measures. For estimation of such nonlinear functionals of the underlying distribution, the standard efficiency analysis is not necessarily applicable. In this paper we therefore study importance sampling algorithms by considering moderate deviations of the associated weighted empirical processes. Using a delta method for large deviations, combined with classical large deviation techniques, the moderate deviation principle is obtained for importance sampling estimators of two of the most common risk measures: value at risk and expected shortfall.

In this paper we deal with an optimal stopping problem whose objective is to maximize the probability of selecting k out of the last ℓ successes, given a sequence of independent Bernoulli trials of length N, where k and ℓ are predetermined integers satisfying 1≤k≤ℓ

In the paper we consider the following modification of a discrete-time branching process with stationary immigration. In each generation a binomially distributed subset of the population will be observed. The number of observed individuals constitute a partially observed branching process. After inspection both observed and unobserved individuals may change their offspring distributions. In the subcritical case we investigate the possibility of using the known estimators for the offspring mean and for the mean of the stationary-limiting distribution of the process when the observation of the population sizes is restricted. We prove that, if both the population and the number of immigrants are partially observed, the estimators are still strongly consistent. We also prove that the `skipped' version of the estimator for the offspring mean is asymptotically normal and the estimator of the stationary distribution's mean is asymptotically normal under additional assumptions.

This paper pioneers a Freidlin–Wentzell approach to stochastic impulse control of exchange rates when the central bank desires to maintain a target zone. Pressure to stimulate the economy forces the bank to implement diffusion monetary policy involving Freidlin–Wentzell perturbations indexed by a parameter ε∈ [0,1]. If ε=0, the policy keeps exchange rates in the target zone for all times t≥0. When ε>0, exchange rates continually exit the target zone almost surely, triggering central bank interventions which force currencies back into the zone or abandonment of all targets. Interventions and target zone deviations are costly, motivating the bank to minimize these joint costs for any ε∈ [0,1]. We prove convergence of the value functions as ε→0 achieving a value function approximation for small ε. Via sample path analysis and cost function bounds, intervention followed by target zone abandonment emerges as the optimal policy.

The contribution of the theory of scan statistics to the study of many real-life applications has been rapidly expanding during the last decades. The multiple scan statistic, defined on a sequence of n Bernoulli trials, enumerates the number of occurrences of k consecutive trials which contain at least r successes among them (r≤k≤n). In this paper we establish some asymptotic results for the distribution of the multiple scan statistic, as n,k,r→∞ and illustrate their accuracy through a simulation study. Our approach is based on an appropriate combination of compound Poisson approximation and random walk theory.

In this paper we aim to apply simple actuarial methods to build an insurance plan protecting against an epidemic risk in a population. The studied model is an extended SIR epidemic in which the removal and infection rates may depend on the number of registered removals. The costs due to the epidemic are measured through the expected epidemic size and infectivity time. The premiums received during the epidemic outbreak are measured through the expected susceptibility time. Using martingale arguments, a method by recursion is developed to calculate the cost components and the corresponding premium levels in this extended epidemic model. Some numerical examples illustrate the effect of removals and the premium calculation in an insurance plan.

In this paper we study the expected rank problem under full information. Our approach uses the planar Poisson approach from Gnedin (2007) to derive the expected rank of a stopping rule that is one of the simplest nontrivial examples combining rank dependent rules with threshold rules. This rule attains an expected rank lower than the best upper bounds obtained in the literature so far, in particular, we obtain an expected rank of 2.326 14.

The sequential stochastic assignment problem (SSAP) allocates distinct workers with deterministic values to sequentially arriving tasks with stochastic parameters to maximize the expected total reward. In this paper we study an extension of the SSAP, in which the worker values are considered to be random variables, taking on new values upon each task arrival. Several SSAP models with different assumptions on the distribution of the worker values and closed-form expressions for optimal assignment policies are presented.

Consider a d-type (dk≥2 individuals in the nth generation, pick k individuals at random from the nth generation and trace their lines of descent backward in time till they meet. In this paper, the limit behaviors of the distributions of the generation number of the most recent common ancestor of any k chosen individuals and of the whole population are studied for both critical and subcritical cases. Also, we investigate the limit distribution of the joint distribution of the generation number and their types.

In this paper we estimate quantile sensitivities for dependent sequences via infinitesimal perturbation analysis, and prove asymptotic unbiasedness, weak consistency, and a central limit theorem for the estimators under some mild conditions. Two common cases, the regenerative setting and ϕ-mixing, are analyzed further, and a new batched estimator is constructed based on regenerative cycles for regenerative processes. Two numerical examples, the G/G/1 queue and the Ornstein–Uhlenbeck process, are given to show the effectiveness of the estimator.

We consider a greedy walk on a Poisson process on the real line. It is known that the walk does not visit all points of the process. In this paper we first obtain some useful independence properties associated with this process which enable us to compute the distribution of the sequence of indices of visited points. Given that the walk tends to +∞, we find the distribution of the number of visited points in the negative half-line, as well as the distribution of the time at which the walk achieves its minimum.

De Bruijn's identity relates two important concepts in information theory: Fisher information and differential entropy. Unlike the common practice in the literature, in this paper we consider general additive non-Gaussian noise channels where more realistically, the input signal and additive noise are not independently distributed. It is shown that, for general dependent signal and noise, the first derivative of the differential entropy is directly related to the conditional mean estimate of the input. Then, by using Gaussian and Farlie–Gumbel–Morgenstern copulas, special versions of the result are given in the respective case of additive normally distributed noise. The previous result on independent Gaussian noise channels is included as a special case. Illustrative examples are also provided.

In this paper we study the asymptotic behaviour near extinction of (sub-)critical continuous-state branching processes. In particular, we establish an analogue of Khintchine's law of the iterated logarithm near extinction time for a continuous-state branching process whose branching mechanism satisfies a given condition.

The coalescent was introduced by Kingman (1982a), (1982b) and Tajima (1983) as a continuous-time Markov chain model describing the genealogical relationship among sampled genes from a panmictic population of a species. The random mating in a population is a strict condition and the genealogical structure of the population has a strong influence on the genetic variability and the evolution of the species. In this paper, starting from a discrete-time Markov chain model, we show the weak convergence to a continuous-time Markov chain, called the structured coalescent model, describing the genealogy of the sampled genes from whole population by means of passing the limit of the population size. Herbots (1997) proved the weak convergence to the structured coalescent on the condition of conservative migration and Wright–Fisher-type reproduction. We will give the proof on the condition of general migration rates and exchangeable reproduction.

We study in this paper a generalized coupon collector problem, which consists in analyzing the time needed to collect a given number of distinct coupons that are drawn from a set of coupons with an arbitrary probability distribution. We suppose that a special coupon called the null coupon can be drawn but never belongs to any collection. In this context, we prove that the almost uniform distribution, for which all the non-null coupons have the same drawing probability, is the distribution which stochastically minimizes the time needed to collect a fixed number of distinct coupons. Moreover, we show that in a given closed subset of probability distributions, the distribution with all its entries, but one, equal to the smallest possible value is the one, which stochastically maximizes the time needed to collect a fixed number of distinct coupons. An computer science application shows the utility of these results., Comment: arXiv admin note: text overlap with arXiv:1402.5245

In this paper we study a special class of size dependent branching processes. We assume that for some positive integer K as long as the population size does not exceed level K, the process evolves as a discrete-time supercritical branching process, and when the population size exceeds level K, it evolves as a subcritical or critical branching process. It is shown that this process does die out in finite time T. The question of when the mean value E(T) is finite or infinite is also addressed.

In this paper we connect Archimedean survival processes (ASPs) with the theory of Markov copulas. ASPs were introduced by Hoyle and Mengütürk (2013) to model the realized variance of two assets. We present some new properties of ASPs related to their dependency structure. We study weak and strong Markovian consistency properties of ASPs. An ASP is weak Markovian consistent, but generally not strong Markovian consistent. Our results contain necessary and sufficient conditions for an ASP to be strong Markovian consistent. These properties are closely related to the concept of Markov copulas, which is very useful in modelling different dependence phenomena. At the end we present possible applications.

In this paper, we investigate the functional central limit theorem for stochastic processes associated to partial sums of additive functionals of reversible Markov chains with general spate space, under the normalization standard deviation of partial sums. For this case, we show that the functional central limit theorem is equivalent to the fact that the variance of partial sums is regularly varying with exponent 1 and the partial sums satisfy the CLT. It is also equivalent to the conditional CLT., 10 pages

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.

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.

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.

Classical Jackson networks are a well-established tool for the analysis of complex systems. In this paper we analyze Jackson networks with the additional features that (i) nodes may have an infinite supply of low priority work and (ii) nodes may be unstable in the sense that the queue length at these nodes grows beyond any bound. We provide the limiting distribution of the queue length distribution at stable nodes, which turns out to be of product form. A key step in establishing this result is the development of a new algorithm based on adjusted traffic equations for detecting unstable nodes. Our results complement the results known in the literature for the subcases of Jackson networks with either infinite supply nodes or unstable nodes by providing an analysis of the significantly more challenging case of networks with both types of nonstandard node present. Building on our product-form results, we provide closed-form solutions for common customer and system oriented performance measures.

In this paper we analyze a tollbooth tandem queueing problem with an infinite number of servers. A customer starts service immediately upon arrival but cannot leave the system before all customers who arrived before him/her have left, i.e. customers depart the system in the same order as they arrive. Distributions of the total number of customers in the system, the number of departure-delayed customers in the system, and the number of customers in service at time t are obtained in closed form. Distributions of the sojourn times and departure delays of customers are also obtained explicitly. Both transient and steady state solutions are derived first for Poisson arrivals, and then extended to cases with batch Poisson and nonstationary Poisson arrival processes. Finally, we report several stochastic ordering results on how system performance measures are affected by arrival and service processes.

In this paper we consider an M/M/c queue modified to allow both mass arrivals when the system is empty and the workload to be removed. Properties of queues which terminate when the server becomes idle are firstly developed. Recurrence properties, equilibrium distribution, and equilibrium queue-size structure are studied for the case of resurrection and no mass exodus. All of these results are then generalized to allow for the removal of the entire workload. In particular, we obtain the Laplace transformation of the transition probability for the absorptive M/M/c queue.

In this paper we give an extension of the results of the generalized waiting time problem given by El-Desouky and Hussen (1990). An urn contains m types of balls of unequal numbers, and balls are drawn with replacement until first duplication. In the case of finite memory of order k, let ni be the number of type i, i = 1, 2, …, m. The probability of success pi = ni/N, i = 1, 2, …, m, where ni is a positive integer and Let Ym,k be the number of drawings required until first duplication. We obtain some new expressions of the probability function, in terms of Stirling numbers, symmetric polynomials, and generalized harmonic numbers. Moreover, some special cases are investigated. Finally, some important new combinatorial identities are obtained.