scholarly journals Exact simulation of generalised Vervaat perpetuities

2019 ◽  
Vol 56 (01) ◽  
pp. 57-75 ◽  
Author(s):  
Angelos Dassios ◽  
Yan Qu ◽  
Jia Wei Lim

AbstractWe consider a generalised Vervaat perpetuity of the form X = Y1W1 +Y2W1W2 + · · ·, where $W_i \sim {\cal U}^{1/t}$ and (Yi)i≥0 is an independent and identically distributed sequence of random variables independent from (Wi)i≥0. Based on a distributional decomposition technique, we propose a novel method for exactly simulating the generalised Vervaat perpetuity. The general framework relies on the exact simulation of the truncated gamma process, which we develop using a marked renewal representation for its paths. Furthermore, a special case arises when Yi = 1, and X has the generalised Dickman distribution, for which we present an exact simulation algorithm using the marked renewal approach. In particular, this new algorithm is much faster than existing algorithms illustrated in Chi (2012), Cloud and Huber (2017), Devroye and Fawzi (2010), and Fill and Huber (2010), as well as being applicable to the general payments case. Examples and numerical analysis are provided to demonstrate the accuracy and effectiveness of our method.

2021 ◽  
Vol 58 (2) ◽  
pp. 347-371
Author(s):  
Yan Qu ◽  
Angelos Dassios ◽  
Hongbiao Zhao

AbstractThere are two types of tempered stable (TS) based Ornstein–Uhlenbeck (OU) processes: (i) the OU-TS process, the OU process driven by a TS subordinator, and (ii) the TS-OU process, the OU process with TS marginal law. They have various applications in financial engineering and econometrics. In the literature, only the second type under the stationary assumption has an exact simulation algorithm. In this paper we develop a unified approach to exactly simulate both types without the stationary assumption. It is mainly based on the distributional decomposition of stochastic processes with the aid of an acceptance–rejection scheme. As the inverse Gaussian distribution is an important special case of TS distribution, we also provide tailored algorithms for the corresponding OU processes. Numerical experiments and tests are reported to demonstrate the accuracy and effectiveness of our algorithms, and some further extensions are also discussed.


1996 ◽  
Vol 33 (01) ◽  
pp. 146-155 ◽  
Author(s):  
K. Borovkov ◽  
D. Pfeifer

In this paper we consider improvements in the rate of approximation for the distribution of sums of independent Bernoulli random variables via convolutions of Poisson measures with signed measures of specific type. As a special case, the distribution of the number of records in an i.i.d. sequence of length n is investigated. For this particular example, it is shown that the usual rate of Poisson approximation of O(1/log n) can be lowered to O(1/n 2). The general case is discussed in terms of operator semigroups.


1980 ◽  
Vol 12 (01) ◽  
pp. 200-221 ◽  
Author(s):  
B. Natvig

In this paper we arrive at a series of bounds for the availability and unavailability in the time interval I = [t A , t B ] ⊂ [0, ∞), for a coherent system of maintained, interdependent components. These generalize the minimal cut lower bound for the availability in [0, t] given in Esary and Proschan (1970) and also most bounds for the reliability at time t given in Bodin (1970) and Barlow and Proschan (1975). In the latter special case also some new improved bounds are given. The bounds arrived at are of great interest when trying to predict the performance process of the system. In particular, Lewis et al. (1978) have revealed the great need for adequate tools to treat the dependence between the random variables of interest when considering the safety of nuclear reactors. Satyanarayana and Prabhakar (1978) give a rapid algorithm for computing exact system reliability at time t. This can also be used in cases where some simpler assumptions on the dependence between the components are made. It seems, however, impossible to extend their approach to obtain exact results for the cases treated in the present paper.


2016 ◽  
Vol 24 (1) ◽  
pp. 29-41 ◽  
Author(s):  
Roman Frič ◽  
Martin Papčo

