Identification and nonparametric estimation of a transformed additively separable model

Let r(x,z)$r(x,z)$ be a function that, along with its derivatives, can be consistently estimated nonparametrically. This paper discusses the identification and consistent estimation of the unknown functions H$H$, M$M$, G$G$ and F$F$, where r(x,z)=H[M(x,z)]$r(x,z)=H\left[M(x,z)\right]$, M(x,z)=G(x)+F(z)$M(x,z)=G\left(x\right)+F\left(z\right)$, and H$H$ is strictly monotonic. An estimation algorithm is proposed for each of the model’s unknown components when r(x,z)$r(x,z)$ represents a conditional mean function. The resulting estimators use marginal integration to separate the components G$G$ and F$F$. Our estimators are shown to have a limiting Normal distribution with a faster rate of convergence than unrestricted nonparametric alternatives. Their small sample performance is studied in a Monte Carlo experiment. We apply our results to estimate generalized homothetic production functions for four industries in the Chinese economy.     

The first-order approach when the cost of effort is money

We provide sufficient conditions for the first-order approach in the principal-agent problem when the agent’s utility has the nonseparable form u(y−c(a))$u(y-c\left(a\right))$ where y$y$ is the contractual payoff and c(a)$c\left(a\right)$ is the money cost of effort. We first consider a decision-maker facing prospects which cost c(a)$c\left(a\right)$ and with distributions of returns y$y$ that depend on a$a$. The decision problem is shown to be concave if the primitive of the cdf of returns is jointly convex in a$a$ and y$y$, a condition we call Concavity of the Cumulative Quantile (CCQ) and which is satisfied by many common distributions. Next we apply CCQ to the distribution of outcomes (or their likelihood-ratio transforms) in the principal-agent problem and derive restrictions on the utility function that validate the first-order approach. We also discuss another condition, log-convexity of the distribution, and show that it allows binding limited liability constraints, which CCQ does not.     

Treatment effects in sample selection models and their nonparametric estimation

In a sample-selection model with the ‘selection’ variable Q$Q$ and the ‘outcome’ variable Y^∗$Y^{\ast }$, Y^∗$Y^{\ast }$ is observed only when Q=1$Q=1$. For a treatment D$D$ affecting both Q$Q$ and Y^∗$Y^{\ast }$, three effects are of interest: ‘participation  ’ (i.e., the selection) effect of D$D$ on Q$Q$, ‘visible performance  ’ (i.e., the observed outcome) effect of D$D$ on Y≡QY^∗$Y\equiv Q{Y}^{\ast }$, and ‘invisible performance  ’ (i.e., the latent outcome) effect of D$D$ on Y^∗$Y^{\ast }$. This paper shows the conditions under which the three effects are identified, respectively, by the three corresponding mean differences of Q$Q$, Y$Y$, and Y|Q=1$Y|Q=1$ (i.e., Y^∗|Q=1$Y^{\ast }|Q=1$) across the control (D=0$D=0$) and treatment (D=1$D=1$) groups. Our nonparametric estimators for those effects adopt a two-sample framework and have several advantages over the usual matching methods. First, there is no need to select the number of matched observations. Second, the asymptotic distribution is easily obtained. Third, over-sampling the control/treatment group is allowed. Fourth, there is a built-in mechanism that takes into account the ‘non-overlapping support problem’, which the usual matching deals with by choosing a ‘caliper’. Fifth, a sensitivity analysis to gauge the presence of unobserved confounders is available. A simulation study is conducted to compare the proposed methods with matching methods, and a real data illustration is provided.     

Quasi ML estimation of the panel AR(1) model with arbitrary initial conditions

In this paper we show that the Quasi ML estimation method yields consistent Random and Fixed Effects estimators for the autoregression parameter ρ$\rho$ in the panel AR(1) model with arbitrary initial conditions and possibly time-series heteroskedasticity even when the error components are drawn from heterogeneous distributions. We investigate both analytically and by means of Monte Carlo simulations the properties of the QML estimators for ρ$\rho$. The RE(Q)MLE for ρ$\rho$ is asymptotically at least as robust to individual heterogeneity and, when the data are i.i.d. and normal, at least as efficient as the FE(Q)MLE for ρ$\rho$. Furthermore, the QML estimators for ρ$\rho$ only suffer from a ‘weak moment conditions’ problem when ρ$\rho$ is close to one if the cross-sectional average of the variances of the errors is (almost) constant over time, e.g. under time-series homoskedasticity. However, in this case the QML estimators for ρ$\rho$ are still consistent when ρ$\rho$ is local to or equal to one although they converge to a non-normal possibly asymmetric distribution at a rate that is lower than N^1/2$N^{1/2}$ but at least N^1/4$N^{1/4}$. Finally, we study the finite sample properties of two types of estimators for the standard errors of the QML estimators for ρ$\rho$, and the bounds of QML based confidence intervals for ρ$\rho$.     

