Next Article in Journal
Market Microstructure Effects on Firm Default Risk Evaluation
Previous Article in Journal
Evaluating Eigenvector Spatial Filter Corrections for Omitted Georeferenced Variables
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of Gini Index within Pre-Specified Error Bound

by
Bhargab Chattopadhyay
1,* and
Shyamal Krishna De
2
1
Department of Mathematical Sciences, The University of Texas at Dallas, Richardson, TX 75080, USA
2
School of Mathematical Sciences, National Institute of Science Education and Research, Jatni 752050, Odisha, India
*
Author to whom correspondence should be addressed.
Econometrics 2016, 4(3), 30; https://doi.org/10.3390/econometrics4030030
Submission received: 19 December 2015 / Revised: 25 April 2016 / Accepted: 3 June 2016 / Published: 24 June 2016

Abstract

:
Gini index is a widely used measure of economic inequality. This article develops a theory and methodology for constructing a confidence interval for Gini index with a specified confidence coefficient and a specified width without assuming any specific distribution of the data. Fixed sample size methods cannot simultaneously achieve both specified confidence coefficient and fixed width. We develop a purely sequential procedure for interval estimation of Gini index with a specified confidence coefficient and a specified margin of error. Optimality properties of the proposed method, namely first order asymptotic efficiency and asymptotic consistency properties are proved under mild moment assumptions of the distribution of the data.

1. Introduction

Economic inequality arises due to inequality in the distribution of income and assets among individuals or groups within a society or region or even between countries. Economic inequality is usually measured to evaluate the effects of economic policies at the micro or macro level. Several inequality indexes that measure the economic inequality are proposed in the economics literature. Among those indexes, Gini inequality index is the most widely used measure. The most celebrated Gini index, as given in [1], is
G F ( X ) = Δ 2 μ , where Δ = E X 1 X 2 , μ = E ( X )
and X 1 & X 2 are two i.i.d. copies of nonnegative random variable X with distribution function F. Gini index compares every individual’s income with other individual’s income. If there are n randomly selected individuals with incomes given by X 1 , , X n , then the estimator of the celebrated Gini index is
G ^ n = Δ ^ n 2 X ¯ n ,
where X ¯ n is the sample mean and Δ ^ n is the sample Gini’s mean difference defined as,
X ¯ n = 1 n i = 1 n X i and Δ ^ n = n 2 1 1 i 1 < i 2 n X i 1 X i 2 .
The Gini index is undefined if X ¯ n = 0 . We ignore this special case.
Inference for inequality measures, including Gini index, has been an area of research interest among many economists in recent years. For the existing literature on inference problems related to inequality index, we refer to [2,3,4,5,6,7]. Even though there exist innovative methods for constructing confidence intervals for G F (e.g., see [7]), due to large standard errors of estimated Gini index as mentioned in [8], we may not get a short confidence interval for Gini index. We know that the confidence interval varies from sample to sample and so does its width. Wider confidence intervals provide less precise information about the true value of the parameter of interest. Since it is desirable to construct shorter confidence intervals, we rather fix the length of the confidence interval, or in other words, the margin of error while achieving the confidence coefficient ( 1 α ) for some specified α in ( 0 , 1 ) . This problem is known as the fixed-width confidence interval estimation problem.
No fixed sample size procedure can provide a solution to the fixed-width confidence interval estimation problem (e.g., see [9]). Thus, one must resort to sampling in stages to construct a 100 ( 1 α ) % confidence interval for G F with a pre-specified width. This problem falls in the domain of sequential analysis. For details about the general theory of fixed-width confidence interval estimation, we refer to [10,11]. Sequential analysis is concerned with studies where sample sizes are not fixed in advance. Instead, the sequential estimation procedure depends on collecting observations until an a-priori specified criterion or stopping rule is satisfied.
We know that Gini’s mean difference is U-statistic with a symmetric kernel of degree 2 and the sample mean is a U-statistic with a symmetric kernel of degree 1 (e.g., see [12]). Under distribution-free scenario, [7] used the central limit theorem for U-statistics to come up with a confidence interval for Gini index. However, this cannot be used to find out a fixed-width confidence interval for Gini index. In this article, we solve the problem of obtaining a fixed-width confidence interval for Gini index using a purely sequential procedure with a stopping rule based on several U-statistics. Apart from being unbiased estimators, U-statistics are also reverse martingales with respect to some non-increasing filtration as proven in [13]. For more literature on reverse martingales, we refer to classical textbooks on probability theory and stochastic processes such as [14,15]. We exploit the reverse martingale property of U-statistics to derive attractive asymptotic properties of our proposed estimation procedure.
In the next section, we formally state the fixed-width confidence interval estimation problem and why a fixed-sample size procedure cannot be used. In Section 3, a purely sequential procedure is proposed to construct a 100 ( 1 α ) % fixed-width confidence interval for unknown population Gini index and implementation and characteristics of the sequential procedure is discussed as well. Section 4 presents simulation study and validate all theoretical results related to our procedure. We conclude this article with some remarks in Section 5.

