Search

A class of games for finding a leader among a group of candidates is studied in detail. This class covers games based on coin tossing and rock-paper-scissors as special cases and its complexity exhibits similar stochastic behaviors: either of logarithmic mean and bounded variance or of exponential mean and exponential variance. Many applications are also discussed.

In this paper, to model cascading failures, a new stochastic failure model is proposed. In a system subject to cascading failures, after each failure of the component, the remaining component suffers from increased load or stress. This results in shortened residual lifetimes of the remaining components. In this paper, to model this effect, the concept of the usual stochastic order is employed along with the accelerated life test model, and a new general class of stochastic failure models is generated.

In this paper, we first use the distribution of the number of records to demonstrate that the right tail probabilities of counts of rare events are generally better approximated by the right tail probabilities of Poisson distribution than {those} of normal distribution. We then show the moderate deviations in Poisson approximation generally require an adjustment and, with suitable adjustment, we establish better error estimates of the moderate deviations in Poisson approximation than those in \cite{CFS}. Our estimates contain no unspecified constants and are easy to apply. We illustrate the use of the theorems in six applications: Poisson-binomial distribution, matching problem, occupancy problem, birthday problem, random graphs and 2-runs. The paper complements the works of \cite{CC92,BCC95,CFS}., 29 pages and 5 figures

This paper generalizes the Kunita–Watanabe decomposition of an $L^2$ space. The generalization comes from using nonlinear stochastic integrals where the integrator is a family of continuous martingales bounded in $L^2$ . This result is also the solution of an optimization problem in $L^2$ . First, martingales are assumed to be stochastic integrals. Then, to get the general result, it is shown that the regularity of the family of martingales with respect to its spatial parameter is inherited by the integrands in the integral representation of the martingales. Finally, an example showing how the results of this paper, with the Clark–Ocone formula, can be applied to polynomial functions of Brownian integrals.

This paper considers two variants of M/G/1 queues with impatient customers, which are denoted by M/G/1+Gw and M/G/1+Gs. In the M/G/1+Gw queue customers have deadlines for their waiting times, and they leave the system immediately if their services do not start before the expiration of their deadlines. On the other hand, in the M/G/1+Gs queue customers have deadlines for their sojourn times, where customers in service also immediately leave the system when their deadlines expire. In this paper we derive comparison results for performance measures of these models. In particular, we show that if the service time distribution is new better than used in expectation, then the loss probability in the M/G/1+Gs queue is greater than that in the M/G/1+Gw queue.

The Ginibre point process (GPP) is one of the main examples of determinantal point processes on the complex plane. It is a recurring distribution of random matrix theory as well as a useful model in applied mathematics. In this paper we briefly overview the usual methods for the simulation of the GPP. Then we introduce a modified version of the GPP which constitutes a determinantal point process more suited for certain applications, and we detail its simulation. This modified GPP has the property of having a fixed number of points and having its support on a compact subset of the plane. See Decreusefond et al. (2013) for an extended version of this paper.

In a defaultable market, an investor trades having only partial information about the behavior of the market. Taking into account the intraday stock movements, the risky asset prices are modelled by marked point processes. Their dynamics depend on an unobservable process, representing the amount of news reaching the market. This is a marked point process, which may have common jump times with the risky asset price processes. The problem of hedging a defaultable claim is studied. In order to discuss all these topics, in this paper we examine stochastic control problems using backward stochastic differential equations (BSDEs) and filtering techniques. The goal of this paper is to construct a sequence of functions converging to the value function, each of these is the unique solution of a suitable BSDE.

In this paper we investigate the limiting behavior of the failure rate for the convolution of two or more life distributions. In a previous paper on mixtures, Block, Mi and Savits (1993) showed that the failure rate behaves like the limiting behavior of the strongest component. We show a similar result here for convolutions. We also show by example that unlike a mixture population, the ultimate direction of monotonicity does not necessarily follow that of the strongest component.

