Estimating Engel curves: a new way to improve the SILC-HBS matching process using GLM methods

Microdata are required to evaluate the distributive impact of the taxation system as a whole (direct and indirect taxes) on individuals or households. However, in European Union countries this information is usually distributed into two separate surveys: the Household Budget Surveys (HBS), including total household expenditure and its composition, and EU Statistics on Income and Living Conditions (EU-SILC), including detailed information about households'' income and direct (but not indirect) taxes paid. We present a parametric statistical matching procedure to merge both surveys. For the first stage of matching, we propose estimating total household expenditure in HBS (Engel curves) using a GLM estimator, instead of the traditionally used OLS method. It is a better alternative, insofar as it can deal with the heteroskedasticity problem of the OLS estimates, while making it unnecessary to retransform the regressors estimated in logarithms. To evaluate these advantages of the GLM estimator, we conducted a computational Monte Carlo simulation. In addition, when an error term is added to the deterministic imputation of expenditure in the EU-SILC, we propose replacing the usual Normal distribution of the error with a Chi-square type, which allows a better approximation to the original expenditures variance in the HBS. An empirical analysis is provided using Spanish surveys for years 2012–2016. In addition, we extend the empirical analysis to the rest of the European Union countries, using the surveys provided by Eurostat (EU-SILC, 2011; HBS, 2010).     

M-quantile models with application to poverty mapping

Over the last decade there has been growing demand for estimates of population characteristics at small area level. Unfortunately, cost constraints in the design of sample surveys lead to small sample sizes within these areas and as a result direct estimation, using only the survey data, is inappropriate since it yields estimates with unacceptable levels of precision. Small area models are designed to tackle the small sample size problem. The most popular class of models for small area estimation is random effects models that include random area effects to account for between area variations. However, such models also depend on strong distributional assumptions, require a formal specification of the random part of the model and do not easily allow for outlier robust inference. An alternative approach to small area estimation that is based on the use of M-quantile models was recently proposed by Chambers and Tzavidis (Biometrika 93(2):255–268, 2006) and Tzavidis and Chambers (Robust prediction of small area means and distributions. Working paper, 2007). Unlike traditional random effects models, M-quantile models do not depend on strong distributional assumption and automatically provide outlier robust inference. In this paper we illustrate for the first time how M-quantile models can be practically employed for deriving small area estimates of poverty and inequality. The methodology we propose improves the traditional poverty mapping methods in the following ways: (a) it enables the estimation of the distribution function of the study variable within the small area of interest both under an M-quantile and a random effects model, (b) it provides analytical, instead of empirical, estimation of the mean squared error of the M-quantile small area mean estimates and (c) it employs a robust to outliers estimation method. The methodology is applied to data from the 2002 Living Standards Measurement Survey (LSMS) in Albania for estimating (a) district level estimates of the incidence of poverty in Albania, (b) district level inequality measures and (c) the distribution function of household per-capita consumption expenditure in each district. Small area estimates of poverty and inequality show that the poorest Albanian districts are in the mountainous regions (north and north east) with the wealthiest districts, which are also linked with high levels of inequality, in the coastal (south west) and southern part of country. We discuss the practical advantages of our methodology and note the consistency of our results with results from previous studies. We further demonstrate the usefulness of the M-quantile estimation framework through design-based simulations based on two realistic survey data sets containing small area information and show that the M-quantile approach may be preferable when the aim is to estimate the small area distribution function.     

Taylor linearization sampling errors and design effects for poverty measures and other complex statistics

A systematic procedure for the derivation of linearized variables for the estimation of sampling errors of complex nonlinear statistics involved in the analysis of poverty and income inequality is developed. The linearized variable extends the use of standard variance estimation formulae, developed for linear statistics such as sample aggregates, to nonlinear statistics. The context is that of cross-sectional samples of complex design and reasonably large size, as typically used in population-based surveys. Results of application of the procedure to a wide range of poverty and inequality measures are presented. A standardized software for the purpose has been developed and can be provided to interested users on request. Procedures are provided for the estimation of the design effect and its decomposition into the contribution of unequal sample weights and of other design complexities such as clustering and stratification. The consequence of treating a complex statistic as a simple ratio in estimating its sampling error is also quantified. The second theme of the paper is to compare the linearization approach with an alternative approach based on the concept of replication, namely the Jackknife repeated replication (JRR) method. The basis and application of the JRR method is described, the exposition paralleling that of the linearization method but in somewhat less detail. Based on data from an actual national survey, estimates of standard errors and design effects from the two methods are analysed and compared. The numerical results confirm that the two alternative approaches generally give very similar results, though notable differences can exist for certain statistics. Relative advantages and limitations of the approaches are identified.     

