Next Article in Journal
The Existence, Transcendence, and Evolution of the Subject—A Method Based on Subject Information
Previous Article in Journal
How Digital Technology Can Reshape the Trust System of Engineering—Taking Beijing Daxing International Airport as an Example
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Proceeding Paper

On Estimation of the Remainder Term in New Asymptotic Expansions in the Central Limit Theorem †

1
Department of Probability Theory and Applied Mathematics, Moscow Technical University of Communications and Informatics, 123423 Moscow, Russia
2
Mechanical Engineering Research Institute RAS, 101990 Moscow, Russia
*
Author to whom correspondence should be addressed.
Presented at the 1st International Online Conference on Mathematics and Applications, 1–15 May 2023; Available online: https://iocma2023.sciforum.net/.
Comput. Sci. Math. Forum 2023, 7(1), 60; https://doi.org/10.3390/IOCMA2023-14423
Published: 28 April 2023

Abstract

:
We offer a new asymptotic expansion with an explicit remainder estimate in the central limit theorem. The results obtained are essentially based on new forms of asymptotic expansions in the central limit theorem. We also present a more accurate estimation of the CLT-expansions remainder which is rigorously proved and backed up numerically. It is shown that our approach can be used for further refinement of allied asymptotic expansions.

1. Introduction

The central limit theorem (CLT) states that under fairly broad conditions the sum of independent (or weakly dependent) identically distributed (i.i.d.) random variables is approximately normally distributed. In this paper, we deal with CLT for i.i.d. variables  X 1 , X 2 ,   . . .  with zero mean and unit variance each. Let  X 1  follow a distribution  P  with cumulative distribution function  F ( x )  and characteristic function  f t .
We assume that  X 1  has a finite absolute moment of order  m + 2  and for some  ν > 0  the function  f t ν  is integrable on  R . Since for any  T > 0  the integral
T + f t ν d t < +
converges and
α T = sup f t : t T < 1
(see [1], p. 43), it follows that for  n ν  the distribution  P n  of the rescaled sum  X 1 + . . . + X n n 1 2  has a continuous density  p n x . In it is its turn, this implies that for this distribution  P n  (with cumulative distribution function  F n x ) the CLT holds [1],  p n x φ x  as  n +  for any  x R ; here  φ x = 1 2 π e x 2 / 2  is the density of the standard normal distribution  Φ x . The density  p n x  can be represented as the inverse Fourier transform [1] (pp. 42, 147) of the characteristic function  f n  of the convolution of n copies of the original distribution.
Here the question naturally arises concerning the accuracy of the CLT-approximation. One of the main results on the subject is the Berry–Esseen theorem (see, for example, [2]). Nevertheless, this theorem is too general (and not attempting to take any specific properties of the original distribution into play) and, therefore, provides rather crude estimates of the convergence rate [3].
One way to improve the accuracy of CLT-approximations is to use asymptotic expansions. Until quite recently, most of such expansions only gave an estimate of the order of approximation at best and thus were of little use for, say, numeric computations.
Yu.V. Prohorov wrote [4] (p. 7) that “it was P.L. Chebyshev’s idea to explore the asymptotic behavior of the difference  F n ( x ) Φ ( x )  and it was them to give a formal expansion of the difference”. A number of such expansions under different restrictions on the original distribution were obtained later by H. Cramer [5] and investigated by C.-G. Esseen [6]. Moreover, H. Cramer [7] claimed that the featured series expansions were introduced by F. Edgeworth [8]. For example, if there exists an integer  m 1 , such that  M X 1 m + 2 <  and Cramer’s (C) condition  lim t sup f ( t ) < 1  is fulfilled, then
F n ( x ) = Φ x + k = 1 m P k ( Φ ) n k / 2 + O 1 n m / 2 , n ,
where  P k ( Φ ) = L 3 k 1 ( x ) φ ( x )  and  L 3 k 1 ( x )  is a polynomial of degree  3 k 1  in x. Explicit formulas for  P k Φ  in terms of semi-invariants were obtained by V.V. Petrov [9] in 1962.
Seeking more efficient forms of such expansions V.V. Senatov took what can be justly characterized a revolutionary step; he offered expansions that allowed for explicit estimation of the remainder (not just barely indicating big-Oh approximation order). At the moment only one type of CLT-asymptotic expansions, namely the Gram–Charlier expansion, was widely known, that is,
p n x = φ x + k = 3 , k 3 m 4 3 m 3 θ k P n H k ( x ) φ ( x ) + O 1 n m / 2 , n ,
where  H k x = 1 k φ ( k ) x φ x  are Chebyshev–Hermite polynomials and
θ k P n = 1 k ! + H k x p n x d x , k 0
are normalized moments [10].
As mentioned by V.V. Senatov [1] (p. 124), already, in 1920, H. Cramer [7] noticed that this accuracy can be attained with only  m + 2  moments at hand (the Gram–Charlier expansion provides the same accuracy only when  3 m 3  moments are used). That was the reason, in Senatov’s opinion, why researchers’ primary focus was on the Edgeworth–Cramer (not Gram–Charlier) expansion.
At the same, it turns out that there exists an expansion of the Gram–Charlier type that makes use of only  m + 2  absolute moment of variable  X 1 . V.V. Senatov to construct this kind of asymptotic expansions proposed together with the moments:
θ k = j = 0 k 2 ( 1 ) j α k 2 j α 2 j ( φ ) ,
use the so-called [10] Senatov incomplete moments:
θ k ( l ) = j = 0 l 2 ( 1 ) j α k 2 j α 2 j ( φ ) , l k ,
where
α j = M X 1 j j ! , β j = M X 1 j j ! , α 2 l φ = 1 2 l l ! .
V.V. Senatov came up [11] with what he dubbed a shortened Gram–Charlier expansion:
p n ( x ) = φ ( x ) + k = 3 m + 1 θ k P n H k x φ x + k = m + 2 , k 3 m 4 3 m 3 θ k ( m + 1 ) P n H k x φ x + R .
He used arbitrary measures (not necessarily probabilistic measures) in their derivation, which imposes additional limitations on the moments of the original distribution.
In attempts to remove the limitations and make the estimation of the remainder more accurate, V.V. Senatov and V.N. Sobolev [12] suggested a novel form of asymptotic expansion that does not impose any additional restrictions on the momenta (as compared to [1]).
The Gram–Charlier expansion differs architecturally from the Edgeworth–Cramer expansion, the former is in powers of the Chebyshev–Hermite polynomials, while the latter is in powers of n. Expansions of the third type [12,13] are obtained as follows, the terms in such an expansion are to be ordered with respect to the number of factors, which are Senatov’s moments.
Before the results of V.V. Senatov and V.N. Sobolev [12,13], only two types of such expansions were widely known, the Gram–Charlier and Edgeworth–Cramer expansions. In the former, the terms are grouped in the order of the Chebyshev–Hermite polynomials, and in the latter, they are grouped in powers of n. In [12], V.V. Senatov and V.N. Sobolev proposed grouping the terms according to the number of factors of Senatov’s moments
p n ( x ) = φ ( x ) + s = 1 m 1 C n s l = 3 s m 1 + 2 s Θ s , l n l / 2 H l x ) φ ( x ) + O 1 n m / 2 , n ,
where
Θ s , l = t 1 + . . . + t s = l θ t 1 . . . θ t s
the summation is carried out over tuples of natural numbers  t 1 , …,  t s , such that  t j 3 j = 1 , , m 1  and  t 1 + + t s = l .

