Next Article in Journal
A Method of Riemann–Hilbert Problem for Zhang’s Conjecture 1 in a Ferromagnetic 3D Ising Model: Trivialization of Topological Structure
Next Article in Special Issue
Limit Theory for Stationary Autoregression with Heavy-Tailed Augmented GARCH Innovations
Previous Article in Journal
Analysis and Computation of Solutions for a Class of Nonlinear SBVPs Arising in Epitaxial Growth
Previous Article in Special Issue
On Small Deviation Asymptotics in the L2-Norm for Certain Gaussian Processes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Chebyshev–Edgeworth-Type Approximations for Statistics Based on Samples with Random Sizes

by
Gerd Christoph
1,*,† and
Vladimir V. Ulyanov
2,†
1
Department of Mathematics, Otto-von-Guericke University Magdeburg, 39016 Magdeburg, Germany
2
Faculty of Computational Mathematics and Cybernetics, Lomonosov Moscow State University, 119991 Moscow, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2021, 9(7), 775; https://doi.org/10.3390/math9070775
Submission received: 3 March 2021 / Revised: 28 March 2021 / Accepted: 29 March 2021 / Published: 2 April 2021
(This article belongs to the Special Issue Analytical Methods and Convergence in Probability with Applications)

Abstract

:
Second-order Chebyshev–Edgeworth expansions are derived for various statistics from samples with random sample sizes, where the asymptotic laws are scale mixtures of the standard normal or chi-square distributions with scale mixing gamma or inverse exponential distributions. A formal construction of asymptotic expansions is developed. Therefore, the results can be applied to a whole family of asymptotically normal or chi-square statistics. The random mean, the normalized Student t-distribution and the Student t-statistic under non-normality with the normal limit law are considered. With the chi-square limit distribution, Hotelling’s generalized T 0 2 statistics and scale mixture of chi-square distributions are used. We present the first Chebyshev–Edgeworth expansions for asymptotically chi-square statistics based on samples with random sample sizes. The statistics allow non-random, random, and mixed normalization factors. Depending on the type of normalization, we can find three different limit distributions for each of the statistics considered. Limit laws are Student t-, standard normal, inverse Pareto, generalized gamma, Laplace and generalized Laplace as well as weighted sums of generalized gamma distributions. The paper continues the authors’ studies on the approximation of statistics for randomly sized samples.

1. Introduction

In classical statistical inference, the number of observations is usually known. If observations are collected in a fixed time span or we lack observations the sample size may be a realization of a random variable. The number of failed devices in the warranty period, the number of new infections each week in a flu season, the number of daily customers in a supermarket or the number of traffic accidents per year are all random numbers.
Interest in studying samples with a random number of observations has grown steadily over the past few years. In medical research, the authors of [1,2,3] examines ANOVA models with unknown sample sizes for the analysis of fixed one-way effects in order to avoid false rejection. Applications of orthogonal mixed models to situations with samples of a random number of observations of a Poisson or binomial distributed random variable are presented. Based a random number of observations [4], Al-Mutairi and Raqab [5] and Barakat et al. [6] examined the mean with known and unknown variances and the variance in the normal model, confidence intervals for quantiles and prediction intervals for the future observations for generalized order statistics. An overview on statistical inference of samples with random sample sizes and some applications are given in [4], see also the references therein.
When the non-random sample size is replaced by a random variable, the asymptotic features of statistics can change radically, as shown by Gnedenko [7]. The monograph by Gnedenko and Korolev [8] deals with below limit distributions for randomly indexed sequences and their applications.
General transfer theorems for asymptotic expansions of the distribution function of statistics based on samples with non-random sample sizes to their analogues for samples of random sizes are proven in [9,10]. In these papers, rates of convergence and first-order expansion are proved for asymptotically normal statistics. The results depend on the rates of convergence with which the distributions of the normalized random sample sizes approach the corresponding limit distribution.
The difficulty of obtaining second-order expansions for the normalized random sample sizes beyond the rates of convergences was overcome by Christoph et al. [11]. Second-order expansions were proved by the authors of [11,12] for the random mean and the median of samples with random sample sizes and the authors of [13,14] for the three geometric statistics of Gaussian vectors, the length of a vector, the distance between two vectors and the angle between two vectors associated with their correlation coefficient when the dimension of the vectors is random.
The classical Chebyshev–Edgeworth expansions strongly influenced the development of asymptotic statistics. The fruitful interactions between Chebyshev Edgeworth expansions and Bootstrap methods are demonstrated in [15]. Detailed reviews of applications of Chebyshev–Edgeworth expansions in statistics were given by, e.g., Bickel [16] and Kolassa [17]. If the arithmetic mean of independent random variables is considered as the statistic, only the expected value and the dispersion are taken into account in the central limit theorem or in the Berry–Esseen inequalities. The two important characteristics of random variables, skewness and kurtosis, have great influence on second order expansions, provided that the corresponding moments exist. The Cornish–Fisher inversion of the Chebyshev–Edgeworth expansion allows the approximation of the quantiles of the test statistics used, for example, in many hypothesis tests. In [11], Theorems 3 and 6, and [12], Corollaries 6.2 and 6.3, Cornish–Fisher expansions for the random mean and median from samples with random sample sizes are obtained. In the same way, Cornish–Fisher expansions for the quantiles of the statistics considered in present paper can be derived from the corresponding Chebyshev–Edgeworth expansions.
In the present paper, we continue our research on approximations if the sample sizes are random. To the best of our knowledge, Chebyshev–Edgeworth-type expansions with asymptotically chi-square statistics have not yet been proven in the literature when the sample sizes are random.
The article is structured as follows. Section 2 describes statistical models with random numbers of observations, the assumptions about statistics and random sample sizes and transfer propositions from samples with non-random to random sample sizes. Section 3 presents statistics with non-random sample sizes with Chebyshev–Edgeworth expansions based on standard normal or chi-square distributions. Corresponding expansions of the negative binomial or discrete Pareto distributions as random sample sizes are considered in Section 4. Section 5 describes the influence of non-random, random or mixed normalization factors on the limit distributions of the examined statistics that are based on samples with random sample sizes. Besides the common Student’s t, normal and Laplace distributions, inverse Pareto, generalized gamma and generalized Laplace as well as weighted sums of generalized gamma distributions also occur as limit laws. The main results for statistic families with different normalization factors and examples are given in Section 6. To prove statements about a family of statistics, formal constructions for the expansions are worked out in Section 7, which are used in Section 8 to prove the theorems. Conclusions are drawn in Section 9. We leave four auxiliary lemmas to Appendix A.

2. Statistical Models with a Random Number of Observations

Let X 1 , X 2 , R = ( ) and N 1 , N 2 , N + = { 1 , 2 , } be random variables defined on a common probability space Ω , A , P . The random variables X 1 , , X m denote the observations and form the random sample with a non-random sample size m N + . Let
T m : = T m X 1 , , X m with m N +
be some statistic obtained from the sample { X 1 , X 2 , , X m } . Consider now the sample X 1 , , X N n . The random variable N n N + denotes the random size of the underlying sample, that is the random number of observations, depending on a parameter n N + . We suppose for each n N + that N n N + is independent of random variables X 1 , X 2 , and N n in probability as n .
Let T N n be a statistic obtained from a random sample X 1 , X 2 , , X N n defined as
T N n ( ω ) : = T N n ( ω ) X 1 ( ω ) , , X N n ( ω ) ( ω ) for   all ω Ω and   every n N + .

2.1. Assumptions on Statistics T m and Random Sample Sizes N n

In further consideration, we restrict ourselves to only those terms in the expansions that are used below.
We assume that the following condition for the statistic T m with E T m = 0 from a sample with non-random sample size m N + is fulfilled:
Assumption 1.
There are differentiable functions for all x 0 distribution function F ( x ) and bounded functions f 1 ( x ) , f 2 ( x ) and real numbers γ { 0 , ± 1 / 2 , ± 1 } , a > 1 and 0 < C 1 < so that for all integers m 1
sup x | P m γ T m x F ( x ) m 1 / 2 f 1 ( x ) m 1 f 2 ( x ) | C 1 m a .
Remark 1.
In contrast to Bening et al. [10], the differentiability of F ( x ) , f 1 ( x ) and f 2 ( x ) is only required for x 0 . In the present article, in addition to the normal distribution, the chi-square distribution with p degrees of freedom is used as F ( x ) , which is not differentiable in x = 0 if p = 1 or p = 2 .
The distribution functions of the normalized random variables N n N + satisfy the following condition:
Assumption 2.
A distribution function H ( y ) with H ( 0 + ) = 0 , a function of bounded variation h 2 ( y ) , a sequence 0 < g n and real numbers b > 0 and C 2 > 0 exist so that for all integers n 1
sup y 0 P g n 1 N n y H ( y ) C 2 n b , f o r 0 < b 1 , sup y 0 P g n 1 N n y H ( y ) n 1 h 2 ( y ) C 2 n b , f o r b > 1 .

2.2. Transfer Proposition from Samples with Non-Random to Random Sample Sizes

Assumptions 1 and 2 allow the construction of expansions for distributions of normalized random-size statistics T N n based on approximate results for fixed-size normalized statistics T m in (1) and for the random size N n in (2).
Proposition 1.
Suppose γ { 0 , ± 1 / 2 , ± 1 } and the statistic T m and the sample size N n satisfy Assumptions 1 and 2. Then, for all n N + , the following inequality applies:
sup x R | P g n γ T N n x G n ( x , 1 / g n ) | C 1 E N n a + ( C 3 D n + C 4 ) n b ,
where
G n ( x , 1 / g n ) = 1 / g n F ( x y γ ) + f 1 ( x y γ ) g n y + f 2 ( x y γ ) g n y d H ( y ) + h 2 ( y ) n ,
D n = sup x 1 / g n y F ( x y γ ) + f 1 ( x y γ ) g n y + f 2 ( x y γ ) y g n d y ,
a > 1 , b > 0 and f 1 ( z ) , f 2 ( z ) , h 2 ( y ) are given in (1) and (2). The constants C 1 , C 3 , C 4 do not depend on n.
General transfer theorems with more terms are proved in [9,10] for γ 0 .
Remark 2.
The approximation function G n ( x , 1 / g n ) is not a polynomial in g n 1 / 2 and n 1 / 2 . The domain [ 1 / g n , ) of integration in (4) depends on g n . Some of the integrals in (4) could tend to infinity with 1 / g n 0 as n .
The following statement clarifies the problem.
Proposition 2.
In addition to the conditions of Proposition 1, let the following conditions be satisfied on the functions H ( . ) and h 2 ( . ) , depending on the rate of convergence b > 0 in (2):
H ( 1 / g n ) c 1 g n b , f o r b > 0 ,
0 1 / g n y 1 / 2 d H ( y ) c 2 g n b + 1 / 2 , f o r b > 1 / 2 ,
i : 0 1 / g n y 1 d H ( y ) c 3 g n b + 1 , i i : h 2 ( 0 ) = 0 a n d | h 2 ( 1 / g n ) | c 4 n g n b , i i i : 0 1 / g n y 1 | h 2 ( y ) | d y c 5 n g n b , f o r b > 1 .
Then, for the function G n ( x , 1 / g n ) defined in (4), one has
sup x | G n ( x , 1 / g n ) G n , 2 ( x ) I 1 ( x , n ) I 2 ( x , n ) I 3 ( x , n ) I 4 ( x , n ) | C g n b
with
G n , 2 ( x ) = 0 F ( x y γ ) d H ( y ) , f o r 0 < b 1 / 2 , 0 F ( x y γ ) + f 1 ( x y γ ) g n y d H ( y ) = : G n , 1 ( x ) , f o r 1 / 2 < b 1 , G n , 1 ( x ) + 0 f 2 ( x y γ ) g n y d H ( y ) + 0 F ( x y γ ) n d h 2 ( y ) , f o r b > 1 ,
I 1 ( x , n ) = 1 / g n f 1 ( x y γ ) g n y d H ( y ) f o r b 1 / 2 , I 2 ( x , n ) = 1 / g n f 2 ( x y γ ) g n y d H ( y ) f o r b 1 ,
I 3 ( x , n ) = 1 / g n f 1 ( x y γ ) n g n y d h 2 ( y ) a n d I 4 ( x , n ) = 1 / g n f 2 ( x y γ ) n g n y d h 2 ( y ) f o r b > 1 .
Remark 3.
The lower limit of integration in I 1 ( x , n ) to I 4 ( x , n ) in (10) and (11) depends on g n . If the sample size N n = N n ( r ) is negative binomial distributed with, e.g., 0 < r < 1 / 2 or 1 < r < 2 and g n = r ( n 1 ) + 1 (see (28) below), then both I 1 ( x , n ) and I 4 ( x , n ) have order n r and not n 1 / 2 or n 2 , as it seems at first glance.
Remark 4.
The additional conditions (6)–(8) guarantee to extend the integration range of the integrals in (9) from [ 1 / g n , ) to ( 0 , ) .
Proof of Propositions 1 and 2:
Evidence of Proposition 1 follows along the similar arguments of the more general Transfer Theorem 3.1 in [10] for γ 0 . The proof was adapted by Christoph and Ulyanov [13] to negative γ < 0 , too. Therefore, the Proposition 1 applies to γ { 0 , ± 1 / 2 , ± 1 } .
The present Propositions 1 and 2 differ from Theorems 1 and 2 in [13] only by the additional term f 1 ( x y γ ) ( g n y ) 1 / 2 and the added condition (7) to estimate this additional term. Therefore, the details are omitted her. □
Remark 5.
In Appendix 2 of the monograph by Gnedenko and Korolev [8], asymptotic expansions for generalized Cox processes are proved (see Theorems A2.6.1–A2.6.3). As random sample size, the authors considered a Cox process N ( t ) controlled by a Poisson process Λ ( t ) (also known as a doubly stochastic Poisson process) and proved asymptotic expansions for the random sum S ( t ) = k = 1 N ( t ) X k , where X 1 , X 2 , are independent identically distributed random variables. For each t 0 , the random variables N ( t ) , X 1 , X 2 , are independent. The above-mentioned theorems are close to Proposition 1. The structure of the functions G 2 ; n ( . ) in (4) and the bounds on the right-hand side of inequality (3) in Proposition 1 differ from the corresponding terms in Theorems A2.6.1–A2.6.3. Thus, the bounds contain little o-terms.