Semiparametric binary regression models under shape constraints with an application to Indian schooling data

We consider estimation of the regression function in a semiparametric binary regression model defined through an appropriate link function (with emphasis on the logistic link) using likelihood-ratio based inversion. The dichotomous response variable Δ$\Delta$ is influenced by a set of covariates that can be partitioned as (X,Z)$(X,Z)$ where Z$Z$ (real valued) is the covariate of primary interest and X$X$ (vector valued) denotes a set of control variables. For any fixed X$X$, the conditional probability of the event of interest (Δ=1$\Delta =1$) is assumed to be a non-decreasing function of Z$Z$. The effect of the control variables is captured by a regression parameter β$\beta$. We show that the baseline conditional probability function (corresponding to X=0$X=0$) can be estimated by isotonic regression procedures and develop a likelihood ratio based method for constructing asymptotic confidence intervals for the conditional probability function (the regression function) that avoids the need to estimate nuisance parameters. Interestingly enough, the calibration of the likelihood ratio based confidence sets for the regression function no longer involves the usual χ²${\chi }^2$ quantiles, but those of the distribution of a new random variable that can be characterized as a functional of convex minorants of Brownian motion with quadratic drift. Confidence sets for the regression parameter β$\beta$ can however be constructed using asymptotically χ²${\chi }^2$ likelihood ratio statistics. The finite sample performance of the methods are assessed via a simulation study. The techniques of the paper are applied to data sets on primary school attendance among children belonging to different socio-economic groups in rural India.     

Panel unit root tests in the presence of a multifactor error structure

This paper extends the cross-sectionally augmented panel unit root test (CIPS$CIPS$) proposed by Pesaran (2007) to the case of a multifactor error structure, and proposes a new panel unit root test based on a simple average of cross-sectionally augmented Sargan–Bhargava statistics (CSB$CSB$). The basic idea is to exploit information regarding the m$m$ unobserved factors that are shared by k$k$ observed time series in addition to the series under consideration. Initially, we develop the tests assuming that m⁰$m^0$, the true number of factors, is known and show that the limit distribution of the tests does not depend on any nuisance parameters, so long as k≥m⁰−1$k\ge m^0-1$. Small sample properties of the tests are investigated by Monte Carlo experiments and are shown to be satisfactory. Particularly, the proposed CIPS$CIPS$ and CSB$CSB$ tests have the correct size for all   combinations of the cross section (N$N$) and time series (T$T$) dimensions considered. The power of both tests rises with N$N$ and T$T$, although the CSB$CSB$ test performs better than the CIPS$CIPS$ test for smaller sample sizes. The various testing procedures are illustrated with empirical applications to real interest rates and real equity prices across countries.     

Estimating deterministic trends with an integrated or stationary noise component

We propose a test for the slope of a trend function when it is a priori unknown whether the series is trend-stationary or contains an autoregressive unit root. The procedure is based on a Feasible Quasi Generalized Least Squares method from an AR(1) specification with parameter α$\alpha$, the sum of the autoregressive coefficients. The estimate of α$\alpha$ is the OLS estimate obtained from an autoregression applied to detrended data and is truncated to take a value 1 whenever the estimate is in a T^−δ$T^{-\delta }$ neighborhood of 1. This makes the estimate “super-efficient” when α=1$\alpha =1$ and implies that inference on the slope parameter can be performed using the standard Normal distribution whether α=1$\alpha =1$ or |α|<1$\left|\alpha \right|<1$. Theoretical arguments and simulation evidence show that δ=1/2$\delta =1/2$ is the appropriate choice. Simulations show that our procedure has better size and power properties than the tests proposed by [Bunzel, H., Vogelsang, T.J., 2005. Powerful trend function tests that are robust to strong serial correlation with an application to the Prebish–Singer hypothesis. Journal of Business and Economic Statistics 23, 381–394] and [Harvey, D.I., Leybourne, S.J., Taylor, A.M.R., 2007. A simple, robust and powerful test of the trend hypothesis. Journal of Econometrics 141, 1302–1330].     