2. Main Result

To improve the accuracy of the expansions from [12] note first that the values
α 0 = 1 , α 1 = 0 , α 2 = 1 2 , α 2 l φ = 1 2 l l !
are known. Then, we introduce the following non-negative quantities:
θ s = β s + 2 + j = 1 s 2 2 α 2 j ( φ ) α s 2 j + j = s 2 1 s 2 ( 1 ) j α 2 j ( φ ) α s 2 j ,
for which the inequalities
θ s β s + 2 + j = 1 s 2 α 2 j ( φ ) α s 2 j
are valid. Therefore, using  θ s  in the estimates instead of the right side of the last inequality increases accuracy of asymptotic expansions.
Let  θ s l  be the abbreviated versions of  θ s : they are calculated by the same formulas in which  α k = 0  for  k > l .  There is also an improvement here.
Another improvement consists of preserving the minus signs between the summands when evaluating the estimates. The following quantities naturally arise
S 0 , 2 ( Θ ) = θ m + 2 , S 0 , 3 ( Θ ) = θ m + 3 m + 1
which for  l 1  read
S l , 2 ( Θ ) = k 1 , k 2 , . . . , k l m θ 3 + k 1 θ 3 + k 2 . . . θ 3 + k l θ m + 2 l + k 1 + k 2 + . . . + k l ,
S l , 3 ( Θ ) = k 1 , k 2 , . . . , k l m θ 3 + k 1 θ 3 + k 2 . . . θ 3 + k l θ m + 3 l k 1 k 2 . . . k l m + 2 l k 1 k 2 . . . k l .
The notation  k 1 , k 2 , . . . , k l m  implies summation over all sets  k 1 ,   k 2 ,   . . . ,   k l  of non-negative numbers, such that  0 k 1 + k 2 + . . . + k l m l 1 .
We will also use the following quantities of the moment type
L l u = 1 2 π t u t l e t 2 2 d t , B l , n k = 1 2 π T n + T n t l μ n k t n d t ,
where the function  μ ( t ) = max f ( t ) , e t 2 / 2  was introduced by V.Yu. Korolev.
Below, we formulate our main theorem that provides a CLT-asymptotic expansion with an improved estimation of the remainder under fairly general conditions.
Theorem 1.
Let identically distributed independent random variables  X 1 ,   X 2 ,   . . .  with zero mean and unit variance each follow the same distribution P. Suppose that P has a finite absolute moment of order  m + 2  and  f t ν d t < .  Where  f t  is the characteristic function of P. Then, for any  n max ν , m + 1  and for all  x R
p n x = φ x l = 0 m 1 C n l n 3 l k 1 , k 2 , . . . , k l m θ 3 + k 1 θ 3 + k 2 . . . θ 3 + k l n k 1 + k 2 + . . . + k l H 3 l + k 1 + . . . + k l x + R n , m x ,
where
R n , m x 1 n m + 2 l = 0 m 2 C n l + 1 n l B m + 2 + 2 l , n 1 S l , 2 ( Θ ) + B m + 3 + 2 l , n 1 S l , 3 ( Θ ) 1 n
+ Λ n T + Λ ¯ n T .
and two last terms in the remainder’s estimate are
Λ n T = n π α n ν T T + f t ν d t ,
Λ ¯ n T = l = 0 m 1 C n l k 1 , k 2 , . . . , k l m θ 3 + k 1 θ 3 + k 2 . . . θ 3 + k l n 3 l + k 1 + k 2 + . . . + k l L 3 l + k 1 + k 2 + . . . + k l T n ,
which decay exponentially fast.