Bootstrapping Spatial Median for Location Problems

In multivariate location problems, the sample mean is most widely used, having various advantages. It is, however, very sensitive to outlying observations and inefficient for data from heavy tailed distributions. In this situation, the spatial median is more robust than the sample mean and could be a reasonable alternative. We reviewed several spatial median based testing methods for multivariate location and compared their significance level and power through Monte Carlo simulations. The results show that bootstrap method is efficient for the estimation of the covariance matrix of the sample spatial median. We also proposed bootstrap simultaneous confidence intervals based on the spatial median for multiple comparisons in the multi-sample case.     

A nonparametric plug-in rule for selecting optimal block lengths for block bootstrap methods

In this paper, we consider the problem of empirical choice of optimal block sizes for block bootstrap estimation of population parameters. We suggest a nonparametric plug-in principle that can be used for estimating ‘mean squared error’-optimal smoothing parameters in general curve estimation problems, and establish its validity for estimating optimal block sizes in various block bootstrap estimation problems. A key feature of the proposed plug-in rule is that it can be applied without explicit analytical expressions for the constants that appear in the leading terms of the optimal block lengths. Furthermore, we also discuss the computational efficacy of the method and explore its finite sample properties through a simulation study.     

Statistical inference for the quintile share ratio

In recent years, the Quintile Share Ratio (or QSR) has become a very popular measure of inequality. In 2001, the European Council decided that income inequality in European Union member states should be described using two indicators: the Gini Index and the QSR. The QSR is generally defined as the ratio of the total income earned by the richest 20% of the population relative to that earned by the poorest 20%. Thus, it can be expressed using quantile shares, where a quantile share is the share of total income earned by all of the units up to a given quantile. The aim of this paper is to propose an improved methodology for the estimation and variance estimation of the QSR in a complex sampling design framework. Because the QSR is a non-linear function of interest, the estimation of its sampling variance requires advanced methodology. Moreover, a non-trivial obstacle in the estimation of quantile shares in finite populations is the non-unique definition of a quantile. Thus, two different conceptions of the quantile share are presented in the paper, leading us to two different estimators of the QSR. Regarding variance estimation, [Osier, 2006] and [Osier, 2009] proposed a variance estimator based on linearization techniques. However, his method involves Gaussian kernel smoothing of cumulative distribution functions. Our approach, also based on linearization, shows that no smoothing is needed. The construction of confidence intervals is discussed and a proposition is made to account for the skewness of the sampling distribution of the QSR. Finally, simulation studies are run to assess the relevance of our theoretical results.     

EMPIRICAL LIKELIHOOD METHODS FOR THE GINI INDEX

The Gini index and its generalizations have been used extensively for measuring inequality and poverty in the social sciences. Recently, interval estimation based on nonparametric statistics has been proposed in the literature, for example the naive bootstrap method, the iterated bootstrap method and the bootstrap method via a pivotal statistic. In this paper, we propose empirical likelihood methods to construct confidence intervals for the Gini index or the difference of two Gini indices. Simulation studies show that the proposed empirical likelihood method performs slightly worse than the bootstrap method based on a pivotal statistic in terms of coverage accuracy, but it requires less computation. However, the bootstrap calibration of the empirical likelihood method performs better than the bootstrap method based on a pivotal statistic.     

Inference for Income Distributions Using Grouped Data

We develop a general approach to estimation and inference for income distributions using grouped or aggregate data that are typically available in the form of population shares and class mean incomes, with unknown group bounds. We derive generic moment conditions and an optimal weight matrix that can be used for generalized method-of-moments (GMM) estimation of any parametric income distribution. Our derivation of the weight matrix and its inverse allows us to express the seemingly complex GMM objective function in a relatively simple form that facilitates estimation. We show that our proposed approach, which incorporates information on class means as well as population proportions, is more efficient than maximum likelihood estimation of the multinomial distribution, which uses only population proportions. In contrast to the earlier work of Chotikapanich, Griffiths, and Rao, and Chotikapanich, Griffiths, Rao, and Valencia, which did not specify a formal GMM framework, did not provide methodology for obtaining standard errors, and restricted the analysis to the beta-2 distribution, we provide standard errors for estimated parameters and relevant functions of them, such as inequality and poverty measures, and we provide methodology for all distributions. A test statistic for testing the adequacy of a distribution is proposed. Using eight countries/regions for the year 2005, we show how the methodology can be applied to estimate the parameters of the generalized beta distribution of the second kind (GB2), and its special-case distributions, the beta-2, Singh–Maddala, Dagum, generalized gamma, and lognormal distributions. We test the adequacy of each distribution and compare predicted and actual income shares, where the number of groups used for prediction can differ from the number used in estimation. Estimates and standard errors for inequality and poverty measures are provided. Supplementary materials for this article are available online.     

