Conditions for finiteness and bounds on moments of record values from iid continuous life distributions

We consider the standard and kth record values arising in sequences of independent identically distributed continuous and positive random variables with finite expectations. We determine necessary and sufficient conditions on the type of record k, its number n and moment order r so that the rth moment of the n value of kth record is finite for every parent distribution. Under the conditions we present the optimal upper bounds on these moments expressed in the scale units being the respective powers of the first population moment. The theoretical results are illustrated by some numerical evaluations.


Introduction
The notion of records in probabilistic setup was defined by Chandler (1952). Given a sequence of random variables X 1 , X 2 , … , we say that the upper record value occurs in the sequence at time j and becomes equal to X j if X j > X i , i < j . Using the notion of order statistics X 1∶j ≤ … ≤ X j∶j arising by ordering X 1 , … , X j in the ascending order, we obtain the record at time j if X j∶j > X j−1∶j−1 . Dziubdziela and Kopociński (1976) generalized the notion introducing so called kth records which appear when the new value of kth maximum becomes greater than the previous one. Formally, for n ≥ 1 the sequences of consecutive kth record times T (k) n and kth record values R (k) n are defined as follows For k = 1 the above definitions describe the classic record times and values of Chandler (1952).
If the joint distribution of X 1 , X 2 , … does not admit the repetitions almost surely (e.g., X i are iid with a continuous parent distribution), then the sequences of kth records are infinite. The kth records with k ≥ 2 , although seemingly less natural and intuitive than the first ones, attract increasing interest of statistical investigators of extreme events. The main reason is that they occur much more often. Indeed, if a new first record appears in the sequence, the all the kth records change: the old kth records become new (k + 1) th records. If a new observation becomes the th maximum, then we get new kth records for k ≥ , but the other kth records k < remain unchanged.
The theory of records is best developed in the model of iid continuously distributed random variables. In particular, the distribution function of the nth value of kth record is the composition G (k) n • F , where and F is the parent distribution function of the original observations. Moreover, if F has a density function f with respect to the Lebesgue measure, then the density function of R (k) n exists and has the form g (k) n • F ⋅ f , where The above establishments can be found, for example, in the monographs by Arnold et al (1998) and Nevzorov (2001), devoted to the record value theory.
In the paper we consider the classic model of iid random variables X 1 , X 2 , … with a continuous distribution function F. We additionally assume that X 1 is positive and has a finite expectation > 0 . Our main results are following. For every type k ∈ ℕ and number n ∈ ℕ of the record value R (k) n we determine the necessary and sufficient conditions on the raw moment order r > 0 such that R (k) n r < ∞ for all parent distribution functions F satisfying the assumptions. Our conditions are following: r ≤ k if n = 1 , and r < k if n ≥ 2 . Moreover, under these conditions we establish the sharp upper bounds on R (k) n r expressed in the scale units being the respective powers r of the mean of the parent distribution which are valid for all continuous F supported on ℝ + . The bounds are presented in Sect. 2. In particular, these results provide the sufficiency proof of our conditions for finiteness of moments of record values. The necessity of the conditions are proven in Sect. 3. In Sect. 4 we present and briefly discuss some numerical examples of the bounds determined in Sect. 2. Some conclusions of our results are presented in Sect. 5. The problem of existence of kth record values has been hitherto solved only for some special cases of k and n. For k = n = 1 it is trivial because R (1) 1 = X 1 . For k = 1 < n the solution was presented by Nagaraja (1978), Sect. 2). It follows that both R (1) n r < ∞ for all baseline distributions with finite mean iff r < 1 . For n = 1 < k , the conclusion can be deduced from Sen (1959) who considered order statistics. This was explicitly stated by Papadatos (2021) Here we focus on the remaining cases k, n ≥ 2 . To the best of our knowledge, the most general sufficient conditions were presented by Cramer et al. (2002), Theorem 2.4 and the following Remark (v). They proved that finiteness of X 1 assures that R (k) n r < ∞ if r < k and r ≤ n.
The first result concerning the bounds on the expectations of record values was presented by Nagaraja (1978) who determined the sharp evaluations of the expectations of the standard record values from the population means expressed in the population standard deviation units. Raqab (1997) derived analogous bounds for the kth records. These results were generalized by Raqab (2000) and Raqab and Rychlik (2002), respectively, by considering more general scale units based on various central absolute moments of the original variables. Klimczak (2007) considered kth records from bounded populations expressing the expectation bounds in terms of the length of the support interval. There are known some optimal mean-variance bounds on the expectations of record values from restricted families of distributions. The most general results are due to Bieniek (2007) and Goroncy (2019) who considered the families with monotone generalized failure rates. We also mention the bounds on differences of record values studied in Danielak (2005), and Danielak and Raqab (2004). Evaluations of record raw moments different from the first ones have not been presented by now in the literature. We can only mention analogous results for the order statistics determined in Papadatos (2021).

