Model‐Assisted Estimation of a Spatial Population Mean

This paper deals with the estimation of the mean of a spatial population. Under a design‐based approach to inference, an estimator assisted by a penalized spline regression model is proposed and studied. Proof that the estimator is design‐consistent and has a normal limiting distribution is provided. A simulation study is carried out to investigate the performance of the new estimator and its variance estimator, in terms of relative bias, efficiency, and confidence interval coverage rate. The results show that gains in efficiency over standard estimators in classical sampling theory may be impressive.     

Cross‐sectional Dependency and Size Distortion in a Small‐sample Homogeneous Panel Data Unit Root Test*

In this paper, we investigate the effects of cross‐sectional disturbance correlation in a homogeneous panel data unit root test. As reported by other authors, the unit root test has incorrect size in the presence of cross‐sectional correlation. We suggest that a previously known estimator can be used to reduce the size distortions. We supply response surface estimates for critical values and study the size characteristics of the proposed test. We find that the suggested estimator performs well in small‐sample homogeneous panel data unit root tests. The reduction in size distortion comes at a small cost of lower power against a stationary alternative.     

On the Generalized Bootstrap for Sample Surveys with Special Attention to Poisson Sampling

We study the generalized bootstrap technique under general sampling designs. We focus mainly on bootstrap variance estimation but we also investigate the empirical properties of bootstrap confidence intervals obtained using the percentile method. Generalized bootstrap consists of randomly generating bootstrap weights so that the first two (or more) design moments of the sampling error are tracked by the corresponding bootstrap moments. Most bootstrap methods in the literature can be viewed as special cases. We discuss issues such as the choice of the distribution used to generate bootstrap weights, the choice of the number of bootstrap replicates, and the potential occurrence of negative bootstrap weights. We first describe the generalized bootstrap for the linear Horvitz‐Thompson estimator and then consider non‐linear estimators such as those defined through estimating equations. We also develop two ways of bootstrapping the generalized regression estimator of a population total. We study in greater depth the case of Poisson sampling, which is often used to select samples in Price Index surveys conducted by national statistical agencies around the world. For Poisson sampling, we consider a pseudo‐population approach and show that the resulting bootstrap weights capture the first three design moments of the sampling error. A simulation study and an example with real survey data are used to illustrate the theory.     

Common correlated effect cross‐sectional dependence corrections for nonlinear conditional mean panel models

This paper provides an approach to estimation and inference for nonlinear conditional mean panel data models, in the presence of cross‐sectional dependence. We modify Pesaran's (Econometrica, 2006, 74(4), 967–1012) common correlated effects correction to filter out the interactive unobserved multifactor structure. The estimation can be carried out using nonlinear least squares, by augmenting the set of explanatory variables with cross‐sectional averages of both linear and nonlinear terms. We propose pooled and mean group estimators, derive their asymptotic distributions, and show the consistency and asymptotic normality of the coefficients of the model. The features of the proposed estimators are investigated through extensive Monte Carlo experiments. We also present two empirical exercises. The first explores the nonlinear relationship between banks' capital ratios and riskiness. The second estimates the nonlinear effect of national savings on national investment in OECD countries depending on countries' openness.     

THE ROLE OF CONDITIONAL HETEROSKEDASTICITY IN IDENTIFYING AND ESTIMATING LINEAR TRIANGULAR SYSTEMS,WITH APPLICATIONS TO ASSET PRICING MODELS THAT INCLUDE A MISMEASURED FACTOR

A new estimator is proposed for linear triangular systems, where identification results from the model errors following a bivariate and diagonal GARCH(1,1) process with potentially time‐varying error covariances. This estimator applies when traditional instruments are unavailable. I demonstrate its usefulness on asset pricing models like the capital asset pricing model and Fama–French three‐factor model. In the context of a standard two‐pass cross‐sectional regression approach, this estimator improves the pricing performance of both models. Set identification bounds and an associated estimator are also provided for cases where the conditions supporting point identification fail. Copyright © 2013 John Wiley & Sons, Ltd.     

Asymptotics and smoothing parameter selection for penalized spline regression with various loss functions