2. Problem Statement and Optimal Sample Size

Consider n randomly selected individuals from some population of interest with incomes denoted by X 1 , X 2 , , X n . Suppose these are nonnegative independent and identically distributed random variables assumed to be drawn from an unknown distribution function F where the support of the distribution is ( 0 , ) . A strongly consistent estimator of population Gini index G F is G ^ n given in (2). For pre-specified α in ( 0 , 1 ) , the goal of this paper is to develop the theory for constructing a 100 ( 1 α ) % fixed-width confidence interval for G F . Formally, we would like to construct a confidence interval J n = ( G ^ n d , G ^ n + d ) such that
P G ^ n d < G F < G ^ n + d 1 α ,
for some prefixed margin of error d > 0 . Using [7], we have
n G ^ n G F N 0 , ξ 2 as n ,
where ξ 2 is the asymptotic variance given by
ξ 2 = Δ 2 4 μ 4 σ 2 Δ τ μ 3 + Δ 2 μ 2 + σ 1 2 μ 2 .
Here,
τ = E ( X 1 X 1 X 2 ) and σ 1 2 = V E X 1 X 2 | X 1 .
Schroder and Yitzhaki [16] proposed a way to come up with the reasonable sample size related to the convergence of the distribution of G ^ n to normality. In this paper, we verify via simulation study that for moderate sample sizes (see Section 4) the distribution of sample Gini index is approximately normal. Based on the asymptotic normality of G ^ n , we observe that the coverage probability is
P G ^ n d < G F < G ^ n + d 2 Φ d n ξ 1 ,
where Φ is the distribution function of standard normal random variable. In order to have 100 ( 1 α ) % confidence interval, sample size n must satisfy
2 Φ d n ξ 1 1 α .
Solving (7) for n, we obtain n d 2 z α / 2 2 ξ 2 , where z α / 2 is the upper α 2 t h quantile of the standard normal distribution. Provided ξ is known, the optimal (minimal) sample size required to construct a fixed-width confidence interval for Gini index with approximately ( 1 α ) coverage probability is
C = d 2 z α / 2 2 ξ 2 ,
where, w the lowest integer which is greater than or equal to w.
The optimal fixed sample size C is unknown since the true value of ξ is unknown in practice. If C were known, one would just draw C observations independently from the population of interest and compute ( G ^ C d , G ^ C + d ) which would satisfy (4) approximately. Since C is unknown, one must draw samples at least in two stages in order to achieve the desired coverage probability at least approximately. In the first stage, one must estimate C by estimating ξ, and then in the subsequent stages one should collect samples until the current sample size is more or equal to the estimated optimal sample size. In this article, we propose a purely sequential sampling procedure to estimate the optimal sample size C and ensure that the fixed-width confidence interval based on the final sample size attains the desired ( 1 α ) coverage probability atleast asymptotically.

3. The Sequential Estimation Procedure