Bounds
In this section we present the sharp bounds on the ratios R (k) n r r for r < k ≥ 1 and n ≥ 2 . We exclude from our considerations the first values of kth records, because these are minima of the first k observations. The bounds for the moments of order statistics, and sample minima in particular were determined by Papadatos (2021). He proved that For r < k the bound is attained by any degenerate parent distribution (see Papadatos (2021), Remark 4). If r = k , we get the equality in Eq. (3) iff the parent distribution is supported on two points 0 and p with probabilities 1 − p and p for arbitrary 0 < p < 1 (see Papadatos (2021), Corollary 3). A classic result by Sen (1959) says that if r > k , then there exist parent life distributions such that for iid X 1 , … , X k we have X 1 < ∞ and X r 1∶k = +∞ . Since in our model continuity of random variables is desired, we can replace X r 1∶k by [R (k) 1 ] r in the above relations, but the attainability conditions should be modified. The bounds are attained in the limit by sequences of continuous parent distributions converging weakly to the corresponding optimal discrete ones.
In the sequel, we use the following lemma.
Lemma 1 Fix k, n ≥ 2 and 1 ≤ r < k . Then the function is maximized by the unique point 0 < = (k) n (r) < 1 − exp − n−1 r+k−2 that is the unique solution to the equation We denote the maximal value of function Eq. (4) by The middle formula follows from Eq. (5). The last one provides the explicit form of it.
The proof of Lemma 1 is postponed to the Appendix. We first take into account the moments of record values of orders 0 < r < 1 for k, n ≥ 2 . To this end we use a simplified version of Theorem 1 in Moriguti (1953) (see Lemma 2) and the Hölder inequality (see Lemma 3).

Lemma 2 Suppose that a real function g defined on [a, b] has a finite integral. Let g denote (the right-continuous version of, say) the derivative of the greatest convex
We use the above construction for the density functions Eq. (2). For k, n ≥ 2 the derivative of the greatest convex minorant of the antiderivative Eq. (1) of Eq. (2) has the form n (1) is the unique solution to Eq. (5) with r = 1 . The greatest convex minorant is equal to the antiderivative G (k) n (see Eq. 1) on the interval (0, (k) n (1)) and less than G (k) n on ( (k) n (1), 1) . These facts were established and applied in a number of papers, see, e.g., Raqab (1997) and Raqab and Rychlik (2002). The Hölder inequality can be found, e.g., in Mitrinović (1970). Theorem 1 Let X 1 , X 2 , … , be an infinite sequence of iid continuous positive random variables with a mean 0 < < ∞ . Then for k, n ≥ 2 and 0 < r < 1 we have

