The ‘coupon collection problem’ refers to a class of occupancy problems in which j identical items are distributed, independently and at random, to n cells, with no restrictions on multiple occupancy. Identifying the cells as coupons, a coupon is ‘collected’ if the cell is occupied by one or more of the distributed items; thus, some coupons may never be collected, whereas others may be collected once or twice or more. We call the number of coupons collected exactly r times coupons of type r. The coupon collection model we consider is general, in that a random number of purchases occurs at each stage of collecting a large number of coupons; the sample sizes at each stage are independent and identically distributed according to a sampling distribution. The joint behavior of the various types is an intricate problem. In fact, there is a variety of joint central limit theorems (and other limit laws) that arise according to the interrelation between the mean, variance, and range of the sampling distribution, and of course the phase (how far we are in the collection processes). According to an appropriate combination of the mean of the sampling distribution and the number of available coupons, the phase is sublinear, linear, or superlinear. In the sublinear phase, the normalization that produces a Gaussian limit law for uncollected coupons can be used to obtain a multivariate central limit law for at most two other types — depending on the rates of growth of the mean and variance of the sampling distribution, we may have a joint central limit theorem between types 0 and 1, or between types 0, 1, and 2. In the linear phase we have a multivariate central limit theorem among the types 0, 1,…, k for any fixed k.

Hazard rates play an important role in various areas, e.g. reliability theory, survival analysis, biostatistics, queueing theory, and actuarial studies. Mixtures of distributions are also of great preeminence in such areas as most populations of components are indeed heterogeneous. In this study we present a sufficient condition for mixtures of two elements of the same natural exponential family (NEF) to have an increasing hazard rate. We then apply this condition to some classical NEFs having either quadratic or cubic variance functions (VFs) and others as well. Particular attention is paid to the hyperbolic cosine NEF having a quadratic VF, the Ressel NEF having a cubic VF, and the NEF generated by Kummer distributions of type 2. The application of such a sufficient condition is quite intricate and cumbersome, in particular when applied to the latter three NEFs. Various lemmas and propositions are needed to verify this condition for such NEFs. It should be pointed out, however, that our results are mainly applied to a mixture of two populations.

We consider regular variation of a Lévy process X := (X t)t≥0 in with Lévy measure Π, emphasizing the dependence between jumps of its components. By transforming the one-dimensional marginal Lévy measures to those of a standard 1-stable Lévy process, we decouple the marginal Lévy measures from the dependence structure. The dependence between the jumps is modeled by a so-called Pareto Lévy measure, which is a natural standardization in the context of regular variation. We characterize multivariate regularly variation of X by its one-dimensional marginal Lévy measures and the Pareto Lévy measure. Moreover, we define upper and lower tail dependence coefficients for the Lévy measure, which also apply to the multivariate distributions of the process. Finally, we present graphical tools to visualize the dependence structure in terms of the spectral density and the tail integral for homogeneous and nonhomogeneous Pareto Lévy measures.

The distributions for continuous, discrete, and conditional discrete scan statistics are studied. The approach of finite Markov chain imbedding, which has been applied to random permutations as well as to runs and patterns, is extended to compute the distribution of the conditional discrete scan statistic, defined from a sequence of Bernoulli trials. It is shown that the distribution of the continuous scan statistic induced by a Poisson process defined on (0, 1] is a limiting distribution of weighted distributions of conditional discrete scan statistics. Comparisons of rates of convergence as well as numerical comparisons of various bounds and approximations are provided to illustrate the theoretical results.

We use a change-of-variable formula in the framework of functions of bounded variation to derive an explicit formula for the Fourier transform of the level crossing function of shot noise processes with jumps. We illustrate the result in some examples and give some applications. In particular, it allows us to study the asymptotic behavior of the mean number of level crossings as the intensity of the Poisson point process of the shot noise process goes to infinity.

We give an efficient method based on minimal deterministic finite automata for computing the exact distribution of the number of occurrences and coverage of clumps (maximal sets of overlapping words) of a collection of words. In addition, we compute probabilities for the number of h-clumps, word groupings where gaps of a maximal length h between occurrences of words are allowed. The method facilitates the computation of p-values for testing procedures. A word is allowed to contain other words of the collection, making the computation more general, but also more difficult. The underlying sequence is assumed to be Markovian of an arbitrary order.

In this paper we provide some sufficient conditions to stochastically compare linear combinations of independent random variables. The main results extend those given in Proschan (1965), Ma (1998), Zhao et al. (2011), and Yu (2011). In particular, we propose a new sufficient condition to compare the peakedness of linear combinations of independent random variables which may have heavy-tailed properties.

Let x denote a vector of length q consisting of 0s and 1s. It can be interpreted as an ‘opinion’ comprised of a particular set of responses to a questionnaire consisting of q questions, each having {0, 1}-valued answers. Suppose that the questionnaire is answered by n individuals, thus providing n ‘opinions’. Probabilities of the answer 1 to each question can be, basically, arbitrary and different for different questions. Out of the 2q different opinions, what number, μn , would one expect to see in the sample? How many of these opinions, μn (k), will occur exactly k times? In this paper we give an asymptotic expression for μn / 2q and the limit for the ratios μn (k)/μn , when the number of questions q increases along with the sample size n so that n = λ2q , where λ is a constant. Let p(x ) denote the probability of opinion x . The key step in proving the asymptotic results as indicated is the asymptotic analysis of the joint behaviour of the intensities np(x ). For example, one of our results states that, under certain natural conditions, for any z > 0, ∑1 {np( x) > z} = d n z −u ,d n = o(2q ).