3. Chebyshev–Edgeworth Expansions Based on Standard Normal and Chi-Square Distributions

We consider two classes of statistics which are asymptotically normal or chi-square distributed.

3.1. Examples for Asymptotically Normally Distributed Statistics

Let X , X 1 , X 2 , be independent identically distributed random variables with
E X 5 < , E ( X ) = μ , 0 < Var ( X ) = σ 2 , skewness λ 3 = σ 3 E X μ 3 and kurtosis λ 4 = σ 4 E X μ 4 .
The random variable X is assumed to satisfy Cramér’s condition
lim sup | t | E e i t X < 1 .
Consider the asymptotically normal sample mean:
X ¯ m = X 1 + + X m / m m = 1 , 2 , ,
It follows from Petrov [18], Theorem 5.18 with k = 5 , that
sup x P ( σ 1 m ( X ¯ m μ ) x ) Φ 2 ; m ( x ) C m 3 / 2 ,
with C being independent of m and second order expansion
Φ 2 ; m ( x ) = Φ ( x ) λ 3 6 m H 2 ( x ) + 1 m λ 4 24 H 3 ( x ) + λ 3 2 72 H 5 ( x ) φ ( x ) ,
where Φ ( x ) and φ ( x ) are standard normal distribution function and its density and H k ( x ) are the Chebyshev–Hermite polynomials
H 2 ( x ) = x 2 1 H 3 ( x ) = x 3 3 x a n d H 5 ( x ) = x 5 10 x 3 + 15 x .
Let the random variable χ d 2 be chi-square distributed with d degrees of freedom having distribution function G d ( x ) and density function g d ( x ) :
G d ( x ) = P ( χ d 2 x ) = 0 x g d ( y ) d y a n d g d ( y ) = 1 2 d / 2 Γ ( d / 2 ) y ( d 2 ) / 2 e y / 2 , y > 0 .
Next, we examine the scale-mixed normalized statistic T m = m Z / χ m 2 , where Z and χ m 2 are independent random variables with the standard normal distribution Φ ( x ) and the chi-square distribution G m ( x ) , respectively. Then, the statistic T m = m Z / χ m 2 follows the Student’s t-distribution with m degrees of freedom. Example 2.1 in [19] indicates
P m Z χ m 2 x Φ ( x ) ( x 3 + x ) φ ( x ) 4 m sup x { | x 5 + 2 x 3 + 3 x | φ ( x ) } 6 m 2 + 6 ( m + 4 ) m 3 30.5 m 2 .
Chebyshev–Edgeworth expansions of Student’s t-statistic under non-normality are well investigated, but only Hall [20] proved these under minimal moment condition. Let conditions (12) and (13) are satisfied for independent identically distributed random variables X 1 , X 2 , Define T m * = m 1 / 2 ( X ¯ m μ ) / σ ^ m with sample mean X ¯ n and biased sample variance σ ^ m 2 = m 1 i = 1 m ( X i X ¯ m ) 2 . It follows from Hall [20] that for Student’s t-statistic T m * :
R m ( x ) = P m 1 / 2 X ¯ m μ σ ^ m x Φ ( x ) φ ( x ) P 1 ( x ) m + P 2 ( x ) m C m 3 / 2 ( 1 + u ( m ) )
uniformly in x, where u ( m ) 0 as m ,
P 1 ( x ) = λ 3 ( 2 x 2 + 1 ) / 6 and P 2 ( x ) = x λ 4 12 ( x 2 3 ) λ 3 2 18 ( x 4 + 2 x 2 3 ) 1 4 ( x 2 + 3 ) .
Remark 6.
The estimate (19) does not satisfy (1) in Assumption 1 because we do not have a computable error bound U with | u ( m ) | U < for all m N + . The estimate (19) does not satisfy (1) in Assumption 1 because we do not have a computable constant C with | u ( m ) | C < for all m N + , if all parameter are given. The remainder in (19) meets order condition R m ( x ) = O ( m 3 / 2 ) as m , but in the equivalent condition s u p x R m ( x ) C m 3 / 2 for all m M the values C > 0 and M > 0 are unknown. About non-asymptotic bounds and order conditions, see the work of Fujikoshi and Ulyanov [19] (Section 1.1).
In [21], an inequality for a first order approximation is proved:
sup x P m 1 / 2 X ¯ m μ σ ^ m x Φ ( x ) P 1 ( x ) φ ( x ) m C m 1 ,
where E | X | 4 + ε < is required for arbitrary ε > 0 and P 1 ( x ) is defined in (20).

3.2. Examples for Asymptotically Chi-Square Distributed Statistics

The baseline distribution of the second order expansions is now the chi-square distribution G d ( x ) occurring as limit distribution in different multivariate tests (see [22], Chapters 5 and 8–10, [19,23]).
At first, we consider statistic T m = T 0 2 = m tr S q S m 1 , where S q and S m are random matrices independently distributed as Wishart distributions W p ( q , I p ) and W p ( m , I p ) , respectively, with identity operator I p in R p . Note that W has Wishart distribution W p ( q , Σ ) if q p and its density is
1 2 p q / 2 Γ p ( q / 2 ) | Σ | q / 2 exp 1 2 tr Σ 1 W | W | ( q p 1 ) / 2 ,
where Γ p ( q / 2 ) = π p ( p 1 ) / 4 Π k = 1 p Γ ( q k + 1 ) / 2 (see [23], Chapter 2, for some basic properties).
Hotelling’s generalized T 0 2 distribution allows approximation
sup x P m tr S q S m 1 x G d ( x ) d 4 m a 0 G d ( x ) + a 1 G d + 2 ( x ) + a 2 G d + 4 ( x ) C m 2
(see [24], Theorem 4.1), where
d = p q , a 0 = q p 1 , a 1 = 2 q and a 2 = q + p + 1 with a 0 + a 1 + a 2 = 0 .
If T m = χ d 2 / χ m 2 is a scale mixture, where χ d 2 and χ m 2 are independent, T m allows asymptotic expansion
sup x P m χ d 2 / χ m 2 x G d ( x ) d 4 m a 0 G d ( x ) + a 1 G d + 2 ( x ) + a 2 G d + 4 ( x ) C m 2
(see [25], Section 5), where now
a 0 = 2 d , a 1 = 2 d and a 2 = ( 2 + d ) with a 0 + a 1 + a 2 = 0 .
Integration by parts gives G k + 2 ( x ) = 2 g k + 2 ( x ) + G k ( x ) . Moreover, g k + 2 ( x ) = ( x / k ) g k ( x ) for k = d and k = d + 2 . Then, it follows for both statistics Z m = tr S q S m 1 in (22) and Z m = χ d 2 / χ m 2 in (24) that
sup x P m Z m x G d ( x ) + g d ( x ) m ( a 1 + a 2 ) x 2 + a 2 x 2 2 ( d + 2 ) C m 2
where the coefficients a 1 and a 2 are defined in (23) and (25).
The scaled mixture T m = m χ 4 2 / χ m 2 is considered in the works by Fujikoshi et al. [23] (Example 13.2.2) and Fujikoshi and Ulyanov [19] (Example 2.2). The estimation given there leads to a computable error bound:
sup x | P m χ 4 2 χ m 2 x G 4 ( x ) + ( 2 x x 2 ) g 4 ( x ) 2 m | x 2 | x 2 4 | e x / 2 12 m 2 + 12 ( m + 4 ) m 3 65.9 m 2 .
Remark 7.
The statistics T m in (15), (18) and (21) satisfy Assumption 1 with the normal limit distribution Φ ( x ) and in (26) and (27) with chi-square distributions G d ( x ) and G 4 ( x ) , respectively.

4. Chebyshev–Edgeworth Expansions for Distributions of Normalized Random Sample Sizes

As in the articles by, e.g., Bening et al. [9,10], Christoph et al. [11,12] and Christoph and Ulyanov [13] and Christoph and Ulyanov [14], we consider as random sample sizes N n the negative binomial random variable N n ( r ) and the maximum of n independent discrete Pareto random variables N n ( s ) where r > 0 and s > 0 are parameters.
“The negative binomial distribution is one of the two leading cases for count models, it accommodates the overdispersion typically observed in count data (which the Poisson model cannot)” [26]. Moreover, E N n ( r ) < and P ( N n ( r ) / E N n ( r ) y ) tends to the gamma distribution G r r ( y ) with identical shape and rate parameters r > 0 .
On the other hand, the mean for the discrete Pareto-like variable N n ( s ) does not exist, yet P ( N n ( s ) / n x ) tends to the inverse exponential distribution W s ( y ) = e s / y with scale parameter s > 0 .
Remark 8.
The authors of [1,2,3,4,27], among others, considered the binomial or Poisson distributions as random number N of observations. If N = N n is binomial (with parameters n and 0 < p < 1 ) or Poisson (with rate λ n , 0 < λ < ) distributed, then P ( N n E N n x ) tends to the degenerated in 1 distribution as n . Therefore, Assumption 2 for the Transfer Proposition 1 is not fulfilled. On the other hand, since binomial or Poisson sample sizes are asymptotically normally distributed and if the statistic T m is also asymptotically normally distributed, so is the statistic T N n , too (see [28]). Chebyshev–Edgeworth expansions for lattice distributed random variables exist so far only with bounds of small-o or large- O order (see [29]). For (2) in Assumption 2, computable error bounds C 2 are required because the constant C 3 in (3) depends on C 2 (see also Remark 6 on large- O -bounds and computable error bounds).

4.1. The Random Sample Size N n = N n ( r ) Has Negative Binomial Distribution with Success Probability 1 / n

The sample size N n ( r ) has a negative binomial distribution shifted by 1 with the parameters 1 / n and r > 0 , the probability mass function
P ( N n ( r ) = k ) = Γ ( k + r 1 ) Γ ( k ) Γ ( r ) 1 n r 1 1 n k 1 , k = 1 , 2 ,
and g n = E ( N n ( r ) ) = r ( n 1 ) + 1 . Bening and Korolev [30] and Schluter and Trede [26] showed
lim n sup y P ( N n ( r ) / g n y ) G r , r ( y ) = 0 ,
where G r , r ( y ) is the gamma distribution function with its density
g r , r ( y ) = r r Γ ( r ) y r 1 e r y I ( 0 ) ( y ) , y R .
In addition to the expansion of N n ( r ) , a bound of the negative moment E ( N n ( r ) ) a in (3) is required, where m a is rate of convergence of the Chebyshev–Edgeworth expansion for T m in (1).
Proposition 3.
Suppose that r > 0 and the discrete random variables N n ( r ) have probability mass function (28) with g n : = E N n ( r ) = r ( n 1 ) + 1 . Then,
sup y 0 P N n ( r ) g n y G r , r ( y ) h 2 ; r ( y ) n C 2 ( r ) n min { r , 2 } ,
for all n N + , where the constant C 2 ( r ) > 0 does not depent on n and
h 2 ; r ( y ) = 0 , f o r r < 1 , ( 2 r ) 1 g r , r ( y ) ( y 1 ) ( 2 r ) + 2 Q 1 g n y , f o r r 1 .
Q 1 ( y ) = 1 / 2 ( y [ y ] ) a n d [ . ] d e n o t e s   t h e   i n t e g e r   p a r t   o f   a   n u m b e r .
Moreover, negative moments E ( N n ( r ) ) a fulfill the estimate for all r > 0 , α > 0
E N n ( r ) α C ( r ) n min { r , α } , r α ln ( n ) n α , r = α
and the convergence rate in case r = α cannot be improved.
Proof. 
In [10] (Formula (21)) and in [31] (Formula (11)), the convergence rate is reported for the case r < 1 . In [11] (Theorem 1), the Chebyshev–Edgeworth expansion for r > 1 is proved. In the case r = 1 , for geometric distributed random variable N n ( 1 ) N + with success probability 1 / n the proof is straightforward:
P N n n y = 1 P ( N n [ n y ] + 1 ) = 1 1 1 n n y τ = 1 e y + e y n y 2 τ + r n ( y ) ,
where sup y | r n ( y ) | C n 2 and τ = n y [ n y ] = 1 / 2 Q 1 ( n y ) [ 0 1 ) . Hence, (31) holds for r = 1 .
In [12] (Corollary 4.2), leading terms for the negative moments of E N n ( r ) p are derived, which lead to (34). □
Remark 9.
The negative binomial random variables N n ( r ) satisfy (2) in Assumption 2 and the additional conditions (6), (7) and (8) in Proposition 2 with H ( y ) = G r , r ( y ) , h 2 ( y ) = h 2 ; r ( y ) , g n = E N n ( r ) = r ( n 1 ) + 1 and b = min { r 2 } . The jumps of the distribution function P ( N n ( r ) g n y ) only affect the function Q 1 ( . ) in the term h 2 ; r ( . ) .