Let 𝓈 be a finite or countable set. Given a matrix F = (F ij ) i,j∈𝓈 of distribution functions on R and a quasistochastic matrix Q = (q ij ) i,j∈𝓈 , i.e. an irreducible nonnegative matrix with maximal eigenvalue 1 and associated unique (modulo scaling) positive left and right eigenvectors u and v, the matrix renewal measure ∑ n≥0 Q n ⊗ F *n associated with Q ⊗ F := (q ij F ij ) i,j∈𝓈 (see below for precise definitions) and a related Markov renewal equation are studied. This was done earlier by de Saporta (2003) and Sgibnev (2006, 2010) by drawing on potential theory, matrix-analytic methods, and Wiener-Hopf techniques. In this paper we describe a probabilistic approach which is quite different and starts from the observation that Q ⊗ F becomes an ordinary semi-Markov matrix after a harmonic transform. This allows us to relate Q ⊗ F to a Markov random walk {(M n , S n )} n≥0 with discrete recurrent driving chain {M n } n≥0. It is then shown that renewal theorems including a Choquet-Deny-type lemma may be easily established by resorting to standard renewal theory for ordinary random walks. The paper concludes with two typical examples.

Suppose that X 1, …, X n are random variables with the same known marginal distribution F but unknown dependence structure. In this paper we study the smallest possible value of P(X 1 + · · · + X n < s) over all possible dependence structures, denoted by m n,F (s). We show that m n,F (ns) → 0 for s no more than the mean of F under weak assumptions. We also derive a limit of m n,F (ns) for any s ∈ R with an error of at most n -1/6 for general continuous distributions. An application of our result to risk management confirms that the worst-case value at risk is asymptotically equivalent to the worst-case expected shortfall for risk aggregation with dependence uncertainty. In the last part of this paper we present a dual presentation of the theory of complete mixability and give dual proofs of theorems in the literature on this concept.

This paper considers risk-sensitive average optimization for denumerable continuous-time Markov decision processes (CTMDPs), in which the transition and cost rates are allowed to be unbounded, and the policies can be randomized history dependent. We first derive the multiplicative dynamic programming principle and some new facts for risk-sensitive finite-horizon CTMDPs. Then, we establish the existence and uniqueness of a solution to the risk-sensitive average optimality equation (RS-AOE) through the results for risk-sensitive finite-horizon CTMDPs developed here, and also prove the existence of an optimal stationary policy via the RS-AOE. Furthermore, for the case of finite actions available at each state, we construct a sequence of models of finite-state CTMDPs with optimal stationary policies which can be obtained by a policy iteration algorithm in a finite number of iterations, and prove that an average optimal policy for the case of infinitely countable states can be approximated by those of the finite-state models. Finally, we illustrate the conditions and the iteration algorithm with an example.

In this paper, we study a large system of $N$ servers each with capacity to process at most $C$ simultaneous jobs and an incoming job is routed to a server if it has the lowest occupancy amongst $d$ (out of N) randomly selected servers. A job that is routed to a server with no vacancy is assumed to be blocked and lost. Such randomized policies are referred to JSQ(d) (Join the Shortest Queue out of $d$) policies. Under the assumption that jobs arrive according to a Poisson process with rate $N\lambda^{(N)}$ where $\lambda^{(N)}=\sigma-\frac{\beta}{\sqrt{N}}$, $\sigma\in\mb{R}_+$ and $\beta\in\mb{R}$, we establish functional central limit theorems (FCLTs) for the fluctuation process in both the transient and stationary regimes when service time distributions are exponential. In particular, we show that the limit is an Ornstein-Uhlenbeck process whose mean and variance depend on the mean-field of the considered model. Using this, we obtain approximations to the blocking probabilities for large $N$, where we can precisely estimate the accuracy of first-order approximations., Comment: 29 pages

Comparison results for Markov processes with respect to function-class-induced (integral) stochastic orders have a long history. The most general results so far for this problem have been obtained based on the theory of evolution systems on Banach spaces. In this paper we transfer the martingale comparison method, known for the comparison of semimartingales to Markovian semimartingales, to general Markov processes. The basic step of this martingale approach is the derivation of the supermartingale property of the linking process, giving a link between the processes to be compared. This property is achieved using the characterization of Markov processes by the associated martingale problem in an essential way. As a result, the martingale comparison method gives a comparison result for Markov processes under a general alternative but related set of regularity conditions compared to the evolution system approach.

