Using Ridge Least Median Squares to Estimate the Parameter by Solving Multicollinearity and Outliers Problems

Save this PDF as:

Size: px
Start display at page:

Download "Using Ridge Least Median Squares to Estimate the Parameter by Solving Multicollinearity and Outliers Problems"

Transcription

1 Modern Applied Science; Vol. 9, No. ; 05 ISSN E-ISSN 9-85 Published by Canadian Center of Science and Education Using Ridge Least Median Squares to Estimate the Parameter by Solving Multicollinearity and Outliers Problems Kafi Dano Pati, Robiah Adnan & Bello Abdulkadir Rasheed Department of Mathematical Sciences, Universiti Teknologi Malaysia, 80 Skudai, Johor, Malaysia Correspondence: Assoc. Prof. Dr. Robiah Adnan, Department of Mathematical Sciences, Faculty of Science, Universiti Teknologi Malaysia, 80 UTM Johor Bahru, Johor Bahru, Johor, Malaysia. Received: October 4, 04 Accepted: November, 04 Online Published: January 0, 05 doi:0.559/mas.v9np9 URL: Abstract In the multiple linear regression analysis, the ridge regression estimator is often used to address the problem of multicollinearity. Besides multicollinearity, outliers also constitute a problem in the multiple linear regression analysis. We propose a new biased estimators of the robust ridge regression called the Ridge Least Median Squares (RLMS) estimator in the presence of multicollinearity and outliers. For this purpose, a simulation study is conducted in order to see the difference between the proposed method and the existing methods in terms of their effectiveness measured by the mean square error. In our simulation studies the performance of the proposed method RLMS is examined for different number of observations and the different percentage of outliers. The results of different illustrative cases are presented. This paper also provides the results of the RLMS on a real-life data set. The results show that RLMS is better than the existing methods of Ordinary Least Squares (OLS), Ridge Least Absolute Value (RLAV) and Ridge Regression (RR) in the presence of multicollinearity and outliers. Keywords: multicollinearity, ridge regression, outliers, robust estimation, robust ridge methods. Introduction It is clear that the classical ordinary least squares (OLS) estimation has many statistical properties, most of the time when errors are normally distributed the variance of OLS estimation is the least of unbiased estimator. However, OLS estimator is often affected by two problems, namely multicollinearity and outliers which leads to the weakness of OLS estimates (Ma and Du, 04). In order to overcome the influence of multicollinearity many biased estimators are introduced. (Hoerl and Kennard, 970 a, b) First proposed a useful tool called ridge regression for improving prediction in regression situations with highly correlated variables. This method has been successfully applied in many fields. However, the ridge regression method is not robust as it is sensitive to observations that are considered as contaminated data, usually called outliers. The robust ridge regression estimation has been investigated, to remedy the effect of multicollinearity and outliers. The parameter estimation becomes more difficult when these problems are present in the data set. Different robust estimation methods have been used in the literature to estimate the parameter of regression model in the presence of outliers and multicollinearity. (Pariente, 977; Askin and Montgomery, 980; Pfaffenberger and Dielman, 984; Pfaffenberger, 985 and Nikolai Krivulin, 99), among others presented several new ideas concerning the LMS estimator so as to provide some theoretical framework for efficient regression algorithms. (Desire and Leonard, 986) introduced a procedure to determine the minimum redundancy of a measurement set and the non-redundant samples of measurement that are robust against outliers and multicollinearity. This paper proposes a new estimation method of a robust ridge estimator that has a high breakdown point called the robust ridge regression estimation method based on least median squares (LMS). The least median squares estimator has 50% breakdown point as proposed by (Rousseeuw, 984, 985), who replaced the least sum of squares (OLS) that has 0% breakdown point with the least median of squares LMS. (Hampel, 975), explains the breakdown point as the smallest percentage of contaminated data that can cause the estimator to take a random large abnormal value. 9