4.2. The Random Sample Size N n = N n ( s ) Is the Maximum of n Independent Discrete Pareto Variables

We consider the continuous Pareto Type II (Lomax) distribution function
F Y * ( x ) = 1 1 + ( x 1 ) / s 1 f o r x 1 .
The discrete Pareto II distribution F Y ( s ) is obtained by discretizing the continuous Pareto distribution F Y * ( x ) , : P ( Y ( s ) = k ) = F Y * ( k ) F Y * ( k 1 ) , k N + . The random variable Y ( s ) is the discrete counterpart on the positive integers to the continuous random variable Y * . Both random variables Y * and Y ( s ) have shape parameter 1 and scale parameter s > 0 (see [32]). The discrete Pareto distributed Y ( s ) has probability mass and distribution functions:
P ( Y ( s ) = k ) = s s + k 1 s s + k a n d P Y ( s ) k = k s + k , f o r k N + .
Let Y 1 ( s ) , Y 2 ( s ) , be a sequence of independent random variables with the common distribution function (35). Define
N n ( s ) = max 1 j n Y j ( s ) w i t h P ( N n ( s ) k ) = k s + k n , n N + , k N + s > 0 .
The random variable N n ( s ) is extremely spread over the positive integers.
Proposition 4.
Consider the discrete random variable N n ( s ) with distribution function (36). Then,
sup y > 0 P N n ( s ) n y W s ( y ) h 2 ; s ( y ) n C 3 ( s ) n 2 f o r a l l n N + a n d f i x e d s > 0 ,
W s ( y ) = e s / y a n d h 2 ; s ( y ) = s e s / y s 1 + 2 Q 1 ( n y ) / 2 y 2 , y > 0
where C 3 ( s ) > 0 does not depend on n and Q 1 ( y ) is defined in (33). Moreover,
E N n ( s ) p C ( p ) n min { p , 2 } ,
where for 0 < p 2 the order of the bound is optimal.
The Chebyshev–Edgeworth expansion (37) is proved in [11] (Theorem 4). In [12] (Corollary 5.2), leading terms for the negative moments E N n ( s ) p are derived for the negative moments that lead to (39).
Remark 10.
Let the random variable V ( s ) is exponentially distributed with rate parameter s > 0 . Then, W ( s ) = 1 / V ( s ) is an inverse exponentially distributed random variable with the continuous distribution function W s ( y ) = e s / y I ( 0 ) ( y ) . Both W s ( y ) and P ( N n ( s ) y ) are heavy tailed with shape parameter 1.
Remark 11.
Since E W ( s ) = and E N n ( s ) = for all n N + , we choose g n = n as normalizing factor for N n ( s ) in (37).
Remark 12.
The random sample sizes N n ( s ) satisfy (2) in Assumption 2 and the additional conditions (6)–(8) in Proposition 2 with W s ( y ) = e s / y , h 2 ( y ) = h 2 ; s ( y ) , g n = n and b = 2 . The jumps of the distribution function P ( N n ( s ) n y ) only affects the function Q 1 ( . ) in the term h 2 ; s ( . ) .
Remark 13.
Lyamin [33] proved a bound | P ( N n ( s ) n y ) W s ( y ) | 0.37 n 1 for integers s 1 .

5. Limit Distributions of Statistics with Random Sample Sizes Using Different Scaling Factors

The statistic T m from a sample with non-random sample size m N + fulfills condition (1) in Assumption 1. Instead of the non-random sample size m, we consider a random sample size N n N + satisfying condition (2) in Assumption 2. Let g n be a sequence with g n as n . Consider the scaling factor g n γ N n γ * γ by the statistics T N n with γ { 0 , ± 1 / 2 } if F ( x ) = Φ ( x ) and γ * = 1 / 2 or γ { 0 , ± 1 } if F ( x ) = G u ( x ) and γ * = 1 . Then, conditioning on N n and using (1) and (2), we have
P g n γ N n γ * γ T N n x = P N n γ T N n x ( N n / g n ) γ = m = 1 P m γ T m x ( m / g n ) γ P ( N n = m ) ( 1 ) E F x ( N n / g n ) γ = 1 / g n F ( x y γ ) d P ( N n / g n y ) ( 2 ) 1 / g n F ( x y γ ) d H ( y ) .
If there exists a limit distribution of P g n γ N n γ * γ T N n x as n , then it has to be a scale mixture of parent distribution F ( x ) and positive mixing parameter H ( y ) : 0 F ( x y γ ) d H ( y ) (see, e.g., [23,34], Chapter 13, and [19] and the references therein).
Remark 14.
Formula (40) shows that different normalization factors at T N n lead to different scale mixtures of the limit distribution of the normalized statistics T N n .

5.1. The Case F ( x ) = Φ ( x ) and H ( y ) = G r , r ( y )

The statistics (15), (18) and (21) considered in Section 3.1 have normal approximations Φ ( x ) . The limit distribution for the normalized random sample size N n ( r ) / E N n ( r ) is the gamma distribution G r , r ( y ) with density (30). We investigate the dependence of the limit distributions in P g n γ N n ( r ) 1 / 2 γ T N n ( r ) x 0 Φ ( x y γ ) d G r , r ( y ) as n for γ { 1 / 2 , 0 , 1 / 2 } .
(i) If γ = 1 / 2 , then the limit distribution is Student-s t distribution S 2 r ( x ) having density
s 2 r ( x ) = Γ ( ( r + 1 / 2 ) 2 r π Γ ( r ) 1 + x 2 2 r r + 1 / 2 , r > 0 , x R .
(ii) If γ = 0 , the standard normal law Φ ( x ) is the limit one with density φ ( x ) .
(iii) For γ = 1 / 2 , the generalized Laplace distributions L r ( x ) occur with density (see [13], Section 5.1.3):
l r ( x ) = r r Γ ( r ) 0 φ ( x y 1 / 2 ) y r 3 / 2 e r y d y = 2 r r Γ ( r ) 2 π | x | 2 r r 1 / 2 K r 1 / 2 ( 2 r | x | ) .
where K α ( u ) is the Macconald function of order α or modified Bessel function of the third kind with index α . The function K α ( u ) is also sometimes called a modified Bessel function of the second kind of order α . For properties of these functions, see, e.g., Chapter 51 in [35] or the Appendix on Bessel functions in [36].
If r N + , the so-called Sargan densities l 1 ( x ) , l 2 ( x ) , and their distribution functions are computable in closed forms (see Formulas (63)–(65) below in Section 7):
l 1 ( x ) = 1 2 e 2 | x | a n d L 1 ( x ) = 1 1 2 e 2 | x | , x > 0 l 2 ( x ) = 1 2 + | x | e 2 | x | a n d L 2 ( x ) = 1 1 2 ( 1 + x ) e 2 | x | , x > 0 l 3 ( x ) = 3 6 16 1 + 6 | x | + 2 x 2 e 6 | x | a n d L 3 ( x ) = 1 1 2 + 5 6 x 16 + 3 x 2 8 e 6 | x | ,
where L r ( x ) = 1 L r ( x ) for x 0 .
The double exponential or standard Laplace density is l 1 ( x ) with variance 1 and distribution function L 1 ( x ) given in (43). The Sargan distributions are therefore a generalisation of the standard Laplace distribution.

5.2. The Case F ( x ) = G d ( x ) and H ( y ) = W s ( y ) = e s / y

The statistics considered in Section 3.2 asymptotically approach chi-square distribution G d ( x ) . The limit distribution for the normalized random sample size N n ( s ) / n is the inverse exponential distribution W s ( y ) = e s / y I ( 0 , ) ( y ) .
(i) If γ = 1 , then the generalized gamma distribution W d ( x ; 2 s ) occurs with density w d ( x ; 2 s ) :
w d ( x ; 2 s ) = s s x 2 d / 4 1 / 2 K d / 2 1 ( 2 s x ) I ( 0 , ) ( x )
where the Macconald function K α ( u ) already appears in Formula (42) with different α and argument. For α = m + 1 / 2 , where m is an integer, the Macconald function K m + 1 / 2 ( u ) has a closed form (see Formulas (63)–(65) below in Section 7). Therefore, if d = 1 , 3 , 5 is an odd number, then the density w d ( x ; 2 s ) may be calculated in closed form. The distribution functions W d ( x ; 2 s ) with density functions w d ( x ; 2 s ) for d = 1 , 3 , 5 and x > 0 are
w 1 ( x ; 2 s ) = s ( 2 s x ) 1 / 2 e 2 s x and W 1 ( x ; 2 s ) = 1 e 2 s x ,
w 3 ( x ; 2 s ) = s e 2 s x and W 3 ( x ; 2 s ) = 1 e 2 s x ( 2 s x + 1 ) ,
w 5 ( x ; 2 s ) = s 3 1 + 2 2 s x e 2 s x and W 5 ( x ; 2 s ) = 1 e 2 s x 2 s x + 2 s x 3 .
Remark 15.
Functions in (45) are Weibull density and distribution functions, in (46) there are density and distribution functions of a generalized gamma distribution, but w 5 ( x ) and W 5 ( x ) are even more general.
The family of generalized gamma distributions contains many absolutely continuous distributions concentrated on the non-negative half-line.
Remark 16.
The generalized gamma distribution G * ( x ; r , α , λ ) corresponds to the density
g * ( x ; r , α , λ ) = | α | λ r Γ ( r ) x α r 1 e λ x α , x 0 , | α | > 0 , r > 0 , λ > 0 ,
where α and r are the two shape parameters and λ the scale parameter. The density representation (48) is suggested in the work of Korolev and Zeifman [37] or Korolev and Gorshenin [38], and many special cases are listed therein. In addition to, e.g., Gamma and Weibull distributions (with a > 0), inverse Gamma, Lévy and Fréche distributions (with a < 0) also belong to that family of generalized gamma distributions.
Remark 17.
The Weibull density in (45) is w 1 ( x ; 2 s ) = g * ( x ; 1 , 1 / 2 , 2 s ) . Moreover, w 3 ( x ; 2 s ) = g * ( x ; 2 , 1 / 2 , 2 s ) . The densities w 5 ( x ) , w 7 ( x ) , w 9 ( x ) , are weighted sums of generalized gamma distribution with different shape parameters r, e.g.,
w 5 ( x ; 2 s ) = 1 3 g * ( x ; 2 , 1 2 , 2 s ) + 2 3 g * ( x ; 3 , 1 2 , 2 s ) w 9 ( x ; 2 s ) = 5 35 g * ( x ; 2 , 1 2 , 2 s ) + 10 35 g * ( x ; 3 , 1 2 , 2 s ) + 12 35 g * ( x ; 4 , 1 2 , 2 s ) + 8 35 g * ( x ; 5 , 1 2 , 2 s ) .
(i) If γ = 1 . For better readability I have introduced for (i) and after (ii).
(ii) If γ = 0 , the standard normal law Φ ( x ) is the limit distribution with density φ ( x ) .
(iii) If γ = 1 , as limit distribution the inverse Pareto distribution occurs V d / 2 ( x ; 2 s ) with shape parameter d / 2 , scale parameter 2 s and density v d / 2 ( x ; 2 s ) :
V d / 2 ( x ; 2 s ) = x 2 s + x d / 2 and v d / 2 ( x ; 2 s ) = s d x d / 2 1 ( x + 2 s ) d / 2 + 1 for x 0
In [39], a robust and efficient estimator for the shape parameter of the inverse Pareto distribution and applications are given.

6. Main Results

We examine asymptotic approximations of P g n γ N n γ * γ T N n x depending on the scaling factor g n γ N n γ * γ for γ { 0 , ± 1 / 2 , ± 1 } , for γ * = 1 / 2 if the statistic T m is asymptotically normal or γ * = 1 if T m is asymptotically chi-square distributed.

6.1. Asymptotically Normal Distributed Statistics with Negative Binomial Distributed Sample Sizes

Consider first the statistics estimated in (15), (18) and (21) with the normal limiting distribution Φ ( x ) . They have the form
P ( m Z m x ) Φ ( x ) m 1 / 2 ( p 0 + p 2 x 2 ) + m 1 ( p 1 x + p 3 x 3 + p 5 x 5 ) I a > 1 ( a ) φ ( x ) m a .
The sample size is negative binomial N n = N n ( r ) with probability mass function (28).
Theorem 1.
Let r > 0 . If inequality (50) for the statistic Z m and inequality (31) for random sample size N n ( r ) with g n = E N n ( r ) = r ( n 1 ) + 1 hold, then, for all n N + , the following expansions apply:
i: 
The non-random scaling factor g n by statistic Z N n ( r ) leads to Student’s t-approximation.
sup x P g n Z N n ( r ) x S 2 r ; n ( x ) C r n r ln ( n ) , r { 1 / 2 , 3 / 2 , 2 } , n min { r , 2 } , r { 1 / 2 , 3 / 2 , 2 }
where
S 2 r ; n ( x ) = S 2 r ( x ) + s 2 r ( x ) g n p 0 x 2 + 2 r 2 r 1 + p 2 x 2 I { r > 1 / 2 } ( r ) + s 2 r ( x ) g n p 1 x ( x 2 + 2 r ) ( 2 r 1 ) + p 3 x 3 + p 5 x 5 ( 2 r + 1 ) x 2 + 2 r + ( 2 r ) x ( x 2 + 1 ) 4 ( 2 r 1 ) I { r > 1 } ( r ) ,
S 2 r ( x ) is Student’s t-distribution having density s 2 r ( x ) , defined in (41), and p k are the coefficients in (50).
ii: 
The standard normal approximation occurs at random scaling factor N n ( r ) by statistic Z N n ( r ) :
sup x P ( N n ( r ) Z N n ( r ) x ) Φ n 2 ( x ) C r n min { r , 2 } , r 2 , ln ( n ) n 2 , r = 2 ,
where
Φ n 2 ( x ) = Φ ( x ) + r Γ ( r 1 / 2 ) Γ ( r ) g n ( p 0 + p 2 x 2 ) φ ( x ) I { r > 1 / 2 } ( r ) + ( p 1 x + p 3 x 3 + p 5 x 5 ) φ ( x ) g n ln n I { r = 1 } ( r ) + r r 1 I { r > 1 } ( r ) ,
iii: 
If r = 2 , the mixed scaling factor g n 1 / 2 N n ( 2 ) by statistic Z N n ( 2 ) leads to generalized Laplace approximation:
sup x P g n 1 / 2 N n ( 2 ) Z N n ( 2 ) x L 2 ( x ) l n ; 2 ( x ) C 2 ln ( n ) n 2
where L 2 ( x ) = 1 1 2 ( 1 + x ) e 2 | x | , L 2 ( x ) = 1 L 2 ( x ) for x 0 and
l n ; 2 ( x ) = e 2 | x | g n p 0 ( | x | + 1 / 2 ) + 2 p 2 x 2 ) e 2 | x | g n p 1 x + 4 p 3 | x | x + 4 p 5 ( 2 x 3 + | x | x ) .
Remark 18.
Analogous to (54) and (55), expansions for all r > 0 can be derived from Formulas (42) and (63)–(66) below in Section 7, whereby closed forms can be presented only for r { 1 , 2 , 3 , } .
The statistics from Section 3.1 are considered with different normalization factors as applications of Theorem 1:
Corollary 1.
Let the conditions of Theorem 1 be satisfied:
i: 
In the case of the Student’s t-statistic Z / χ m 2 estimated in (18), one has (51) with p 0 = p 2 = p 5 = 0 and p 1 = p 3 = 1 / 4 using non-random scaling factor g n :
sup x P g n Z χ N n ( r ) 2 x S 2 r ( x ; n ) C r n min { r , 2 } , r 2 , ln ( n ) n 2 , r = 2 ,
where
S 2 r ( x ; n ) = S 2 r ( x ) s 2 r ( x ) 2 r ( x + x 3 ) ( 2 r ) x ( x 2 + 1 ) 4 ( 2 r 1 ) g n I { r > 1 } ( r )
ii: 
In the case of Student’s one-sample t-test statistic under non-normality T m = ( X ¯ m μ ) / σ ^ m estimated in (21) with a = 1 , the first-order approximation defined in (52) for 0 < r 1 and (53) with p 0 = λ 3 / 6 and p 2 = λ 3 / 3 using random scaling factor N n ( r ) leads uniformly in x to:
P ( N n ( r ) T N n ( r ) x ) Φ ( x ) + r Γ ( r 1 / 2 ) λ 3 ( 2 x 2 + 1 ) 6 Γ ( r ) g n φ ( x ) C r n min { r , 1 } , r 1 , ln ( n ) n 1 , r = 1 .
iii: 
Considering sample mean X ¯ m estimated in (15), one has (55) with p 0 = p 2 = λ 3 / 6 , p 1 = λ 4 / 8 5 λ 3 2 / 24 , p 3 = λ 4 / 24 + 5 λ 3 2 / 36 , and p 5 = λ 3 2 / 72 using mixed scaling factor g n 1 / 2 N n ( 2 ) :
sup x P g n 1 / 2 N n ( 2 ) X ¯ N n ( 2 ) x L 2 ( x ) l 2 ; n ( x ) C 2 ln ( n ) n 2 ,
where the generalized Laplace distributions L 2 ( x ) is defined in (43) and
l 2 ; n ( x ) = λ 3 e 2 | x | 6 g n 2 x 2 | x | 1 / 2 + e 2 | x | 36 g n λ 4 ( 9 x 6 x | x | ) + 12 λ 3 2 ( 4 x 3 + 18 x | x | 15 x ) .
Remark 19.
The approximating functions in the expansions for P g n γ N n ( r ) 1 / 2 γ T N n ( r ) x with the statistics estimated in (15), (18) and (21) can only be given in closed form for all r > 0 in the case of non-random ( γ = 1 / 2 ) or random ( γ = 0 ) normalization factors. In the case of the mixed ( γ = 1 / 2 ) normalization factor, only for positive integer r closed forms are available, while in the other cases Macconald functions are involved.