Let η = (η1,…,ηn ) be a positive random vector. If its coordinates ηi and ηj are exchangeable, i.e. the distribution of η is invariant with respect to the swap πij of its ith and jth coordinates, then Ef(η) = Ef(πij η) for all integrable functions f. In this paper we study integrable random vectors that satisfy this identity for a particular family of functions f, namely those which can be written as the positive part of the scalar product 〈u, η〉 with varying weights u. In finance such functions represent payoffs from exchange options with η being the random part of price changes, while from the geometric point of view they determine the support function of the so-called zonoid of η. If the expected values of such payoffs are πij -invariant, we say that η is ij-swap-invariant. A full characterisation of the swap-invariance property and its relationship to the symmetries of expected payoffs of basket options are obtained. The first of these results relies on a characterisation theorem for integrable positive random vectors with equal zonoids. Particular attention is devoted to the case of asset prices driven by Lévy processes. Based on this, concrete semi-static hedging techniques for multi-asset barrier options, such as weighted barrier swap options, weighted barrier quanto-swap options, or certain weighted barrier spread options, are suggested.

Signature-based representations of the reliability functions of coherent systems with independent and identically distributed component lifetimes have proven very useful in studying the ageing characteristics of such systems and in comparing the performance of different systems under varied criteria. In this paper we consider extensions of these results to systems with heterogeneous components. New representation theorems are established for both the case of components with independent lifetimes and the case of component lifetimes under specific forms of dependence. These representations may be used to compare the performance of systems with homogeneous and heterogeneous components.

In this paper, some ordering properties of convolutions of heterogeneous Bernoulli random variables are discussed. It is shown that, under some suitable conditions, the likelihood ratio order and the reversed hazard rate order hold between convolutions of two heterogeneous Bernoulli sequences. The results established here extend and strengthen the previous results of Pledger and Proschan (1971) and Boland, Singh and Cukic (2002).

We use the properties of the Matuszewska indices to show asymptotic inequalities for hazard rates. We discuss the relation between membership in the classes of dominatedly or extended rapidly varying tail distributions and corresponding hazard rate conditions. Convolution closure is established for the class of distributions with extended rapidly varying tails.

Exact lower bounds on the exponential moments of min(y, X) and X 1{X < y} are provided given the first two moments of a random variable X. These bounds are useful in work on large deviation probabilities and nonuniform Berry-Esseen bounds, when the Cramér tilt transform may be employed. Asymptotic properties of these lower bounds are presented. Comparative advantages of the so-called Winsorization min(y, X) over the truncation X 1{X < y} are demonstrated. An application to option pricing is given.

Improved bounds on the copula of a bivariate random vector are computed when partial information is available, such as the values of the copula on a given subset of [0, 1]2, or the value of a functional of the copula, monotone with respect to the concordance order. These results are then used to compute model-free bounds on the prices of two-asset options which make use of extra information about the dependence structure, such as the price of another two-asset option.

This paper is motivated by relations between association and independence of random variables. It is well known that, for real random variables, independence implies association in the sense of Esary, Proschan and Walkup (1967), while, for random vectors, this simple relationship breaks. We modify the notion of association in such a way that any vector-valued process with independent increments also has associated increments in the new sense - association between blocks. The new notion is quite natural and admits nice characterization for some classes of processes. In particular, using the covariance interpolation formula due to Houdré, Pérez-Abreu and Surgailis (1998), we show that within the class of multidimensional Gaussian processes, block association of increments is equivalent to supermodularity (in time) of the covariance functions. We also define corresponding versions of weak association, positive association, and negative association. It turns out that the central limit theorem for weakly associated random vectors due to Burton, Dabrowski and Dehling (1986) remains valid, if the weak association is relaxed to the weak association between blocks.

In this work we provide sufficient conditions under which a general counting process stopped at a random time independent from the process belongs to the reliability decreasing reversed hazard rate (DRHR) or increasing failure rate (IFR) class. We also give some applications of these results in generalized renewal and trend renewal processes stopped at a random time.

We present a method for computing the probability density function (PDF) and the cumulative distribution function (CDF) of a nonnegative infinitely divisible random variable X. Our method uses the Lévy-Khintchine representation of the Laplace transform Ee-λX = e-ϕ(λ), where ϕ is the Laplace exponent. We apply the Post-Widder method for Laplace transform inversion combined with a sequence convergence accelerator to obtain accurate results. We demonstrate this technique on several examples, including the stable distribution, mixtures thereof, and integrals with respect to nonnegative Lévy processes.

Kochar and Xu (2009) proved that a parallel system with heterogeneous exponential component lifetimes is more skewed (according to the convex transform order) than the system with independent and identically distributed exponential components. In this paper we extend this study to the general k-out-of-n systems for the case when there are only two types of component in the system. An open problem proposed in Pǎltǎnea (2008) is partially solved.

We establish an upper bound on the tails of a random variable that arises as a solution of a stochastic difference equation. In the nonnegative case our bound is similar to a lower bound obtained in Goldie and Grübel (1996).

Consider a continuous-time renewal risk model with a constant force of interest. We assume that claim sizes and interarrival times correspondingly form a sequence of independent and identically distributed random pairs and that each pair obeys a dependence structure described via the conditional tail probability of a claim size given the interarrival time before the claim. We focus on determining the impact of this dependence structure on the asymptotic tail probability of discounted aggregate claims. Assuming that the claim size distribution is subexponential, we derive an exact locally uniform asymptotic formula, which quantitatively captures the impact of the dependence structure. When the claim size distribution is extended regularly varying tailed, we show that this asymptotic formula is globally uniform.