Lemma 3 Let g and h be non-negative non-zero elements of the Banach spaces
The bound in Eq. (8) is attained in the limit by continuous parent distributions tending weakly to the distribution function n −1 denotes here the inverse of the increasing part of Eq. 2) with the atom of size 1 − (k) n (1) at the right end-point of the support, and absolutely continuous part between 0 and the atom.
The proof of Theorem 1 can be found in the Appendix. Now we consider the cases k = 1 with n ≥ 2 and 0 < r < 1 . The idea of getting the bounds is analogous to that of the previous theorem. However, for k = 1 we can obtain explicit formulae and we decided to present them below. The reason is that in this case the method presented in Lemma 1 generates g (1) (2) for k = 1 is then increasing, its antiderivative is convex, and so it coincides with its greatest convex minorant. This substantially simplifies the evaluations.
Theorem 2 For X 1 , X 2 , … being iid random variables with an expectation 0 < < ∞ , and the respective classic records with k = 1 , n ≥ 2 and 0 < r < 1 we get We postpone the proof of Theorem 2 to the Appendix. If n−r 1−r is an integer number m, say, i.e. r = i+1 n+i for some i = 0, 1, … , then the right-hand side of Eq. (11) has yet nicer form [(m−1)!] n−1 m−1 (n−1)! . We now concentrate on evaluations of rth moments of kth record values for r ≥ 1 . At the beginning, we exclude from our investigations the classic moments R (1) n , n ≥ 2 , because there exist parent distributions such that X 1 = < ∞ and R (1) 2 = +∞ (see Nagaraja (1978)). Otherwise we use the following lemma which was presented in Papadatos (2021), Corollary 4.

and the equality is attained if F is a two-point distribution function supported on 0 and a positive number.
Observe that for r = 1 relation Eq. (12) becomes a trivial equality with no assumptions on F.
Theorem 3 Let X 1 , X 2 , … satisfy the assumptions of Theorem 1. For n, k ≥ 2 and 1 ≤ r < k we have where with (k) n (r) determined by the Eq. (5).

The equality in Eq. (13) is attained in the limit by continuous parent distributions tending weakly to the two-point distribution
The proof of Theorem 3 is presented in the Appendix.

Conditions for moment finiteness
Here we present the necessary and sufficient conditions for finiteness of rth moments of nth values of kth records for arbitrary continuous life distributions of the baseline sequence with a finite mean. As we mentioned in the introduction, the conditions are known in some special cases. If n = 1 then R (k) 1 = X 1∶k , and the result is immediately concluded from the classic paper by Sen (1959) on order statistics (see also Papadatos (2021)). If r ≤ k , then for every parent distribution with X 1 < ∞ we have [R (k) 1 ] r < ∞ as well. For r > k there exist parent distribution functions such that X 1 < ∞ and [R (k) 1 ] r = ∞ . An example of such distribution function is (see Papadatos (2021), Remark 3). If n ≥ 2 and k = 1 , the respective condition r < 1 was established in Nagaraja (1978). In Lemma 2.1 of the paper he proved that finiteness of X 1 and [R (1) n ] r < ∞ for all n ≥ 2 and 0 < r < 1 . On the other hand, for the distribution function Eq. (14) we have X 1 = 2e < ∞ and R (1) n = ∞ for all n ≥ 2.
x > e, In Theorem 4 below we treat all the remaining cases with n, k ≥ 2 . Cramer et al. (2002) proved that finiteness of X 1 = assures that R (k) n r < ∞ if r < k and r ≤ n . Our conditions are essentially weaker.
Theorem 4 Let X 1 , X 2 , … be positive iid random variables with a continuous marginal distribution and finite mean. The necessary and sufficient condition for finiteness of R (k) n r for all parent distributions and some n, k ≥ 2 is r < k.
We present the proof of Theorem 4 in the Appendix. It consists in delivering an example of a single absolutely continuous life marginal distribution of the elements of the sequence X 1 , X 2 , … , such that X 1 < +∞ and R (k) n r = +∞ for all k ≥ 2 , n ≥ 2 , and r ≥ k . We finally present a simple generalization of the results of this section.
Corollary 1 Consider a sequence of positive and continuous iid random variables X 1 , X 2 , … . Then for arbitrary parent distribution with finite pth moment we have R (k) n r < ∞ iff either n = 1 and r ≤ pk or n ≥ 2 and r < pk.