6.2. Asymptotically Chi-Square Distributed Statistics with Pareto-Like Distributed Sample Sizes

Consider now the statistics, estimated in (26) and (27) with limit chi-square distributions. They have the form
P ( m T m x ) G d ( x ) m 1 ( q 1 x + q 2 x 2 ) g d ( x ) C ( s ) m 2 .
The sample size is the Pareto-like random variable N n = N n ( s ) with probability mass function (35).
Theorem 2.
Let s > 0 and (36) be the distribution function of the random sample size N n = N n ( s ) . If for the statistic T m the inequality (56) with limiting chi-square distribution G d ( x ) and the inequality (37) with g n = n for the random sample size N n ( s ) hold, then for all n N + one has the following approximation:
i: 
The non-random scaling factor n by T N n ( s ) leads to the limiting generalized gamma distributions.
sup x > 0 P n T N n ( s ) x W d ; n ( x ; 2 s ) C ( s ) n 2 ln n f o r d = 1 a n d d = 3 ,
W 1 ; n ( x ; 2 s ) = W 1 ( x ; 2 s ) + n 1 w 1 ( x ; 2 s ) q 1 x 1 + 2 s x + q 2 x 2 ( s 1 ) x ( 2 s x + 1 ) 4 s
and
W 3 ; n ( x ; 2 s ) = W 3 ( x ; 2 s ) + n 1 w 3 ( x ; 2 s ) q 1 x 2 s x 2 s + q 2 x 2 ( s 1 ) x 2 2 2 s x .
where the limit law W d ( x ; 2 s ) with density w d ( x ; 2 s ) for d = 1 and d = 3 are given in (45) and (46).
ii: 
The random scaling factor N n ( s ) by T N n ( s ) induces the limiting chi-square distribution.
sup x P ( N n ( s ) T N n ( s ) x ) G d ( x ) g d ( x ) s n q 1 x + q 2 x 2 C ( s ) ln n n 2 ,
iii: 
Limiting inverse Pareto distributions occur at mixed scaling factor n 1 N n 2 ( s ) by T N n ( s ) .
sup x > 0 P N n 2 ( s ) n T N n ( s ) x V d / 2 ( x ; 2 s ) 1 n v d / 2 ; n ( x ; 2 s ) C ( s ) ln n n 2 ,
where
v d / 2 ; n ( x ; 2 s ) = v d / 2 ( x ; 2 s ) q 1 x ( d + 2 ) x + 2 s + q 2 x 2 ( d + 4 ) ( d + 2 ) ( x + 2 s ) 2 + ( s 1 ) x ( 2 + d ) 2 ( x + 2 s )
with inverse Pareto distribution V d / 2 ( x ; 2 s ) having shape parameter d / 2 , scale parameter 2 s and density v d / 2 ( x ; 2 s ) defined in (49).
Remark 20.
Analogous to (57), expansions for all d N + can be derived from Formulas (44), (63)–(65) and (69) below in Section 7, whereby closed forms can be given for d { 1 , 3 , 5 , } .
The statistics from Section 3.2 are considered with different normalization factors as applications of Theorem 2.
Corollary 2.
Let the conditions of Theorem 2 be satisfied.
i: 
Let χ d 2 / χ m 2 be scale mixture, estimated in (24), where χ d 2 and χ m 2 are independent. Then, using non-random scaling factor, n limiting generalized gamma distributions occur with q 1 = ( d 2 ) / 2 and q 2 = 1 / 2 in (58) and (59):
sup x > 0 P n tr χ d 2 / χ N n ( s ) 2 x W d ; n ( x ; 2 s ) C ( s ) n 2 ln n , f o r d = 1 a n d d = 3 ,
W 1 ; n ( x ; 2 s ) = W 1 ( x ; 2 s ) + n 1 w 1 ( x ; 2 s ) x 1 + 2 s x x 2 2 ( s 1 ) x ( 2 s x + 1 ) 4 s
x > 0 and
W 3 ; n ( x ; 2 s ) = W 3 ( x ; 2 s ) + n 1 w 3 ( x ; 2 s ) x 4 s 2 s x x 2 2 ( s 1 ) x 2 2 2 s x , x > 0 ,
where the limit law W d ( x ; 2 s ) with density w d ( x ; 2 s ) for d = 1 and d = 3 are given in (45) and (46).
ii: 
For the scaled mixture χ 4 2 / χ m 2 estimated in (27), one gets the limiting chi-square distribution with a random scaling factor N n ( s ) in (60) with q 1 = 1 and q 2 = 1 / 2 :
sup x P ( N n ( s ) χ 4 2 / χ N n ( s ) 2 x ) G d ( x ) g d ( x ) s n x x 2 / 2 ) C ( s ) ln n n 2 ,
iii: 
In the case of the Hotelling’s generalized T 0 2 statistic T 0 2 = m tr S q S m 1 estimated in (22), one has the limiting inverse Pareto distributions with mixed scaling factor n 1 N n 2 ( s ) by tr S q S N n ( s ) 1 . Here, (61) holds with q 1 = ( p + 1 q ) / 2 and q 2 = ( p + 1 + q ) / ( 2 d + 4 ) .
sup x > 0 P N n 2 ( s ) n tr S q S N n ( s ) 1 x V d / 2 ( x ; 2 s ) 1 n v d / 2 ; n ( x ; 2 s ) C ( s ) ln n n 2 ,
and
v d / 2 ; n ( x ; 2 s ) = v d / 2 ( x ; 2 s ) ( p + 1 q ) x ( d + 2 ) 2 ( x + 2 s ) + ( p + 1 + q ) x 2 ( d + 4 ) 2 ( x + 2 s ) 2 + ( s 1 ) x ( 2 + d ) 2 ( x + 2 s )
where the inverse Pareto distribution V d / 2 ( x ; 2 s ) with shape parameter d / 2 , scale parameter 2 s and density v d / 2 ( x ; 2 s ) is defined in (49).
Remark 21.
For the statistics estimated in (26) and (27), the approximating functions in the expansions for P g n γ N n ( s ) 1 γ T N n ( s ) x can only be given in closed form for all integer d in the case of non-random ( γ = 1 ) or random ( γ = 0 ) normalization factors. In the case of the mixed ( γ = 1 ) normalization factor, only for odd integer d in closed form can be presented; for even integer d, the Macconald functions are involved.

7. Formal Construction of the Expansions

Expansions of the statistics considered in (15), (18), (21), (26) and (27) have the structure:
G ( x ) + g ( x ) m 1 / 2 P 1 ( x ; j 1 * ) + m 1 P 2 ( x ; j 2 * )
with g ( x ) = G ( x ) and polynomials P 1 ( x ; j 1 * ) , P 2 ( x ; j 2 * ) of degrees j 1 * and j 2 * , respectively. Here, G ( x ) = Φ ( x ) or G ( x ) = P ( χ d 2 x ) .
We calculate the integrals with k = 1 , 2 and j = 0 , 1 , , j k * :
J 1 ( x ; γ ) = 0 G ( x y γ ) d H ( y ) and J 2 ( x ; γ , k , j ) = x j 0 y γ j k / 2 g ( x y γ ) d H ( y ) .
The limit distributions of the random sizes N n are H ( y ) = G r , r ( y ) and H ( y ) = W s ( y ) = e s / y with corresponding second approximation h 2 ( y ) .
We use the following formulas several times: Formula 2.3.3.1 in [40]
M α ( p ) = 0 y α 1 e p y d y = y = 1 / z 0 z α 1 e p / z d z = Γ ( α ) p α α > 0 , p > 0 .
and Formula 2.3.16.1 in [40] with real α and p , q > 0 :
K α * ( p , q ) = 0 y α 1 e p y q / y d y = 2 q p α / 2 K α ( 2 p q ) ,
where the Macconald function K α ( u ) already appears in Formula (42) with different α and argument.
For α = m + 1 / 2 , where m is an integer, the Macdonald function K m + 1 / 2 ( u ) has a closed form (see Formulas 2.3.16.2 and 2.3.16.3 in [40] with p , q > 0 ):
K m + 1 / 2 * ( p , q ) = K m * * ( p , q ) if m is   an   integer ,
where
K m * * ( p , q ) = 0 y m 1 / 2 e p y q / y d y = ( 1 ) m π m p m p 1 / 2 e 2 p q , m = 0 , 1 , 2 , , ( 1 ) m π p m q m e 2 p q , m = 0 , 1 , 2 ,