3. Conclusions

We obtain new explicit estimates for accuracy of approximation in the CLT-expansions. Our approach can be used for further refinement of allied asymptotic expansions.

Author Contributions

Conceptualization, V.S. and S.R.; methodology, V.S. and S.R.; software, V.S. and S.R.; validation, V.S. and S.R.; formal analysis, V.S. and S.R.; investigation, V.S. and S.R.; resources, V.S. and S.R.; data curation, V.S. and S.R.; writing—original draft preparation, V.S. and S.R.; writing—review and editing, V.S. and S.R.; visualization, V.S. and S.R.; supervision, V.S. and S.R.; project administration, V.S. and S.R.; funding acquisition, V.S. and S.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CLTCentral Limit Theorem
i.i.d.independent and identically distributed

References

  1. Senatov, V.V. The Central Limit Theorem: Approximation Accuracy and Asymptotic Expansions; Librokom: Moscow, Russia, 2009. (In Russian) [Google Scholar]
  2. Shevtsova, I.G. Optimizaciya struktury momentnyh ocenok tochnosti normal’noj approksimacii dlya raspredelenij summ nezavisimyh sluchajnyh velichin [Optimization of the Structure of Moment Estimates for the Accuracy of Normal Approximation of the Distribution of Sums of Independent Random Variables]. Ph.D. Dissertation, Moscow State Univiversity, Moscow, Russia, 2013; 354p. (In Russian). [Google Scholar]
  3. Sobolev, V.N. On the estimation of the accuracy of asymptotic expansions in the central limit theorem. Surv. Appl. Ind. Math. 2011, 18, 45–51. (In Russian) [Google Scholar]
  4. Prohorov, Y.V. Predel’nye teoremy dlya summ nezavisimyh sluchajnyh velichin [Limit Theorems for Sums of Independent Random Variables]. Ph.D. Thesis, Moscow State Univiversity, Moscow, Russia, 1952. (In Russian). [Google Scholar]
  5. Cramér, H. On an asymptotic expansion occurring in the theory of probability. J. Lond. Math. Soc. 1927, 2, 262–265. [Google Scholar] [CrossRef]
  6. Esseen, C.-G. Fourier analysis of distribution functions: A mathematical study of the Laplace–Gaussian law. Acta Math. 1945, 106, 1–125. [Google Scholar] [CrossRef]
  7. Cramér, H. Mathematical Methods of Statistics; Princeton University Press: Princeton, NJ, USA, 1946. [Google Scholar]
  8. Edgeworth, F.Y. The law of error. Camb. Phil. Soc. Proc. 1905, 20, 36–141. [Google Scholar] [CrossRef]
  9. Petrov, V.V. On some polynomials encountered in probability. Vestnik Leningrad Univ. Math. 1962, 19, 150–153. [Google Scholar]
  10. Sobolev, V.N.; Kondratenko, A.E. On Senatov Moments in Asymptotic Expansions in the Central Limit Theorem. Theory Probab. Appl. 2022, 67, 154–157. [Google Scholar] [CrossRef]
  11. Senatov, V.V. On Asymptotic Expansions in the Central Limit Theorem with Explicit Estimates of Remainder Terms. Theory Probab. Appl. 2006, 51, 810–816. [Google Scholar] [CrossRef]
  12. Senatov, V.V.; Sobolev, V.N. New forms of asymptotic expansions in the central limit theorem. Theory Probab. Appl. 2013, 57, 82–96. [Google Scholar] [CrossRef]
  13. Sobolev, V.N. On asymptotic expansions in CLT. Vestn. TVGU. Seriya Prikl. Mat. [Herald Tver State Univ. Ser. Appl. Math.] 2010, 18, 35–48. (In Russian) [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sobolev, V.; Ramodanov, S. On Estimation of the Remainder Term in New Asymptotic Expansions in the Central Limit Theorem. Comput. Sci. Math. Forum 2023, 7, 60. https://doi.org/10.3390/IOCMA2023-14423

AMA Style

Sobolev V, Ramodanov S. On Estimation of the Remainder Term in New Asymptotic Expansions in the Central Limit Theorem. Computer Sciences & Mathematics Forum. 2023; 7(1):60. https://doi.org/10.3390/IOCMA2023-14423

Chicago/Turabian Style

Sobolev, Vitaly, and Sergey Ramodanov. 2023. "On Estimation of the Remainder Term in New Asymptotic Expansions in the Central Limit Theorem" Computer Sciences & Mathematics Forum 7, no. 1: 60. https://doi.org/10.3390/IOCMA2023-14423

Article Metrics

Back to TopTop