District-level poverty estimation: a proposed method

This paper develops a method of estimating micro-level poverty in cases where data are scarce. The method is applied to estimate district-level poverty using the household level Indian national sample survey data for two states, viz., West Bengal and Madhya Pradesh. The method involves estimation of state-level poverty indices from the data formed by pooling data of all the districts (each time excluding one district) and multiplying this poverty vector with a known weight matrix to obtain the unknown district-level poverty vector. The proposed method is expected to yield reliable estimates at the district level, because the district-level estimate is now based on a much larger sample size obtained by pooling data of several districts. This method can be an alternative to the “small area estimation technique” for estimating poverty at sub-state levels in developing countries.     

Bootstrapping mean‐squared errors of robust small‐area estimators: Application to the method‐of‐payments surveys data

This paper proposes a new bootstrap procedure for mean‐squared errors of robust small‐area estimators. We formally prove the asymptotic validity of the proposed bootstrap method and examine its finite‐sample performance through Monte Carlo simulations. The results show that our procedure performs well and competes with existing ones. We also provide an application to the estimation of the total volume and value of cash, debit card, and credit card transactions in Canada as well as in its provinces and subgroups of households. In particular, we found that there is a significant average annual decline rate of 3.1% in the volume of cash transactions and that this decline is relatively higher among high‐income households living in heavily populated provinces. Our bootstrap estimator also provides indicators of quality useful in selecting the best small‐area predictor among several alternatives in practice.     

A general class of poverty measures

Measurement of poverty in the spirit of Sen essentially consists of two steps, namely (1) the identification of the poor and (2) the aggregation of the income distribution of the poor into an indicator called poverty measure. We contribute to the second step from a mathematical and statistical point of view by defining a class of poverty measures and investigating its properties including sensitivity with respect to changes in the income distribution. Various poverty orderings are defined and their relation to stochastic dominance is investigated. Finally estimation of the class of poverty measures from individual and grouped data is considered.     

U-Statistics and Their Asymptotic Results for Some Inequality and Poverty Measures

U-statistics form a general class of statistics that have certain important features in common. This class arises as a generalization of the sample mean and the sample variance, and typically members of the class are asymptotically normal with good consistency properties. The class encompasses some widely used income inequality and poverty measures, in particular the variance, the Gini index, the poverty rate, the average poverty gap ratios, the Foster-Greer-Thorbecke index, and the Sen index and its modified form. This paper illustrates how these measures come together within the class of U-statistics, and thereby why U-statistics are useful in econometrics.     

U-Statistics and Their Asymptotic Results for Some Inequality and Poverty Measures

U-statistics form a general class of statistics that have certain important features in common. This class arises as a generalization of the sample mean and the sample variance, and typically members of the class are asymptotically normal with good consistency properties. The class encompasses some widely used income inequality and poverty measures, in particular the variance, the Gini index, the poverty rate, the average poverty gap ratios, the Foster-Greer-Thorbecke index, and the Sen index and its modified form. This paper illustrates how these measures come together within the class of U-statistics, and thereby why U-statistics are useful in econometrics.     

收入增长、分配公平与贫困减少

 利用CHNS调查数据,分两个阶段：1990-1999、1999-2005,依据收入分布函数定量讨论了收入增长、分配公平与贫困减少三者之间的关系。均值回归与O-B分解结果显示：收入决定机制是影响人均收入增长的主导因素,然而需要关注在1999-2005年间教育的系数效应为负。分位数回归与M-M分解显示：两个阶段,收入决定机制是导致居民收入分布向右偏移的主要原因;第一阶段的居民收入差距扩大主要源于劳动力特征的变化,第二阶段的居民收入差距扩大难以解释。Shapley分解结果显示：收入增长能显著减少贫困,而收入差距扩大会累及减贫效果;在两个不同阶段,贫困减少的增长效应贡献始终最大,而分配效应与贫困线变动效应贡献顺序发生了交替变化。     

Comparison of bootstrap estimation intervals to forecast arithmetic mean and median air passenger demand

