Iterated distributions were introduced by Averous and Meste , in order to construct a classification of lifetime distributions, and were initially studied in a more systematic way by Fagiuoli and Pellerey . The iteration procedure described below produces what is also known as equilibrium distributions in economics or actuarial activities, since they describe the distribution of the first drop below the initial reserve. Moreover, equilibrium distributions play an important role in ageing relations (see, for example, Chatterjee and Mukherjee ) or in renewal theory (see Cox 
). A suitable representation of the iterated distributions describes their tails as normalized moments of stop-loss risk premiums with a given deductible, a common type of contract in actuarial activity, with interest to the characterization of ruin probabilities and insolvency. We obtain a characterization of the asymptotics of these moments, with respect to their order, for initial random variables with distribution in the Gamma or Weibull families. These results provide simple numerical approximations for the stop-loss premiums.
Let us describe our framework and introduce the basic notation and transformations. We will be assuming throughout that is a nonnegative random variable with density function , distribution function , and tail function .
For each , define
For each , define the iterated distribution induced by , by its tail as follows:
assuming the integrals above are finite.
The distribution is known as the equilibrium distribution of . Hence, the iteration process above defines, for each , the equilibrium distribution of a random variable with tail . Thus, the iteration procedure given in Definition 1 may be restated in terms of the iterated equilibrium distribution.
The iterated equilibrium distribution of has tail given by , that is, the iterated distribution induced by described in Definition 1.
Taking into account this identification, we shall refer throughout to iterated distributions instead of iterated equilibrium, as these only differ on the count of the iteration steps.
It is easily verified that the iterated distribution induced by an exponential random variable
is exactly the same exponential distribution as. That is, the exponential distributions are fixed points with respect to the iteration procedure introduced in Definition 1. Moreover, it is also easily verified that explicit identification of the iterated distributions may become in general quite complex once we leave the case of being exponentially distributed.
As already mentioned, iterated distributions were used by Averous and Meste 
to classify the distribution ofwith respect to its tail behavior. General properties of the iterated distributions were studied in Fagiuoli and Pellerey  and Nanda et al. , but to the best of our knowledge, the behavior of the iterated distributions with respect to the iteration step has not been studied before. In this paper, we will be mainly interested on this behavior of iterated distributions as the iteration step grows.
Although the iterated distributions are defined in a recursive way, the following theorem gives a useful closed form representation.
Theorem 3 (Lemma 2 and Remark 3 in Arab and Oliveira )
Assume is an absolutely continuous nonnegative random variable with finite moment of order for . Then the iterated tail may be represented as
where is the residual lifetime at age .
It follows from (3) that the iterated distribution may be interpreted as a normalized survival moment of order . It is also worth mentioning that (3) means that, up to a normalizing factor, the iterated distribution is the stop-loss transform of order , , a common quantity of interest in actuarial models (see, for example, Cheng and Pai , Nair et al. , Tsai  or Rachev and Rüschendorf , among many other references).
Although Theorem 3 provides a closed representation for , it does not seem to be very helpful for actual calculations. An illustrative example, that we will be exploring in detail later, is obtained by assuming that has distribution in the Gamma or Weibull families, which are important classes of distributions in many different research fields such as reliability theory.
The paper is organized as follows: In Section 2, we provide a formula for the higher order moments of the iterated distributions. In Section 3, a recursive representation for the high order iterated distributions of a convolution is obtained and this result is used to find an expression for the iterated distributions of a random variable that is Gamma distributed with integer shape parameter. In Section 4, we derive an explicit expression for the iterated distribution of a Gamma distributed random variable, still assuming the shape parameter is an integer and furthermore, we study the limit behavior of the iterated distribution as the iteration step tends to infinity, now for general shape parameter. Finally in Section 5, we identify the limit behavior for the iterated distributions induced from the Weibull family of distributions.
2 Moments of iterated distributions
We shall start by a characterization of the moments of iterated distributions. First, remark that the normalizing constants in (2) are obviously the mathematical expectation of the iterated distributions. These first order moments have been characterized, in terms of moments of the inducing random variable, in Corollary 2.1 by Nanda et al. :
The following result extends (4), characterizing higher order moments of the iterated distributions.
Assume is an absolutely continuous nonnegative random variable with finite moment of order , where and . Then the iterated distribution induced by has a finite moment of order given by
By inverting the integration order, we easily find that
where is the Euler function, which leads to the desired result.
The following description for the variance of iterated distributions is now straightforward.
Assume is an absolutely continuous nonnegative random variable with finite moment of order . The variance of the iterated distribution induced by X is
As mentioned earlier, the exponential distributions are fixed points for the iteration procedure, so one could expect to have iterated distributions converging to the exponential. The above corollary shows that this may not happen. Indeed, it follows immediately that if the quotient is bounded with respect to , then , implying that the equilibrium distributions are, in such cases, converging to a degenerate distribution. We will prove later that the Weibull distributions satisfy this asymptotic degeneracy.
3 Iterated distribution of convolutions
Summing independent random variables is a common way to introduce new families of distributions, appearing as convolution powers built upon some initial distribution. With this in mind, we shall derive a characterization for the iterated distribution for the -th convolution power based on expressions for the -th convolution power. For this purpose, it is convenient to describe the distribution in terms of their densities rather than using distribution functions. Of course, as follows from Definition 1, the density for the iterated distribution is, up to multiplication by a constant, the tail for the iterated distribution:
In the sequel, we will be representing the convolution of density functions by defined as follows:
Moreover, the expression will represent the -th convolution power of a density function .
Let be nonnegative i.i.d random variables with the same distribution as and define . For every , the density of the iterated distribution induced by is given by
where and .
Proof. Using (5) for the -th convolution power and reversing the integration order, we find that
We rewrite as
As what regards , we have that
The conclusion (6) now follows by combining (7), (8) and (9).
The result above provides a recursive formula for the high order iterated distributions of the convolution . A simple application is given next, identifying explicitly the iterated distributions induced by integer shape parameter Gamma distributions.
The is the 2-nd convolution power of the exponential distribution with hazard rate , whose density function is , for . With the notation introduced above, the density function of the distribution is represented as , hence the iterated distribution induced by the distribution has density , and we may use (6) to obtain a recursive representation.
As the exponential is a fixed point for the iterative procedure introduced in Definition 1, we have , therefore is the moment of order of an exponential random variable with hazard rate , that is, . Moreover, as the 2-nd convolution power is the , we also know that . Furthermore,
Hence, the iterated distribution induced by the distribution has density
This approach may be extended to general Gamma distributions with integer valued shape parameter. The argument goes along the same line, needing some extra effort on the manipulation of some combinatorial sums. We first state an auxiliary combinatorial lemma that is useful for obtaining the results that follow.
For every , we have that
Proof. This follows easily by induction on .
Let be a random variable with distribution , and be the density of an exponential distribution with hazard rate . For every , the density function of the iterated distribution induced by is
Proof. As , we will use the recursive representation from Theorem 7. Note that the moments of convolutions are just the moments of order of the distribution. Hence
The representation just proved in Proposition 10 allows for a first characterization of the behavior of the iterated distributions as the iteration step goes to .
Let be a random variable with distribution . Then, for each fixed, , where is the density of an exponential distribution with hazard rate .
Proof. The proof follows easily by observing that the convolution that appears in the first term on the right hand side of (11) is bounded by . Allowing to tend to we have the desired result.
One could prefer to have a recursive characterization for the density of the iterated distribution involving only elementary operations, that is, avoiding convolutions on the recursive expression. This can be obtained by iteratively using (6) to describe the -th convolution power that appears in the right-hand side of (6). We apply this technique for random variables that are Gamma distributed and this leads to a long algebraic manipulation proving the representation given below. Note that the iterated distributions induced by the Gamma family will be studied with a more explicit approach in the next section.
Let be a random variable with distribution , with , and let . Then,
We may rewrite this for the -fold convolution, to get
The first term on the right in (15) may be written as , so it can be replaced by the previous expression to find, after a rearrangement of the terms:
Again, we apply (15) to the first term on the right-hand side above. So, iterating this substitution, we finally get the representation:
We may, of course, rewrite the representation above for -fold convolution:
from which follows that
Since is the exponential density, the iterated density , hence, we may rewrite:
We may now compute each of these three terms. Recalling that the convolution of exponentials is a gamma density, and the expressions for the moments, it follows easily that
So, we need to compute two summations:
and, using Lemma 9,
We now compute
Note first that , so the summation ranges only from up to . Further, for ,
Replacing these expressions and inverting the summations, we find
using Lemma 9 for the final equality.
4 Iterated distributions induced by Gamma distributions
Given a random variable , the iterated distributions induced by , where , are easily related to the iterated distributions induced by . Indeed, it follows immediately from (3) that, for every and , . Therefore, in order to characterize the iterated distributions induced by a Gamma distributed random variable, it is enough to treat the case where has distribution . We first derive an explicit expression for the iterated distribution when the shape parameter is an integer.
Assume is distributed with integer shape parameter . For every , the iterated distribution induced by has tail given by
Proof. We start by calculating the stop-loss transform of order for the random variable .
where the last equality follows by applying the binomial expansion. The result follows by writing