In sequential estimation procedures, the parameter estimates are updated as the data is observed. In the first step, a small sample, called the pilot sample, is observed to gather preliminary information about the parameter of interest. Then, in each successive step, one or more additional observations are collected and the estimates of the parameters are updated. After each and every step a decision is taken whether to continue or to terminate the sampling process. This decision is based on a pre-defined stopping rule.
From (8) we note that the optimal sample size needed to find a fixed-width confidence interval depends on unknown parameter ξ 2 . So, let us first find a good estimator of the unknown parameter ξ 2 . Following [7,17], we consider the following strongly consistent estimator of ξ 2 based on U-statistics. Let us define a U-statistic, for each j = 1 , 2 , , n ,
Δ ^ n ( j ) = n 1 2 1 T j X i 1 X i 2 ,
where T j = { ( i 1 , i 2 ) : 1 i 1 < i 2 n and i 1 , i 2 j } . Define W j n = n Δ ^ n ( n 2 ) Δ ^ n ( j ) for j = 1 , , n , and W ¯ n = n 1 j = 1 n W j n . According to [17], a strongly consistent estimator of 4 σ 1 2 is
s w n 2 = ( n 1 ) 1 i = 1 n ( W j n W ¯ n ) 2 .
Using [7],
τ ^ n = 2 n ( n 1 ) ( n , 2 ) 1 2 ( X i 1 + X i 2 ) X i 1 X i 2
is an estimator of τ. Let S n 2 be the sample variance. Thus, the estimator of ξ 2 is
V n 2 = m a x 0 , Δ ^ n 2 S n 2 4 X ¯ n 4 Δ ^ n X ¯ n 3 τ ^ n + Δ ^ n 2 X ¯ n 2 + s w n 2 4 X ¯ n 2
similar to [7]. Note that τ ^ n , S n 2 , Δ ^ n are U-statistics of degree 2 (e.g., [7,18]) and the sample mean X ¯ n is a U-statistic of degree 1. Using continuous mapping theorem [17], and Theorem 3.2.1 in [11], we observe that a strongly consistent estimator of ξ 2 is V n 2 which will be used in our proposed sequential procedure to estimate C.
Several plug-in estimators of the asymptotic variance parameter of Gini index are proposed in the economics and statistics literature. To find the details about several plug-in estimators of the asymptotic variance of Gini index under different sampling schemes, we refer to [4,19]. The proposed plug-in estimator in [4] is simpler than V n 2 . However, it is not known whether the estimator enjoys the almost sure convergence property which is very important for us as we need this property to prove the asymptotic optimality properties of the proposed sequential procedure. Moreover, with the high-end computing facilities available these days, V n 2 can be computed in seconds.
Using V n 2 as the estimator of ξ 2 , we define the stopping rule N d , for every d > 0 , as
N d = the smallest integer n ( m ) such that n z α / 2 d 2 V n 2 + n 1 .
Here, m is called the initial or pilot sample size, and the term n 1 is known as a correction term. Note that V n can be very close to zero with positive probability. Without the correction term, the inequality (12) may be satisfied for very small n terminating the sampling process too early. Thus the correction term n 1 ensures that the sampling process for estimating the optimal sample size does not stop too early. For details about the correction term, we refer to [11].
From (12), we note that, N d z α / 2 d 2 N d 1 , i.e., the final sample size must be at least z α / 2 / d . Therefore, we consider the pilot sample size to be m = max { 4 , z α / 2 / d } . This technique of estimating pilot sample size can also be found in [10].
Recall that the optimal sample size required to achieve 100 ( 1 α ) % confidence interval for Gini index is C which is unknown in practice. The stopping variable N d defined in (12) serves as an estimator of C. Below, we develop a purely sequential procedure to estimate the optimal sample size C.

Implementation and Characteristics

We propose the following purely sequential procedure to estimate the optimal sample size C.
Stage 1: Compute the pilot sample size m = max { 4 , z α / 2 / d } and draw a random sample of size m from the population of interest. Based on this pilot sample of size m, obtain an estimate of ξ 2 by finding V m 2 as given in (11) and check whether m ( z α / 2 / d ) 2 V m 2 + m 1 . If m < ( z α / 2 / d ) 2 V m 2 + m 1 then go to the next step. Otherwise, set the final sample size N d = m .
Stage 2: Draw an additional observation independent of the pilot sample and update the estimate of ξ 2 by computing V m + 1 2 . Check if m + 1 ( z α / 2 / d ) 2 V m + 1 2 + ( m + 1 ) 1 . If m + 1 < ( z α / 2 / d ) 2 V m + 1 2 + ( m + 1 ) 1 then go to the next step. Otherwise, stop sampling and report the final sample size as N d = m + 1 .
The process of collecting observations one by one is continued until there are N d observations such that N d ( z α / 2 / d ) 2 V N d 2 + N d 1 . At this stage, we stop sampling and report the final sample size as N d .
Based on the above algorithm, the sampling process will stop at some stage. This is proved in Lemma A1 which states that if observations are collected using (12), under appropriate conditions, P ( N d < ) = 1 . This is a very important property of any sequential procedure since it mathematically ensures that the sampling will be terminated eventually.
Next, we establish some desirable asymptotic properties of our proposed sequential procedure. First, we prove that the final sample size N d required by our sampling strategy is close to the optimal sample size C at least asymptotically. We also prove the asymptotic efficiency property of sequential procedure which ensures that, on average, we collect only the minimum number of samples to achieve certain accuracy of estimation. Second, we show that the fixed-width confidence interval G ^ N d d , G ^ N d + d contains the true value of Gini index G F nearly with probability 1 α . We formally state these results in Theorems 1 and 2.
Theorem 1. 
If the parent distribution F of X is such that E [ X 4 ] exists, then the stopping rule in (12) yields the following asymptotic optimality properties:
(i) 
N d / C a . s . 1 as d 0 .
(ii) 
P G ^ N d d < G F < G ^ N d + d 1 α a s d 0 .
Theorem 2. 
If the parent distribution F of X is such that the support of the distribution being ( t , ) with t > 0 and E [ X 4 ] exists, then the stopping rule in (12) yields
E N d / C 1 a s d 0 .
Theorems 1 and 2 are proved in the appendix. Part (i) of Theorem 1 implies that the ratio of final sample size of our procedure and the optimal sample size, C asymptotically converges to 1. Part (ii) of Theorem 1 implies that the coverage probability produced by the fixed-width confidence interval G ^ N d d , G ^ N d + d attains the desired level 1 α asymptotically. Theorem 2 implies that the ratio of the average final sample size of our procedure asymptotically converges to the optimal sample size, C.