Abstract The influence of “Grundbegriffe” by A. N. Kolmogorov (published in 1933) on education in the area of probability and its impact on research in stochastics cannot be overestimated. We would like to point out three aspects of the classical probability theory “calling for” an upgrade: (i) classical random events are black-and-white (Boolean); (ii) classical random variables do not model quantum phenomena; (iii) basic maps (probability measures and observables { dual maps to random variables) have very different “mathematical nature”. Accordingly, we propose an upgraded probability theory based on Łukasiewicz operations (multivalued logic) on events, elementary category theory, and covering the classical probability theory as a special case. The upgrade can be compared to replacing calculations with integers by calculations with rational (and real) numbers. Namely, to avoid the three objections, we embed the classical (Boolean) random events (represented by the f0; 1g-valued indicator functions of sets) into upgraded random events (represented by measurable {0; 1}-valued functions), the minimal domain of probability containing “fractions” of classical random events, and we upgrade the notions of probability measure and random variable.


Mathematics ◽  
2019 ◽  
Vol 7 (8) ◽  
pp. 727 ◽  
Author(s):  
Dongming Nie ◽  
Saima Rashid ◽  
Ahmet Ocak Akdemir ◽  
Dumitru Baleanu ◽  
Jia-Bao Liu

In this article, we aim to establish several inequalities for differentiable exponentially convex and exponentially quasi-convex mapping, which are connected with the famous Hermite–Hadamard (HH) integral inequality. Moreover, we have provided applications of our findings to error estimations in numerical analysis and higher moments of random variables.


2008 ◽  
Vol 24 (5) ◽  
pp. 1443-1455 ◽  
Author(s):  
James Davidson ◽  
Jan R. Magnus ◽  
Jan Wiegerinck

We consider the Breitung (2002, Journal of Econometrics 108, 343–363) statistic ξn, which provides a nonparametric test of the I(1) hypothesis. If ξ denotes the limit in distribution of ξn as n → ∞, we prove (Theorem 1) that 0 ≤ ξ ≤ 1/π2, a result that holds under any assumption on the underlying random variables. The result is a special case of a more general result (Theorem 3), which we prove using the so-called cotangent method associated with Cauchy's residue theorem.


Author(s):  
T. P. Speed

AbstractEarlier work of the author exploiting the role of partition lattices and their Mbius functions in the theory of cumulants, k-statistics and their generalisations is extended to multiply-indexed arrays of random variables. The natural generalisations of cumulants and k-statistics to this context are shown to include components of variance and the associated linear combinations of mean-squares which are used to estimate them. Expressions for the generalised cumulants of arrays built up as sums of independent arrays of effects as in anova models are derived in terms of the generalized cumulants of the effects. The special case of degree two, covering the unbiased estimation of components of variance, is discussed in some detail.


1980 ◽  
Vol 17 (01) ◽  
pp. 102-111 ◽  
Author(s):  
Arthur V. Peterson ◽  
Richard A. Kronmal

We obtain a representation of an arbitrary discrete distribution with n mass points by an equiprobable mixture of r distributions, each of which has no more than a (≧2) mass points, where r is the smallest integer greater than or equal to (n – 1)/(a – 1). An application to the generation of discrete random variables on a computer is described, which has as an important special case Walker's (1977) alias method.


Author(s):  
JANUSZ WYSOCZAŃSKI

We define a deformation of free creations (and annihilations), given by operators on the full Fock space, acting nontrivially only between the vacuum subspace ℂΩ and the twofold tensor product [Formula: see text]. Then we study the distribution of the deformed free gaussian operators, with the deformation containing also a real parameter d. The recurrence formula for moments is shown, and the Cauchy transform of the distribution measure is computed. This yields the description of the measure: absolutely continuous part and the atomic part. The existence of atoms depends on the parameter d. The special case d =1 is studied with all details, with the formula for moments is given as values of the hypergeometric series. Finally we show the formula for computing the mixed moments of the deformed operators.


Sign in / Sign up

Export Citation Format

Share Document