A CONSISTENT MODEL SPECIFICATION TEST FOR A REGRESSION FUNCTION BASED ON NONPARAMETRIC WAVELET ESTIMATION

In this paper we present a consistent specification test of a parametric regression function against a general nonparametric alternative. The proposed test is based on wavelet estimation and it is shown to have similar rates of convergence to the more commonly used kernel based tests. Monte Carlo simulations show that this test statistic has adequate size and high power and that it compares favorably with its kernel based counterparts in small samples.     

ESTIMATING THE ROOT OF A NONPARAMETRIC REGRESSION FUNCTION IN A ROBUST FASHION

For a nonparametric regression model y = m(x)+e with n independent observations, we analyze a robust method of finding the root of m(x) based on an M-estimation first discussed by Härdle & Gasser (1984). It is shown here that the robustness properties (minimaxity and breakdown function) of such an estimate are quite analogous to those of an M -estimator in the simple location model, but the rate of convergence is somewhat limited due to the nonparametric nature of the problem.     

VARIANCE ESTIMATION IN THE ERROR COMPONENTS REGRESSION MODEL

ABSTRACT

In a regression model with a random individual and a random time effect explicit representations of the nonnegative quadratic minimum biased estimators of the corresponding variances are deduced. These estimators always exist and are unique. Moreover, under normality assumption of the dependent variable unbiased estimators of the mean squared errors of the variance estimates are derived. Finally, confidence intervals on the variance components are considered.     

DO NOT WEIGHT FOR HETEROSCEDASTICITY IN NONPARAMETRIC REGRESSION

The potential role of weighting in kernel regression is examined. The concept that weighting has something to do with heteroscedastic errors is shown to be false. However, weighting does affect bias, and ways in which this might be exploited are indicated.     

SADDLEPOINT METHODS FOR BOOTSTRAP CONFIDENCE BANDS IN NONPARAMETRIC REGRESSION

Bootstrap techniques have been used to construct confidence bands in nonparametric regression problems (Härdle & Bowman, 1988). Yet the required simulation is generally computationally intensive and therefore makes it difficult to conduct further investigations. In this paper, two saddlepoint methods are considered as alternatives to the naive simulation procedure. Some improvements to Härdle & Bowman's bootstrap method are suggested. The improvements are numerically verified using these efficient and accurate analytic methods.     

RANDOM PERMUTATION TESTING IN MULTIPLE LINEAR REGRESSION

We developed an alternative random permutation testing method for multiple linear regression, which is an improvement over the existing one proposed by [1] Kennedy, P. E. 1995. Randomization tests in econometrics. Journal of Business and Economic Statistics, 13: 85–94. [Taylor & Francis Online], [Web of Science ®] , [Google Scholar] or [2] Freedman, D. and Lane, D. 1983. A nonstochastic interpretation of reported significance levels. Journal of Business and Economic Statistics, 1: 292–298. [Taylor & Francis Online] , [Google Scholar].     

NONPARAMETRIC ESTIMATION OF ANGULAR DISTRIBUTIONS FROM CARTESIAN OBSERVATIONS WITH ERROR

If angular data are obtained from Cartesian observations, then any measurement error in these observations will produce a particular error structure in the angular data. The paper shows how non-parametric density estimation by orthogonal series may be performed in this case.     

MODEL-BASED VARIANCE ESTIMATION IN SURVEYS WITH STRATIFIED CLUSTERED DESIGN

A model-based method for estimating the sampling variances of estimators of (sub-)population means, proportions, quantiles, and regression parameters in surveys with stratified clustered design is described and applied to a survey of US secondary education. The method is compared with the jackknife by a simulation study. The model-based estimators of the sampling variances have much smaller mean squared errors than their jackknife counterparts. In addition, they can be improved by incorporating information about the unknown parameters (variances) from external sources. A regression-based smoothing method for estimating the sampling variances of the estimators for a large number of subpopulation means is proposed. Such smoothing may be invaluable when subpopulations are represented in the sample by only few subjects.     