4. Simulation Study

In this section, we validate the asymptotic properties of our method stated in Theorems 1 and 2 through Monte Carlo study. To implement the sequential procedure, we fix d ( = 0 . 01 , 0 . 02 ) and α ( = 0 . 1 , 0 . 05 ) . Using the pilot sample size formula m = max { 4 , z α 2 / d } , the pilot sample size considered here is 165. Then, we implement the sequential procedure described in Subsection 3.1 and estimate the average sample size ( N ¯ ), the maximum sample size ( max ( N ) ), the standard error ( s ( N ¯ ) ) of N ¯ , the coverage probability (p), and its standard error ( s p ) based on 2000 replications by drawing random samples from gamma (shape = 2 . 649 , rate = 0 . 84 ) distribution truncated at t = 0 . 001 ( G F = 0 . 3308 , ξ 2 = 0 . 0468 ) , log-normal distribution (mean = 2 . 185 , sd = 0 . 562 ) truncated at t = 0 . 001 ( G F = 0 . 3089 , ξ 2 = 0 . 0532 ) , and Pareto (20,000, 5).
Table 1 summarizes the numerical results obtained from the simulation study. The parameters of log-normal and gamma distributions are the same as used by [20]. From the fourth column of Table 1, we find that the ratio of the average final sample size and C is close to 1. Moreover, column 6 of Table 1 illustrates that the attained coverage probability is very close to the desired level of 90 % . Thus, we find that the simulation results validate all theoretical results mentioned in the previous section, and the performance of the procedure is satisfactory for the above mentioned distributions.
To verify whether the distribution of sample Gini index converges to normality, we do the following simulation study. We draw samples of sizes n = 470 , 509, and 245 from gamma (shape = 2 . 649 , rate = 0 . 84 ), log-normal (mean = 2 . 185 , sd = 0 . 562 ), and Pareto (20,000, 5) distributions respectively and compute the sample Gini indexes. Please note that the choice of the parameters of lognormal and gamma distributions are same as in [20]. The sample size n is chosen according to the smallest N ¯ in Table 1 for different scenarios. We compute 200 replications for the sample Gini index (using set.seed(123) in “R” language) and observe that Shapiro-Wilk’s test for normality returns p-values 0.1938, 0.5066, and 0.3984 for gamma (shape = 2 . 649 , rate = 0 . 84 ), log-normal (mean = 2 . 185 , sd = 0 . 562 ), and Pareto (20,000, 5) distributions respectively. Thus, we observe that for the above scenarios sample Gini index converges to normality for moderate sample sizes.

5. Discussion and Concluding Remarks

Gini index is a widely used measure of economic inequality index. In order to evaluate the economic policies adopted by a government, it is important to estimate Gini index at any specific time period. If the income data for all households in the region of interest is not available, one needs to estimate Gini index by drawing a sample of households from that region. Typically, large income surveys are associated with different sampling schemes. To review several sampling techniques, we refer to [19,21,22,23,24]. The sampling technique chosen to collect data usually depends on the socio-economic diversity and size of the country or region. For regions or smaller countries with lesser socio-economic diversity, the simple random sampling technique can be used to collect income or expenditure data to estimate Gini index. Several research articles (e.g., [4,5,7,25,26]) are devoted to drawing statistical inference on inequality indexes which are computed from household income or expenditure by means of simple random sampling from the population of interest. In this paper, we also use simple random sampling technique to collect income or expenditure data in order to estimate Gini index accurately. Even though the sequential methodology is introduced under i.i.d. framework, which may be considered as a practical limitation of our work, sequential methodologies may be adopted to different sampling schemes (e.g., see [27]). In the following Subsection 5.1, we discuss the possibility of extending our work to stratified sampling.
Without assuming any specific distribution of the data, we show that the ratio of the final sample size and the optimal sample size approaches 1. We also show that the confidence interval constructed using our proposed sequential method attains the required coverage probability. Thus, based on these results, we conclude that the proposed sequential estimation strategy can efficiently construct a 100 ( 1 α ) % fixed-width confidence interval for Gini index. In this article, we consider that after pilot sample, one additional observation is collected in each step. If instead, a group of r ( 1 ) observations are collected in each step after the pilot sample stage, the same properties will hold. The proofs will be similar to the ones in Appendix.
The theory of the sequential procedure revolves around the idea of “learn-as-you-go”. In our proposed method of estimation, the final sample size is not fixed in advance, and the observations are collected until the estimated optimal sample size is obtained achieving the required 100 ( 1 α ) % fixed-width confidence interval for Gini index. We hope that number of sequential procedures to estimate income inequalities will be developed following this article, whereupon the idea of sequential procedure can be applied to other sampling schemes used in economic surveys, taking into account of the cost considerations as well.

Possible Extension to Stratified Sampling