7.1. The Case G ( x ) = Φ ( x ) and H ( y ) = G r , r ( y )

Consider statistics that meet the condition (50).
Let J 1 ( x ; γ ) = 0 Φ ( x y γ ) d G r , r ( y ) with γ { 0 , ± 1 / 2 } . Then, J 1 ( x ; γ ) = Φ ( x ) for γ = 0 and
x J 1 ( x ; γ ) = r r Γ ( r ) 2 π 0 y r + γ 1 e ( x 2 y 2 γ / 2 + r y ) d y = r r Γ ( r ) 2 π M r + 1 / 2 r 1 + x 2 / ( 2 r ) , for γ = 1 / 2 , K r 1 / 2 * ( r , x 2 / 2 ) = 2 x 2 2 r r / 2 1 / 4 K r 1 / 2 ( 2 r x ) , for γ = 1 / 2 .
If r > 0 is an integer number then using (64) with m = r 1 , the density of J 1 ( x ; 1 / 2 ) can be calculated with (65) in a closed form.
Let γ { 0 , ± 1 / 2 } . Let k = 1 , 2 and j = 0 , 1 , , 5 be the exponents at m k / 2 and x j in (50), respectively.
J 2 ( x ; γ , k , j ) = r r x j Γ ( r ) 2 π 0 y j γ + r 1 k / 2 e ( x 2 y 2 γ / 2 + r y ) d y = r r x j Γ ( r ) 2 π M j / 2 + r k / 2 r 1 + x 2 / ( 2 r ) , for γ = 1 / 2 , e x 2 / 2 M r k / 2 ( r ) , for γ = 0 , K r ( j + k ) / 2 * ( r , x 2 / 2 ) , for γ = 1 / 2 .
In (50) k + j are odd integers. If r > 0 is an integer, then K r ( j + k ) / 2 * ( r , x 2 / 2 ) = K r ( j + k + 1 ) / 2 * * ( r , x 2 / 2 ) .
Define p j * = p j I 2 ( x ; γ , k , j ) with coefficient p j from (50) and calculate the terms in (67):
γ = 1 / 2 , k = 1 : p 0 * = p 0 ( 2 r + x 2 ) 2 r 1 s 2 r ( x ) , p 2 * = p 2 x 2 s 2 r ( x ) , k = 2 : p 1 * = p 1 x ( x 2 + 2 r ) 2 r 1 s 2 r ( x ) , p 3 * = p 3 x 3 s 2 r ( x ) , p 5 * = p 5 x 5 ( 2 r + 1 ) x 2 + 2 r s 2 r ( x ) , γ = 1 / 2 , k = 1 , r = 2 : p 0 * = p 0 ( | x | + 1 / 2 ) e 2 | x | , p 2 * = p 2 2 x 2 e 2 | x | , k = 2 , r = 2 : p 1 * = p 1 2 x e 2 | x | , p 3 * = p 3 4 | x | x e 2 | x | , p 5 * = p 5 4 ( 2 x 3 + x | x | ) e 2 | x | , γ = 0 , k = 1 , 2 , j = 0 , 1 , 2 , 3 , 5 : p j * = p j x j φ ( x ) r k / 2 Γ ( r k / 2 ) / Γ ( r ) ,

7.2. The Case G ( x ) = G d ( x ) and H ( y ) = W s ( y ) = e s / y

Consider statistics that meet the condition (56). Let J 1 ( x ; γ ) = 0 G u ( x y γ ) s y 2 e s / y d y , γ { 0 , ± 1 } . Then, J 1 ( x ; 0 ) = G u ( x ) and
x J 1 ( x ; γ ) = s 0 y γ 2 g u ( x y γ ) e s / y d y = s x d / 2 1 2 d / 2 Γ ( d / 2 ) 0 y γ d / 2 2 e x y γ / 2 s / y d y , γ = ± 1 .
Let γ = 1 . Using (63) with α = d / 2 1 , p = x / 2 and q = s , we find
x J 1 ( x ; 1 ) = s 0 y g u ( x y ) e s / y d y = s x d / 2 1 2 d / 2 Γ ( d / 2 ) 0 y d / 2 2 e x y / 2 s / y d y = s x d / 2 2 d / 2 Γ ( d / 2 ) K d / 2 1 * ( x / 2 , s ) = s s x 2 d / 4 1 / 2 K d / 2 1 ( 2 s x ) .
If d = 1 , 3 , 5 , is an odd number, using the closed form K m * * ( p , q ) in (65) with m = ( d 3 ) / 2 , p = x / 2 and q = s , then J 1 ( x ; 1 ) = W d ( x , 2 s ) and its density w d ( x ; 2 s ) may be calculated in closed form:
w d ( x ; 2 s ) = x J 1 ( x ; 1 ) = s x d / 2 1 2 d / 2 1 Γ ( d / 2 ) K ( d 3 ) / 2 * * ( x / 2 , s ) for d = 1 , 3 , 5 ,
The distribution functions W d ( x ; 2 s ) and their densities w d ( x ; 2 s ) for d = 1 , 3 , 5 are given in (45)–(47).
If γ = 1 , we use (62) with α = d / 2 + 1 , p = ( x + 2 s ) / 2 , q = s and the substitution y = 1 / z :
x J 1 ( x ; 1 ) = s x d / 2 1 2 d / 2 Γ ( d / 2 ) 0 y d / 2 2 e ( x / 2 + s ) / y d y = y = 1 / z s x d / 2 1 2 d / 2 Γ ( d / 2 ) 0 z d / 2 e ( x + 2 s ) z / 2 d z = s x d / 2 1 Γ ( d / 2 + 1 ) 2 d / 2 + 1 2 d / 2 Γ ( d / 2 ) ( x + 2 s ) d / 2 + 1 = s d x d / 2 1 ( x + 2 s ) d / 2 + 1 = s d x 2 1 + 2 s x d / 2 1 = v d / 2 ( x ; 2 s ) .
where v d / 2 ( x ; 2 s ) is the density of the inverse Pareto distribution defined in (49).
Suppose γ { 0 , ± 1 } . Let j = 1 , 2 be the exponent at x j in (56). Then, by (65) for positive odd numbers d with α = j + ( d 7 ) / 2 , p = x / 2 , q = s if γ = 1 , by (62) with α = 2 , p = s for γ = 0 and with α = j + d / 2 1 , p = ( x + 2 s ) / 2 for γ = 1 :
J 2 ( x ; γ , 2 , j ) = s x j 2 d / 2 Γ ( d / 2 ) 0 y j γ 3 ( x y γ ) d / 2 1 e ( x y γ / 2 + s / y ) d y = s x j + d / 2 1 2 d / 2 Γ ( d / 2 ) 0 y γ ( j + d / 2 1 ) 3 e x y γ / 2 + s / y d y = s x j + d / 2 1 2 d / 2 Γ ( d / 2 ) K j + ( d 7 ) / 2 * * ( x / 2 , s ) , γ = 1 , e x / 2 M 2 ( s ) , γ = 0 , M j + 1 + d / 2 ( x + 2 s ) / 2 , γ = 1 .
If d is not an odd number, K j + ( d 7 ) / 2 * * ( x / 2 , s ) in (71) has to be replaced by K j + d / 2 3 * ( x / 2 , s ) , which may be calculated with (63) where the Macdonald functions K j + d / 2 3 ( 2 s x ) are involved.
Define q j * = q j I 2 ( x ; γ , 2 , j ) for j = 1 , 2 with the coefficient q j from (56). Calculating the corresponding terms in (71) we find
γ = 1 , d = 1 : q 1 * = q 1 x ( 2 s ) 1 1 + 2 s x w 1 ( x ) , q 2 * = q 2 x 2 w 1 ( x ) , γ = 1 , d = 3 : q 1 * = q 1 ( 2 s ) 1 x 2 s x w 3 ( x ) , q 2 * = q 2 x 2 w 3 ( x ) , γ = 1 : q 1 * = q 1 x ( d + 2 ) x + 2 s v d / 2 ( x ; 2 s ) , q 2 * = q 2 x 2 ( d + 4 ) ( d + 2 ) ( x + 2 s ) 2 v d / 2 ( x ; 2 s ) , γ = 0 : q 1 * = q 1 x s 1 g d / 2 ( x ) q * 21 = q 2 x 2 s 1 g d / 2 ( x )

8. Proof of Theorems

We find from Lemmas A1 and A2 that D n in (5) in Proposition 1 is bounded and the integrals in (10) and (11) in Proposition 2 have the necessary convergence rates. It remains to calculate the integrals in (9).
Proof of Theorem 1.
Let F ( x ) = Φ ( x ) , H ( y ) = G r , r ( y ) and h 2 ( y ) = h 2 ; r ( y ) defined in (32).
Suppose J 1 ( x ; γ ) = 0 Φ ( x y γ ) d G r , r ( y ) with γ { 0 , ± 1 / 2 } , which are the limit distributions in (9) for P g n γ N n ( r ) 1 / 2 γ Z N n ( r ) x under the condition of Theorem 1. Then, J 1 ( x ; γ ) = Φ ( x ) for γ = 0 . It follows from (66), (62) for γ = 1 / 2 and (65) for γ = 1 / 2 that
x J 1 ( x ; γ ) = s 2 r ( x ) = Γ ( r + 1 / 2 ) 2 r π Γ ( r ) 1 + x 2 2 r ( r + 1 / 2 ) , γ = 1 / 2 with J 1 ( x ; 1 / 2 ) = S 2 r ( x ) , φ ( x ) = 1 2 π e x 2 / 2 , γ = 0 with J 1 ( x ; 0 ) = Φ ( x ) , l 2 ( x ) = 1 2 + | x | e 2 | x | , r = 2 , γ = 1 / 2 , with J 1 ( x ; 1 / 2 ) = L 2 ( x ) ,
where s 2 r ( x ) is the density of Student’s t-distribution with 2 r degrees of freedom and l 2 ( x ) is the density of a generalized Laplace distribution.
Integral J 2 ( x ; γ ) = 0 y 1 / 2 p 0 + p 2 x 2 y 2 γ φ ( x y γ ) d G r , r ( y ) is the integral by g n 1 / 2 in the expansion (9). Then, using (67) and (68) with k = 1 , we obtain
J 2 ( x ; γ ) = p 0 J 2 ( x ; γ , 1 , 0 ) + p 2 J 2 ( x ; γ , 1 , 2 ) = p 0 * + p 2 * ,
Integral J 3 ( x ; γ ) = 0 y 1 p 1 x y γ + p 3 x 3 y 3 γ + p 5 x 5 y 5 γ φ ( x y γ ) d G r , r ( y ) is the integral by g n 1 in the expansion (9). Then, using again (67) and (68) with k = 2 , we obtain
J 3 ( x ; γ ) = p 1 J 2 ( x ; γ , 2 , 1 ) + p 3 J 2 ( x ; γ , 2 , 3 ) + p 5 J 2 ( x ; γ , 2 , 5 ) = p 1 * + p 3 * + p 5 *
Integration by parts in the last integral by n 1 in (9) for γ = ± 1 / 2 and r > 1 leads to
J 4 ( x ; γ ) = 0 Φ ( x y γ ) d h 2 , r ( y ) = γ x r r 2 r 2 π Γ ( r ) 0 y γ + r 2 e x 2 y 2 γ / 2 + r y ( y 1 ) ( 2 r ) + 2 Q 1 g n y d y .
Suppose γ = 1 / 2 . We find from (62)
J 4 ( x ; 1 / 2 ) = x r r ( 2 r ) 4 r 2 π Γ ( r ) M r 1 / 2 ( r + x 2 / 2 ) M r + 1 / 2 ( r + x 2 / 2 ) J 4 * ( x ; 1 / 2 ) = ( 2 r ) x ( x 2 + 1 ) 4 r ( 2 r 1 ) s 2 r ( x ) J 4 * ( x ; 1 / 2 ) .
with
J 4 * ( x ; 1 / 2 ) = x r r 1 2 2 π Γ ( r ) 0 y r 3 / 2 e ( r + x 2 / 2 ) y Q 1 g n y d y ,
where Q 1 ( y ) is defined in (33). It follows from Lemma A3 that for r > 1
sup x n 1 | J 4 * ( x ; 1 / 2 ) | c ( r ) n r .
Hence, because of 0 g n 1 ( r n ) 1 ( n g n ) 1 for r 1 , we obtain
1 n 0 Φ ( x y ) d h 2 ( y ) ( 2 r ) x ( x 2 + 1 ) 4 ( 2 r 1 ) g n s 2 r ( x ) 1 n | J 4 * | + C ( r ) n g n c 1 ( r ) n min { r 2 } .
For γ = 1 / 2 , we only consider the case r = 2 , which results in J 4 ( x ; 1 / 2 ) = 0 and
J 4 * ( x ; 1 / 2 ) = x 2 2 π 0 y 1 / 2 Q 1 g n y e ( 2 y + x 2 / ( 2 y ) ) d y ,
where sup x n 1 J 4 * ( x ; 1 / 2 ) C n 2 is proved in Lemma A3.
If γ = 0 , then J 4 ( x ; 0 ) = Φ ( x ) ( h 2 ; r ( ) h 2 ; r ( 0 ) ) = 0 since Q 1 ( 0 ) = 1 / 2 .
The proof of Theorem 1 follows from (73)–(75) and (77) and Lemma A3. □
Proof of Theorem 2.
Let F ( x ) = G d ( x ) , H ( y ) = W s ( y ) = e s / y and h 2 ( y ) = h 2 ; s ( y ) defined in (32).
Suppose J 1 ( x ; γ ) = 0 G d ( x y γ ) s y 2 e s / y d y with γ { 0 , ± 1 } which are the limit distributions in (9) for P g n γ N n ( s ) 1 / 2 γ Z N n ( s ) x under the condition of Theorem 2.
Then, J 1 ( x ; γ ) = G d ( x ) for γ = 0 . It follows from (69) and (65) for γ = 1 and (70) for γ = 1 that
x J 1 ( x ; γ ) = w 1 ( x ; 2 s ) = s ( 2 s x ) 1 / 2 e 2 s x γ = 1 with J 1 ( x ; 1 ) = W 1 ( x ) , w 3 ( x ; 2 s ) = s e 2 s x γ = 1 with J 1 ( x ; 1 ) = W 3 ( x ) , φ ( x ) = 1 2 π e x 2 / 2 , γ = 0 with J 1 ( x ; 0 ) = Φ ( x ) , v d / 2 ( x ; 2 s ) = s d x d / 2 1 ( x + 2 s ) d / 2 + 1 γ = 1 , with J 1 ( x ; 1 / 2 ) = V d / 2 ( x ; 2 s ) ,
where w 1 ( x ; 2 s ) is the Weibull density (see (45)), w 3 ( x ; 2 s ) is the generalized gamma density0 (see (46)) and v d / 2 ( x ; 2 s ) is the density of the inverse Pareto distribution V d / 2 ( x ; 2 s ) defined in (49).
Integral J 2 ( x ; γ ) = 0 y 1 q 1 x y γ + q 2 x 2 y 2 γ g d / 2 ( x y γ ) s y 2 e s / y d y is the integral by g n 1 in the expansion (9). Then, the use of (71) and (72) leads to
J 2 ( x ; γ ) = q 1 J 2 ( x ; γ , 2 , 1 ) + q 3 J 2 ( x ; γ , 2 , 3 ) = q 1 * + q 3 * .
Integration by parts in the last integral by n 1 in (9) for γ = ± 1 leads to -4.6cm0cm
J 3 ( x ; γ , d ) = 0 G d ( x y γ ) d h 2 ; s ( y ) = J 4 ( x ; γ , d ) + J 4 * ( x ; γ , d ) , J 4 ( x ; γ , d ) = s ( s 1 ) γ x d / 2 2 d / 2 + 1 Γ ( d / 2 ) 0 y γ d / 2 3 e x y γ / 2 s / y d y = s ( s 1 ) γ x d / 2 2 d / 2 + 1 Γ ( d / 2 ) K ( γ d 5 ) / 2 * * ( x / 2 , s ) , J 4 * ( x ; γ , d ) = s γ x d / 2 2 d / 2 Γ ( d / 2 ) 0 y γ d / 2 3 e x y γ / 2 s / y Q 1 n y d y ,
where Q 1 ( y ) is defined in (33). Suppose γ = 1 . We get with (65)
J 4 ( x ; 1 , 1 ) = ( s 1 ) x ( 2 s x + 1 ) 4 s w 1 ( x ; 2 s ) and J 4 ( x ; 1 , 3 ) = ( s 1 ) x 2 2 2 s x w 3 ( x ; 2 s )
For γ = 1 using (65), we see that
J 4 ( x ; 1 , d ) = s ( s 1 ) x d / 2 2 s x 2 1 + d / 2 Γ ( d / 2 ) M ( d + 4 ) / 2 ( x + 2 s ) / 2 = ( s 1 ) x ( 2 + d ) 2 ( x + 2 s ) v d / 2 ( x ; 2 s ) .
In Lemma A4, sup x n 1 J 4 * ( x ; γ , d ) c ( s ) n 2 for γ = ± 1 is proved.
If γ = 0 , then J 3 ( x ; 0 , d ) = G d ( x ) ( h 2 ; s ( ) lim y 0 h 2 ; s ( y ) ) = 0 .
Combining the above estimates proves Theorem 2. □