BOOTSTRAP TESTS FOR THE ERROR DISTRIBUTION IN LINEAR AND NONPARAMETRIC REGRESSION MODELS

In this paper we investigate several tests for the hypothesis of a parametric form of the error distribution in the common linear and non‐parametric regression model, which are based on empirical processes of residuals. It is well known that tests in this context are not asymptotically distribution‐free and the parametric bootstrap is applied to deal with this problem. The performance of the resulting bootstrap test is investigated from an asymptotic point of view and by means of a simulation study. The results demonstrate that even for moderate sample sizes the parametric bootstrap provides a reliable and easy accessible solution to the problem of goodness‐of‐fit testing of assumptions regarding the error distribution in linear and non‐parametric regression models.     

THE CONSTRUCTION OF EQUILEVERAGE DESIGNS FOR MULTIPLE LINEAR REGRESSION

The hat matrix is widely used as a diagnostic tool in linear regression because it contains the leverages which the independent variables exert on the fitted values. In some experiments, cases with high leverage may be avoided by judicious choice of design for the independent variables. A variety of methods for constructing equileverage designs for linear regression are discussed. Such designs remove one of the factors, namely large leverage points, which can lead to nonrobust estimators and tests. In addition, a method is given for combining equileverage designs to test for lack of fit of the linear model.     

EXPLICIT NONPARAMETRIC CONFIDENCE INTERVALS FOR THE VARIANCE WITH GUARANTEED COVERAGE

ABSTRACT

In this paper, we provide a method for constructing confidence intervals for the variance which exhibits guaranteed coverage probability for any sample size, uniformly over a wide class of probability distributions. In contrast, standard methods achieve guaranteed coverage only in the limit for a fixed distribution or for any sample size over a very restrictive (parametric) class of probability distributions. Of course, it is impossible to construct effective confidence intervals for the variance without some restriction, due to a result of Bahadur and Savage.[1] Bahadur, R. and Savage, L. 1956. The Nonexistence of Certain Statistical Procedures in Nonparametric Problems. Annals of Mathematical Statistics, 25: 1115–1122.  [Google Scholar] However, it is possible if the observations lie in a fixed compact set. We also consider the case of lower confidence bounds without any support restriction. Our method is based on the behavior of the variance over distributions that lie within a Kolmogorov–Smirnov confidence band for the underlying distribution. The method is a generalization of an idea of Anderson,[2] Anderson, T. 1967. Confidence Limits for the Expected Value of an Arbitrary Bounded Random Variable with a Continuous Distribution Function. Bull. ISI, 43: 249–251.  [Google Scholar] who considered only the case of the mean; it applies to very general parameters, and particularly the variance. While typically it is not clear how to compute these intervals explicitly, for the special case of the variance we provide an algorithm to do so. Asymptotically, the length of the intervals is of order n ^?/2 (in probability), so that, while providing guaranteed coverage, they are not overly conservative. A small simulation study examines the finite sample behavior of the proposed intervals.     

EFFICIENT ESTIMATION IN REGRESSION FRAILTY OR TRANSFORMATION MODELS BASED ON AN ALGORITHM

This paper considers a regression frailty or transformation model in which the structural parameter is the vector of regression coefficients and the nuisance parameter is a vector of arbitrarily high dimension. It proposes jointly (implicitly) defined parameter estimators which have been proved to be consistent and asymptotically efficient, and develops an algorithmic procedure that provides these estimators. The behaviour of the algorithm is illustrated by analysing simulated and real data.     

ESTIMATION OF PARAMETERS IN THE MARGINAL FISHER DISTRIBUTION1

The Fisher distribution is frequently used as a model for the probability distribution of directional data, which may be specified either in terms of unit vectors or angular co-ordinates (co-latitude and azimuth). If, in practical situations, only the co-latitudes can be observed, the available data must be regarded as a sample from the corresponding marginal distribution. This paper discusses the estimation by Maximum Likelihood (ML) and the Method of Moments of the two parameters of this marginal Fisher distribution. The moment estimators are generally simpler to compute than the ML estimators, and have high asymptotic efficiency.     