The sequential procedure that is proposed under i.i.d. framework may be extended to non i.i.d framework where stratified sampling is used. Suppose we divide the population into S strata. In the population, stratum s contains a mass of H s households. The total number of households in the population is H = s = 1 S H s . The density of the household income or expenditure, X in the s th stratum is denoted by d F ( x | s ) .
Now, a sample of n s households (indexed by h s ) is drawn by simple random sample with replacement from every strata so that the total number of households in the sample is n = s = 1 S n s and n s = n a s with s = 1 S a s = 1 , where a s = H s / H . Let x s h s be the total income of the h s th household belonging to the s th stratum. If w s h s = H s n s is the weight of h s th household in the s th stratum, then following [21,22], Gini index can be estimated by the estimator
G ^ = 1 2 μ ^ s = 1 S h s = 1 n s w s h s x s h s 1 F ^ ( x s h s ) ,
where
μ ^ = s = 1 S h s = 1 n s w s h s x s h s and F ^ ( x s h s ) = i = 1 S j = 1 n s w i j I [ x i j x h c s ] / i = 1 S j = 1 n s w i j .
Now, following [21,22] we have
n G ^ G F d N 0 , V * ,
where V * is the asymptotic variance given in [22], modified to take into account of the stratified sampling only. Now,
P G ^ n G F d 2 Φ n V * d 1
The confidence coefficient will be approximately 1 α provided n d / V * z α / 2 . In order to have a fixed-width confidence interval, we need sample size n satisfying
n d 2 z α / 2 2 V * C , say .
Since C is unknown it must be estimated in the first stage and continue sampling until the sample size n is bigger than corresponding estimated value of C. Note that C is the optimum (i.e., minimum) household size to be sampled to achieve ( 1 α ) confidence level provided V * were known. The optimal number of households to be sampled in the s th stratum ( s = 1 , , S ) will be C s = C a s which is also unknown since C is unknown. Bhattacharya [22] proposed an estimator of the asymptotic variance V * of the Gini index under complex household survey which can be used in Equation (18). Then the stopping rule developed in Equation (12) may be modified taking into account of the stratification and finite sampling scenario to find out an estimate of the optimum number of households in order to find a fixed-width confidence interval for Gini index under stratified sampling. However, we do not intend to explore this possibility in this article, and we believe that this could be a good topic of future research.

Acknowledgments

We thank the three anonymous referees and the editor whose insightful comments helped us improve the paper. We remain deeply indebted to Professor Gautam Tripathi and Professor Nitis Mukhopadhyay for their comments and suggestions.

Author Contributions

The authors contributed equally.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Lemma A1. 
Under the assumption that ξ < , for any d > 0 , the stopping time N d is finite, that is, P ( N d < ) = 1 .
Proof. 
The Lemma A1 is proved by using (12) and the fact that V n 2 is strongly consistent estimator of ξ 2 and N d as d 0 almost surely.  ☐
Lemma A2. 
The value of sample Gini index lies between 0 and 1.
Proof. 
Let Y 1 , , Y n be the ordered incomes of n persons where Y 1 represents the income of the poorest person and Y n represents the income of the richest person. Using [28], Gini index can be rewritten as
0 G ^ n = 2 i = 1 n i Y i n i = 1 n Y i n + 1 n 2 n i = 1 n Y i n i = 1 n Y i n + 1 n = n 1 n 1 .
This proves the lemma.  ☐

Appendix A.1 Proof of Theorem 1