2 Modern Applied Science Vol. 9, No. ; 05 The aim of evaluating these alternatives biased robust methods is to find the estimator of the parameters of the model that is highly efficient and effective in the presence of multicollinearity and outliers.. The performance of the robust ridge estimators is examined by using the Standard Errors (SE) on a real data set presented by (Hald, 95). The remainder of the paper is organized as follows; after a brief review of ridge regression in Section, we define a robust regression and ridge robust estimators in Section. A numerical example is given in section 4, and the results of simulation studies are presented in Section 5. In order to compare the proposed with the existing methods. We shall illustrate that the new estimation method is better than the existing methods when multicollinearity and outliers occur simultaneously. Finally, in section 6 conclusions are presented.. Method. Ridge Regression When the least-squares method is applied to non-orthogonal data, we obtain weak estimates of the regression coefficients (Hoerl and Kennard, 976). The assumptions of the least squares method that ˆβ is an unbiased estimator of β. The Gauss-Markov property assures that the least-squares estimator has a minimum variance in the class of unbiased linear estimators in the presence of multicollinearity (Marquardt and Snee, 975). One way to reduce this problem is to descend the requirement that the estimator of β be unbiased. Assume that a biased estimator of β is found to be say ˆβ, that has smaller variance than the unbiased estimator. The mean square error of β is defined as ˆ ˆ ˆ MSE( β)=var( β) + E( β β ) () Assuming a small amount of bias in ˆβ, the variance of ˆβ can be made small such that the Mean Square Error (MSE) of ˆβ is less than the variance of the unbiased estimator β. A number of procedures have been employed in obtaining biased estimates of regression coefficients. One of these procedures is ridge regression, primarily proposed by (Hoerl and Kennard, 970). More specifically, this estimator is defined as the solution to βˆ = ( X'X+ KI) X'y where K 0 is a constant selected by the analyst. The constant K is usually referred to the biasing parameter, it is obvious that if K = 0 then the ridge estimator is the least squares estimator. Hence, on the other hand, when K > 0 the bias increases. However, the variance ridge decreases as K increases. In this ridge regression, we would like to select a value of K such that the decrease in the variance term is greater than the increase in the squared bias. If this can be done, the mean square error of the ridge estimator will ridge be less than the variance of the least square estimator ˆβ. (Hoerl and Kennard, 976), proved that there exists a non- zero value of K for which the Mean Square Error (MSE) is less than the variance of the least squares estimator ˆβ. We can select a reasonably small value of K at which the ridge estimates are steady. Mostly, this will produce a set of estimates with smaller MSE than ridge the least-squares estimates. (Hoerl and Kennard, 976), suggested that an appropriate choice of K is where ˆβ and ˆ σ are found by the least squares solution.. Robust Regression Estimators (RRE) Ordinary least squares (OLS) regression is very useful, easily interpretable statistical method. However, it is not complete. When we apply an OLS regression, one needs to be aware of its sensitivity to outliers. By "sensitivity to outliers", we mean that an OLS regression model can be highly affected by a few records in the dataset and can then yield results that do not reflect the relationship between the outcome variable and the explanatory variables seen in the rest of the records. Robust regression gives rise to an alternative to OLS regression that is less sensitive to outliers and still defines a linear relationship between the outcome and the explanatory variables. 9

3 Modern Applied Science Vol. 9, No. ; 05 There are several different classifications of robust estimates that exist to handle these violations. One important estimator is called the least median squares estimator (Rousseeuw and Leroy, 987), which has the advantage of minimizing the influence of the residuals. According to (Venables and Ripley, 999), this algorithm minimizes the median of ordered squares of residuals in order to get the regression coefficients β and can be written as equation. Least Median of Squares = Min median i LMS y-xβ () (Martin, 00) describes the median squared residuals lack a smooth squared residual function and the asymptotic rate is n to converge efficiently under normality it means that LMS converges less smoothly and weaker than the LTS algorithm. And also, takes a long time to converge. Least Median of Squares (LMS) estimator is one of the true high breakdown point estimators that reached the above mentioned upper boundary of the breakdown point. (Rousseeuw, 984)... Robust Ridge Regression Estimators (RRRE) (Pfaffenberger and Dielman 985, 990), proposed robust ridge regression by extending the development of their technique by performing Monte Carlo simulation studies to compare various approaches. The proposed method in this paper combines the LMS estimator with the ridge regression estimator which is referred to as the RLMS estimator. So, RLMS robust ridge estimators will be resistant to multicollinearity problems and less affected by outliers. The RLMS estimator can be written as equation. - βˆ RLMS = ( X'X + KLMSI ) X'Y () The value of K is determined from data using equation 4. ps K = LMS LMS βˆ ' ˆ LMSβLMS and n εilms S = i= LMS n p where p is the number of independent variables, n is the number of observations in the data, is the LMS estimates of β and ε is the residuals from LMS method. LMS. Results. Numerical Result A real data set of (Hald, 95) is considered to assess the effectiveness of the proposed robust ridge regression method. This data set consists of 4 variables and observations with outliers. The response variable (y) is the heat evolved from a particular mixture of cement, and the covariates are tricalcium aluminate (x), tricalcium silicate (x), tetracalcium alumina ferrite (x), and dicalcium silicate (x4). (Table ) contains the parameters estimates, the standard error of the parameters estimates and Variance Inflation Factor (VIF) analysis of the Hald Data. The VIF for each of these predictors are all extremely high, indicating troublesome multicollinearity with the presence of outliers in the data, the use of robust method provides more stable parameter estimates. With this aim, initial robust regression estimates was first calculated to obtain robust ridge estimates in the presence of both multicollinearity and outliers; these estimates are given in (Table ) (4) Table. Estimated parameter and SE of,, and 4 for the different methods with the variance inflation factor. Coef. βˆ 4 Estimate VIF OLS RIDGE RLMS RLAV parameter s.e. parameter s.e. parameter s.e. parameter s.e

4 Modern Applied Science Vol. 9, No. ; 05 The standard error of the proposed RLMS method is lower than the existing RLAV in the presence of outliers. (Table ) shows that the data have a high VIF for all variables with three outliers. Also, the standard errors of the parameters estimates for RLMS are less than all except for the ridge regression whose difference is smaller.. Simulation Study We carry out a simulation study to compare the performance of the different methods OLS, RR and RLAV with the proposed estimator RLMS. The simulation is designed to allow both multicollinearity and non-normality to exist simultaneously. The non-normal distributions are used to generate outliers. Suppose, we have the following linear regression model. (Siti Meriam et al., 0) y = β + β x +β x +β x +e i o i i i i, where i=,, (5) The parameter values β o, β, β and β are set equal to one. The explanatory variables x i, x i and x i are generated using equation (6) x ij ( ρ ) = + i=,,, n, j =,, (6) ρz ij where, z are independent standard normal random numbers generated by the normal distribution. ij The explanatory variables values were generated for a given sample size n. The sample sizes used were 50 and 00. The value of ρ representing the correlation between any two explanatory variables, and its values were chosen as: 0.0, 0.5 and The percentage of outliers present in this data set is 0%. The number of replications used is 500. The statistics computed are the bias, root of mean squared error (RMSE), standard error (SE), and 6 pairwise MSE ratios of the estimators. The bias and MSE are given as: Bias = β β where β i= 500 MSE = ( ˆ ) 500 β -β i= i i, therefore, the RMSE is given by i ˆ / [ ( )] The VIF for the simulated data are shown in (Table ) i k βi = i k k = 500 MSE β j where j= 0,,,, and the mean squared error ( MSE ) is Table. The VIF for the simulated data ρ = 0.99 Var X X X VIF (N=50) VIF (N=00) Here the maximum VIF is when the correlation between independent variables was very high with different sample size. So it is clear that the multicollinearity problem exists. The measure of convergence was computed as the number of times estimator or to the true parameter while the value in (Table and Table 4) show the summary statistics such as bias, root mean squared error (RMSE) and standard error (SE) of the estimators of the normal distributions for sample size 50 and 00 with 0% and 0% of outliers and different value of ρ. On the other hand, when we apply these methods to the simulated data with different sample size in the presence of different percentage of outliers and different levels of multicollineaity, we obtained that the standard errors for proposed method are less than the standard errors of the aforementioned existing methods. 94

5 Modern Applied Science Vol. 9, No. ; 05 (Table 5) shows the efficiency of the estimators by looking at the MSE ratios of the estimators written as follows. RMSE (proposed method) MSE ratios = (7) RMSE (existing method) This Table supported the result obtained from (Table and Table 4). These ratios show the efficiency of RLMS relative to other estimators then MSE ratios less than denotes that the estimator is more efficient, however, values more than denote that the other estimator is more efficient. From (Table and Table 4) we can see that the RMSE of the OLS is relatively smaller than the other estimators when the errors are normally distributed that is, without multicollinearity and no outliers. As expected, the OLS gives the best result in the ideal situation. Also, the result in (Table 5) is in favor of LS. However, we see in (Table 5) that the ratio of MSE of RLMS to OLS is greater than.00 denoting that the OLS is more efficient than the RLMS when no multicollinearity and no outliers. On the other hand, we can see from the (Table and Table 4) that the RMSE of the RIDGE is relatively smaller than the RLAV. Also, the values of MSE ratios for the estimators ridge regression are less than indicates that the estimator is less efficient than RLAV and RLMS when the errors are normally distributed without outliers and no multicollinearity. While, for non-normal error distribution and when correlation and outliers are present in the data, RLMS is better than the OLS, RIDGE and RLAV. RLAV is almost as good as RIDGE and OLS. Consequently, we can see that the RMSE of the RLMS is relatively smaller than the other estimators when the errors are normally distributed in the presence of multicollinearity and outliers. It obviously shows that RLMS is more efficient than RLAV and RIDGE but certainly much more efficient than OLS when multicollinearity and outliers are present. Table. Bias, RMSE and SE of, and with error normal (0,) distribution of the sample size n=50 correlation 0.0, 0.5, 0.99 outliers 0% and 0% 0% outliers Value of ρ 0.0 ˆ 0.5 ˆ Bias RMSE S.e β Bias RMSE S.e Bias RMSE S.e % outiers Bias RMSE S.e β Bias RMSE S.e Bias RMSE S.e % outiers Bias RMSE S.e β Bias ˆ 95

6 Modern Applied Science Vol. 9, No. ; RMSE S.e Bias RMSE S.e Table 4. Bias, RMSE and SE of, and with error normal (0,) distribution of the sample size n=00 correlation 0.0, 0.5, 0.99 outliers 0% and 0% 0% outliers Value of Bias ρ ˆ RMSE β S.e Bias RMSE S.e Bias RMSE S.e % outiers Bias RMSE S.e Bias RMSE S.e Bias RMSE S.e % outiers Bias RMSE S.e Bias RMSE S.e Bias RMSE S.e Table 5. MSE ratios of 6 pairwise estimators of,, with errors normal (0,) distribution and 0% and 0% of outliers ˆ β Estimator vs Estimator vs Estimator Values of ρ RLMS OLS RIDGE RLAV

7 Modern Applied Science Vol. 9, No. ; RLAV RIDGE OLS RIDGE OLS The simulation results for larger samples, that is for n=00 are consistent with the results of the smaller sized samples. The results also show that the estimator for larger samples are more efficient than those of smaller samples since the values of RMSE are smaller. 4. Discussion Multicollinearity data sets with outliers are very common in real life situation. In order to remedy both problems, robust biased estimation methods are applied. The best model is chosen by looking at the RMSE value. The simulation study in section 5 provides RMSE values, bias and standard errors of the OLS, ridge, RLAV and RLMS estimators. It can be observed that the RMSE obtained from RLMS is the minimum. Thus, RLMS is the best method. When there is no multicollinearity and outliers, ridge regression has the least RMSE value, thus the best method. However, when there is multicollinearity and outliers in the data, then RLMS has the least RMSE value, thus it is considered the best method. We also use a real-life data set to study the effect of multicollinearity and outliers. The results obtained from our proposed method of RLMS are better than the OLS, Ridge and RLAV in terms of their RMSE values. Consequently, in this study, it is shown that the results obtained from the numerical data set and a simulation study with both multicollinearity and outliers, RLMS gives better results followed by RLAV method. This result is true for ρ =0.0, 0.5 and 0.99 with the sample size of n=50 and 00. Acknowledgements This work was supported by the Universiti Teknologi Malaysia, under the Research University Grant (QJ H68) and Ministry of Higher Education (MOHE) of Malayisa. References Askin, R. G., & Montgomery, D. C. (980). Augmented robust estimators technometrics,, -4. Insightful Corp Presentation, New York Desire, L. M., & Leonard, K. (986). A robust method for outlier and model error detection in regression and calibration. Analytica chimica Acta., 87, Haled. (95). Statistical Theory with Engineering Application. Wiley, New York. p Hampel, F. R. (975). Beyond location parameters: Robust concepts and methods proceedings of the 40 th session of the ISI., 46(), Hoerl, A. E., & Baldwin, K. F. (976). Ridge Regression: Some Simulations. Communications in Statistics, (4), 04-. Hoerl, A. E., & Kennard, R. W. (970a). Ridge Regression: Biased Estimation for Non-orthogonal Problems. Technometrics, (55). Hoerl, A. E., & Kennard, R. W. (970b). Ridge Regression: Applications for Non- orthogonal Problems. Technometrics, (69). Hoerl, A. E., & Kennard, R. W. (976). Ridge regression: Iterative estimation of the biasing parameter. Commun. Statist. Theory Meth., 5(), Mal, C. Z., & Dul, Y. L. (04). Generalized shrunken type-gm estimator and its application. International Conference on Applied Sciences (ICAS0). Marquardt, D. W., & Snee, R. D. (975). Ridge regression in practice. Amer. Statist., (9), -9. Martin, R. D. (00). Robust Statistics with the S-Plus Robust Library and Financial Applications. 97

IDENTIFYING MULTIPLE OUTLIERS IN LINEAR REGRESSION : ROBUST FIT AND CLUSTERING APPROACH

SESSION X : THEORY OF DEFORMATION ANALYSIS II IDENTIFYING MULTIPLE OUTLIERS IN LINEAR REGRESSION : ROBUST FIT AND CLUSTERING APPROACH Robiah Adnan 2 Halim Setan 3 Mohd Nor Mohamad Faculty of Science, Universiti

MULTICOLLINEARITY DIAGNOSTIC MEASURES BASED ON MINIMUM COVARIANCE DETERMINATION APPROACH

Professor Habshah MIDI, PhD Department of Mathematics, Faculty of Science / Laboratory of Computational Statistics and Operations Research, Institute for Mathematical Research University Putra, Malaysia

PAijpam.eu M ESTIMATION, S ESTIMATION, AND MM ESTIMATION IN ROBUST REGRESSION

International Journal of Pure and Applied Mathematics Volume 91 No. 3 2014, 349-360 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu doi: http://dx.doi.org/10.12732/ijpam.v91i3.7

Least Absolute Value vs. Least Squares Estimation and Inference Procedures in Regression Models with Asymmetric Error Distributions

Journal of Modern Applied Statistical Methods Volume 8 Issue 1 Article 13 5-1-2009 Least Absolute Value vs. Least Squares Estimation and Inference Procedures in Regression Models with Asymmetric Error

Ridge Estimator in Logistic Regression under Stochastic Linear Restrictions

British Journal of Mathematics & Computer Science 15(3): 1-14, 2016, Article no.bjmcs.24585 ISSN: 2231-0851 SCIENCEDOMAIN international www.sciencedomain.org Ridge Estimator in Logistic Regression under

Least Squares Estimation-Finite-Sample Properties

Least Squares Estimation-Finite-Sample Properties Ping Yu School of Economics and Finance The University of Hong Kong Ping Yu (HKU) Finite-Sample 1 / 29 Terminology and Assumptions 1 Terminology and Assumptions

MOVING WINDOW REGRESSION (MWR) IN MASS APPRAISAL FOR PROPERTY RATING. Universiti Putra Malaysia UPM Serdang, Malaysia

MOVING WINDOW REGRESSION (MWR IN MASS APPRAISAL FOR PROPERTY RATING 1 Taher Buyong, Suriatini Ismail, Ibrahim Sipan, Mohamad Ghazali Hashim and 1 Mohammad Firdaus Azhar 1 Institute of Advanced Technology

Homoskedasticity. Var (u X) = σ 2. (23)

Homoskedasticity How big is the difference between the OLS estimator and the true parameter? To answer this question, we make an additional assumption called homoskedasticity: Var (u X) = σ 2. (23) This

Study of the Relationship between Dependent and Independent Variable Groups by Using Canonical Correlation Analysis with Application

Modern Applied Science; Vol. 9, No. 8; 2015 ISSN 1913-1844 E-ISSN 1913-1852 Published by Canadian Center of Science and Education Study of the Relationship between Dependent and Independent Variable Groups

ISyE 691 Data mining and analytics

ISyE 691 Data mining and analytics Regression Instructor: Prof. Kaibo Liu Department of Industrial and Systems Engineering UW-Madison Email: kliu8@wisc.edu Office: Room 3017 (Mechanical Engineering Building)

Estimating σ 2. We can do simple prediction of Y and estimation of the mean of Y at any value of X.

Estimating σ 2 We can do simple prediction of Y and estimation of the mean of Y at any value of X. To perform inferences about our regression line, we must estimate σ 2, the variance of the error term.

Multiple Linear Regression

Multiple Linear Regression Asymptotics Asymptotics Multiple Linear Regression: Assumptions Assumption MLR. (Linearity in parameters) Assumption MLR. (Random Sampling from the population) We have a random

Variable Selection and Model Building

LINEAR REGRESSION ANALYSIS MODULE XIII Lecture - 37 Variable Selection and Model Building Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur The complete regression

HETEROSKEDASTICITY, TEMPORAL AND SPATIAL CORRELATION MATTER

ACTA UNIVERSITATIS AGRICULTURAE ET SILVICULTURAE MENDELIANAE BRUNENSIS Volume LXI 239 Number 7, 2013 http://dx.doi.org/10.11118/actaun201361072151 HETEROSKEDASTICITY, TEMPORAL AND SPATIAL CORRELATION MATTER

A Short Introduction to Curve Fitting and Regression by Brad Morantz

A Short Introduction to Curve Fitting and Regression by Brad Morantz bradscientist@machine-cognition.com Overview What can regression do for me? Example Model building Error Metrics OLS Regression Robust

A Modified M-estimator for the Detection of Outliers

A Modified M-estimator for the Detection of Outliers Asad Ali Department of Statistics, University of Peshawar NWFP, Pakistan Email: asad_yousafzay@yahoo.com Muhammad F. Qadir Department of Statistics,

Multiple Regression Analysis: Heteroskedasticity

Multiple Regression Analysis: Heteroskedasticity y = β 0 + β 1 x 1 + β x +... β k x k + u Read chapter 8. EE45 -Chaiyuth Punyasavatsut 1 topics 8.1 Heteroskedasticity and OLS 8. Robust estimation 8.3 Testing

Multiple Regression Analysis. Part III. Multiple Regression Analysis

Part III Multiple Regression Analysis As of Sep 26, 2017 1 Multiple Regression Analysis Estimation Matrix form Goodness-of-Fit R-square Adjusted R-square Expected values of the OLS estimators Irrelevant

MULTICOLLINEARITY AND VARIANCE INFLATION FACTORS. F. Chiaromonte 1

MULTICOLLINEARITY AND VARIANCE INFLATION FACTORS F. Chiaromonte 1 Pool of available predictors/terms from them in the data set. Related to model selection, are the questions: What is the relative importance

Multiple Regression Methods

Chapter 1: Multiple Regression Methods Hildebrand, Ott and Gray Basic Statistical Ideas for Managers Second Edition 1 Learning Objectives for Ch. 1 The Multiple Linear Regression Model How to interpret

Statistical Inference with Regression Analysis

Introductory Applied Econometrics EEP/IAS 118 Spring 2015 Steven Buck Lecture #13 Statistical Inference with Regression Analysis Next we turn to calculating confidence intervals and hypothesis testing

Comparing the performance of modified F t statistic with ANOVA and Kruskal Wallis test

Appl. Math. Inf. Sci. 7, No. 2L, 403-408 (2013) 403 Applied Mathematics & Information Sciences An International ournal http://dx.doi.org/10.12785/amis/072l04 Comparing the performance of modified F t statistic

Robustness of location estimators under t- distributions: a literature review

IOP Conference Series: Earth and Environmental Science PAPER OPEN ACCESS Robustness of location estimators under t- distributions: a literature review o cite this article: C Sumarni et al 07 IOP Conf.

Eric Shou Stat 598B / CSE 598D METHODS FOR MICRODATA PROTECTION

Eric Shou Stat 598B / CSE 598D METHODS FOR MICRODATA PROTECTION INTRODUCTION Statistical disclosure control part of preparations for disseminating microdata. Data perturbation techniques: Methods assuring

CHAPTER 6: SPECIFICATION VARIABLES

Recall, we had the following six assumptions required for the Gauss-Markov Theorem: 1. The regression model is linear, correctly specified, and has an additive error term. 2. The error term has a zero

Small sample corrections for LTS and MCD

Metrika (2002) 55: 111 123 > Springer-Verlag 2002 Small sample corrections for LTS and MCD G. Pison, S. Van Aelst*, and G. Willems Department of Mathematics and Computer Science, Universitaire Instelling

Econometrics Summary Algebraic and Statistical Preliminaries

Econometrics Summary Algebraic and Statistical Preliminaries Elasticity: The point elasticity of Y with respect to L is given by α = ( Y/ L)/(Y/L). The arc elasticity is given by ( Y/ L)/(Y/L), when L

The Precise Effect of Multicollinearity on Classification Prediction

Multicollinearity and Classification Prediction The Precise Effect of Multicollinearity on Classification Prediction Mary G. Lieberman John D. Morris Florida Atlantic University The results of Morris and

A Brief Overview of Robust Statistics

A Brief Overview of Robust Statistics Olfa Nasraoui Department of Computer Engineering & Computer Science University of Louisville, olfa.nasraoui_at_louisville.edu Robust Statistical Estimators Robust

The multiple regression model; Indicator variables as regressors

The multiple regression model; Indicator variables as regressors Ragnar Nymoen University of Oslo 28 February 2013 1 / 21 This lecture (#12): Based on the econometric model specification from Lecture 9

Regression #3: Properties of OLS Estimator

Regression #3: Properties of OLS Estimator Econ 671 Purdue University Justin L. Tobias (Purdue) Regression #3 1 / 20 Introduction In this lecture, we establish some desirable properties associated with

Multilevel modeling and panel data analysis in educational research (Case study: National examination data senior high school in West Java)

Multilevel modeling and panel data analysis in educational research (Case study: National examination data senior high school in West Java) Pepi Zulvia, Anang Kurnia, and Agus M. Soleh Citation: AIP Conference

Response Surface Methodology

Response Surface Methodology Process and Product Optimization Using Designed Experiments Second Edition RAYMOND H. MYERS Virginia Polytechnic Institute and State University DOUGLAS C. MONTGOMERY Arizona

Some Monte Carlo Evidence for Adaptive Estimation of Unit-Time Varying Heteroscedastic Panel Data Models

Some Monte Carlo Evidence for Adaptive Estimation of Unit-Time Varying Heteroscedastic Panel Data Models G. R. Pasha Department of Statistics, Bahauddin Zakariya University Multan, Pakistan E-mail: drpasha@bzu.edu.pk

A NOTE ON ROBUST ESTIMATION IN LOGISTIC REGRESSION MODEL

Discussiones Mathematicae Probability and Statistics 36 206 43 5 doi:0.75/dmps.80 A NOTE ON ROBUST ESTIMATION IN LOGISTIC REGRESSION MODEL Tadeusz Bednarski Wroclaw University e-mail: t.bednarski@prawo.uni.wroc.pl

The Simple Regression Model. Part II. The Simple Regression Model

Part II The Simple Regression Model As of Sep 22, 2015 Definition 1 The Simple Regression Model Definition Estimation of the model, OLS OLS Statistics Algebraic properties Goodness-of-Fit, the R-square

On Monitoring Shift in the Mean Processes with. Vector Autoregressive Residual Control Charts of. Individual Observation

Applied Mathematical Sciences, Vol. 8, 14, no. 7, 3491-3499 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.12988/ams.14.44298 On Monitoring Shift in the Mean Processes with Vector Autoregressive Residual

Multiple Regression Analysis

Multiple Regression Analysis y = 0 + 1 x 1 + x +... k x k + u 6. Heteroskedasticity What is Heteroskedasticity?! Recall the assumption of homoskedasticity implied that conditional on the explanatory variables,

Linear Regression Models. Based on Chapter 3 of Hastie, Tibshirani and Friedman

Linear Regression Models Based on Chapter 3 of Hastie, ibshirani and Friedman Linear Regression Models Here the X s might be: p f ( X = " + " 0 j= 1 X j Raw predictor variables (continuous or coded-categorical

Using Estimating Equations for Spatially Correlated A

Using Estimating Equations for Spatially Correlated Areal Data December 8, 2009 Introduction GEEs Spatial Estimating Equations Implementation Simulation Conclusion Typical Problem Assess the relationship

CHAPTER 4: Forecasting by Regression

CHAPTER 4: Forecasting by Regression Prof. Alan Wan 1 / 57 Table of contents 1. Revision of Linear Regression 3.1 First-order Autocorrelation and the Durbin-Watson Test 3.2 Correction for Autocorrelation

Multivariate Regression Analysis

Matrices and vectors The model from the sample is: Y = Xβ +u with n individuals, l response variable, k regressors Y is a n 1 vector or a n l matrix with the notation Y T = (y 1,y 2,...,y n ) 1 x 11 x

ECON3150/4150 Spring 2016

ECON3150/4150 Spring 2016 Lecture 6 Multiple regression model Siv-Elisabeth Skjelbred University of Oslo February 5th Last updated: February 3, 2016 1 / 49 Outline Multiple linear regression model and

Linear Models 1. Isfahan University of Technology Fall Semester, 2014

Linear Models 1 Isfahan University of Technology Fall Semester, 2014 References: [1] G. A. F., Seber and A. J. Lee (2003). Linear Regression Analysis (2nd ed.). Hoboken, NJ: Wiley. [2] A. C. Rencher and

Linear models. Linear models are computationally convenient and remain widely used in. applied econometric research

Linear models Linear models are computationally convenient and remain widely used in applied econometric research Our main focus in these lectures will be on single equation linear models of the form y

Prediction Intervals in the Presence of Outliers

Prediction Intervals in the Presence of Outliers David J. Olive Southern Illinois University July 21, 2003 Abstract This paper presents a simple procedure for computing prediction intervals when the data

Useful Numerical Statistics of Some Response Surface Methodology Designs

Journal of Mathematics Research; Vol. 8, No. 4; August 20 ISSN 19-9795 E-ISSN 19-9809 Published by Canadian Center of Science and Education Useful Numerical Statistics of Some Response Surface Methodology

Measurement Error and Linear Regression of Astronomical Data. Brandon Kelly Penn State Summer School in Astrostatistics, June 2007

Measurement Error and Linear Regression of Astronomical Data Brandon Kelly Penn State Summer School in Astrostatistics, June 2007 Classical Regression Model Collect n data points, denote i th pair as (η

Matematické Metody v Ekonometrii 7.

Matematické Metody v Ekonometrii 7. Multicollinearity Blanka Šedivá KMA zimní semestr 2016/2017 Blanka Šedivá (KMA) Matematické Metody v Ekonometrii 7. zimní semestr 2016/2017 1 / 15 One of the assumptions

PBAF 528 Week 8. B. Regression Residuals These properties have implications for the residuals of the regression.

PBAF 528 Week 8 What are some problems with our model? Regression models are used to represent relationships between a dependent variable and one or more predictors. In order to make inference from the

Prevention of the quality of cement from chemical and mineralogical parameters of the clinker

International Journal of Innovation and Applied Studies ISSN 2028-9324 Vol. 10 No. 4 Mar. 2015, pp. 1291-1298 2015 Innovative Space of Scientific Research Journals http://www.ijias.issr-journals.org/ Prevention

Statistics GIDP Ph.D. Qualifying Exam Methodology

Statistics GIDP Ph.D. Qualifying Exam Methodology January 9, 2018, 9:00am 1:00pm Instructions: Put your ID (not your name) on each sheet. Complete exactly 5 of 6 problems; turn in only those sheets you

INFLUENCE OF USING ALTERNATIVE MEANS ON TYPE-I ERROR RATE IN THE COMPARISON OF INDEPENDENT GROUPS ABSTRACT

Mirtagioğlu et al., The Journal of Animal & Plant Sciences, 4(): 04, Page: J. 344-349 Anim. Plant Sci. 4():04 ISSN: 08-708 INFLUENCE OF USING ALTERNATIVE MEANS ON TYPE-I ERROR RATE IN THE COMPARISON OF

Linear Regression. Junhui Qian. October 27, 2014

Linear Regression Junhui Qian October 27, 2014 Outline The Model Estimation Ordinary Least Square Method of Moments Maximum Likelihood Estimation Properties of OLS Estimator Unbiasedness Consistency Efficiency

Chapter 26 Multiple Regression, Logistic Regression, and Indicator Variables

Chapter 26 Multiple Regression, Logistic Regression, and Indicator Variables 26.1 S 4 /IEE Application Examples: Multiple Regression An S 4 /IEE project was created to improve the 30,000-footlevel metric

Multiple Regression. Midterm results: AVG = 26.5 (88%) A = 27+ B = C =

Economics 130 Lecture 6 Midterm Review Next Steps for the Class Multiple Regression Review & Issues Model Specification Issues Launching the Projects!!!!! Midterm results: AVG = 26.5 (88%) A = 27+ B =

Applied Regression Analysis

Applied Regression Analysis Chapter 3 Multiple Linear Regression Hongcheng Li April, 6, 2013 Recall simple linear regression 1 Recall simple linear regression 2 Parameter Estimation 3 Interpretations of

Lecture Notes Autumn 2010 Dr. Getinet Haile, University of Mannheim 1. Introduction Introduction & CLRM, Autumn Term 2010 1 What is econometrics? Econometrics = economic statistics economic theory mathematics

Lecture 2 Linear Regression: A Model for the Mean. Sharyn O Halloran

Lecture 2 Linear Regression: A Model for the Mean Sharyn O Halloran Closer Look at: Linear Regression Model Least squares procedure Inferential tools Confidence and Prediction Intervals Assumptions Robustness

A New Class of Generalized Exponential Ratio and Product Type Estimators for Population Mean using Variance of an Auxiliary Variable

ISSN 1684-8403 Journal of Statistics Volume 21, 2014. pp. 206-214 A New Class of Generalized Exponential Ratio and Product Type Abstract Hina Khan 1 and Ahmad Faisal Siddiqi 2 In this paper, a general

Lecture 13: Simple Linear Regression in Matrix Format

See updates and corrections at http://www.stat.cmu.edu/~cshalizi/mreg/ Lecture 13: Simple Linear Regression in Matrix Format 36-401, Section B, Fall 2015 13 October 2015 Contents 1 Least Squares in Matrix

OSU Economics 444: Elementary Econometrics. Ch.10 Heteroskedasticity

OSU Economics 444: Elementary Econometrics Ch.0 Heteroskedasticity (Pure) heteroskedasticity is caused by the error term of a correctly speciþed equation: Var(² i )=σ 2 i, i =, 2,,n, i.e., the variance

Multiple Regression. Peerapat Wongchaiwat, Ph.D.

Peerapat Wongchaiwat, Ph.D. wongchaiwat@hotmail.com The Multiple Regression Model Examine the linear relationship between 1 dependent (Y) & 2 or more independent variables (X i ) Multiple Regression Model

Estimation of the Parameters of Bivariate Normal Distribution with Equal Coefficient of Variation Using Concomitants of Order Statistics

International Journal of Statistics and Probability; Vol., No. 3; 013 ISSN 197-703 E-ISSN 197-7040 Published by Canadian Center of Science and Education Estimation of the Parameters of Bivariate Normal

Bias-Correction in Vector Autoregressive Models: A Simulation Study

Econometrics 2014, 2, 45-71; doi:10.3390/econometrics2010045 OPEN ACCESS econometrics ISSN 2225-1146 www.mdpi.com/journal/econometrics Article Bias-Correction in Vector Autoregressive Models: A Simulation

Sizing Mixture (RSM) Designs for Adequate Precision via Fraction of Design Space (FDS)

for Adequate Precision via Fraction of Design Space (FDS) Pat Whitcomb Stat-Ease, Inc. 61.746.036 036 fax 61.746.056 pat@statease.com Gary W. Oehlert School of Statistics University it of Minnesota 61.65.1557

Model Mis-specification

Model Mis-specification Carlo Favero Favero () Model Mis-specification 1 / 28 Model Mis-specification Each specification can be interpreted of the result of a reduction process, what happens if the reduction

Combining data from two independent surveys: model-assisted approach

Combining data from two independent surveys: model-assisted approach Jae Kwang Kim 1 Iowa State University January 20, 2012 1 Joint work with J.N.K. Rao, Carleton University Reference Kim, J.K. and Rao,

Multicollinearity occurs when two or more predictors in the model are correlated and provide redundant information about the response.

Multicollinearity Read Section 7.5 in textbook. Multicollinearity occurs when two or more predictors in the model are correlated and provide redundant information about the response. Example of multicollinear

Basic econometrics. Tutorial 3. Dipl.Kfm. Johannes Metzler

Basic econometrics Tutorial 3 Dipl.Kfm. Introduction Some of you were asking about material to revise/prepare econometrics fundamentals. First of all, be aware that I will not be too technical, only as

Robust statistics. Michael Love 7/10/2016

Robust statistics Michael Love 7/10/2016 Robust topics Median MAD Spearman Wilcoxon rank test Weighted least squares Cook's distance M-estimators Robust topics Median => middle MAD => spread Spearman =>

ECON2228 Notes 7. Christopher F Baum. Boston College Economics. cfb (BC Econ) ECON2228 Notes / 41

ECON2228 Notes 7 Christopher F Baum Boston College Economics 2014 2015 cfb (BC Econ) ECON2228 Notes 6 2014 2015 1 / 41 Chapter 8: Heteroskedasticity In laying out the standard regression model, we made

Small Domain Estimation for a Brazilian Service Sector Survey

Proceedings 59th ISI World Statistics Congress, 5-30 August 013, Hong Kong (Session CPS003) p.334 Small Domain Estimation for a Brazilian Service Sector Survey André Neves 1, Denise Silva and Solange Correa

ECON3150/4150 Spring 2015

ECON3150/4150 Spring 2015 Lecture 3&4 - The linear regression model Siv-Elisabeth Skjelbred University of Oslo January 29, 2015 1 / 67 Chapter 4 in S&W Section 17.1 in S&W (extended OLS assumptions) 2

3 Multiple Linear Regression

3 Multiple Linear Regression 3.1 The Model Essentially, all models are wrong, but some are useful. Quote by George E.P. Box. Models are supposed to be exact descriptions of the population, but that is

x 21 x 22 x 23 f X 1 X 2 X 3 ε

Chapter 2 Estimation 2.1 Example Let s start with an example. Suppose that Y is the fuel consumption of a particular model of car in m.p.g. Suppose that the predictors are 1. X 1 the weight of the car

The regression model with one stochastic regressor.

The regression model with one stochastic regressor. 3150/4150 Lecture 6 Ragnar Nymoen 30 January 2012 We are now on Lecture topic 4 The main goal in this lecture is to extend the results of the regression

PENALIZING YOUR MODELS AN OVERVIEW OF THE GENERALIZED REGRESSION PLATFORM Michael Crotty & Clay Barker Research Statisticians JMP Division, SAS Institute Copyr i g ht 2012, SAS Ins titut e Inc. All rights

Iterative Stepwise Selection and Threshold for Learning Causes in Time Series

Iterative Stepwise Selection and Threshold for Learning Causes in Time Series Jianxin Yin Shaopeng Wang Wanlu Deng Ya Hu Zhi Geng School of Mathematical Sciences Peking University Beijing 100871, China

The general linear regression with k explanatory variables is just an extension of the simple regression as follows

3. Multiple Regression Analysis The general linear regression with k explanatory variables is just an extension of the simple regression as follows (1) y i = β 0 + β 1 x i1 + + β k x ik + u i. Because

IEOR 165 Lecture 7 1 Bias-Variance Tradeoff

IEOR 165 Lecture 7 Bias-Variance Tradeoff 1 Bias-Variance Tradeoff Consider the case of parametric regression with β R, and suppose we would like to analyze the error of the estimate ˆβ in comparison to

Ma 3/103: Lecture 24 Linear Regression I: Estimation

Ma 3/103: Lecture 24 Linear Regression I: Estimation March 3, 2017 KC Border Linear Regression I March 3, 2017 1 / 32 Regression analysis Regression analysis Estimate and test E(Y X) = f (X). f is the

Multilevel Statistical Models: 3 rd edition, 2003 Contents

Multilevel Statistical Models: 3 rd edition, 2003 Contents Preface Acknowledgements Notation Two and three level models. A general classification notation and diagram Glossary Chapter 1 An introduction

Analysis of Cross-Sectional Data

Analysis of Cross-Sectional Data Kevin Sheppard http://www.kevinsheppard.com Oxford MFE This version: November 8, 2017 November 13 14, 2017 Outline Econometric models Specification that can be analyzed

Iris Wang.

Chapter 10: Multicollinearity Iris Wang iris.wang@kau.se Econometric problems Multicollinearity What does it mean? A high degree of correlation amongst the explanatory variables What are its consequences?

Response Surface Methodology:

Response Surface Methodology: Process and Product Optimization Using Designed Experiments RAYMOND H. MYERS Virginia Polytechnic Institute and State University DOUGLAS C. MONTGOMERY Arizona State University

Applications of Mathematics

Applications of Mathematics Tomáš Cipra Exponential smoothing for irregular data Applications of Mathematics Vol. 5 2006) No. 6 597--604 Persistent URL: http://dml.cz/dmlcz/34655 Terms of use: Institute

The Multiple Regression Model

Multiple Regression The Multiple Regression Model Idea: Examine the linear relationship between 1 dependent (Y) & or more independent variables (X i ) Multiple Regression Model with k Independent Variables:

ASA Section on Survey Research Methods

REGRESSION-BASED STATISTICAL MATCHING: RECENT DEVELOPMENTS Chris Moriarity, Fritz Scheuren Chris Moriarity, U.S. Government Accountability Office, 411 G Street NW, Washington, DC 20548 KEY WORDS: data

Economic modelling and forecasting

Economic modelling and forecasting 2-6 February 2015 Bank of England he generalised method of moments Ole Rummel Adviser, CCBS at the Bank of England ole.rummel@bankofengland.co.uk Outline Classical estimation

FinQuiz Notes

Reading 10 Multiple Regression and Issues in Regression Analysis 2. MULTIPLE LINEAR REGRESSION Multiple linear regression is a method used to model the linear relationship between a dependent variable

Mean squared error matrix comparison of least aquares and Stein-rule estimators for regression coefficients under non-normal disturbances

METRON - International Journal of Statistics 2008, vol. LXVI, n. 3, pp. 285-298 SHALABH HELGE TOUTENBURG CHRISTIAN HEUMANN Mean squared error matrix comparison of least aquares and Stein-rule estimators

Multiple Linear Regression

Andrew Lonardelli December 20, 2013 Multiple Linear Regression 1 Table Of Contents Introduction: p.3 Multiple Linear Regression Model: p.3 Least Squares Estimation of the Parameters: p.4-5 The matrix approach

INFERENCE FOR MULTIPLE LINEAR REGRESSION MODEL WITH EXTENDED SKEW NORMAL ERRORS

Pak. J. Statist. 2016 Vol. 32(2), 81-96 INFERENCE FOR MULTIPLE LINEAR REGRESSION MODEL WITH EXTENDED SKEW NORMAL ERRORS A.A. Alhamide 1, K. Ibrahim 1 M.T. Alodat 2 1 Statistics Program, School of Mathematical

Multiple Regression Analysis

1 OUTLINE Basic Concept: Multiple Regression MULTICOLLINEARITY AUTOCORRELATION HETEROSCEDASTICITY REASEARCH IN FINANCE 2 BASIC CONCEPTS: Multiple Regression Y i = β 1 + β 2 X 1i + β 3 X 2i + β 4 X 3i +

Ch3. TRENDS. Time Series Analysis

3.1 Deterministic Versus Stochastic Trends The simulated random walk in Exhibit 2.1 shows a upward trend. However, it is caused by a strong correlation between the series at nearby time points. The true

Regularization: Ridge Regression and the LASSO

Agenda Wednesday, November 29, 2006 Agenda Agenda 1 The Bias-Variance Tradeoff 2 Ridge Regression Solution to the l 2 problem Data Augmentation Approach Bayesian Interpretation The SVD and Ridge Regression

Apart from this page, you are not permitted to read the contents of this question paper until instructed to do so by an invigilator.

B. Sc. Examination by course unit 2014 MTH5120 Statistical Modelling I Duration: 2 hours Date and time: 16 May 2014, 1000h 1200h Apart from this page, you are not permitted to read the contents of this