Numerical results
Here we illustrate the results of Sect. 2 by some numerical examples. In Tables 1, 2 and 3 we present the bounds on the ratios of rth moments of kth records and the rth powers of the population means for some values of the first, second and third records, respectively. In all the cases we consider the second, third, fifth and tenth values of kth records. Table 1 is based on Theorem 2, and contains the evaluations of moments of orders r = 0.25 , 0.5, and 0.9. In Tables 2 and 3 we present the bounds for the moments of orders r = 0.5 i , i = 1, … , 2k − 1, and r = k − 0.1 , where k = 2 and 3, respectively. The results of the first columns of Tables 2 and  3 are concluded from Theorem 1, and the other ones are calculated with use of Theorem 3. All the bounds are greater than 1, which is a consequence of the sharp bounds

3
Conditions for finiteness and bounds on moments of record values… proved by Papadatos (2021). We also note that the bounds are increasing in n and decreasing in k, which is also obvious, because we have the relations R (k) n < R (k) n+1 and R (k+1) n > R (k) n following from the definitions. Moreover, they increase with respect to r and tend to infinity as r approaches k. In the case k = 1 we can apply Eq. (11) with the factorial approximation of the Gamma function and Stirling formula in order to obtain The right-hand side expression converges to 1 as n increases. We cannot provide such approximations for with r → k ≥ 2 , because (k) n (r) do not have analytic representations, and the bounds, dependent on them, do not have explicit formulae.  Anyway, we can see that for k = 2, 3 and r = k − 0.1 the bounds B (k) 10 (r) are greater than 10 9 .

Conclusions
We consider the moments of kth upper record values in the classic model of sequences of independent and identically continuously distributed positive random variables. Our purpose is to determine the necessary and sufficient conditions for the moment orders r > 0 such that the rth moment R (k) n r of the nth value, n ≥ 1 , of kth record, k ≥ 1 , is finite for all parent distributions with a finite mean. Since the solution is known in the literature for the particular cases k = 1 ≤ n and n = 1 ≤ k , we focus on the remaining cases k, n ≥ 2 proving that the necessary and sufficient conditions are r < k for arbitrary n ≥ 2 . The necessity proof consists in constructing the parent distribution with a finite expectation such that R (k) n r = +∞ for all n ≥ 2 and r ≥ k ≥ 2 . Instead of the sufficiency proof, we provide a stronger result: for every n ≥ 2 and r < k ≥ 1 we determine the sharp upper bounds on R (k) n r over all the parent distributions with finite means, expressed in the scale units being the rth powers of the mean of the single observation. Numerical exemplary evaluations show that the bounds are extremely large when the moment orders r are close to the border value k even for moderate n.
Our findings allow the researchers not to bother about existence of moments of kth records of orders less than pk if they consider sequences of random variables with arbitrary parent life distributions with finite pth moments. We hope that the tools presented in the paper shall be useful in determining the necessary and sufficient conditions for finiteness of moments in other models of ordered random variables, e.g., for progressively type II censored order statistics and generalized order statistics.

Appendix
The Appendix contains the proofs of Lemma 1 and Theorems 1-4.
Proof of Lemma 1. The function Eq. (4) is positive and satisfies h(0) = 1. By the de l'Hospital rule and The numerator N(x), say, has the derivative

