Berry-Esséen rate in asymptotic normality for perturbed sample quantiles

In estimating quantiles with a sample of sizeN obtained from a distributionF, the perturbed sample quantiles based on a kernel functionk have been investigated by many authors. It is well known that their behaviour depends on the choices of “window-width”, sayw _N. Under suitable and reasonably mild assumptions onF andk, Ralescu and Sun (1993) have recently proven that lim _N→∞ N ^1/4w_N=0 is the necessary and sufficient condition for the asymptotic normality of the perturbed sample quantiles. In this paper, their rate of convergence is investigated. It turns out that the optimal Berry-Esséen rate ofO(N^?1/2) can be achieved by choosing the window-width suitably, sayw _N=O(N^?1/2). The obtained results, in addition to being explicit enough to verify the sufficient condition for the asymptotic normality, improve Ralescu's (1992) result of which the rate is of order (logN)N ^?1/2.     

Two kinds of variance/covariance estimates in linear mixed models

For longitudinal data, the within-subject covariance matrix plays an important role in statistical inference and it is of great interest to investigate this. In the paper, two kinds of estimators are investigated for the random effect covariance matrix D ₁ and the error variance σ ² in linear mixed models. One is to estimate D ₁ first and then to estimate σ ²; the other kind is to estimate σ ² first and then for D ₁. Both kinds of estimators are consistent. The covariance matrices of these covariance estimators and the variances of these two error variance estimators are calculated. In particular, the mean square errors of these estimators are also derived for one dimensional random effects. Besides, a simulation study is conducted to investigate the performances of these estimators.     

Generalised residuals and heterogeneous duration models: With applications to the Weilbull model

This paper deals with models for the duration of an event that are misspecified by the neglect of random multiplicative heterogeneity in the hazard function. This type of misspecification has been widely discussed in the literature [e.g., Heckman and Singer (1982), Lancaster and Nickell (1980)], but no study of its effect on maximum likelihood estimators has been given. This paper aims to provide such a study with particular reference to the Weibull regression model which is by far the most frequently used parametric model [e.g., Heckman and Borjas (1980), Lancaster (1979)]. In this paper we define generalised errors and residuals in the sense of Cox and Snell (1968, 1971) and show how their use materially simplifies the analysis of both true and misspecified duration models. We show that multiplicative heterogeneity in the hazard of the Weibull model has two errors in variables interpretations. We give the exact asymptotic inconsistency of M.L. estimation in the Weibull model and give a general expression for the inconsistency of M.L. estimators due to neglected heterogeneity for any duration model to O(σ²), where σ² is the variance of the error term. We also discuss the information matrix test for neglected heterogeneity in duration models and consider its behaviour when σ²>0.     

Robustness of statistical inferences using linear models with meta-analytic correlation matrices

To examine complex relationships among variables, researchers in human resource management, industrial-organizational psychology, organizational behavior, and related fields have increasingly used meta-analytic procedures to aggregate effect sizes across primary studies to form meta-analytic correlation matrices, which are then subjected to further analyses using linear models (e.g., multiple linear regression). Because missing effect sizes (i.e., correlation coefficients) and different sample sizes across primary studies can occur when constructing meta-analytic correlation matrices, the present study examined the effects of missingness under realistic conditions and various methods for estimating sample size (e.g., minimum sample size, arithmetic mean, harmonic mean, and geometric mean) on the estimated squared multiple correlation coefficient (R²) and the power of the significance test on the overall R² in linear regression. Simulation results suggest that missing data had a more detrimental effect as the number of primary studies decreased and the number of predictor variables increased. It appears that using second-order sample sizes of at least 10 (i.e., independent effect sizes) can improve both statistical power and estimation of the overall R² considerably. Results also suggest that although the minimum sample size should not be used to estimate sample size, the other sample size estimates appear to perform similarly.     

A theory of coalition formation in committees

Let G = (N,W) be a strong weighted majority game and let A be a set of alternatives. Denote by L the set of linear orders on A. A social choice function F:L^N → A is a representation of G if the simple game G¹(F) associated with F equals G. A coalition S is determining in G if it satisfies the following condition. Let F be a representation of G and let R^N ? L^N. Then, if a simple majority of the members of S consider an alternative x to be their best choice, then S can ‘enforce’ x to be a Nash equilibrium payoff in the resulting non-cooperative voting game g(F,R^N). In this paper we generalize the above notion of a determining coalition to committees (i.e., proper and monotonic simple games), and give a complete characterization of the set of determining coalitions of a committee. Furthermore, we discuss our notion of a determining coalition in the light of some real-life data on formation of coalitions in town councils in Israel.     