(i)
The definition of stopping rule N d in (12) yields
z α / 2 d 2 V N d 2 N d m I ( N d = m ) + z α / 2 d 2 V N d 1 2 + ( N d 1 ) 1 .
Since N d a.s. as d 0 and V n ξ a.s. as n , by Theorem 2.1 of [29], V N d 2 ξ 2 a.s. Hence, dividing all sides of (A1) by C and letting d 0 , we prove N d / C 1 a.s. as d 0 .
(ii)
In order to show that our procedure satisfies the asymptotic consistency property, we will derive an Anscombe-type random central limit theorem for Gini index. This requires the existence of usual central limit theorem of Gini index and uniform continuity in probability (u.c.i.p.) condition. For details about the u.c.i.p. condition, we refer to [17,30,31,32] etc.
First of all, let us define n 1 = ( 1 ρ ) C and n 2 = ( 1 + ρ ) C for 0 < ρ < 1 . Now, we know from [7] that Y n = n ( Δ ^ n Δ ) , n ( X ¯ n μ ) L N 2 ( 0 , Σ ) , where
Σ = 4 σ 1 2 2 ( τ μ Δ ) 2 ( τ μ Δ ) σ 2 .
First, let us prove that Y N d L N 2 ( 0 , Σ ) . Define D = ( a 0 a 1 ) . Note that D Y N d = D Y C + ( D Y N d D Y C ) . Thus, it is enough to show that ( D Y N d D Y C ) P 0 as d 0 . We can write
( D Y N d D Y C ) = a 0 N d ( Δ ^ N Δ ^ C ) + a 1 N d ( X ¯ N d X ¯ C ) + ( N d / C 1 ) D Y C .
Fix some ϵ > 0 and note that
P | a 0 N d ( Δ ^ N d Δ ^ C ) + a 1 N d ( X ¯ N d X ¯ C ) | > ϵ P | a 0 N d ( Δ ^ N d Δ ^ C ) + a 1 N d ( X ¯ N d X ¯ C ) | > ϵ , | N d C | < ρ C + P [ | N d C | > ρ C ] P m a x n 1 < n < n 2 | n | Δ ^ n Δ ^ C | > ϵ 2 | a 0 | + P m a x n 1 < n < n 2 | n | X ¯ n X ¯ C | > ϵ 2 | a 1 | + P [ | N d C | > ρ C ]
Here, Δ ^ n and X ¯ n are both U-statistics which satisfy Anscombe’s u.c.i.p. condition (for e.g., see [17]). Using u.c.i.p. condition and the fact that N d / C a . s . 1 , we conclude that for given ϵ > 0 , there exist η > 0 and d 0 > 0 such that
P { | a 0 N d ( Δ ^ N d Δ ^ C ) + a 1 N d ( X ¯ N d X ¯ C ) | > ϵ } < η for all d d 0 .
This implies a 0 N d ( Δ ^ N d Δ ^ C ) + a 1 N d ( X ¯ N d X ¯ C ) P 0 as d 0 . Also, note that N d / C 1 D Y C P 0 as d 0 since N d / C 1 almost surely and D Y C L N 2 ( 0 , Σ ) . Thus, from (A2), we conclude ( D Y N d D Y C ) P 0 , that is, Y N d L N 2 ( 0 , Σ ) . Now, define G ( u , v ) = u 2 v , if v 0 . Using Taylor’s expansion, we can write
N d ( G ( Δ ^ N d , X ¯ N d ) G ( Δ , μ ) ) = N d Δ ^ N d Δ 2 μ Δ 2 μ 2 ( X ¯ N d μ ) + R N d ,
where R N d = 2 ( Δ ^ N d Δ ) ( X ¯ N d μ ) / b 2 + 4 a ( X ¯ N d μ ) 2 / b 3 , a = Δ + p ( Δ ^ N d Δ ) , b = 2 μ + p ( 2 X ¯ N d 2 μ ) , and p ( 0 , 1 ) . Rewriting (A3) in the vector-matrix form, we get
N d ( G ( Δ ^ N d , X ¯ N d ) G ( Δ , μ ) ) = D Y N d + N d R N d ,
where D = 1 2 μ , Δ 2 μ 2 . Note that N d ( X ¯ N d μ ) converges in distribution to a normal distribution by Anscombe’s CLT and both Δ ^ N d Δ and X ¯ N d μ converges to 0 almost surely. This yields N d R N d P 0 as d 0 . Hence, N d ( G ^ N d G F ) L N ( 0 , D Σ D ) as d 0 . This completes the proof of Theorem 1.   ☐

Appendix A.2 Proof of Theorem 2

In this subsection, we prove a lemma that is essential to establish Theorem 2. Note from (12) that N d z α / 2 2 d 2 N d 1 , i.e., N d z α / 2 d ( = m ) with probability 1. Suppose X ( n ) = ( X ( 1 ) , , X ( n ) ) denotes the n dimensional vector of order statistics from the sample X 1 , , X n , and F n is the σ-algebra generated by ( X ( n ) , X n + 1 , X n + 2 , ) . By [13], X ¯ n , F n , S n 2 , F n , τ ^ n , F n , Δ ^ n , F n , and their convex functions are all reverse submartingales. Using reverse submartingale properties, let us prove the following lemma along the lines of [33].
Lemma A3. 
If E ( X 1 2 p ) is finite for some p > 1 , then E sup n m V n 2 < for m 4 .
Proof. 
To prove Lemma A3, it is enough to show that: E sup n m s w n 2 X ¯ n 2 , E sup n m Δ ^ n X ¯ n 3 τ ^ n , E sup n m Δ ^ n 2 X ¯ n 2 , and E sup n m Δ ^ n 2 X ¯ n 4 S n 2 are finite.
We note that, 0 Δ ^ n 2 X ¯ 1 . So, it is enough to show that E sup n m s w n 2 X ¯ n 2 , E sup n m τ ^ n X ¯ n 2 and E sup n m S n 2 X ¯ n 2 are finite. Following [34] (p. 338), we have E sup n m s w n 2 < if E [ X 1 α ] < for α > 2 and m 4 . Therefore,
E sup n m s w n 2 X ¯ n 2 t 2 E sup n m s w n 2 < .
We note that τ ^ n and S n 2 are U-statistics. Using Lemma 9.2.4 of [35], for p > 1 ,
E sup n m τ ^ n p p 1 p p E τ ^ m p and E sup n m S n 2 p p 1 p p E S m 2 p .
Since E ( X 1 2 p ) if finite for some p > 1 , E τ ^ m p and E S m 2 p are finite which yield
E sup n m τ ^ n X ¯ n 2 t 2 E sup n m τ ^ n < ,
and
E sup n m S n 2 X ¯ n 2 t 2 E sup n m S n 2 < .
This completes the proof of Lemma A3.  ☐
Below, we prove Theorem 2 by using Lemma A3.
Since N d m a.s., dividing (A1) by C yields
N d / C m I ( N d = m ) / C 1 ξ 2 sup d > 0 V N d 1 2 + ( m 1 ) 1 almost surely .
Since E sup d > 0 V N d 1 2 < by Lemma A3 and N d / C 1 a.s. as d 0 , by the dominated convergence theorem, we conclude that lim d 0 E ( N d / C ) = 1 .
This completes the proof of Theorem 2. ☐