3
Conditions for finiteness and bounds on moments of record values… Since we have which means that it is negative for 0 < x < = (k) n (r) = 1 − exp − n−1 r+k−2 < 1 , and positive for < x < 1 . It follows that N(x) is first decreasing and then increasing. Since N(0) = r ≥ 1 , and N(1) = 0 , the function is positive on some interval [0, ) , = (k) n (r) < (k) n (r) , and negative elsewhere. The same concerns Eq. (16), because its denominator is always positive. Since h(0) = 1 and Eq. (15) holds, we conclude that Eq. (4) is increasing on [0, ) , and decreasing on ( , 1) . Obviously its unique maximum = (k) n (r) is the only point that satisfies Eq. (5) (cf Eq. 16). It is also less than (k) n (r) = 1 − exp − n−1 r+k−2 . ▪ Proof of Theorem 1. We have The first inequality holds by Lemma 2 due to the facts that [F −1 (x)] r is nondecreasing, and Eq. (7) holds. The latter the an application of the Hölder inequality with parameters p = 1 r > 1 and q = p p−1 = 1 1−r > 1 . Note that changing the variables twice y = − ln(1 − x) and z = k−r 1−r y , we obtain Conditions for finiteness and bounds on moments of record values… By Lemma 1, the equality in the former inequality occurs if the only value of F on (0, ∞) different from 1 is (k) n (r) . By Lemma 4, for r > 1 , we get the equality in the latter inequality if F is supported on 0 and a positive number. Combining the conditions, we deduce that the only distribution that attains the equality in Eq. (13) is one that assigns probability (k) n (r) to 0 and 1 − (k) n (r) to a positive number. This distribution satisfies the moment condition if this positive number amounts to . This ends the proof. ▪ Proof of Theorem 4. The sufficiency of the condition is deduced from Theorem 3. To verify its necessity we should present examples of the parent distributions with finite means such that R (k) n r = ∞ for given n ≥ 2 and r ≥ k ≥ 2 . In fact, it suffices to construct distributions such that R (k) 2 k = ∞ for k ≥ 2 . Our construction works for all k ≥ 2.
We take a sequence of positive numbers e j , j = 0, 1, … , tending to ∞ by the recursive relation and a family of disjoint intervals Consider a mixture of distributions ∑ ∞ j=0 j U j , where U j is the uniform distribution over the interval I j , j = 0, 1, … , and the weights are defined as The mixture is a proper probability measure, because e 0 = 1, e j = e e j−1 , j = 1, 2, … , I j = e j , e j + 1 , j = 0, 1, … .
The measure is obviously absolutely continuous with respect to the Lebesgue measure, and its density function has the form Suppose that the sequence of iid life random variables X 1 , X 2 , … has the common density function Eq. (22). It has a finite expectation For proving that R (k) 2 k = ∞ we need to calculate for j = 1, 2, … . Remind that the corresponding distribution function of the second value of kth record is The distribution function of Eq. (23) satisfies (23) , j = 1, 2, … we have and We are now in a position to prove that It suffices to show that the summands do not tend to 0 as j → ∞ . We have Note that for all p ∈ ℝ , because Therefore 1 − F e j k = 1 j 2k (e j − e j−1 ) k , − ln 1 − F e j = 2 ln j + ln(e j − e j−1 ), = e k j [1 + 2k ln j + k ln(e j − e j−1 )] j 2k (e j − e j−1 ) k − e k j [1 + 2k ln(j + 1) + k ln(e j+1 − e j )] (j + 1) 2k (e j+1 − e j ) k .  x (see Eq. 22) for arbitrary > 0 we obtain X 1 = and R (k) 2 k = ∞ .

3
Conditions for finiteness and bounds on moments of record values… The integer part of e 5 has approximately 10 1.012558⋅10 1656520 digits. It follows that the support intervals I j , j = 0, 1, … , are very dispersed on the positive half-axis, and the corresponding distribution function has a very heavy tail. One could decrease the dispersion by replacing e in the construction of above proof by some 1 < a < e , especially one close to 1 but this would generate more sophisticated formulae. Anyway, for any a > 1 the tails are heavier than these of the distribution functions with ln m x denoting m-fold composition of logarithm functions and any integer m ≥ 2 and > 0 . It can be checked that these distribution functions have finite expectations and satisfy m, R (k) n k < ∞ iff n ≤ k and m, R (k) n r = ∞ for all r > k with n ≥ 2 and r=k with n ≥ k + 1.