ON VARIANCE ESTIMATION IN UNEQUAL PROBABILITY SAMPLING1

The purpose of this article is to propose a model-based estimator of the variance of the Horvitz-Thompson estimator. Empirical investigations reveal that the estimator is seldom greatly biased and is quite satisfactory from the stability point of view.     

SEPARATE VERSUS SYSTEM METHODS OF STEIN-RULE ESTIMATION IN SEEMINGLY UNRELATED REGRESSION MODELS

ABSTRACT

Despite the sizeable literature associated with the seemingly unrelated regression models, not much is known about the use of Stein-rule estimators in these models. This gap is remedied in this paper, in which two families of Stein-rule estimators in seemingly unrelated regression equations are presented and their large sample asymptotic properties explored and evaluated. One family of estimators uses a shrinkage factor obtained solely from the equation under study while the other has a shrinkage factor based on all the equations of the model. Using a quadratic loss measure and Monte-Carlo sampling experiments, the finite sample risk performance of these estimators is also evaluated and compared with the traditional feasible generalized least squares estimator.     

A MORE GENERAL CRITERION FOR SUBSET SELECTION IN MULTIPLE LINEAR REGRESSION

ABSTRACT

In this article, we propose a more general criterion called S_p -criterion, for subset selection in the multiple linear regression Model. Many subset selection methods are based on the Least Squares (LS) estimator of β, but whenever the data contain an influential observation or the distribution of the error variable deviates from normality, the LS estimator performs ‘poorly’ and hence a method based on this estimator (for example, Mallows’ C_p -criterion) tends to select a ‘wrong’ subset. The proposed method overcomes this drawback and its main feature is that it can be used with any type of estimator (either the LS estimator or any robust estimator) of β without any need for modification of the proposed criterion. Moreover, this technique is operationally simple to implement as compared to other existing criteria. The method is illustrated with examples.     

Weighted local linear composite quantile estimation for the case of general error distributions

It is known that for nonparametric regression, local linear composite quantile regression (local linear CQR) is a more competitive technique than classical local linear regression since it can significantly improve estimation efficiency under a class of non-normal and symmetric error distributions. However, this method only applies to symmetric errors because, without symmetric condition, the estimation bias is non-negligible and therefore the resulting estimator is inconsistent. In this paper, we propose a weighted local linear CQR method for general error conditions. This method applies to both symmetric and asymmetric random errors. Because of the use of weights, the estimation bias is eliminated asymptotically and the asymptotic normality is established. Furthermore, by minimizing asymptotic variance, the optimal weights are computed and consequently the optimal estimate (the most efficient estimate) is obtained. By comparing relative efficiency theoretically or numerically, we can ensure that the new estimation outperforms the local linear CQR estimation. Finite sample behaviors conducted by simulation studies further illustrate the theoretical findings.     

A NONPARAMETRIC TEST OF CHANGING CONDITIONAL VARIANCES IN AUTOREGRESSIVE TIME SERIES

A nonparametric test for detecting changing conditional variances in stationary AR(p) time series is proposed in this paper. For AR(1) models, the test statistic is a Kolmogorov-Smirnov type statistic and the asymptotic theory is developed under both the null and the alternative hypotheses. For AR(p) models (p ≥ 2), an approximate test procedure is proposed. The empirical upper percentage points for our test are tabulated for both p = 1 and p = 2 cases and a bootstrap procedure is suggested for the p ≥ 3 case. Monte Carlo simulations demonstrate that the test has very good powers for finite samples under both normal and non-normal errors.     

LOGISTIC REGRESSION COMPARED TO NORMAL DISCRIMINATION FOR NON-NORMAL POPULATIONS'

An investigation is undertaken of the logistic regression procedure for estimating the posterior probability of an object belonging to one of two populations. The asymptotic bias and mean square error associated with the procedure are derived for univariate populations whose distributions satisfy the general Day-Kerridge model for which the logistic form is valid for the posterior probability. These properties are compared with those of the normal discrimination method based on the classical assumption of normal populations with common variances. The asymptotic relative efficiency of logistic regression is considered on the basis of asymptotic mean square error.