Tankov (2011) improved the Fréchet bounds for a bivariate copula when its values on a compact subset of [0, 1]2 are given. He showed that the best possible bounds are quasi-copulas and gave a sufficient condition for these bounds to be copulas. In this note we give weaker sufficient conditions to ensure that the bounds are copulas. We also show how this can be useful in portfolio selection. It turns out that finding a copula as a lower bound plays a key role in determining optimal investment strategies explicitly for investors with some type of state-dependent constraints.

Nonnegative infinitely divisible (i.d.) random variables form an important class of random variables. However, when this type of random variable is specified via Lévy densities that have infinite integrals on (0, ∞), except for some special cases, exact sampling is unknown. We present a method that can sample a rather wide range of such i.d. random variables. A basic result is that, for any nonnegative i.d. random variable X with its Lévy density explicitly specified, if its distribution conditional on X ≤ r can be sampled exactly, where r > 0 is any fixed number, then X can be sampled exactly using rejection sampling, without knowing the explicit expression of the density of X. We show that variations of the result can be used to sample various nonnegative i.d. random variables.

In this paper we consider three-dimensional random tessellations that are stable under iteration (STIT tessellations). STIT tessellations arise as a result of subsequent cell division, which implies that their cells are not face-to-face. The edges of the cell-dividing polygons are the so-called I-segments of the tessellation. The main result is an explicit formula for the distribution of the number of vertices in the relative interior of the typical I-segment. In preparation for its proof, we obtain other distributional identities for the typical I-segment and the length-weighted typical I-segment, which provide new insight into the spatiotemporal construction process.

Given a probability density, we estimate the rate of decay of the measure of the level sets of its evolutes by the Ornstein–Uhlenbeck semigroup. The rate is faster than what follows from the preservation of mass and Markov’s inequality.

Almost all populations existing in the real world are finite populations. Specifically, in the areas relevant to lifetime modeling and analysis, finite populations are frequently encountered. However, descriptions of failure/survival patterns of elements in the finite population have not yet been properly established. In particular, it is questionable whether the ordinary failure rate can be defined for finite populations in the same way and whether the corresponding interpretations are still valid. In this paper we consider two kinds of finite mixed population and provide new definitions for their failure rates. Then we clarify the notion of failure rate in finite populations.

We consider a multivariate distributional recursion of sum type, as arises in the probabilistic analysis of algorithms and random trees. We prove an upper tail bound for the solution using Chernoff's bounding technique by estimating the Laplace transform. The problem is traced back to the corresponding problem for binary search trees by stochastic domination. The result obtained is applied to the internal path length and Wiener index of random b-ary recursive trees with weighted edges and random linear recursive trees. Finally, lower tail bounds for the Wiener index of these trees are given.

We study a family of distributions that satisfy the stability-under-addition property, provided that the number ν of random variables in a sum is also a random variable. We call the corresponding property ν-stability and investigate the situation when the semigroup generated by the generating function of ν is commutative. Using results from the theory of iterations of analytic functions, we describe ν-stable distributions generated by summations with rational generating functions. A new case in this class of distributions arises when generating functions are linked with Chebyshev polynomials. The analogue of normal distribution corresponds to the hyperbolic secant distribution.

In the study of the reliability of technical systems in reliability engineering, coherent systems play a key role. In this paper we consider a coherent system consisting of n components with independent and identically distributed components and propose two time-dependent criteria. The first criterion is a measure of the residual lifetime of live components of a coherent system having some of the components alive when the system fails at time t. The second criterion is a time-dependent measure which enables us to investigate the inactivity times of the failed components of a coherent system still functioning though some of its components have failed. Several ageing and stochastic properties of the proposed measures are then established.

The concept of complete mixability is relevant to some problems of optimal couplings with important applications in quantitative risk management. In this paper we prove new properties of the set of completely mixable distributions, including a completeness and a decomposition theorem. We also show that distributions with a concave density and radially symmetric distributions are completely mixable.

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 := (Xt)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).