This paper investigates a financial market where stock returns depend on an unobservable Gaussian mean reverting drift process. Information on the drift is obtained from returns and randomly arriving discrete-time expert opinions. Drift estimates are based on Kalman filter techniques. We study the asymptotic behavior of the filter for high-frequency experts with variances that grow linearly with the arrival intensity. The derived limit theorems state that the information provided by discrete-time expert opinions is asymptotically the same as that from observing a certain diffusion process. These diffusion approximations are extremely helpful for deriving simplified approximate solutions of utility maximization problems.

Block-structured Markov chains model a large variety of queueing problems and have many important applications in various areas. Stability properties have been well investigated for these Markov chains. In this paper we will present transient properties for two specific types of block-structured Markov chains, including M/G/1 type and GI/M/1 type. Necessary and sufficient conditions in terms of system parameters are obtained for geometric transience and algebraic transience. Possible extensions of the results to continuous-time Markov chains are also included.

In this paper we address the question of finding the point which maximizes the pth moment of the exit time of planar Brownian motion from a given domain. We present a geometrical method for excluding parts of the domain from consideration which makes use of a coupling argument and the conformal invariance of Brownian motion. In many cases the maximizing point can be localized to a relatively small region. Several illustrative examples are presented.

In this paper we study first passage percolation on a random graph model, the configuration model. We first introduce the notions of weighted diameter, which is the maximum of the weighted lengths of all optimal paths between any two vertices in the graph, and the flooding time, which represents the time (weighted length) needed to reach all the vertices in the graph starting from a uniformly chosen vertex. Our result consists in describing the asymptotic behavior of the diameter and the flooding time, as the number of vertices n tends to infinity, in the case where the weight distribution G has an exponential tail behavior, and proving that this category of distributions is the largest possible for which the asymptotic behavior holds.

In this paper the behaviour of the failure rate and reversed failure rate of an n-component coherent system is studied, where it is assumed that the lifetimes of the components are independent and have a common cumulative distribution function F. Sufficient conditions are provided under which the system failure rate is increasing and the corresponding reversed failure rate is decreasing. We also study the stochastic and ageing properties of doubly truncated random variables for coherent systems.

The signature representation shows that the reliability of the system is a mixture of the reliability functions of the k-out-of-n systems. The first representation was obtained for systems with independent and identically distributed (IID) components and after it was extended to exchangeable (EXC) components. The purpose of the present paper is to extend it to the class of systems with identically distributed (ID) components which have a diagonal-dependent copula. We prove that this class is much larger than the class with EXC components. This extension is used to compare systems with non-EXC components.

This paper investigates the random horizon optimal stopping problem for measure-valued piecewise deterministic Markov processes (PDMPs). This is motivated by population dynamics applications, when one wants to monitor some characteristics of the individuals in a small population. The population and its individual characteristics can be represented by a point measure. We first define a PDMP on a space of locally finite measures. Then we define a sequence of random horizon optimal stopping problems for such processes. We prove that the value function of the problems can be obtained by iterating some dynamic programming operator. Finally we prove via a simple counter-example that controlling the whole population is not equivalent to controlling a random lineage.

This paper presents a flexible family which we call the $\alpha$ -mixture of survival functions. This family includes the survival mixture, failure rate mixture, models that are stochastically closer to each of these conventional mixtures, and many other models. The $\alpha$ -mixture is endowed by the stochastic order and uniquely possesses a mathematical property known in economics as the constant elasticity of substitution, which provides an interpretation for $\alpha$ . We study failure rate properties of this family and establish closures under monotone failure rates of the mixture’s components. Examples include potential applications for comparing systems.

In this paper, we obtain an asymptotic formula for the persistence probability in the positive real line of a random polynomial arising from evolutionary game theory. It corresponds to the probability that a multi-player two-strategy random evolutionary game has no internal equilibria. The key ingredient is to approximate the sequence of random polynomials indexed by their degrees by an appropriate centered stationary Gaussian process., revised version

In this paper, we analyze the set of all possible aggregate distributions of the sum of standard uniform random variables, a simply stated yet challenging problem in the literature of distributions with given margins. Our main results are obtained for two distinct cases. In the case of dimension two, we obtain four partial characterization results. An analytical result with full generality is not available in this case, and it seems to be out of reach with existing techniques. For dimension greater or equal to three, we obtain a full characterization of the set of aggregate distributions, which is the first complete characterization result of this type in the literature for any choice of continuous marginal distributions.