References

  1. B.C. Arnold. “Inequality measures for multivariate distributions.” Metron 63 (2005): 317–327. [Google Scholar]
  2. R. Andres, and C. Samuel. Inference on Income Inequality and Tax Progressivity Indices: U-Statistics and Bootstrap Methods. ECINEQ working paper 2005-9; Palma, Spain: ECINEQ, 2005. [Google Scholar]
  3. J.A. Bishop, J.P. Formby, and B. Zheng. “Statistical inference and the sen index of poverty.” Int. Econ. Rev. 38 (1997): 381–387. [Google Scholar] [CrossRef]
  4. R. Davidson. “Reliable inference for the Gini index.” J. Econom. 150 (2009): 30–40. [Google Scholar] [CrossRef]
  5. J.L. Gastwirth. “The estimation of the Lorenz curve and Gini index.” Rev. Econ. Stat. 54 (1972): 306–316. [Google Scholar] [CrossRef]
  6. P. Palmitesta, P. Corrado, and S. Cosimo. “Confidence interval estimation for inequality indices of the Gini family.” Comput. Econ. 16 (2000): 137–147. [Google Scholar] [CrossRef]
  7. K. Xu. “U-statistics and their asymptotic results for some inequality and poverty measures.” Econom. Rev. 26 (2007): 567–577. [Google Scholar] [CrossRef]
  8. E. Maasoumi. “Empirical analysis of inequality and welfare.” In Handbook of Applied Microeconomics. Edited by S. Schmidt and H. Pesaran. Malden, MA, USA: Blackwell Publishers Inc., 1997. [Google Scholar]
  9. G.B. Dantzig. “On the non-existence of tests of “student’s” hypothesis having power functions independent of σ.” Ann. Math. Stat. 11 (1940): 186–192. [Google Scholar] [CrossRef]
  10. N. Mukhopadhyay, and B.M. de Silva. Sequential Methods and Their Applications. Boca Raton, FL, USA: CRC Press, 2009. [Google Scholar]
  11. P.K. Sen. Sequential Nonparametrics: Invariance Principles and Statistical Inference. New York, NY, USA: Wiley, 1981. [Google Scholar]
  12. W. Hoeffding. “A class of statistics with asymptotically normal distribution.” Ann. Math. Stat. 19 (1948): 293–325. [Google Scholar] [CrossRef]
  13. A.J. Lee. U-Statistics: Theory and Practice. New York, NY, USA: CRC Press, 1990. [Google Scholar]
  14. M. Loève. Probability Theory. Princeton, NJ, USA: Van Nostrand, 1963. [Google Scholar]
  15. J.L. Doob. Stochastic Processes. New York, NY, USA: Wiley, 1953. [Google Scholar]
  16. C. Schröder, and S. Yitzhaki. Reasonable Sample Sizes for Convergence to Normality. No. 714, SOEP Papers on Multidisciplinary Panel Data Research. 2014. Available online: http://papers.ssrn.com/sol3/papers.cfm?abstractid=2539096 (accessed on 5 March 2016).
  17. R. Sproule. “A Sequential Fixed-Width Confidence Interval for the Mean of a U-Statistic.” Ph.D. Thesis, University of North Carolina, Chapel Hill, NC, USA, 1969. [Google Scholar]
  18. B. Chattopadhyay, and N. Mukhopadhyay. “Two-stage fixed-width confidence intervals for a normal mean in the presence of suspect outliers.” Seq. Anal. 32 (2013): 134–157. [Google Scholar] [CrossRef]
  19. M. Langel, and Y. Tillè. “Variance estimation of the Gini index: Revisiting a result several times published.” J. R. Stat. Soc. Ser. A Stat. Soc. 176 (2013): 521–540. [Google Scholar] [CrossRef]
  20. M.R. Ransom, and J.S. Cramer. “Income distribution functions with disturbances.” Eur. Econ. Rev. 22 (1983): 363–372. [Google Scholar] [CrossRef]
  21. D. Bhattacharya. “Asymptotic inference from multi-stage samples.” J. Econom. 126 (2005): 145–171. [Google Scholar] [CrossRef]
  22. D. Bhattacharya. “Inference on inequality from household survey data.” J. Econom. 137 (2007): 674–707. [Google Scholar] [CrossRef]
  23. D.A. Binder, and M.S. Kovacevic. “Estimating some measures of income inequality from survey data: An application of the estimating equations approach.” Surv. Methodol. 21 (1995): 137–146. [Google Scholar]
  24. W.G. Cochran. Sampling Techniques. New York, NY, USA: Wiley & Sons, 1977. [Google Scholar]
  25. C.M. Beach, and R. Davidson. “Distribution-free statistical inference with Lorenz curves and income shares.” Rev. Econ. Stud. 50 (1983): 723–735. [Google Scholar] [CrossRef]
  26. R. Davidson, and J. Duclos. “Statistical inference for stochastic dominance and for the measurement of poverty and inequality.” Econometrica 68 (2000): 1435–1464. [Google Scholar] [CrossRef]
  27. S. Zacks. Stage-Wise Adaptive Designs. New York, NY, USA: Wiley, 2009. [Google Scholar]
  28. C. Damgaard, and J. Weiner. “Describing inequality in plant size or fecundity.” Ecology 81 (2000): 1139–1142. [Google Scholar] [CrossRef]
  29. A. Gut. Stopped random walks: Limit theorems and applications. New York, NY, USA: Springer, 2009. [Google Scholar]
  30. F.J. Anscombe. “Sequential estimation.” J. R. Stat. Soc. Ser. B 15 (1953): 1–29. [Google Scholar]
  31. E. Isogai. “Asymptotic consistency of fixed-width sequential confidence intervals for a multiple regression function.” Ann. Inst. Stat. Math. 38 (1986): 69–83. [Google Scholar] [CrossRef]
  32. N. Mukhopadhyay, and B. Chattopadhyay. “A tribute to Frank Anscombe and random central limit theorem from 1952.” Seq. Anal. 31 (2012): 265–277. [Google Scholar]
  33. S.K. De, and B. Chattopadhyay. “Minimum Risk Point Estimation of Gini Index.” Available online: http://arxiv.org/abs/1503.08148 (accessed on 27 March 2015).
  34. P.K. Sen, and M. Ghosh. “Sequential point estimation of estimable parameters based on U-statistics.” Sankhyā Indian J. Stat. Ser. A 43 (1981): 331–344. [Google Scholar]
  35. M. Ghosh, N. Mukhopadhyay, and P.K. Sen. Sequential Estimation. New York, NY, USA: Wiley, 1997. [Google Scholar]