Supplier evaluation and demand allocation among suppliers in a supply chain

This paper presents a hybrid algorithm that prioritizes the suppliers and then allocates the demand among the suppliers. The objective here is to maximize the total purchase value of the items taking into consideration budget constraint, demand condition, delivery lead-time and supplier capacity. Since the problem is multi-criteria decision making, we solve this problem by integrating the supplier rating with mixed linear integer programming method. The customer demand is allocated by using a hybrid algorithm based on the technique for order preference by similarity to ideal solution (TOPSIS) and the mixed linear integer programming (MILP) approaches. The effectiveness of the proposed algorithm is validated with computational results. Drawing to a case, a supplier S₃ is identified as the best supplier by using the TOPSIS method for demand allocation under no restrictions. On the contrary, under constrained scenario, supplier S₂ is selected as the best supplier by using the hybrid algorithm for demand allocation and maximum units are allocated to S₂.     

Blocked-Error-R 2: A Conceptually Improved Definition of the Proportion of Explained Variance in Models Containing Loops or Correlated Residuals

Bentler and Raykov (2000, Journal of Applied Psychology 85: 125–131), and Jöreskog (1999a, http://www.ssicentral.com/lisrel/column3.htm, 1999b http://www.ssicentral. com/lisrel/column5.htm) proposed procedures for calculating R ² for dependent variables involved in loops or possessing correlated errors. This article demonstrates that Bentler and Raykov’s procedure can not be routinely interpreted as a “proportion” of explained variance, while Jöreskog’s reduced-form calculation is unnecessarily restrictive. The new blocked-error-R ² (beR ²) uses a minimal hypothetical causal intervention to resolve the variance-partitioning ambiguities created by loops and correlated errors. Hayduk (1996) discussed how stabilising feedback models – models capable of counteracting external perturbations – can result in an acceptable error variance which exceeds the variance of the dependent variable to which that error is attached. For variables included within loops, whether stabilising or not, beR ² provides the same value as Hayduk’s (1996) loop-adjusted-R ². For variables not involved in loops and not displaying correlated residuals, beR ² reports the same value as the traditional regression R ². Thus, beR ² provides a conceptualisation of the proportion of explained variance that spans both recursive and nonrecursive structural equation models. A procedure for calculating beR ² in any SEM program is provided.     

An exact non-cooperative support for the sequential Raiffa solution

This article provides an exact non-cooperative foundation of the sequential Raiffa solution for two-person bargaining games. Based on an approximate foundation due to Myerson (1991) for any two-person bargaining game (S, d) an extensive form game G^S,d is defined that has an infinity of weakly subgame perfect equilibria whose payoff vectors coincide with that of the sequential Raiffa solution of (S, d). Moreover all those equilibria share the same equilibrium path consisting of proposing the Raiffa solution and accepting it in the first stage of the game.By a modification of G^S,d the analogous result is provided for subgame perfect equilibria. These results immediately extend to implementation of a sequential Raiffa (solution based) social choice rule in subgame perfect equilibrium.     

Asymptotic inference for a one-dimensional simultaneous autoregressive model

A nonstationary simultaneous autoregressive model \({X^{(n)}_k=\alpha \Big(X^{(n)}_{k-1}+X^{(n)}_{k+1}\Big)+\varepsilon_k, k=1, 2, \ldots , n-1}\), is investigated, where \({X^{(n)}_0}\) and \({X^{(n)}_n}\) are given random variables. It is shown that in the unstable case α = 1/2 the least squares estimator of the autoregressive parameter converges to a functional of a standard Wiener process with a rate of convergence n ², while in the stable situation |α| < 1/2 the estimator is biased but asymptotically normal with a rate n ^1/2.     

Quadratic semiparametric Von Mises calculus

We discuss a new method of estimation of parameters in semiparametric and nonparametric models. The method is based on U-statistics constructed from quadratic influence functions. The latter extend ordinary linear influence functions of the parameter of interest as defined in semiparametric theory, and represent second order derivatives of this parameter. For parameters for which the matching cannot be perfect the method leads to a bias-variance trade-off, and results in estimators that converge at a slower than n ^?1/2-rate. In a number of examples the resulting rate can be shown to be optimal. We are particularly interested in estimating parameters in models with a nuisance parameter of high dimension or low regularity, where the parameter of interest cannot be estimated at n ^?1/2-rate.     

Heteroskedasticity, autocorrelation, and spatial correlation robust inference in linear panel models with fixed-effects

This paper develops an asymptotic theory for test statistics in linear panel models that are robust to heteroskedasticity, autocorrelation and/or spatial correlation. Two classes of standard errors are analyzed. Both are based on nonparametric heteroskedasticity autocorrelation (HAC) covariance matrix estimators. The first class is based on averages of HAC estimators across individuals in the cross-section, i.e. “averages of HACs”. This class includes the well known cluster standard errors analyzed by Arellano (1987) as a special case. The second class is based on the HAC of cross-section averages and was proposed by Driscoll and Kraay (1998). The ”HAC of averages” standard errors are robust to heteroskedasticity, serial correlation and spatial correlation but weak dependence in the time dimension is required. The “averages of HACs” standard errors are robust to heteroskedasticity and serial correlation including the nonstationary case but they are not valid in the presence of spatial correlation. The main contribution of the paper is to develop a fixed-b asymptotic theory for statistics based on both classes of standard errors in models with individual and possibly time fixed-effects dummy variables. The asymptotics is carried out for large time sample sizes for both fixed and large cross-section sample sizes. Extensive simulations show that the fixed-b approximation is usually much better than the traditional normal or chi-square approximation especially for the Driscoll-Kraay standard errors. The use of fixed-b critical values will lead to more reliable inference in practice especially for tests of joint hypotheses.     

Nonparametric transformation to white noise

We consider a semiparametric distributed lag model in which the “news impact curve” m is nonparametric but the response is dynamic through some linear filters. A special case of this is a nonparametric regression with serially correlated errors. We propose an estimator of the news impact curve based on a dynamic transformation that produces white noise errors. This yields an estimating equation for m that is a type two linear integral equation. We investigate both the stationary case and the case where the error has a unit root. In the stationary case we establish the pointwise asymptotic normality. In the special case of a nonparametric regression subject to time series errors our estimator achieves efficiency improvements over the usual estimators, see Xiao et al. [2003. More efficient local polynomial estimation in nonparametric regression with autocorrelated errors. Journal of the American Statistical Association 98, 980–992]. In the unit root case our procedure is consistent and asymptotically normal unlike the standard regression smoother. We also present the distribution theory for the parameter estimates, which is nonstandard in the unit root case. We also investigate its finite sample performance through simulation experiments.     

Some Decompositions of OLSEs and BLUEs Under a Partitioned Linear Model

We consider in this paper a partitioned linear model { y , X ₁ β ₁ + X ₂ β ₂ , σ ² σ } and two corresponding small models { y , X ₁ β ₁ , σ ² σ } and { y , X ₂ β ₂ , σ ² σ } . We derive necessary and sufficient conditions for (i) the ordinary least squares estimator under the full model to be the sum of the ordinary least squares estimators under the two small models; (ii) the best linear unbiased estimator under the full model to be the sum of the best linear unbiased estimators under the two small models; (iii) the best linear unbiased estimator under the full model to be the sum of the ordinary least squares estimators under the two small models. The proofs of the main results in this paper also demonstrate how to use the matrix rank method for characterizing various equalities of estimators under general linear models.     

Estimation of the location parameter under LINEX loss function: multivariate case

The Baysian estimation of the mean vector θ of a p-variate normal distribution under linear exponential (LINEX) loss function is studied when as a special restricted model, it is suspected that for a p × r known matrix Z the hypothesis θ = Zβ, ${\beta\in\Re^r}The Baysian estimation of the mean vector θ of a p-variate normal distribution under linear exponential (LINEX) loss function is studied when as a special restricted model, it is suspected that for a p × r known matrix Z the hypothesis θ = Zβ, b ? ?^r{\beta\in\Re^r} may hold. In this area we show that the Bayes and empirical Bayes estimators dominate the unrestricted estimator (when nothing is known about the mean vector θ).     

Computation and analysis of multiple structural change models

In a recent paper, Bai and Perron ( 1998 ) considered theoretical issues related to the limiting distribution of estimators and test statistics in the linear model with multiple structural changes. In this companion paper, we consider practical issues for the empirical applications of the procedures. We first address the problem of estimation of the break dates and present an efficient algorithm to obtain global minimizers of the sum of squared residuals. This algorithm is based on the principle of dynamic programming and requires at most least‐squares operations of order O(T²) for any number of breaks. Our method can be applied to both pure and partial structural change models. Second, we consider the problem of forming confidence intervals for the break dates under various hypotheses about the structure of the data and the errors across segments. Third, we address the issue of testing for structural changes under very general conditions on the data and the errors. Fourth, we address the issue of estimating the number of breaks. Finally, a few empirical applications are presented to illustrate the usefulness of the procedures. All methods discussed are implemented in a GAUSS program. Copyright © 2002 John Wiley & Sons, Ltd.     

Estimating the autocorrelated error model with trended data

A Monte Carlo study of the small sample properties of various estimators of the linear regression model with first-order autocorrelated errors. When independent variables are trended, estimators using Ttransformed observations (Prais-Winsten) are much more efficient than those using T–1 (Cochrane–Orcutt). The best of the feasible estimators isiterated Prais-Winsten using a sum-of-squared-error minimizing estimate of the autocorrelation coefficient ?. None of the feasible estimators performs well in hypothesis testing; all seriously underestimate standard errors, making estimated coefficients appear to be much more significant than they actually are.     

Estimating quadratic variation consistently in the presence of endogenous and diurnal measurement error

We propose an econometric model that captures the effects of market microstructure on a latent price process. In particular, we allow for correlation between the measurement error and the return process and we allow the measurement error process to have a diurnal heteroskedasticity. We propose a modification of the TSRV estimator of quadratic variation. We show that this estimator is consistent, with a rate of convergence that depends on the size of the measurement error, but is no worse than n^−1/6$n^{-1/6}$. We investigate in simulation experiments the finite sample performance of various proposed implementations.     

Generalized variance-ratio tests for serial correlation in multivariate regression models

Exact tests for rth order serial correlation in the multivariate linear regression model are devised which are based on a multivariate generalization of the F-distribution. The tests require the computation of two multivariate regressions. In the special case of a single-equation regression model the procedures reduce to simple always-conclusive F-tests. The tests are illustrated by applications to the Rotterdam Model of consumer demand.     

Asymptotic distributions of some scale estimators in nonlinear models

Often in the robust analysis of regression and time series models there is a need for having a robust scale estimator of a scale parameter of the errors. One often used scale estimator is the median of the absolute residuals s ₁. It is of interest to know its limiting distribution and the consistency rate. Its limiting distribution generally depends on the estimator of the regression and/or autoregressive parameter vector unless the errors are symmetrically distributed around zero. To overcome this difficulty it is then natural to use the median of the absolute differences of pairwise residuals, s ₂, as a scale estimator. This paper derives the asymptotic distributions of these two estimators for a large class of nonlinear regression and autoregressive models when the errors are independent and identically distributed. It is found that the asymptotic distribution of a suitably standardizes s ₂ is free of the initial estimator of the regression/autoregressive parameters. A similar conclusion also holds for s ₁ in linear regression models through the origin and with centered designs, and in linear autoregressive models with zero mean errors.  This paper also investigates the limiting distributions of these estimators in nonlinear regression models with long memory moving average errors. An interesting finding is that if the errors are symmetric around zero, then not only is the limiting distribution of a suitably standardized s ₁ free of the regression estimator, but it is degenerate at zero. On the other hand a similarly standardized s ₂ converges in distribution to a normal distribution, regardless of the errors being symmetric or not. One clear conclusion is that under the symmetry of the long memory moving average errors, the rate of consistency for s ₁ is faster than that of s ₂.     

Properties of technical efficiency estimators in the stochastic frontier model

This paper considers the disturbance specification ε = v ? u of the stochastic frontier model. For v distributed zero-mean normal and u half normal or exponential, we evaluate the population correlation coefficients between u and three estimators of u, E(u|ε) and two linear estimators, for various values of the signal-to-noise ratio.