The aim of this paper is to compare passenger (pax) demand between airports based on the arithmetic mean (MPD) and the median pax demand (MePD). A three phases approach is applied. First phase, we use bootstrap procedures to estimate the distribution of the arithmetic MPD and the MePD for each block of routes distance; second phase, we use percentile, standard, bias corrected, and bias corrected accelerated methods to calculate bootstrap confidence bands for the MPD and the MePD; and third phase, we implement Monte Carlo (MC) experiments to analyse the finite sample performance of the applied bootstrap. Our results conclude that it is more meaningful to use the estimation of MePD rather than the estimation of MPD in the air transport industry. By carrying out MC experiments, we demonstrate that the bootstrap methods produce coverages close to the nominal for the MPD and the MePD.     

Importance resampling for the smoothed bootstrap

Alternative methods of estimating properties of unknown distributions include the bootstrap and the smoothed bootstrap. In the standard bootstrap setting, Johns (1988) introduced an importance resam¬pling procedure that results in more accurate approximation to the bootstrap estimate of a distribution function or a quantile. With a suitable “exponential tilting” similar to that used by Johns, we derived a smoothed version of importance resampling in the framework of the smoothed bootstrap. Smoothed importance resampling procedures were developed for the estimation of distribution functions of the Studentized mean, the Studentized variance, and the correlation coefficient. Implementation of these procedures are presented via simulation results which concentrate on the problem of estimation of distribution functions of the Studentized mean and Studentized variance for different sample sizes and various pre-specified smoothing bandwidths for the normal data; additional simulations were conducted for the estimation of quantiles of the distribution of the Studentized mean under an optimal smoothing bandwidth when the original data were simulated from three different parent populations: lognormal, t(3) and t(10). These results suggest that in cases where it is advantageous to use the smoothed bootstrap rather than the standard bootstrap, the amount of resampling necessary might be substantially reduced by the use of importance resampling methods and the efficiency gains depend on the bandwidth used in the kernel density estimation.     

Kernel estimations for multivariate density functional with bootstrap

In this article the bootstrap method is discussed for the kernel estimation of the multivariate density function. We have considered sample mean functional and constructed its consistency and asymptotic normality by bootstrap estimator. It has been shown that the bootstrap works for kernel estimates of multivariate density functional. The convergence rate with bootstrap for density has been proved. Finally, two simulations of application are given.     

A simulation study of balanced and antithetic bootstrap resampling methods

An extensive simulation study is conducted to compare the performance between balanced and antithetic resampling for the bootstrap in estimation of bias, variance, and percentiles when the statistic of interest is the median, the square root of the absolute value of the mean, or the median absolute deviations from the median. Simulation results reveal that balanced resampling provide better efficiencies in most cases; however, antithetic resampling is superior in estimating bias of the median. We also investigate the possibility of combining an existing efficient bootstrap computation of Efron (1990) with balanced or antithetic resampling for percentile estimation. Results indicate that the combination method does indeed offer gains in performance though the gains are much more dramatic for the bootstrap t statistic than for any of the three statistics of interest as described above.     

中国养老金体系待遇充足性的多维度评估

养老金制度的首要目标是防止老年贫困和实现一定水平的收入替代,从而养老金的“待遇充足”至关重要。本文借鉴国际上关于养老金充足性的多维度评估思路,结合中国实际,构建了中国养老金体系待遇充足性的三维度、六指标评价体系。在此基础上,基于国家统计局公开的数据以及全国综合社会调查（CGSS）提供的微观数据进行测算分析,发现我国养老金体系在收入水平、防止贫困和性别差异三个维度上的表现并不理想,各地区养老金充足性在不同时间上存在差异。进而提出我国应在不同层次养老金制度设计中强化收入替代功能、防止老年贫困和缩小性别差异等政策建议。     

Robust, distribution-free inference for income share ratios under complex sampling

The quintile share ratio of disposable income is the primary inequality indicator of the European Union. As an inequality indicator, it must be sensitive to extreme large observations. Therefore, outliers have a strong impact on the bias and the variance of the classical quintile share ratio estimator. This may mislead the interpretation of income inequality. A class of estimators which are robust against outliers is introduced. They have a bounded influence function, they may reduce the bias incurred by the robustification and they reduce variability. Based on an asymptotic framework which respects the design-based, non-parametric approach, inference for these robust estimators is developed. A large simulation study with close to reality universes derived from the Statistics of Living Conditions Surveys of the EU allows to study the performance of the proposed estimators.