Table 1. Performance of the proposed sequential procedure when the data is from gamma, log-normal, and Pareto distribution.
Table 1. Performance of the proposed sequential procedure when the data is from gamma, log-normal, and Pareto distribution.
dDistribution N ¯ C N ¯ / C p
α s ( N ¯ ) s p
d = 0 . 01 Gamma1283.745012671.01320.9090
α = 0 . 1 1.7561 0.0064
d = 0 . 02 Gamma469.16504501.04260.9535
α = 0 . 05 1.0485 0.0047
d = 0 . 01 Lognormal1435.364014400.99680.8965
α = 0 . 1 4.091604 0.0068
d = 0 . 02 Lognormal509.00205110.99610.9430
α = 0 . 05 2.0538 0.0052
d = 0 . 01 Pareto654.53646860.95410.9018
α = 0 . 1 4.2151 0.0063
d = 0 . 02 Pareto244.33302441.00140.9470
α = 0 . 05 2.0099 0.0050

Share and Cite

MDPI and ACS Style

Chattopadhyay, B.; De, S.K. Estimation of Gini Index within Pre-Specified Error Bound. Econometrics 2016, 4, 30. https://doi.org/10.3390/econometrics4030030

AMA Style

Chattopadhyay B, De SK. Estimation of Gini Index within Pre-Specified Error Bound. Econometrics. 2016; 4(3):30. https://doi.org/10.3390/econometrics4030030

Chicago/Turabian Style

Chattopadhyay, Bhargab, and Shyamal Krishna De. 2016. "Estimation of Gini Index within Pre-Specified Error Bound" Econometrics 4, no. 3: 30. https://doi.org/10.3390/econometrics4030030

APA Style

Chattopadhyay, B., & De, S. K. (2016). Estimation of Gini Index within Pre-Specified Error Bound. Econometrics, 4(3), 30. https://doi.org/10.3390/econometrics4030030

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