9. Conclusions

Chebyshev–Edgeworth expansions are derived for the distributions of various statistics from samples with random sample sizes. The construction of these asymptotic expansions is based on the given asymptotic expansions for the distributions of statistics of samples with a fixed sample sizes as well as those of the distributions of the random sample sizes.
The asymptotic laws are scale mixtures of the underlying standard normal or chi-square distributions with gamma or inverse exponential mixing distributions. The results hold for a whole family of asymptotically normal or chi-squared statistics since a formal construction of asymptotic expansions are developed. In addition to the random sample size, a normalization factor for the examined statistics also has a significant influence on the limit distribution. As limit laws, Student, standard normal, Laplace, inverse Pareto, generalized gamma, generalized Laplace and weighted sums of generalized gamma distributions occur. As statistica the random mean, the scale-mixed normalized Student t-distribution and the Student’s t-statistic under non-normality with normal limit law, as well as Hotelling’s generalized T 0 2 and scale mixture of chi-squared statistics with chi-square limit laws, are considered. The bounds for the corresponding residuals are presented in terms of inequalities.

Author Contributions

Conceptualization, G.C. and V.V.U.; methodology, V.V.U. and G.C.; writing—original draft, G.C. and V.V.U.; and writing—review and editing, V.V.U. and G.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was done within the framework of the Moscow Center for Fundamental and Applied Mathematics, Lomonosov Moscow State University and University Basic Research Programs.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the managing editor for the assistance and the reviewers for their careful reading of the manuscript and the relevant comments. Their constructive feedback helped to improve the quality of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Auxiliary Statements and Lemmas

In Section 3, we consider statistics satisfying (1) in Assumption 1 and in Section 4 random sample sizes satisfying (2) in Assumption 2. The statistics T m in (15), (18) and (21) satisfy Assumption 1 with the normal limit distribution Φ ( x ) and in (26) and (27) with chi-square limit distributions G d ( x ) and G 4 ( x ) , defined in (17), respectively.
Further, we estimate the functions f k ( x y γ ) and y f k ( x y γ ) y k / 2 , k = 1 , 2 that appear in ( 1 ) of Assumption 1 and in the term D n in ( 5 ) . Since the functions f k ( z ) are products of a polynomial P k ( z ) and a density function p ( z ) with p ( z ) = φ ( z ) or p ( z ) = g r , r ( z ) , it follows for γ { ± 1 / 2 , ± 1 } that, if
f k ( x y γ ) = P k ( x y γ ) p ( x y γ ) then y P k ( x y γ ) p ( x y γ ) y k / 2 = Q k ( x y γ ) p ( x y γ ) y 1 + k / 2 ,
with some polynomial Q k ( z ) . For γ = 0 , we have f k ( x y γ ) = f k ( x ) and Q k ( x ) = ( k / 2 ) f k ( x ) . Hence, (A1) also holds for γ = 0 .
For example, with f 1 ( z ) = λ 3 6 ( z 2 1 ) φ ( z ) occurring in (16) for the sample mean X ¯ m and f 2 ( z ) = ( A z B z 2 ) g d ( z ) occurring in (27) with d = 4 for scale mixture of chi-square statistics, we obtain
y f 1 ( x y γ ) y 1 / 2 = Q 1 ( x y γ ) φ ( x y γ ) y 3 / 2 and y f 2 ( x y γ ) y = Q 2 ( x y γ ) g d ( x y γ ) y 2
with Q 1 ( z ) = λ 3 ( 1 / 2 + ( 2 γ 1 / 2 ) z 2 γ z 4 ) / 6 and Q 2 ( z ) = ( B z 3 ( B d 2 ) z 2 + A ( d 2 ) z ) / 4 .
Remark A1.
If P k ( 0 ) 0 , i.e., the absolute term of the polynomial P k ( z ) is not equal to zero, then it is also the absolute term of Q k ( z ) , i.e., Q k ( 0 ) 0 .
The functions φ ( z ) and e r z in g r , r ( z ) allow obtaining the estimates for
c k * = sup z | f k ( z ) | < and c k * * = sup z | Q k ( z ) | p ( z ) < , k = 1 , 2 .

Appendix A.1. Lemmas A1 and A2