The De Vylder and Goovaerts conjecture is an open problem in risk theory, stating that the finite-time ruin probability in a standard risk model is greater than or equal to the corresponding ruin probability evaluated in an associated model with equalized claim amounts. Equalized means here that the jump sizes of the associated model are equal to the average jump in the initial model between 0 and a terminal time T.In this paper, we consider the diffusion approximations of both the standard risk model and its associated risk model. We prove that the associated model, when conveniently renormalized, converges in distribution to a Gaussian process satisfying a simple SDE. We then compute the probability that this diffusion hits the level 0 before time T and compare it with the same probability for the diffusion approximation for the standard risk model. We conclude that the De Vylder and Goovaerts conjecture holds for the diffusion limits.

A new network evolution model is introduced in this paper. The model is based on cooperations of N units. The units are the nodes of the network and the cooperations are indicated by directed links. At each evolution step N units cooperate, which formally means that they form a directed N-star subgraph. At each step either a new unit joins the network and it cooperates with N − 1 old units, or N old units cooperate. During the evolution both preferential attachment and uniform choice are applied. Asymptotic power law distributions are obtained both for in-degrees and for out-degrees.

The signature of a coherent system has been studied extensively in the recent literature. Signatures are particularly useful in the comparison of coherent or mixed systems under a variety of stochastic orderings. Also, certain signature-based closure and preservation theorems have been established. For example, it is now well known that certain stochastic orderings are preserved from signatures to system lifetimes when components have independent and identical distributions. This applies to the likelihood ratio order, the hazard rate order, and the stochastic order. The point of departure of the present paper is the question of whether or not a similar preservation result will hold for the mean residual life order. A counterexample is provided which shows that the answer is negative. Classes of distributions for the component lifetimes for which the latter implication holds are then derived. Connections to the theory of order statistics are also considered.

In this paper we integrate two strands of the literature on stability of general state Markov chains: conventional, total-variation-based results and more recent order-theoretic results. First we introduce a complete metric over Borel probability measures based on ‘partial’ stochastic dominance. We then show that many conventional results framed in the setting of total variation distance have natural generalizations to the partially ordered setting when this metric is adopted.

In this paper we prove asymptotic upper bounds on the variance of the number of vertices and the missed area of inscribed random disc-polygons in smooth convex discs whose boundary isC+2. We also consider a circumscribed variant of this probability model in which the convex disc is approximated by the intersection of random circles.

Consider a sample 𝒳n={X1,…,Xn} of independent and identically distributed variables drawn with a probability distribution ℙX supported on a compact set M⊂ℝd. In this paper we mainly deal with the study of a natural estimator for the geodesic distance on M. Under rather general geometric assumptions on M, we prove a general convergence result. Assuming M to be a compact manifold of known dimension d′≤d, and under regularity assumptions on ℙX, we give an explicit convergence rate. In the case when M has no boundary, knowledge of the dimension d′ is not needed to obtain this convergence rate. The second part of the work consists in building an estimator for the Fréchet expectations on M, and proving its convergence under regularity conditions, applying the previous results.

De Finetti’s optimal dividend problem has recently been extended to the case when dividend payments can be made only at Poisson arrival times. In this paper we consider the version with bail-outs where the surplus must be nonnegative uniformly in time. For a general spectrally negative Lévy model, we show the optimality of a Parisian-classical reflection strategy that pays the excess above a given barrier at each Poisson arrival time and also reflects from below at 0 in the classical sense.

The study of finite approximations of probability measures has a long history. In Xu and Berger (2017), the authors focused on constrained finite approximations and, in particular, uniform ones in dimensiond=1. In the present paper we give an elementary construction of a uniform decomposition of probability measures in dimensiond≥1. We then use this decomposition to obtain upper bounds on the rate of convergence of the optimal uniform approximation error. These bounds appear to be the generalization of the ones obtained by Xu and Berger (2017) and to be sharp for generic probability measures.

In this paper we consider the degree-wise effect of a second step for a random walk on a graph. We prove that under the configuration model, for any fixed degree sequence the probability of exceeding a given degree threshold is smaller after two steps than after one. This builds on recent work of Krameret al.(2016) regarding the friendship paradox under random walks.

The study of high-dimensional distributions is of interest in probability theory, statistics, and asymptotic convex geometry, where the object of interest is the uniform distribution on a convex set in high dimensions. The ℓp-spaces and norms are of particular interest in this setting. In this paper we establish a limit theorem for distributions on ℓp-spheres, conditioned on a rare event, in a high-dimensional geometric setting. As part of our proof, we establish a certain large deviation principle that is also relevant to the study of the tail behavior of random projections of ℓp-balls in a high-dimensional Euclidean space.

In this paper we consider an infinite system of instantaneously coalescing rate 1 simple symmetric random walks on ℤ2, started from the initial condition with all sites in ℤ2 occupied. Two-dimensional coalescing random walks are a `critical' model of interacting particle systems: unlike coalescence models in dimension three or higher, the fluctuation effects are important for the description of large-time statistics in two dimensions, manifesting themselves through the logarithmic corrections to the `mean field' answers. Yet the fluctuation effects are not as strong as for the one-dimensional coalescence, in which case the fluctuation effects modify the large time statistics at the leading order. Unfortunately, unlike its one-dimensional counterpart, the two-dimensional model is not exactly solvable, which explains a relative scarcity of rigorous analytic answers for the statistics of fluctuations at large times. Our contribution is to find, for any N≥2, the leading asymptotics for the correlation functions ρN(x1,…,xN) as t→∞. This generalises the results for N=1 due to Bramson and Griffeath (1980) and confirms a prediction in the physics literature for N>1. An analogous statement holds for instantaneously annihilating random walks. The key tools are the known asymptotic ρ1(t)∼logt∕πt due to Bramson and Griffeath (1980), and the noncollision probability 𝒑NC(t), that no pair of a finite collection of N two-dimensional simple random walks meets by time t, whose asymptotic 𝒑NC(t)∼c0(logt)-(N2) was found by Cox et al. (2010). We re-derive the asymptotics, and establish new error bounds, both for ρ1(t) and 𝒑NC(t) by proving that these quantities satisfy effective rate equations; that is, approximate differential equations at large times. This approach can be regarded as a generalisation of the Smoluchowski theory of renormalised rate equations to multi-point statistics.

The mean time to failure (MTTF) function in age replacement is used to evaluate the performance and effectiveness of the age replacement policy. In this paper, based on the MTTF function, we introduce two new nonparametric classes of lifetime distributions with nonmonotonic mean time to failure in age replacement; increasing then decreasing MTTF (IDMTTF) and decreasing then increasing MTTF (DIMTTF). The implications between these classes of distributions and some existing classes of nonmonotonic ageing classes are studied. The characterizations of IDMTTF and DIMTTF in terms of the scaled total time on test transform are also obtained.

In this paper we are concerned with the reliability properties of two coherent systems having shared components. We assume that the components of the systems are two overlapping subsets of a set of n components with lifetimes X1,...,Xn. Further, we assume that the components of the systems fail according to the model of sequential order statistics (which is equivalent, under some mild conditions, to the failure model corresponding to a nonhomogeneous pure-birth process). The joint reliability function of the system lifetimes is expressed as a mixture of the joint reliability functions of the sequential order statistics, where the mixing probabilities are the bivariate signature matrix associated to the structures of systems. We investigate some stochastic orderings and dependency properties of the system lifetimes. We also study conditions under which the joint reliability function of systems with shared components of order m can be equivalently written as the joint reliability function of systems of order n (n>m). In order to illustrate the results, we provide several examples.

Goldman (2010) proved that the distribution of a stationary determinantal point process (DPP) Φ can be coupled with its reduced Palm version Φ0,! such that there exists a point process η where Φ=Φ0,!∪η in distribution and Φ0,!∩η=∅. The points of η characterize the repulsive nature of a typical point of Φ. In this paper we use the first-moment measure of η to study the repulsive behavior of DPPs in high dimensions. We show that many families of DPPs have the property that the total number of points in η converges in probability to 0 as the space dimension n→∞. We also prove that for some DPPs, there exists an R∗ such that the decay of the first-moment measure of η is slowest in a small annulus around the sphere of radius √nR∗. This R∗ can be interpreted as the asymptotic reach of repulsion of the DPP. Examples of classes of DPP models exhibiting this behavior are presented and an application to high-dimensional Boolean models is given.

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.