Penalized splines are used in various types of regression analyses, including non‐parametric quantile, robust and the usual mean regression. In this paper, we focus on the penalized spline estimator with general convex loss functions. By specifying the loss function, we can obtain the mean estimator, quantile estimator and robust estimator. We will first study the asymptotic properties of penalized splines. Specifically, we will show the asymptotic bias and variance as well as the asymptotic normality of the estimator. Next, we will discuss smoothing parameter selection for the minimization of the mean integrated squares error. The new smoothing parameter can be expressed uniquely using the asymptotic bias and variance of the penalized spline estimator. To validate the new smoothing parameter selection method, we will provide a simulation. The simulation results show that the consistency of the estimator with the proposed smoothing parameter selection method can be confirmed and that the proposed estimator has better behavior than the estimator with generalized approximate cross‐validation. A real data example is also addressed.     

Survey Estimates by Calibration on Complex Auxiliary Information

In the last decade, calibration estimation has developed into an important field of research in survey sampling. Calibration is now an important methodological instrument in the production of statistics. Several national statistical agencies have developed software designed to compute calibrated weights based on auxiliary information available in population registers and other sources. This paper reviews some recent progress and offers some new perspectives. Calibration estimation can be used to advantage in a range of different survey conditions. This paper examines several situations, including estimation for domains in one‐phase sampling, estimation for two‐phase sampling, and estimation for two‐stage sampling with integrated weighting. Typical of those situations is complex auxiliary information, a term that we use for information made up of several components. An example occurs when a two‐stage sample survey has information both for units and for clusters of units, or when estimation for domains relies on information from different parts of the population. Complex auxiliary information opens up more than one way of computing the final calibrated weights to be used in estimation. They may be computed in a single step or in two or more successive steps. Depending on the approach, the resulting estimates do differ to some degree. All significant parts of the total information should be reflected in the final weights. The effectiveness of the complex information is mirrored by the variance of the resulting calibration estimator. Its exact variance is not presentable in simple form. Close approximation is possible via the corresponding linearized statistic. We define and use automated linearization as a shortcut in finding the linearized statistic. Its variance is easy to state, to interpret and to estimate. The variance components are expressed in terms of residuals, similar to those of standard regression theory. Visual inspection of the residuals reveals how the different components of the complex auxiliary information interact and work together toward reducing the variance.     

On solving bias‐corrected non‐linear estimation equations with an application to the dynamic linear model

In a seminal paper, Mak, Journal of the Royal Statistical Society B, 55, 1993, 945, derived an efficient algorithm for solving non‐linear unbiased estimation equations. In this paper, we show that when Mak's algorithm is applied to biased estimation equations, it results in the estimates that would come from solving a bias‐corrected estimation equation, making it a consistent estimator if regularity conditions hold. In addition, the properties that Mak established for his algorithm also apply in the case of biased estimation equations but for estimates from the bias‐corrected equations. The marginal likelihood estimator is obtained when the approach is applied to both maximum likelihood and least squares estimation of the covariance matrix parameters in the general linear regression model. The new approach results in two new estimators when applied to the profile and marginal likelihood functions for estimating the lagged dependent variable coefficient in the dynamic linear regression model. Monte Carlo simulation results show the new approach leads to a better estimator when applied to the standard profile likelihood. It is therefore recommended for situations in which standard estimators are known to be biased.     

Efficient estimation of factor models with time and cross‐sectional dependence

This paper studies the efficient estimation of large‐dimensional factor models with both time and cross‐sectional dependence assuming (N,T) separability of the covariance matrix. The asymptotic distribution of the estimator of the factor and factor‐loading space under factor stationarity is derived and compared to that of the principal component (PC) estimator. The paper also considers the case when factors exhibit a unit root. We provide feasible estimators and show in a simulation study that they are more efficient than the PC estimator in finite samples. In application, the estimation procedure is employed to estimate the Lee–Carter model and life expectancy is forecast. The Dutch gender gap is explored and the relationship between life expectancy and the level of economic development is examined in a cross‐country comparison.     

Small‐sample bias in synthetic cohort models of labor supply

This paper investigates small‐sample biases in synthetic cohort models (repeated cross‐sectional data grouped at the cohort and year level) in the context of a female labor supply model. I use the Current Population Survey to compare estimates when group sizes are extremely large to those that arise from randomly drawing subsamples of observations from the large groups. I augment this approach with Monte Carlo analysis so as to precisely quantify biases and coverage rates. In this particular application, thousands of observations per group are required before small‐sample issues can be ignored in estimation and sampling error leads to large downward biases in the estimated income elasticity. Copyright © 2007 John Wiley & Sons, Ltd.     

Differences Between Classical and Bayesian Estimates for Mixed Logit Models: A Replication Study

The mixed logit model is widely used in applied econometrics. Researchers typically rely on the free choice between the classical and Bayesian estimation approach. However, empirical evidence of the similarity of their parameter estimates is sparse. The presumed similarity is mainly based on one empirical study that analyzes a single dataset (Huber J, Train KE. 2001. On the similarity of classical and Bayesian estimates of individual mean partworths. Marketing Letters 12 (3): 259–269). Our replication study offers a generalization of their results by comparing classical and Bayesian parameter estimates from six additional datasets and specifically for panel versus cross‐sectional data. In general, our results suggest that the two methods provide similar results, with less similarity for cross‐sectional data than for panel data. Copyright © 2016 John Wiley & Sons, Ltd.     

Successive Sampling to Estimate Quantiles with P-Auxiliary Variables

The successive sampling is a known technique that can be used in longitudinal surveys to estimate population parameters and measurements of difference or change of a study variable. The paper discusses the estimation of quantiles for the current occasion based on sampling in two successive occasions and using p-auxiliary variables obtained of the previous occasion. A multivariate ratio estimator from the matched portion is used to provide the optimum estimate of a quantile by weighting the estimates inversely to derived optimum weights. Its properties are studied under large–sample approximation and the expressions of the variances are established. The behavior of these asymptotic variances is analyzed on the basis of data from natural populations. A simulation study is also used to measure the precision of the proposed estimator.     

On the Construction of Imputation Classes in Surveys

This paper explores the problem of the construction of imputation classes using the score method, sometimes called predictive mean stratification or response propensity stratification, depending on the context. This method was studied in Thomsen (1973) , Little (1986) and Eltinge & Yansaneh (1997) . We use a different framework to evaluate the properties of the resulting imputed estimator of a population mean. In our framework, we condition on the realized sample. This enables us to considerably simplify our theoretical developments in the frequent situation where the boundaries and the number of classes are sample‐dependent. We find that the key factor for reducing the non‐response bias is to form classes homogeneous with respect to the response probabilities and/or the conditional expectation of the variable of interest. In the latter case, the non‐response/imputation variance is also reduced. Finally, we performed a simulation study to fully evaluate various versions of the score method and to compare them with a cross‐classification method, which is frequently used in practice. The results showed the superiority of the score method in general.     

Multiple imputation of incomplete zero‐inflated count data

Empirical count data are often zero‐inflated and overdispersed. Currently, there is no software package that allows adequate imputation of these data. We present multiple‐imputation routines for these kinds of count data based on a Bayesian regression approach or alternatively based on a bootstrap approach that work as add‐ons for the popular multiple imputation by chained equations (mice ) software in R (van Buuren and Groothuis‐Oudshoorn , Journal of Statistical Software, vol. 45, 2011, p. 1). We demonstrate in a Monte Carlo simulation that our procedures are superior to currently available count data procedures. It is emphasized that thorough modeling is essential to obtain plausible imputations and that model mis‐specifications can bias parameter estimates and standard errors quite noticeably. Finally, the strengths and limitations of our procedures are discussed, and fruitful avenues for future theory and software development are outlined.     

Variable selection via composite quantile regression with dependent errors

We propose composite quantile regression for dependent data, in which the errors are from short‐range dependent and strictly stationary linear processes. Under some regularity conditions, we show that composite quantile estimator enjoys root‐n consistency and asymptotic normality. We investigate the asymptotic relative efficiency of composite quantile estimator to both single‐level quantile regression and least‐squares regression. When the errors have finite variance, the relative efficiency of composite quantile estimator with respect to the least‐squares estimator has a universal lower bound. Under some regularity conditions, the adaptive least absolute shrinkage and selection operator penalty leads to consistent variable selection, and the asymptotic distribution of the non‐zero coefficient is the same as that of the counterparts obtained when the true model is known. We conduct a simulation study and a real data analysis to evaluate the performance of the proposed approach.     

Common correlated effects estimation of heterogeneous dynamic panel quantile regression models

This paper proposes a quantile regression estimator for a heterogeneous panel model with lagged dependent variables and interactive effects. The paper adopts the Common Correlated Effects (CCE) approach proposed in the literature and demonstrates that the extension to the estimation of dynamic quantile regression models is feasible under similar conditions to the ones used in the literature. The new quantile regression estimator is shown to be consistent and its asymptotic distribution is derived. Monte Carlo studies are carried out to study the small sample behavior of the proposed approach. The evidence shows that the estimator can significantly improve on the performance of existing estimators as long as the time series dimension of the panel is large. We present an application to the evaluation of Time-of-Use pricing using a large randomized control trial.     

Bias Reduction in Dynamic Panel Data Models by Common Recursive Mean Adjustment*

The within‐group estimator (same as the least squares dummy variable estimator) of the dominant root in dynamic panel regression is known to be biased downwards. This article studies recursive mean adjustment (RMA) as a strategy to reduce this bias for AR(p) processes that may exhibit cross‐sectional dependence. Asymptotic properties for N,T→∞ jointly are developed. When ( log ²T)(N/T)→ζ, where ζ is a non‐zero constant, the estimator exhibits nearly negligible inconsistency. Simulation experiments demonstrate that the RMA estimator performs well in terms of reducing bias, variance and mean square error both when error terms are cross‐sectionally independent and when they are not. RMA dominates comparable estimators when T is small and/or when the underlying process is persistent.     

The Impact of Short‐ and Long‐run Exchange Rate Uncertainty on Investment: A Panel Study of Industrial Countries*

We examine the relationship between aggregate investment and exchange rate uncertainty in the G7, using panel estimation and decomposition of volatility derived from the components generalized autoregressive conditionally heteroscedastic (GARCH) model. Our dynamic panel approach takes account of potential cross‐sectional heterogeneity, which can lead to bias in estimation. We find that for a poolable subsample of European countries, it is the transitory and not the permanent component of volatility which adversely affects investment. To the extent that short‐run uncertainty in the CGARCH model characterizes higher frequency shocks generated by volatile short‐term capital flows, these are most deleterious for investment.     

Spline-based quasi-likelihood estimation of mixed Poisson regression with single-index models

We consider spline-based quasi-likelihood estimation for mixed Poisson regression with single-index models. The unknown smooth function is approximated by B-splines, and a modified Fisher scoring algorithm is employed to compute the estimates. The spline estimate of the nonparametric component is shown to achieve the optimal rate of convergence, and the asymptotic normality of the regression parameter estimates is still valid even if the variance function is misspecified. The semiparametric efficiency of the model can be established if the variance function is correctly specified. The variance of the regression parameter estimates can be consistently estimated by a simple procedure based on the least-squares estimation. The proposed method is evaluated via an extensive Monte Carlo study, and the methodology is illustrated on an air pollution study.     

Three Sides of Smoothing: Categorical Data Smoothing, Nonparametric Regression, and Density Estimation

The past forty years have seen a great deal of research into the construction and properties of nonparametric estimates of smooth functions. This research has focused primarily on two sides of the smoothing problem: nonparametric regression and density estimation. Theoretical results for these two situations are similar, and multivariate density estimation was an early justification for the Nadaraya-Watson kernel regression estimator.
A third, less well-explored, strand of applications of smoothing is to the estimation of probabilities in categorical data. In this paper the position of categorical data smoothing as a bridge between nonparametric regression and density estimation is explored. Nonparametric regression provides a paradigm for the construction of effective categorical smoothing estimates, and use of an appropriate likelihood function yields cell probability estimates with many desirable properties. Such estimates can be used to construct regression estimates when one or more of the categorical variables are viewed as response variables. They also lead naturally to the construction of well-behaved density estimates using local or penalized likelihood estimation, which can then be used in a regression context. Several real data sets are used to illustrate these points.