Journal of Mathematical Sciences, Vol. 121, No. 5, 2004
A NOTE ON BLACKWELL’S RENEWAL THEOREM
E. A. M. Omey (Brussels, Belgium) and J. L. Teugels (Heverlee, Belgium)
UDC 519.2
1. Introduction Our point of departure is a renewal process for which we define a number of crucial quantities. All of these concepts can be found in [1, 5, 7, 10]. Definition 1. Let {Xi , i ∈ N} be a sequence of independent identically distributed random variables with common distribution F of X, where X > 0. The sequence {Xi , i ∈ N} is called a renewal process. Let S0 = 0 and Sn = Xn + Sn−1 , n ≥ 1. The sequence {Sn , n ∈ N} constitutes the set of renewal (time) points. Let also t ≥ 0, and define N (t) = sup[n: Sn ≤ t]; then the processes {N (t), t ≥ 0} and N0 (t) := 1 + N (t) are called the renewal counting processes. The renewal functions are defined and denoted by U (t) =: EN (t) and U0 (t) =: I(t) + U (t). Finally, we write µ := EX ≤ ∞. We call the processes of the associated functions generated by F or by X. Definition 2. Let {Si , i ∈ N} and {N (t), t ≥ 0} be as defined above. For every t ≥ 0 we define (i) the age of the renewal process by Y (t) = t − SN(t) ; (ii) the residual life (time) of the renewal process by Z(t) = SN(t)+1 − t = SN0 (t) − t. We discuss the precise relationship between the asymptotic behavior of the age or residual life and that of the difference process D(t, y) := N (t + y) − N (t) = N0 (t + y) − N0 (t). For µ < ∞ we show that while t ↑ ∞ Blackwell’s theorem ED(t, y) → y/µ is equivalent to the convergence in distribution of Y (t) or of Z(t) to a nondegenerate limit. We also give a result for the case where µ = ∞. The technical link between age and residual life time is expressed by the equation (t ≥ y) t−y Z (1 − F (t + z − u)) dU0 (u), P{Y (t) ≥ y, Z(t) ≥ z} =
(1)
0
which implies in particular that P{Y (t) ≥ y, Z(t) ≥ z} = P{Y (t + z) ≥ y + z} = P{Z(t − y) ≥ y + z}.
(2)
2. The Finite Mean Case D
From (1) and the key renewal theorem one infers that for F nonlattice and µ ≤ ∞, Z(t) −→ Z ∗ . More explicitly, for all u ≥ 0, Zu 1 ∗ (1 − F (v)) dv = µ−1 m(u), (3) lim P{Z(t) ≤ u} = P{Z ≤ u} = t↑∞ µ 0
where
Zx (1 − F (y)) dy
m(x) := 0
is the integrated tail of F . D D From (2) it follows that (3) holds if and only if Y (t) −→ Y ∗ (or even (Y (t), Z(t)) −→ (Y ∗ , Z ∗ )), where P{Y ∗ ≥ y, Z ∗ ≥ z} = P{Y ∗ ≥ y + z} = P{Z ∗ ≥ y + z}. For the finite mean case we have the following result. THEOREM 1. Suppose F is nonlattice and has finite mean µ. Then the following statements are equivalent: (2.1) lim {U (t + y) − U (t)} = y/µ; t↑∞
Proceedings of the Seminar on Stability Problems for Stochastic Models, Varna, Bulgaria, 2002, Part I. 2688
c 2004 Plenum Publishing Corporation 1072-3374/04/1215-2688
D
(2.2) Z(t) −→ Z with Z nondegenerate. Either statement implies that Z and Z ∗ have the same distribution. Proof. Fix y ≥ 0 during the following proof. Together with the original renewal process we define a shadow renewal process {Xn0 , n ≥ 1} generated by X 0 with the same distribution F as the original process. The renewal quantities Sn0 , N 0 (t), N00 (t), U 0 (t), U00 (t), Y 0 (t), and Z 0 (t) are defined as for the first process. We first prove two auxiliary results. LEMMA 1. Difference processes D(t, y) and N00 (y − Z(t)) have the same distribution. Proof. We need to show that for all integral values of k ≥ 0 P{D(t, y) = k} = P{N00 (y − Z(t)) = k}. For convenience, write q(n, k) = P{N0 (t) = n, D(t, y) = k}. For k = 0 q(n, 0) = P{N0 (t) = n, y − Z(t) < 0}. Further, q(n, 1) = P{N0 (t) = n, Sn ≤ t + y < Sn+1 } = P{N0 (t) = n, 0 ≤ y − Z(t) < Xn+1 }. Since the event {N (t) = n, Z(t) ≤ z} is independent of Xn+1 , we have q(n, 1) = P{N (t) = n, 0 ≤ y − Z(t) < X10 } = P{N0 (t) = n, N00 (y − Z(t)) = 1}. For k ≥ 2 we have similarly q(n, k) = P{N0 (t) = n, Sn+k−1 ≤ t + y < Sn+k } 0 ≤ y − Z(t) < Sk0 } = P{N0 (t) = n, N00 (y − Z(t)) = k}. = P{N0 (t) = n, Sk−1 Summing over n gives the required result. D
LEMMA 2. If F is nonlattice and µ < ∞, then D(t, y) −→ N00 (y − Z ∗ ). Proof. As in the previous proof, first take k = 0. Then, by the previous lemma, P{D(t, y) = 0} = P{N00 (y − Z(t)) = 0} = P{y − Z(t) < 0} −→ P{y − Z ∗ ≤ 0} = P{N00 (y − Z ∗ ) = 0}, where we used (3). For k > 0 we show that
P{D(t, y) = k} −→ P{N00 (y − Z ∗ ) = k}.
Using the lemma again and (3), we have for k > 0 P{D(t, y) > k} = P{N00 (y − Z(t)) > k} = P{Z(t) ≤ y − Sk0 } −→ P{Z ∗ ≤ y − Sk0 } = P{N00 (y − Z ∗ ) > k}. This proves the lemma. Proof of Theorem 1. We only need to prove that (2.2) implies (2.1), since the other implication follows from (3). Now, Lemma 1 implies the following explicit expression: Zy P{Z(t) ≤ y − v} dU0 (v) =: Gt ∗ U0 (y),
ED(t, y) =
(4)
0
where we use the notation for a convolution and the abbreviation Gt (z) = P{Z(t) ≤ z}.
(5)
On the right-hand side, we can take the limit for t ↑ ∞. Indeed, Gt (z) → G(z) := P{Z ≤ z} by assumption and hence dominated convergence applies. Therefore, lim ED(t, y) =: k(y) t↑∞
2689
exists. Obviously, 0 ≤ k(y) is nondecreasing and satisfies the functional equation k(u + v) = k(u) + k(v). Hence, k(u) = βu for some nonnegative quantity β. That β = µ−1 follows from the elementary renewal theorem and Ces`aro’s ˆ ˆ Uˆ0 (s) for the Laplace–Stieltjes transforms. Since Uˆ0 (s) = lemma. Alternatively, k(y) = G ∗ U (y) leads to k(s) = G(s) −1 ˆ ˆ ˆ ˆ − Fˆ (s)}−1 , G(s) = β(1 − Fˆ (s))/s. {1 − F (s)} with F , the transform of the generating distribution β/s = G(s){1 −1 ∗ ˆ Since Z is nondegenerate, β > 0. Furthermore, since G(0) = 1, β = µ . That Z and Z have the same distribution is obvious. While Blackwell’s theorem shows that lim ED(t, y) = y/µ, Lemma 2 shows that D(t, y) itself converges weakly to N00 (y − Z ∗ ). 3. The Infinite Mean Case Since the derivation of Dynkin [3] and Lamperti [8] (see also [2]), it is traditional to assume that the renewal process is generated by a distribution with regularly varying tail, i.e., for some slowly varying function `(x) and an index of regular variation α ∈ (0, 1) 1 − F (x) ∼ x−α `(x) for x ↑ ∞. A slightly more appropriate formulation goes in terms of the integrated tail m(x). Denote the class of functions of regular variation with index β by Rβ . Let Z ∗ (α) denote a positive random variable with density pα (x) = sin(πα)(πxα (1 + x))−1 on x > 0. Then the following result is almost classical and of elementary renewal type. LEMMA 3. Let `(t) be slowly varying and 0 < α < 1. Then as t ↑ ∞ the following conditions are equivalent: (i) m(t) ∈ R1−α ; (ii) U (t) ∈ Rα ; D (iii) Z(t)/t −→ Zα∗ . The proof can be found by a combination of results in [2, Sec. 8.6.2]. A more general class than the class of regularly varying functions is the class L of positive, measurable functions a on <+ for which limt↑∞ a(t + y)/a(t) = 1 locally uniformly in y. We link the quantity D(t, y) to Z(t) through the following result. THEOREM 2. Suppose a(t) ∈ L. Then the following statements are equivalent: (3.1) lim a(t)(U (t + y) − U (t)) = k(y), y > 0; t↑∞
(3.2) lim a(t)P{Z(t) ≤ y} = v(y), t↑∞
y > 0.
Moreover, if (3.2) holds, then k(y) = βy for some β ≥ 0 and v(y) = βm(y). Proof. That (3.2) implies (3.1) follows from (4). Indeed, Zy a(t)P{Z(t) ≤ y − v} dU0 (v)
a(t)E(D(t, y)) = 0
and (3.1) follows from Lebesgue’s theorem. The limit equals k(y) = v ∗ U0 (y). Before embarking on the converse, note that the existence of the limit in (3.1) implies its specific form. By the L-condition on a, one immediately sees that the limit k satisfies the functional equation k(u + v) = k(u) + k(v). Since, furthermore, k(u) ≥ 0, k(y) = βy for some β ≥ 0. For the converse we start from the distribution of Z(t). Take y ≤ t. It is not difficult to derive from (1) that Zy P{Z(t) ≤ y} = (1 − F (y))(U0 (t) − U0 (t − y)) +
(U (t) − U (t − u)) dF (u). 0
But then (3.1) easily implies that a(t)P(Z(t) ≤ y) converges. The limit v(y) is given by Zy v(y) = (1 − F (y))k(y) +
Zy (1 − F (u)) dk(u).
k(u) dF (u) = 0
0
Since k(u) = βu, it follows that v(y) = βm(y). We make a number of remarks with respect to the last result. The L-condition is not used for the implication (3.2) to (3.1). Further, k(y) and v(y) can be replaced by O(1). For general results in such a direction, see [9]. 2690
An immediate application of the above theorem results from a Blackwell-type theorem on [4], where the function a(x) = m(x) is regularly varying with index α ∈ (1/2, 1]. Of course, when a(x) is constant, then we recover part of Theorem 1. 4. The Lifespan of a Renewal Process It would be an interesting adaptation of the results if we could prove the following results. Let L(t) := Y (t) + Z(t) be the lifespan of a renewal process. Alternatively, L(t) = SN0 (t) − SN(t) = XN0 (t). It is clear that in the finite mean Ry D case the conditions of Theorem 1 imply that also L(t) −→ L∗ , where P(L∗ ≤ y) = µ−1 0 v dF (v). D
The reciprocal implication that L(t) −→ L∗ (not constant) implies Blackwell’s result is less obvious. We had no success in trying to prove this. Relations that might help and that link L(t) with the quantity D(t, y) are only in integral form and go the wrong way. For example, the following facts are known: E{exp −θL(t)} = Kθ ∗ U0 (t), where
Z∞ Kθ (x) =
e−θv dF (v).
x
This, together with the key renewal theorem, also proves the result on L∗ . A further transform with respect to t is possible, but then it does not seem obvious how to use the fact that t ↑ ∞. For results on L(t), see [6, 11]. REFERENCES 1. G. Alsmeyer, Erneuerungstheorie, Teubner, Stuttgart (1991). 2. N. H. Bingham, C. M. Goldie, and J. L. Teugels, Regular Variation, Cambridge Univ. Press, Cambridge (1987). 3. E. B. Dynkin, “Some limit theorems for sums of independent random variables with infinite mathematical expectations,” Select. Trans. Math. Statist. Probab., 1, 171–189 (1961). 4. K. B. Erickson, “Strong renewal theorems with infinite mean,” Trans. Amer. Math. Soc., 185, 371–381 (1970). 5. W. Feller, An Introduction to Probability Theory and Its Applications, Vol. 2, Wiley, New York (1971). 6. K. Hinderer, “A unifying method for some computations in renewal theory,” Z. angew. Math. Mech., 65, 199–206 (1985). 7. S. Karlin and H. M. Taylor, A First Course in Stochastic Processes, Academic Press, New York (1975). 8. J. Lamperti, “An invariance principle in renewal theory,” Ann. Math. Statist., 33, 685–696 (1962). 9. E. Omey, “On a subclass of regularly varying functions,” J. Statist. Planning Inf., 45, 275–290 (1995). 10. S. M. Ross, Stochastic Processes, Wiley, New York (1983). 11. J. L. Teugels, “The lifespan of a renewal,” in: Stochastic Musings, Festschrift for C. Kevork and P. Tzortopoulos, Institute of Statistical Documentation Research and Analysis, Athens, Greece (2000) (to appear).
EHSAL, Stormstraat, 2, 1000 Brussels, Belgium e-mail:
[email protected] University Centre for Statistics, K.U.Leuven, De Croylaan, 52B, 3001 Heverlee, Belgium e-mail:
[email protected].
2691