Efficiency bounds for estimating linear functionals of nonparametric regression models with endogenous regressors

Let Y=μ^∗(X)+ε$Y={\mu }^{\ast }\left(X\right)+\epsilon$, where μ^∗${\mu }^{\ast }$ is unknown and E[ε|X]≠0$\mathbb{E}\left[\epsilon \right|X]\ne 0$ with positive probability but there exist instrumental variables W$W$ such that E[ε|W]=0$\mathbb{E}\left[\epsilon \right|W]=0$ w.p.1. It is well known that such nonparametric regression models are generally “ill-posed” in the sense that the map from the data to μ^∗${\mu }^{\ast }$ is not continuous. In this paper, we derive the efficiency bounds for estimating certain linear functionals of μ^∗${\mu }^{\ast }$ without assuming μ^∗${\mu }^{\ast }$ itself to be identified.     

On the properties of the coefficient of determination in regression models with infinite variance variables

We examine the asymptotic properties of the coefficient of determination, R²$R^2$, in models with α-stable$\alpha \text{-stable}$   random variables. If the regressor and error term share the same index of stability α<2$\alpha <2$, we show that the R²$R^2$  statistic does not converge to a constant but has a nondegenerate distribution on the entire [0,1]$[0,1]$ interval. We provide closed-form expressions for the cumulative distribution function and probability density function of this limit random variable, and we show that the density function is unbounded at 0 and 1. If the indices of stability of the regressor and error term are unequal, we show that the coefficient of determination converges in probability to either 0 or 1, depending on which variable has the smaller index of stability, irrespective of the value of the slope coefficient. In an empirical application, we revisit the Fama and MacBeth (1973) two-stage regression and demonstrate that in the infinite-variance case the R²$R^2$  statistic of the second-stage regression converges to 0 in probability even if the slope coefficient is nonzero. We deduce that a small value of the R²$R^2$  statistic should not, in itself, be used to reject the usefulness of a regression model.     

Beyond panel unit root tests: Using multiple testing to determine the nonstationarity properties of individual series in a panel

Most panel unit root tests are designed to test the joint null hypothesis of a unit root for each individual series in a panel. After a rejection, it will often be of interest to identify which series can be deemed to be stationary and which series can be deemed nonstationary. Researchers will sometimes carry out this classification on the basis of n$n$ individual (univariate) unit root tests based on some ad hoc significance level. In this paper, we suggest and demonstrate how to use the false discovery rate (FDR)$\left(FDR\right)$ in evaluating I(1)/I(0)$I\left(1\right)/I\left(0\right)$ classifications.     

A simple,robust and powerful test of the trend hypothesis

In this paper we develop a simple test procedure for a linear trend which does not require knowledge of the form of serial correlation in the data, is robust to strong serial correlation, and has a standard normal limiting null distribution under either I(0)$I(0)$ or I(1)$I(1)$ shocks. In contrast to other available robust linear trend tests, our proposed test achieves the Gaussian asymptotic local power envelope in both the I(0)$I(0)$ and I(1)$I(1)$ cases. For near-I(1)$I(1)$ errors our proposed procedure is conservative and a modification for this situation is suggested. An estimator of the trend parameter, together with an associated confidence interval, which is asymptotically efficient, again regardless of whether the shocks are I(0)$I(0)$ or I(1)$I(1)$, is also provided.     

One-step R-estimation in linear models with stable errors

Classical estimation techniques for linear models either are inconsistent, or perform rather poorly, under α$\alpha$-stable error densities; most of them are not even rate-optimal. In this paper, we propose an original one-step R-estimation method and investigate its asymptotic performances under stable densities. Contrary to traditional least squares, the proposed R-estimators remain root-n$n$ consistent (the optimal rate) under the whole family of stable distributions, irrespective of their asymmetry and tail index. While parametric stable-likelihood estimation, due to the absence of a closed form for stable densities, is quite cumbersome, our method allows us to construct estimators reaching the parametric efficiency bounds associated with any prescribed values (_α0,_b0)$({\alpha }_0,b_0)$ of the tail index α$\alpha$ and skewness parameter b$b$, while preserving root-n$n$ consistency under any (α,b)$(\alpha ,b)$ as well as under usual light-tailed densities. The method furthermore avoids all forms of multidimensional argmin computation. Simulations confirm its excellent finite-sample performances.