Lemma A1.
Consider the statistics estimated in (15), (18), (21), (26) and (27). Let g n be a sequence with 0 < g n as n and γ { 1 / 2 , 0 , 1 / 2 , 1 } . Then, with some computable constant 0 < C * ( γ ) < , we obtain
D n = sup x 1 / g n y F ( x y γ ) + f 1 ( x y γ ) g n y + f 2 ( x y γ ) y g n d y C * ( γ ) ,
where F ( x ) , f 1 ( x ) and f 2 ( x ) are defined in the approximation estimates (15), (18), (21), (26) and (27).
Proof of Lemma A1.
The statistics in (15), (18) and (21) satisfy Assumption 1 with the normal limit distribution Φ ( x ) . To estimate D n = sup x | D n ( x ) | , we consider the cases x 0 and x = 0 .
Let x 0 . Since y Φ ( x y γ ) = 0 for γ = 0 and y Φ ( x y γ ) = γ x y γ 1 φ ( x y γ ) d y has constant sign ( γ x ) for γ = ± 1 / 2 and y > 0 , we find
1 / g n y Φ ( x y γ ) d y = 1 / g n γ x y γ 1 φ ( x y γ ) d y = 1 Φ ( x g n γ ) 1 / 2 for x > 0 . Φ ( x g n γ ) 1 / 2 for x < 0 .
From (A1) and (A2), it follows that y f k ( x y γ ) ( g n y ) k / 2 2 c k * * / k , which proves the first case. Moreover, D n ( 0 ) = | Q 1 ( 0 ) | since Q 2 ( 0 ) = 0 for the considered statistics.
Consider now the statistics estimated in (26) and (27) with limit chi-square distributions. We only need to examine x > 0 and γ { 0 , ± 1 } . In the cases now under review, we have f 1 ( x ) = 0 and f 2 ( z ) = ( A z + B z 2 ) g d ( z ) with some real constants A and B. The proof is completed with (A1) and (A2),
1 / g n y G d ( x y γ ) d y = 1 / g n y G d ( x y γ ) d y 1 and 1 / g n y f 2 ( x y γ ) g n y d y c 2 * * . □
Next, the integrals in (10) and (11) in Proposition 2 for the gamma limit distributions H ( y ) = G r , r ( y ) and the inverse exponential limit distribution H ( y ) = exp { s / y } are estimated.
Lemma A2.
(i) The conditions (6), (7) and (8) in Proposition 2 are satisfied for G r , r ( y ) and W s ( y ) = e s / y .
(ii) Let γ { 0 , ± 1 / 2 } . Consider f 1 ( x ) given in (16) and f 2 ( z ) is given in (16) or (18) for statistics occurring in (15), (18) and (21) with limiting distribution Φ ( x ) .
(iia) Let the mixing distribution be H ( y ) = G r , r ( y ) with g n = r ( n 1 ) + 1 , b = min { r , 2 } and h 2 ( y ) = h 2 ; r ( y ) = g r , r ( y ) ( y 1 ) ( 2 r ) + 2 Q 1 g n y , y > 0 . Then, we obtain with k = 1 , 2
sup x | I 1 ( x , n ) | sup x 1 / g n f 1 ( x y γ ) ( g n y ) 1 / 2 d G r , r ( y ) c 1 * g n r f o r 0 < r < 1 / 2 , γ { 0 , ± 1 / 2 } ,
sup x | I 1 ( x , n ) | c 1 * g n 1 / 2 ln g n f o r r = 1 / 2 , γ = ± 1 2 ,
sup x | I 2 ( x , n ) | sup x 1 / g n f 2 ( x y γ ) g n y d G r , r ( y ) c 2 * g n r f o r 0 < r < 1 , γ { 0 , ± 1 / 2 } ,
sup x | I 2 ( x , n ) | c 3 g n 1 f o r r = 1 , γ = ± 1 / 2 ,
sup x | I k ( x , n ) g n k / 2 f k ( x ) ln g n | c 4 g n k / 2 f o r r = k / 2 , γ = 0 ,
sup x | I 2 + k ( x , n ) | = sup x 1 / g n f k ( x y γ ) n g n k / 2 y d h 2 ; r ( y ) c 5 ( r ) g n r , r > 1 , r 1 + k / 2 , c 6 ( r ) g n 1 k / 2 ln n , r = 1 + k / 2 .
(iib) Now, consider the mixing distribution H ( y ) = W s ( y ) = e s / y with g n = n , b = 2 and h 2 ( y ) = h 2 ; s ( y ) = s e s / y s 1 + 2 Q 1 ( n y ) / 2 y 2 , y > 0 . Then, apply it to I 4 ( x , n ) in (11)
sup x | I 4 ( x , n ) | = sup x 1 / n f 2 ( x y γ ) n 2 y d h 2 ; s ( y ) c 6 ( s ) n 2 , f o r s > 0 , γ { 0 , ± 1 / 2 } .
(iii) Put γ { 0 , ± 1 } . Consider f 2 ( x ) = ( A x + B x 2 ) g d ( x ) for statistics occurring in (26) and (27) satisfying Assumption 1 with the chi-square distribution G d ( x ) .
(iiia) The mixing distribution is H ( y ) = G r , r ( y ) as in Case (ia) above. Then, for γ { 0 , ± 1 } ,
sup x | I 2 ( x , n ) | sup x 1 / g n f 2 ( x y γ ) g n y d G r , r ( y ) c 2 * r r ( r 1 ) Γ ( r ) g n r f o r r < 1 ,
sup x | I 2 ( x , n ) | c 2 * n 1 w i t h γ = ± 1 / 2 , sup x | I 2 ( x , n ) n 1 f 2 ( x ) ln n | c 2 * n 1 w i t h γ = 0 , f o r r = 1 ,
sup x | I 4 ( x , n ) | = sup x 1 / g n f 2 ( x y γ ) n g n y d h 2 ; r ( y ) c 3 ( r ) g n min { r , 2 } , f o r r > 1 , r 2 , c 4 ( r ) g n 2 ln n , f o r r = 2 .
(iiib) The mixing distribution is H ( y ) = W s ( y ) = e s / y with g n = n and b = 2 as in Case (iib). Then,
sup x | I 4 ( x , n ) | = sup x 1 / g n f 2 ( x y γ ) n 2 y d h 2 ; s ( y ) C 5 ( r ) n 2 . f o r s > 0 , γ { 0 , ± 1 / 2 } .
Proof of Lemma A2.
(i) Insertion of G r , r ( y ) with h 2 , r ( y ) and W s ( y ) = e s / y with h 2 , s ( y ) and simple calculation result in the necessary estimates in (6)–(8). In the case of W s ( y ) = e s / y , one even gets for all terms exponentially fast decrease.
(ii) The limit distribution of the considered statistics is standard normal Φ ( x ) .
(iia) Let H ( x ) = G r , r ( x ) . Using (A2), the estimations (A3) and (A5) for r < k / 2 , with k = 1 , 2 , are
sup x | I k ( x , n ) | c k * r r g n k / 2 Γ ( r ) 1 / g n y r 1 k / 2 d y c k * r r ( k / 2 r ) Γ ( r ) g n r .
Taking into account
0 ln g n 1 / g n 1 e r y y d y = 1 / g n 1 1 e r y y d y r and 1 e r y y d y e r / r f o r r > 0 ,
the bound (A4) follows from
| I 1 ( x , n ) | c 1 * ( 2 g n ) 1 / 2 Γ ( 1 / 2 ) 1 e y / 2 y d y + 1 / g n 1 e y / 2 y d y c 1 * ( 2 e 1 / 2 + ln g n ) Γ ( 1 / 2 ) ( 2 g n ) 1 / 2 .
If r = 1 with d 2 * = sup z { | z 1 f 2 ( z ) | φ ( z / 2 ) } , we find | f 2 ( z ) | d 2 * | z | φ ( z / 2 ) and the bound (A6) follows from
| I 2 ( x , n ) | d 2 * | x | 2 π n 1 / n y γ 1 e ( y + x 2 y 2 γ / 4 ) d y with γ = ± 1 / 2 ,
where for γ = 1 / 2 using | x | ( 1 + x 2 / 4 ) 1 / 2 2 , we obtain
| I 2 ( x , n ) | d 2 * | x | 2 π n 1 / n y 1 / 2 1 e ( 1 + x 2 / 4 ) y d y d 2 * | x | Γ ( 1 / 2 ) 2 π ( 1 + x 2 / 4 ) 1 / 2 n 1 2 d 2 * n
and, in the case of γ = 1 / 2 , the substitution z = x 2 / ( 4 y ) for x 0 leads to
I 2 ( x , n ) c 2 * | x | 2 π n 1 / n y 1 1 / 2 e ( y + x 2 / ( 4 y ) ) d y 2 c 2 * 2 π n 0 z 1 / 2 e z d z 2 d 2 * n .
Finally, if γ = 0 , then f k ( x y γ ) = f k ( x ) does not depend on y. Then, (A7) follows from (A1), (A2), and (A14) for r = k / 2 .
Let r > 1 . Integration by parts for Lebesgue–Stieltjes integrals I k + 2 ( x , n ) , k = 1 , 2 , in (11) leads to
sup x I k + 2 ( x , n ) 1 n g n k / 2 c k * | h 2 ; r ( 1 / g n ) | + sup x 1 / g n | Q k ( x y γ ) | y 1 + k / 2 | h 2 ; r ( y ) | d y
with bound c k * given in (A2). Defining C r * = r r 2 r Γ ( r ) sup y { e r y ( | y 1 | | 2 r | + 1 ) } < , we find
1 / g n | h 2 ; r ( y ) | y 1 + k / 2 d y C r * 1 / g n y r 2 k / 2 d y = C r * ( 1 + k / 2 r ) g n r + 1 + k / 2 for 1 < r < 1 + k / 2
and, with C r * * = r r 1 2 Γ ( r ) sup y { ( e r y / 2 ( | y 1 | | 2 r | + 1 ) } < , we obtain
1 / g n | h 2 ; r ( y ) | y 1 + k / 2 d y C r * * 1 / g n y r 2 k / 2 e r y / 2 d y C r * * Γ ( r 1 k / 2 ) ( r / 2 ) r 1 k / 2 for r > 1 + k / 2 .
Hence, using g n r n for r > 1 , we obtain (A8) and, for r > 1 , r 1 + k / 2 .
For r = 1 + k / 2 , the second integral in the line above is an exponential integral. Therefore, with (A14), we find (A8) for r = 1 + k / 2 , too.
(iib) The mixing distribution is H ( y ) = W s ( y ) = e s / y with g n = n . Since b = 2 , only I 4 ( x , n ) has to be estimated. Integration by parts for I 4 ( x , n ) in (11) leads to (A15) with k = 2 , g n = n and h 2 , s ( y ) instead of h 2 , r ( y ) . Hence, (A9) follows from
1 / n | h 2 ; s ( y ) | y 2 d y s ( s + 2 ) 1 / n y 4 e s / y d y s + 2 s 2 0 s n z 2 e z d z ( s + 2 ) Γ ( 3 ) s 2 .
(iii) The limit distribution of statistics in (26) and (27) is chi-square distribution G u ( x ) defined in (17). In the considered cases f 1 ( x ) = 0 . Let γ { 0 , ± 1 } . Consider f 2 ( x ) = ( A x + B x 2 ) g d ( x ) with chi-square density g d ( x ) .
(iiia) Let H ( y ) = G r , r ( y ) . We have to estimate I 2 ( x , n ) for r 1 and I 4 ( x , n ) for r > 1 . The bound (A10) for 0 < r < 1 follows from (A2) and
sup x | I 2 ( x , n ) | c 2 * r r g n Γ ( r ) 1 / g n y r 2 e r y d y c 2 * r r ( r 1 ) Γ ( r ) g n r for γ { 0 , ± 1 } .
If r = 1 with C 2 * = sup z | A + B z | 1 2 d / 2 Γ ( d / 2 ) e z / 4 < , we find | f 2 ( z ) | C 2 * z d / 2 e z / 4 and
| I 2 ( x , n ) | C 2 * x d / 2 g n 1 / g n y 1 + d / 2 e ( 1 + x / 4 ) y d y C 2 * x d / 2 g n ( 1 + x / 4 ) d / 2 4 d / 2 C 2 * g n for γ = 1 .
in the case of γ = 1 using variable transformation z = x / ( 4 y ) for x > 0 one has
| I 2 ( x , n ) | C 2 * x d / 2 g n 1 / g n e x / ( 4 y ) y 1 + d / 2 d y C 2 * x d / 2 g n ( x / 4 ) d / 2 0 z 1 + d / 2 e z d z C 2 * 4 d / 2 Γ ( d / 2 ) g n .
If γ = 0 then f 2 ( x y γ ) = f 2 ( x ) , noting (A14) and g n = n for r = 1 , we prove (A11).
Let now r > 1 . It remains to estimate I 4 ( x , n ) . Using (A15) with k = 2 , remembering (A2), we obtain (A12) in the same way as for r > 1 with k = 2 in case (iia) above.
(iiib) The limit distribution H ( y ) = W s ( y ) = e s / y with g n = n and b = 2 . As in Case (iib), taking into consideration (A16), we obtain (A13). □

Appendix A.2. Lemmas A3 and A4

We show that the integrals J 4 * ( x ; γ ) and J 4 * ( x ; γ , d ) in the proofs of Theorems 1 and 2 have the order of the remaining terms. Therefore, the involved jump correcting function Q 1 ( y ) = 1 / 2 ( y [ y ] ) occurring in (32) and (38) has no effect on the second approximation. The function Q 1 ( y ) is periodic with period 1. The Fourier series expansion of Q 1 ( y ) at all non-integer points y is
Q 1 ( y ) = 1 / 2 ( y [ y ] ) = k = 1 sin ( 2 π k y ) k π y [ y ]
(see formula 5.4.2.9 in [40] with a = 0 ).
Lemma A3.
Let J 4 * ( x ; ± 1 ) be defined by (76) and (78), respectively. Then, n 1 J 4 * ( x ; 1 / 2 ) C n r for r > 1 and n 1 J 4 * ( x ; 1 / 2 ) C n 2 .
Proof of Lemma A3.
We begin by considering Q 1 ( y ) in J 4 * ( x ; 1 / 2 ) defined in (A17) following the estimate of J 4 * ( x ) in the proof of Theorem 2 in [11]. Inserting Fourier series expansion of Q 1 ( y ) into the integral J 4 * ( x ; 1 / 2 ) , interchanging the integral and sum and applying formula (2.5.31.4) in [40] with α = r 1 / 2 , p = ( r + x 2 / 2 ) and b = 2 π k g n , then
J 4 * ( x ; 1 / 2 ) = x r r 1 ( 2 π ) 3 / 2 Γ ( r ) k = 1 1 k 0 y r 3 / 2 e ( r + x 2 / 2 ) y sin 2 π k g n y d y = x r r 1 Γ ( r 1 / 2 ) ( 2 π ) 3 / 2 Γ ( r ) k = 1 sin ( r 1 / 2 ) arctan 4 π k g n / ( x 2 + 2 r ) k 2 π k g n 2 + ( r + x 2 / 2 ) 2 ( r 1 / 2 ) / 2 = r r 1 Γ ( r 1 / 2 ) 2 π 2 π Γ ( r ) k = 1 a k ( x ; n ) k .
Now, we split the exponent ( r 1 / 2 ) / 2 = ( r 1 ) / 2 + 1 / 4 and obtain
| a k ( x ; n ) | | x | 2 π k g n 2 + r + x 2 / 2 2 ( r 1 ) / 2 + 1 / 4 | x | ( 2 π k g n ) r 1 ( r + x 2 / 2 ) 1 / 2 2 ( 2 π k ( n 1 ) ) r 1 .
Since r > 1 and n 2 , the first statement in Lemma A3 follows:
sup x n 1 | J 4 * ( x ; 1 / 2 ) | c ( r ) n r k = 1 k r = c 1 ( r ) n r .
To prove the second statement about J 4 * ( x ; 1 / 2 ) , we insert again the Fourier series expansion of Q 1 ( y ) given in (A17) into J 4 * ( x ; 1 / 2 ) and interchange the integral and sum
J 4 * ( x ; 1 / 2 ) = x 2 2 π k = 1 1 π k 0 y 1 / 2 e ( 2 y + x 2 / ( 2 y ) sin ( 2 π k g n y ) d y .
Further, we use formula 2.5.37.4 in [40]
0 y 1 / 2 e p y q / y sin ( b y ) d y = π p 2 + b 2 e 2 q z + z + sin ( 2 q z ) + z cos ( 2 q z )
with 2 z ± 2 = p 2 + b 2 ± p , p = 2 > 0 , q = x 2 / 2 > 0 and b = 2 π g n k > 0 . Use of the estimates
0 < z z + , | x | z + e | x | z + e 1 , p 2 + b 2 b = 2 π g n k a n d k = 1 k 2 = π 2 / 6
leads to the inequalities
sup x 1 n | J 4 * ( x ; 1 / 2 ) | sup x 1 2 2 π k = 1 π 2 | x | z + 2 π 2 g n n k 2 e 2 | x | 1 e 2 12 g n n = C n 2
and Lemma A3 is proven. □
Lemma A4.
Let J 4 * ( x ; γ , d ) be defined by (79), then n 1 J 4 * ( x ; γ , d ) C n 2 for γ = ± 1 .
Proof of Lemma A4.
Using the Fourier series expansion (A17) of the periodic function Q 1 ( y ) , given in (33), and interchange integral and sum, we find
J 4 * ( x ; γ , d ) = s γ x d / 2 2 d / 2 Γ ( d / 2 ) k = 1 1 k π 0 y ( γ d 6 ) / 2 e x y γ / 2 s / y sin ( 2 π k n y ) d y .
We begin by estimating J 4 * ( x ; 1 , 3 ) , i.e., the exponent by y in (A19) is 3 / 2 . Thus, we can use formula 2.5.37.3 in [40]
0 y 3 / 2 e p y s / y sin ( b y ) d y = π s e 2 s z + sin ( 2 s z )
where p = x / 2 > 0 , s > 0 , b = 2 π k n > 0 and 2 z ± 2 = p 2 + b 2 ± p = x 2 / 4 + ( 2 π k n ) 2 ± x / 2 . Since
z + = 1 2 x 2 / 4 + 2 π k n 1 / 2 + x / 2 1 / 2 1 4 x 1 / 2 + 1 8 ( 2 π ) 1 / 4 k 1 / 4 + n 1 / 4
it results in
1 n J 4 * ( x : 1 , 3 ( s x ) 3 / 2 n s 2 1 / 2 π k = 1 exp ( s x ) 1 / 2 / 2 ( 2 π ) 1 / 4 k 1 / 4 + n 1 / 4 / 4 C ( s ) n 2 .
Let now γ = 1 and d = 1 . The main difference compared with the previous estimate of J 4 * ( x ; 1 , 3 ) is that we are facing more technical trouble in order to estimate J 4 * ( x ; 1 , 1 ) . The exponent by y in (A19) is 5 / 2 and we cannot find a closed formula similar to (A19) for this case. To estimate J 4 * in the proof of Theorem 5 in [11], we show that differentiation with respect to s under the integral sign in (A20) is allowed. Hence,
0 y 5 / 2 e p y s / y sin ( b y ) d y = ( π / 2 ) e 2 s z + ( s 3 / 2 sin ( 2 s z ) + 2 s 1 z + sin ( 2 s z ) 2 s 1 z cos ( 2 s z ) ) .
with the same coefficients p, s, b and z ± as in (A20). The use of (A21) and the obvious inequalities 0 < z z + and z + 1 2 z + + 1 8 x 1 / 2 + 1 16 ( 2 π ) 1 / 4 ( k 1 / 4 leads to
1 n J 4 * ( x ; 1 , 1 ) ( s x ) 1 / 2 n 2 1 / 2 π e ( s x ) 1 / 2 / 4 ( 2 π ) 1 / 4 n 1 / 4 / 8 k = 1 1 k e ( 2 π ) 1 / 4 k 1 / 4 / 8 C ( s ) n 2 .
Finally, let now γ = 1 .
J 4 * ( x ; 1 , d ) = s x d / 2 2 d / 2 Γ ( d / 2 ) k = 1 1 k π 0 y ( d / 2 + 3 ) e ( x / 2 + s ) / y sin ( 2 π k n y ) d y .
Partial integration in the integral with A = d / 2 + 3 , B = x / 2 + s and C = 2 π k n leads to
0 y A e B / y sin ( C y ) d y = 0 1 C A y ( A + 1 ) + B y ( A + 2 ) e B / y cos ( C y ) d y
and using (62) to
0 1 C A y ( A + 1 ) + B y ( A + 2 ) e B / y d y 1 C A Γ ( A ) B A + B Γ ( A + 1 ) B A + 1 = Γ ( A + 1 ) C B A .
Therefore,
1 n | J 4 * ( x ; 1 , d ) | s x d / 2 n 2 d / 2 Γ ( d / 2 ) k = 1 Γ ( d / 2 + 4 ) k 2 2 π 2 n ( x / 2 + s ) d / 2 + 3 C ( s ) n 2 ,
and Lemma A4 is proven. □

References

  1. Nunes, C.; Capistrano, G.; Ferreira, D.; Ferreira, S.S.; Mexia, J.T. Exact critical values for one-way fixed effects models with random sample sizes. J. Comput. Appl. Math. 2019, 354, 112–122. [Google Scholar] [CrossRef]
  2. Nunes, C.; Capistrano, G.; Ferreira, D.; Ferreira, S.S.; Mexia, J.T. Random sample sizes in orthogonal mixed models with stability. Comp. Math. Methods 2019, 1, e1050. [Google Scholar] [CrossRef] [Green Version]
  3. Nunes, C.; Mário, A.; Ferreira, D.; Moreira, E.M.; Ferreira, S.S.; Mexia, J.T. An algorithm for simulation in mixed models with crossed factors considering the sample sizes as random. J. Comput. Appl. Math. 2021. [Google Scholar] [CrossRef]
  4. Esquível, M.L.; Mota, P.P.; Mexia, J.T. On some statistical models with a random number of observations. J. Stat. Theory Pract. 2016, 10, 805–823. [Google Scholar] [CrossRef]
  5. Al-Mutairi, J.S.; Raqab, M.Z. Confidence intervals for quantiles based on samples of random sizes. Statist. Pap. 2020, 61, 261–277. [Google Scholar] [CrossRef]
  6. Barakat, H.M.; Nigm, E.M.; El-Adll, M.E.; Yusuf, M. Prediction of future generalized order statistics based on exponential distribution with random sample size. Statist. Pap. 2018, 59, 605–631. [Google Scholar] [CrossRef]
  7. Gnedenko, B.V. Estimating the unknown parameters of a distribution with a random number of independent observations. (Probability theory and mathematical statistics (Russian)). Trudy Tbiliss. Mat. Inst. Razmadze Akad. Nauk Gruzin. SSR 1989, 92, 146–150. [Google Scholar]
  8. Gnedenko, B.V.; Korolev, V.Y. Random Summation. Limit Theorems and Applications; CRC Press: Boca Raton, FL, USA, 1996. [Google Scholar]
  9. Bening, V.E.; Galieva, N.K.; Korolev, V.Y. On rate of convergence in distribution of asymptotically normal statistics based on samples of random size. Ann. Math. Inform. 2012, 39, 17–28. (In Russian) [Google Scholar]
  10. Bening, V.E.; Galieva, N.K.; Korolev, V.Y. Asymptotic expansions for the distribution functions of statistics constructed from samples with random sizes. Inform. Appl. 2013, 7, 75–83. (In Russian) [Google Scholar]
  11. Christoph, G.; Monakhov, M.M.; Ulyanov, V.V. Second-order Chebyshev-Edgeworth and Cornish-Fisher expansions for distributions of statistics constructed with respect to samples of random size. J. Math. Sci. 2020, 244, 811–839. [Google Scholar] [CrossRef]
  12. Christoph, G.; Ulyanov, V.V.; Bening, V.E. Second Order Expansions for Sample Median with Random Sample Size. arXiv 2020, arXiv:1905.07765v2. [Google Scholar]
  13. Christoph, G.; Ulyanov, V.V. Second order expansions for high-dimension low-sample-size data statistics in random setting. Mathematics 2020, 8, 1151. [Google Scholar] [CrossRef]
  14. Christoph, G.; Ulyanov, V.V. Short Expansions for High-Dimension Low-Sample-Size Data Statistics in a Random Setting. Recent Developments in Stochastic Methods and Applications. In Proceedings in Mathematics & Statistics; Shiryaev, A.N., Samouylov, K.E., Kozyrev, D.V., Eds.; Springer: Cham, Switzerland, 2021; (to appear). [Google Scholar]
  15. Hall, P. The Bootstrap and Edgeworth Expansion; Springer Series in Statistics; Springer: New York, NY, USA, 1992. [Google Scholar]
  16. Bickel, P.J. Edgeworth expansions in nonparametric statistics. Ann. Statist. 1974, 2, 1–20. [Google Scholar] [CrossRef]
  17. Kolassa, J.E. Series Approximation Methods in Statistics, 3rd ed.; Lecture Notes in Statistics 88; Springer: New York, NY, USA, 2006. [Google Scholar]
  18. Petrov, V.V. Limit Theorems of Probability Theory, Sequences of Independent Random Variables; Clarendon Press: Oxford, UK, 1995. [Google Scholar]
  19. Fujikoshi, Y.; Ulyanov, V.V. Non-Asymptotic Analysis of Approximations for Multivariate Statistics; Springer: Singapore, 2020. [Google Scholar]
  20. Hall, P. Edgeworth Expansion for Student’s t-statistic under minimal moment conditions. Ann. Probab. 1987, 15, 920–931. [Google Scholar] [CrossRef]
  21. Bentkus, V.; Götze, F.; van Zwet, W.R. An Edgeworth expansion for symmetric statistics. Ann. Statist. 1997, 25, 851–896. [Google Scholar] [CrossRef]
  22. Anderson, T.W. An Introduction to Multivariate Statistical Analysis, 3rd ed.; Wiley Series in Probability and Mathematical Statistics; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2003. [Google Scholar]
  23. Fujikoshi, Y.; Ulyanov, V.V.; Shimizu, R. Multivariate Statistics. High-Dimensional and Large-Sample Approximations; Wiley Series in Probability and Statistics; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2010. [Google Scholar]
  24. Fujikoshi, Y.; Ulyanov, V.V.; Shimizu, R. L1-norm error bounds for asymptotic expansions of multivariate scale mixtures and their applications to Hotelling’s generalized T 0 2 . J. Multivar. Anal. 2005, 96, 1–19. [Google Scholar] [CrossRef] [Green Version]
  25. Ulyanov, V.V.; Fujikoshi, Y. On accuracy of improved χ2-approximation. Georgian Math. J. 2001, 8, 401–414. [Google Scholar] [CrossRef]
  26. Schluter, C.; Trede, M. Weak convergence to the Student and Laplace distributions. J. Appl. Probab. 2016, 53, 121–129. [Google Scholar] [CrossRef]
  27. Döbler, C. New Berry-Esseen and Wasserstein bounds in the CLT for non-randomly centered random sums by probabilistic methods. ALEA Lat. Am. J. Probab. Math. Stat. 2015, 12, 863–902. [Google Scholar]
  28. Robbins, H. The asymptotic distribution of the sum of a random number of random variables. Bull. Am. Math. Soc. 1948, 54, 1151–1161. [Google Scholar] [CrossRef] [Green Version]
  29. Kolassa, J.E.; McCullagh, P. Edgeworth Series for Lattice Distributions. Ann. Statist. 1990, 18, 981–985. [Google Scholar] [CrossRef]
  30. Bening, V.E.; Korolev, V.Y. On the use of Student’s distribution in problems of probability theory and mathematical statistics. Theory Probab. Appl. 2005, 49, 377–391. [Google Scholar] [CrossRef]
  31. Gavrilenko, S.V.; Zubov, V.N.; Korolev, V.Y. The rate of convergence of the distributions of regular statistics constructed from samples with negatively binomially distributed random sizes to the Student distribution. J. Math. Sci. 2017, 220, 701–713. [Google Scholar] [CrossRef]
  32. Buddana, A.; Kozubowski, T.J. Discrete Pareto distributions. Econ. Qual. Control 2014, 29, 143–156. [Google Scholar] [CrossRef]
  33. Lyamin, O.O. On the rate of convergence of the distributions of certain statistics to the Laplace distribution. Mosc. Univ. Comput. Math. Cybern. 2010, 34, 126–134. [Google Scholar] [CrossRef]
  34. Choy, T.B.; Chan, J.E. Scale mixtures distributions in statistical modellings. Aust. N. Z. J. Stat. 2008, 50, 135–146. [Google Scholar] [CrossRef]
  35. Oldham, K.B.; Myland, J.C.; Spanier, J. An Atlas of Functions, 2nd ed.; Springer Science+Business Media: New York, NY, USA, 2009. [Google Scholar]
  36. Kotz, S.; Kozubowski, T.J.; Podgórski, K. The Laplace Distribution and Generalizations: A Revisit with Applications to Communications, Economics, Engineering, and Finance; Birkhäuser: Boston, MA, USA, 2001. [Google Scholar]
  37. Korolev, V.Y.; Zeifman, A.I. Generalized negative binomial distributions as mixed geometric laws and related limit theorems. Lith. Math. J. 2019, 59, 366–388. [Google Scholar] [CrossRef] [Green Version]
  38. Korolev, V.Y.; Gorshenin, A. Probability models and statistical tests for extreme precipitation based on generalized negative binomial distributions. Mathematics 2020, 8, 604. [Google Scholar] [CrossRef] [Green Version]
  39. Safari, M.A.M.; Masseran, N.; Ibrahim, K.; Hussain, S.I. A robust and efficient estimator for the tail index of inverse Pareto distribution. Phys. A Stat. Mech. Its Appl. 2019, 517, 431–439. [Google Scholar] [CrossRef]
  40. Prudnikov, A.P.; Brychkov, Y.A.; Marichev, O.I. Integrals and Series, Vol. 1: Elementary Functions, 3rd ed.; Gordon & Breach Science Publishers: New York, NY, USA, 1992. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Christoph, G.; Ulyanov, V.V. Chebyshev–Edgeworth-Type Approximations for Statistics Based on Samples with Random Sizes. Mathematics 2021, 9, 775. https://doi.org/10.3390/math9070775

AMA Style

Christoph G, Ulyanov VV. Chebyshev–Edgeworth-Type Approximations for Statistics Based on Samples with Random Sizes. Mathematics. 2021; 9(7):775. https://doi.org/10.3390/math9070775

Chicago/Turabian Style

Christoph, Gerd, and Vladimir V. Ulyanov. 2021. "Chebyshev–Edgeworth-Type Approximations for Statistics Based on Samples with Random Sizes" Mathematics 9, no. 7: 775. https://doi.org/10.3390/math9070775

APA Style

Christoph, G., & Ulyanov, V. V. (2021). Chebyshev–Edgeworth-Type Approximations for Statistics Based on Samples with Random Sizes. Mathematics, 9(7), 775. https://doi.org/10.3390/math9070775

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop