 Linear regression

In statistics, linear regression is an approach to modeling the relationship between a scalar variable y and one or more explanatory variables denoted X. The case of one explanatory variable is called simple regression. More than one explanatory variable is multiple regression.
In linear regression, data are modeled using linear functions, and unknown model parameters are estimated from the data. Such models are called linear models. Most commonly, linear regression refers to a model in which the conditional mean of y given the value of X is an affine function of X. Less commonly, linear regression could refer to a model in which the median, or some other quantile of the conditional distribution of y given X is expressed as a linear function of X. Like all forms of regression analysis, linear regression focuses on the conditional probability distribution of y given X, rather than on the joint probability distribution of y and X, which is the domain of multivariate analysis.
Linear regression was the first type of regression analysis to be studied rigorously, and to be used extensively in practical applications. This is because models which depend linearly on their unknown parameters are easier to fit than models which are nonlinearly related to their parameters and because the statistical properties of the resulting estimators are easier to determine.
Linear regression has many practical uses. Most applications of linear regression fall into one of the following two broad categories:
 If the goal is prediction, or forecasting, linear regression can be used to fit a predictive model to an observed data set of y and X values. After developing such a model, if an additional value of X is then given without its accompanying value of y, the fitted model can be used to make a prediction of the value of y.
 Given a variable y and a number of variables X_{1}, ..., X_{p} that may be related to y, linear regression analysis can be applied to quantify the strength of the relationship between y and the X_{j}, to assess which X_{j} may have no relationship with y at all, and to identify which subsets of the X_{j} contain redundant information about y.
Linear regression models are often fitted using the least squares approach, but they may also be fitted in other ways, such as by minimizing the “lack of fit” in some other norm (as with least absolute deviations regression), or by minimizing a penalized version of the least squares loss function as in ridge regression. Conversely, the least squares approach can be used to fit models that are not linear models. Thus, while the terms “least squares” and linear model are closely linked, they are not synonymous.
Contents
Introduction to linear regression
Given a data set of n statistical units, a linear regression model assumes that the relationship between the dependent variable y_{i} and the pvector of regressors x_{i} is linear. This relationship is modeled through a socalled “disturbance term” ε_{i} — an unobserved random variable that adds noise to the linear relationship between the dependent variable and regressors. Thus the model takes the form
where ′ denotes the transpose, so that x_{i}′β is the inner product between vectors x_{i} and β.
Often these n equations are stacked together and written in vector form as
where
Some remarks on terminology and general use:
 is called the regressand, endogenous variable, response variable, measured variable, or dependent variable (see dependent and independent variables.) The decision as to which variable in a data set is modeled as the dependent variable and which are modeled as the independent variables may be based on a presumption that the value of one of the variables is caused by, or directly influenced by the other variables. Alternatively, there may be an operational reason to model one of the variables in terms of the others, in which case there need be no presumption of causality.
 are called regressors, exogenous variables, explanatory variables, covariates, input variables, predictor variables, or independent variables (see dependent and independent variables, but not to be confused with independent random variables). The matrix X is sometimes called the design matrix.
 Usually a constant is included as one of the regressors. For example we can take x_{i1} = 1 for i = 1, ..., n. The corresponding element of β is called the intercept. Many statistical inference procedures for linear models require an intercept to be present, so it is often included even if theoretical considerations suggest that its value should be zero.
 Sometimes one of the regressors can be a nonlinear function of another regressor or of the data, as in polynomial regression and segmented regression. The model remains linear as long as it is linear in the parameter vector β.
 The regressors x_{i} may be viewed either as random variables, which we simply observe, or they can be considered as predetermined fixed values which we can choose. Both interpretations may be appropriate in different cases, and they generally lead to the same estimation procedures; however different approaches to asymptotic analysis are used in these two situations.
 is a pdimensional parameter vector. Its elements are also called effects, or regression coefficients. Statistical estimation and inference in linear regression focuses on β.
 is called the error term, disturbance term, or noise. This variable captures all other factors which influence the dependent variable y_{i} other than the regressors x_{i}. The relationship between the error term and the regressors, for example whether they are correlated, is a crucial step in formulating a linear regression model, as it will determine the method to use for estimation.
Example. Consider a situation where a small ball is being tossed up in the air and then we measure its heights of ascent h_{i} at various moments in time t_{i}. Physics tells us that, ignoring the drag, the relationship can be modeled as
where β_{1} determines the initial velocity of the ball, β_{2} is proportional to the standard gravity, and ε_{i} is due to measurement errors. Linear regression can be used to estimate the values of β_{1} and β_{2} from the measured data. This model is nonlinear in the time variable, but it is linear in the parameters β_{1} and β_{2}; if we take regressors x_{i} = (x_{i1}, x_{i2}) = (t_{i}, t_{i}^{2}), the model takes on the standard form
 h_{i} = x'_{i}β + ε_{i}.
Assumptions
Two key assumptions are common to all estimation methods used in linear regression analysis:
 The design matrix X must have full column rank p. For this property to hold, we must have n > p, where n is the sample size (this is a necessary but not a sufficient condition). If this condition fails this is called the multicollinearity in the regressors. In this case the parameter vector β will be not identifiable — at most we will be able to narrow down its value to some linear subspace of R^{p}.
Methods for fitting linear models with multicollinearity have been developed,^{[1]}^{[2]}^{[3]}^{[4]} but require additional assumptions such as “effect sparsity” — that a large fraction of the effects are exactly zero.
 A simpler statement of this is that there must be enough data available compared to the number of parameters to be estimated. If there is too little data, then you end up with a system of equations with no unique solution. See partial least squares regression.
 The regressors x_{i} are assumed to be errorfree, that is they are not contaminated with measurement errors. Although not realistic in many settings, dropping this assumption leads to significantly more difficult errorsinvariables models.
Beyond these two assumptions, several other statistical properties of the data strongly influence the performance of different estimation methods:
 Some estimation methods are based on a lack of correlation, among the n observations . Statistical independence of the observations is not needed, although it can be exploited if it is known to hold.
 The statistical relationship between the error terms and the regressors plays an important role in determining whether an estimation procedure has desirable sampling properties such as being unbiased and consistent.
 The variances of the error terms may be equal across the n units (termed homoscedasticity) or not (termed heteroscedasticity). Some linear regression estimation methods give less precise parameter estimates and misleading inferential quantities such as standard errors when substantial heteroscedasticity is present.
 The arrangement, or probability distribution of the predictor variables x has a major influence on the precision of estimates of β. Sampling and design of experiments are highlydeveloped subfields of statistics that provide guidance for collecting data in such a way to achieve a precise estimate of β.
Interpretation
A fitted linear regression model can be used to identify the relationship between a single predictor variable x_{j} and the response variable y when all the other predictor variables in the model are “held fixed”. Specifically, the interpretation of β_{j} is the expected change in y for a oneunit change in x_{j} when the other covariates are held fixed. This is sometimes called the unique effect of x_{j} on y. In contrast, the marginal effect of x_{j} on y can be assessed using a correlation coefficient or simple linear regression model relating x_{j} to y.
Care must be taken when interpreting regression results, as some of the regressors may not allow for marginal changes (such as dummy variables, or the intercept term), while others cannot be held fixed (recall the example from the introduction: it would be impossible to “hold t_{i} fixed” and at the same time change the value of t_{i}^{2}).
It is possible that the unique effect can be nearly zero even when the marginal effect is large. This may imply that some other covariate captures all the information in x_{j}, so that once that variable is in the model, there is no contribution of x_{j} to the variation in y. Conversely, the unique effect of x_{j} can be large while its marginal effect is nearly zero. This would happen if the other covariates explained a great deal of the variation of y, but they mainly explain variation in a way that is complementary to what is captured by x_{j}. In this case, including the other variables in the model reduces the part of the variability of y that is unrelated to x_{j}, thereby strengthening the apparent relationship with x_{j}.
The meaning of the expression “held fixed” may depend on how the values of the predictor variables arise. If the experimenter directly sets the values of the predictor variables according to a study design, the comparisons of interest may literally correspond to comparisons among units whose predictor variables have been “held fixed” by the experimenter. Alternatively, the expression “held fixed” can refer to a selection that takes place in the context of data analysis. In this case, we “hold a variable fixed” by restricting our attention to the subsets of the data that happen to have a common value for the given predictor variable. This is the only interpretation of “held fixed” that can be used in an observational study.
The notion of a “unique effect” is appealing when studying a complex system where multiple interrelated components influence the response variable. In some cases, it can literally be interpreted as the causal effect of an intervention that is linked to the value of a predictor variable. However, it has been argued that in many cases multiple regression analysis fails to clarify the relationships between the predictor variables and the response variable when the predictors are correlated with each other and are not assigned following a study design.^{[5]}
Estimation methods
Numerous procedures have been developed for parameter estimation and inference in linear regression. These methods differ in computational simplicity of algorithms, presence of a closedform solution, robustness with respect to heavytailed distributions, and theoretical assumptions needed to validate desirable statistical properties such as consistency and asymptotic efficiency.
Some of the more common estimation techniques for linear regression are summarized below.
 Ordinary least squares (OLS) is the simplest and thus most common estimator. It is conceptually simple and computationally straightforward. OLS estimates are commonly used to analyze both experimental and observational data.
The OLS method minimizes the sum of squared residuals, and leads to a closedform expression for the estimated value of the unknown parameter β:
In simple linear regression, where there is only one regressor (with a constant), the OLS coefficient estimates have a simple form that is closely related to the correlation coefficient between the covariate and the response.  Generalized least squares (GLS) is an extension of the OLS method, that allows efficient estimation of β when either heteroscedasticity, or correlations, or both are present among the error terms of the model, as long as the form of heteroscedasticity and correlation is known independently of the data. To handle heteroscedasticity when the error terms are uncorrelated with each other, GLS minimizes a weighted analogue to the sum of squared residuals from OLS regression, where the weight for the i^{th} case is inversely proportional to var(ε_{i}). This special case of GLS is called “weighted least squares”. The GLS solution to estimation problem is
 Percentage least squares focuses on reducing percentage errors, which is useful in the field of forecasting or time series analysis. It is also useful in situations where the dependent variable has a wide range without constant variance, as here the larger residuals at the upper end of the range would dominate if OLS were used. When the percentage or relative error is normally distributed, least squares percentage regression provides maximum likelihood estimates. Percentage regression is linked to a multiplicative error model, whereas OLS is linked to models containing an additive error term.^{[7]}
 Iteratively reweighted least squares (IRLS) is used when heteroscedasticity, or correlations, or both are present among the error terms of the model, but where little is known about the covariance structure of the errors independently of the data.^{[8]} In the first iteration, OLS, or GLS with a provisional covariance structure is carried out, and the residuals are obtained from the fit. Based on the residuals, an improved estimate of the covariance structure of the errors can usually be obtained. A subsequent GLS iteration is then performed using this estimate of the error structure to define the weights. The process can be iterated to convergence, but in many cases, only one iteration is sufficient to achieve an efficient estimate of β.^{[9]}^{[10]}
 Instrumental variables regression (IV) can be performed when the regressors are correlated with the errors. In this case, we need the existence of some auxiliary instrumental variables z_{i} such that E[z_{i}ε_{i}] = 0. If Z is the matrix of instruments, then the estimator can be given in closed form as
 Optimal instruments regression is an extension of classical IV regression to the situation where E[ε_{i}z_{i}] = 0.
 Least absolute deviation (LAD) regression is a robust estimation technique in that it is less sensitive to the presence of outliers than OLS (but is less efficient than OLS when no outliers are present). It is equivalent to maximum likelihood estimation under a Laplace distribution model for ε.^{[11]}
 Quantile regression focuses on the conditional quantiles of y given X rather than the conditional mean of y given X. Linear quantile regression models a particular conditional quantile, often the conditional median, as a linear function β′x of the predictors.
 Maximum likelihood estimation can be performed when the distribution of the error terms is known to belong to a certain parametric family ƒ_{θ} of probability distributions.^{[12]} When f_{θ} is a normal distribution with mean zero and variance θ, the resulting estimate is identical to the OLS estimate. GLS estimates are maximum likelihood estimates when ε follows a multivariate normal distribution with a known covariance matrix.
 Adaptive estimation. If we assume that error terms are independent from the regressors , the optimal estimator is the 2step MLE, where the first step is used to nonparametrically estimate the distribution of the error term.^{[13]}
 Mixed models are widely used to analyze linear regression relationships involving dependent data when the dependencies have a known structure. Common applications of mixed models include analysis of data involving repeated measurements, such as longitudinal data, or data obtained from cluster sampling. They are generally fit as parametric models, using maximum likelihood or Bayesian estimation. In the case where the errors are modeled as normal random variables, there is a close connection between mixed models and generalized least squares.^{[14]} Fixed effects estimation is an alternative approach to analyzing this type of data.
 Principal component regression (PCR) ^{[3]}^{[4]} is used when the number of predictor variables is large, or when strong correlations exist among the predictor variables. This twostage procedure first reduces the predictor variables using principal component analysis then uses the reduced variables in an OLS regression fit. While it often works well in practice, there is no general theoretical reason that the most informative linear function of the predictor variables should lie among the dominant principal components of the multivariate distribution of the predictor variables. The partial least squares regression is the extension of the PCR method which does not suffer from the mentioned deficiency.
 Total least squares (TLS) ^{[15]} is an approach to least squares estimation of the linear regression model that treats the covariates and response variable in a more geometrically symmetric manner than OLS. It is one approach to handling the "errors in variables" problem, and is sometimes used when the covariates are assumed to be errorfree.
 Ridge regression,^{[16]}^{[17]}^{[18]} and other forms of penalized estimation such as the Lasso,^{[1]} deliberately introduce bias into the estimation of β in order to reduce the variability of the estimate. The resulting estimators generally have lower mean squared error than the OLS estimates, particularly when multicollinearity is present. They are generally used when the goal is to predict the value of the response variable y for values of the predictors x that have not yet been observed. These methods are not as commonly used when the goal is inference, since it is difficult to account for the bias.
 Leastangle regression ^{[2]} is an estimation procedure for linear regression models that was developed to handle highdimensional covariate vectors, potentially with more covariates than observations.
 The Theil–Sen estimator is a simple robust estimation technique that choose the slope of the fit line to be the median of the slopes of the lines through pairs of sample points. It has similar statistical efficiency properties to simple linear regression but is much less sensitive to outliers.^{[19]}
 Other robust estimation techniques, including the αtrimmed mean approach, and L, M, S, and Restimators have been introduced.
Further discussion
In statistics, the problem of numerical methods for linear least squares is an important one because linear regression models are one of the most important types of model, both as formal statistical models and for exploration of data sets. The majority of statistical computer packages contain facilities for regression analysis that make use of linear least squares computations. Hence it is appropriate that considerable effort has been devoted to the task of ensuring that these computations are undertaken efficiently and with due regard to numerical precision.
Individual statistical analyses are seldom undertaken in isolation, but rather are part of a sequence of investigatory steps. Some of the topics involved in considering numerical methods for linear least squares relate to this point. Thus important topics can be
 Computations where a number of similar, and often nested, models are considered for the same data set. That is, where models with the same dependent variable but different sets of independent variables are to be considered, for essentially the same set of data points.
 Computations for analyses that occur in a sequence, as the number of data points increases.
 Special considerations for very extensive data sets.
Fitting of linear models by least squares often, but not always, arises in the context of statistical analysis. It can therefore be important that considerations of computational efficiency for such problems extend to all of the auxiliary quantities required for such analyses, and are not restricted to the formal solution of the linear least squares problem.
Matrix calculations, like any others, are affected by rounding errors. An early summary of these effects, regarding the choice of computational methods for matrix inversion, was provided by Wilkinson.^{[20]}
Extensions
 General linear model considers the situation when the response variable y is not a scalar but a vector. Conditional linearity of E(yx) = Bx is still assumed, with a matrix B replacing the vector β of the classical linear regression model. Multivariate analogues of OLS and GLS have been developed.
 Generalized linear models are a framework for modeling a response variable y in the form g(β′x) + ε, where g is an arbitrary link function. Single index models allow some degree of nonlinearity in the relationship between x and y, while preserving the central role of the linear predictor β′x as in the classical linear regression model. Under certain conditions, simply applying OLS to data from a singleindex model will consistently estimate β up to a proportionality constant.^{[21]}
 Hierarchical linear models (or multilevel regression) organizes the data into a hierarchy of regressions, for example where A is regressed on B, and B is regressed on C. It is often used where the data have a natural hierarchical structure such as in educational statistics, where students are nested in classrooms, classrooms are nested in schools, and schools are nested in some administrative grouping such as a school district. The response variable might be a measure of student achievement such as a test score, and different covariates would be collected at the classroom, school, and school district levels.
 Errorsinvariables models (or “measurement error models”) extend the traditional linear regression model to allow the predictor variables X to be observed with error. This error causes standard estimators of β to become biased. Generally, the form of bias is an attenuation, meaning that the effects are biased toward zero.
 In Dempster–Shafer theory, or a linear belief function in particular, a linear regression model may be represented as a partially swept matrix, which can be combined with similar matrices representing observations and other assumed normal distributions and state equations. The combination of swept or unswept matrices provides an alternative method for estimating linear regression models.
Applications of linear regression
Linear regression is widely used in biological, behavioral and social sciences to describe possible relationships between variables. It ranks as one of the most important tools used in these disciplines.
Trend line
Main article: Trend estimationA trend line represents a trend, the longterm movement in time series data after other components have been accounted for. It tells whether a particular data set (say GDP, oil prices or stock prices) have increased or decreased over the period of time. A trend line could simply be drawn by eye through a set of data points, but more properly their position and slope is calculated using statistical techniques like linear regression. Trend lines typically are straight lines, although some variations use higher degree polynomials depending on the degree of curvature desired in the line.
Trend lines are sometimes used in business analytics to show changes in data over time. This has the advantage of being simple. Trend lines are often used to argue that a particular action or event (such as training, or an advertising campaign) caused observed changes at a point in time. This is a simple technique, and does not require a control group, experimental design, or a sophisticated analysis technique. However, it suffers from a lack of scientific validity in cases where other potential changes can affect the data.
Epidemiology
Early evidence relating tobacco smoking to mortality and morbidity came from observational studies employing regression analysis. In order to reduce spurious correlations when analyzing observational data, researchers usually include several variables in their regression models in addition to the variable of primary interest. For example, suppose we have a regression model in which cigarette smoking is the independent variable of interest, and the dependent variable is lifespan measured in years. Researchers might include socioeconomic status as an additional independent variable, to ensure that any observed effect of smoking on lifespan is not due to some effect of education or income. However, it is never possible to include all possible confounding variables in an empirical analysis. For example, a hypothetical gene might increase mortality and also cause people to smoke more. For this reason, randomized controlled trials are often able to generate more compelling evidence of causal relationships than can be obtained using regression analyses of observational data. When controlled experiments are not feasible, variants of regression analysis such as instrumental variables regression may be used to attempt to estimate causal relationships from observational data.
Finance
The capital asset pricing model uses linear regression as well as the concept of Beta for analyzing and quantifying the systematic risk of an investment. This comes directly from the Beta coefficient of the linear regression model that relates the return on the investment to the return on all risky assets.
Economics
Linear regression is the predominant empirical tool in economics. For example, it is used to predict consumption spending,^{[22]} fixed investment spending, inventory investment, purchases of a country's exports,^{[23]} spending on imports,^{[23]} the demand to hold liquid assets,^{[24]} labor demand,^{[25]} and labor supply.^{[25]}
Environmental science
Linear regression finds application in a wide range of environmental science applications. In Canada, the Environmental Effects Monitoring Program uses statistical analyses on fish and benthic surveys to measure the effects of pulp mill or metal mine effluent on the aquatic ecosystem.^{[26]}
See also
 Analysis of variance
 Anscombe's quartet
 Crosssectional regression
 Curve fitting
 Empirical Bayes methods
 Logistic regression
 Mestimator
 Nonlinear regression
 Nonparametric regression
 Multivariate adaptive regression splines
 Lackoffit sum of squares
 Truncated regression model
 Censored regression model
 Simple linear regression
 Segmented linear regression
Further reading
 Pedhazur, Elazar J (1982). Multiple regression in behavioral research: Explanation and prediction (2nd ed.). New York: Holt, Rinehart and Winston. ISBN 0030417600
 Barlow, Jesse L. (1993). "Chapter 9: Numerical aspects of Solving Linear Least Squares Problems". In Rao, C.R.. Computational Statistics. Handbook of Statistics. 9. NorthHolland. ISBN 0444880968
 Björck, Åke (1996). Numerical methods for least squares problems. Philadelphia: SIAM. ISBN 0898713609.
 Goodall, Colin R. (1993). "Chapter 13: Computation using the QR decomposition". In Rao, C.R.. Computational Statistics. Handbook of Statistics. 9. NorthHolland. ISBN 0444880968
 National Physical Laboratory (1961). "Chapter 1: Linear Equations and Matrices: Direct Methods". Modern Computing Methods. Notes on Applied Science. 16 (2nd ed.). Her Majesty's Stationery Office
 National Physical Laboratory (1961). "Chapter 2: Linear Equations and Matrices: Direct Methods on Automatic Computers". Modern Computing Methods. Notes on Applied Science. 16 (2nd ed.). Her Majesty's Stationery Office
Notes
 ^ ^{a} ^{b} Tibshirani, Robert (1996). "Regression Shrinkage and Selection via the Lasso". Journal of the Royal Statistical Society. Series B (Methodological) 58 (1): 267–288. JSTOR 2346178.
 ^ ^{a} ^{b} Efron, Bradley; Hastie,Trevor; Johnstone,Iain Johnstone;Tibshirani,Robert (2004). "Least Angle Regression". The Annals of Statistics 32 (2): 407–451. doi:10.1214/009053604000000067. JSTOR 3448465.
 ^ ^{a} ^{b} Hawkins, Douglas M. (1973). "On the Investigation of Alternative Regressions by Principal Component Analysis". Journal of the Royal Statistical Society. Series C (Applied Statistics) 22 (3): 275–286. JSTOR 2346776.
 ^ ^{a} ^{b} Jolliffe, Ian T. (1982). "A Note on the Use of Principal Components in Regression". Journal of the Royal Statistical Society. Series C (Applied Statistics) 31 (3): 300–303. JSTOR 2348005.
 ^ Berk, Richard A.. Regression Analysis: A Constructive Critique. Sage. doi:10.1177/0734016807304871.
 ^ Lai, T.L.; Robbins,H; Wei, C.Z. (1978). "Strong consistency of least squares estimates in multiple regression". Proceedings of the National Academy of Sciences USA 75 (7).
 ^ Tofallis, C (2009). "Least Squares Percentage Regression". Journal of Modern Applied Statistical Methods 7: 526–534. http://papers.ssrn.com/sol3/papers.cfm?abstract_id=1406472.
 ^ del Pino, Guido (1989). "The Unifying Role of Iterative Generalized Least Squares in Statistical Algorithms". Statistical Science 4 (4): 394–403. doi:10.1214/ss/1177012408. JSTOR 2245853.
 ^ Carroll, Raymond J. (1982). "Adapting for Heteroscedasticity in Linear Models". The Annals of Statistics 10 (4): 1224–1233. doi:10.1214/aos/1176345987. JSTOR 2240725.
 ^ Cohen, Michael; Dalal, Siddhartha R.; Tukey,John W. (1993). "Robust, Smoothly Heterogeneous Variance Regression". Journal of the Royal Statistical Society. Series C (Applied Statistics) 42 (2): 339–353. JSTOR 2986237.
 ^ Narula, Subhash C.; Wellington, John F. (1982). "The Minimum Sum of Absolute Errors Regression: A State of the Art Survey". International Statistical Review 50 (3): 317–326. doi:10.2307/1402501. JSTOR 1402501.
 ^ Lange, Kenneth L.; Little, Roderick J. A.; Taylor,Jeremy M. G. (1989). "Robust Statistical Modeling Using the t Distribution". Journal of the American Statistical Association 84 (408): 881–896. doi:10.2307/2290063. JSTOR 2290063.
 ^ Stone, C. J. (1975). "Adaptive maximum likelihood estimators of a location parameter". The Annals of Statistics 3 (2): 267–284. doi:10.1214/aos/1176343056. JSTOR 2958945.
 ^ Goldstein, H. (1986). "Multilevel Mixed Linear Model Analysis Using Iterative Generalized Least Squares". Biometrika 73 (1): 43–56. doi:10.1093/biomet/73.1.43. JSTOR 2336270.
 ^ Nievergelt, Yves (1994). "Total Least Squares: StateoftheArt Regression in Numerical Analysis". SIAM Review 36 (2): 258–264. doi:10.1137/1036055. JSTOR 2132463.
 ^ Swindel, Benee F. (1981). "Geometry of Ridge Regression Illustrated". The American Statistician 35 (1): 12–15. doi:10.2307/2683577. JSTOR 2683577.
 ^ Draper, Norman R.; van Nostrand,R. Craig (1979). "Ridge Regression and JamesStein Estimation: Review and Comments". Technometrics 21 (4): 451–466. doi:10.2307/1268284. JSTOR 1268284.
 ^ Hoerl, Arthur E.; Kennard,Robert W.; Hoerl,Roger W. (1985). "Practical Use of Ridge Regression: A Challenge Met". Journal of the Royal Statistical Society. Series C (Applied Statistics) 34 (2): 114–120. JSTOR 2347363.
 ^ Theil, H. (1950). "A rankinvariant method of linear and polynomial regression analysis. I, II, III". Nederl. Akad. Wetensch., Proc. 53: 386–392, 521–525, 1397–1412. MR0036489; Sen, Pranab Kumar (1968). "Estimates of the regression coefficient based on Kendall's tau". Journal of the American Statistical Association 63 (324): 1379–1389. doi:10.2307/2285891. JSTOR 2285891. MR0258201.
 ^ Wilkinson, J.H. (1963) "Chapter 3: Matrix Computations", Rounding Errors in Algebraic Processes, London: Her Majesty's Stationery Office (National Physical Laboratory, Notes in Applied Science, No.32)
 ^ Brillinger, David R. (1977). "The Identification of a Particular Nonlinear Time Series System". Biometrika 64 (3): 509–515. doi:10.1093/biomet/64.3.509. JSTOR 2345326.
 ^ Deaton, Angus (1992). Understanding Consumption. Oxford University Press. ISBN 0198288247.
 ^ ^{a} ^{b} Krugman, P.R.; M. Obstfeld (1988). International Economics: Theory and Policy. Glenview: Scott, Foresman.
 ^ Laidler, David E.W. (1993). The Demand for Money: Theories, Evidence, and Problems, 4th ed. Description.
 ^ ^{a} ^{b} Ehrenberg and Smith, Modern Labor Economics, HarperCollins, 2008.
 ^ EEMP webpage
References
 Cohen, J., Cohen P., West, S.G., & Aiken, L.S. (2003). Applied multiple regression/correlation analysis for the behavioral sciences. (2nd ed.) Hillsdale, NJ: Lawrence Erlbaum Associates
 Charles Darwin. The Variation of Animals and Plants under Domestication. (1868) (Chapter XIII describes what was known about reversion in Galton's time. Darwin uses the term "reversion".)
 Draper, N.R. and Smith, H. Applied Regression Analysis Wiley Series in Probability and Statistics (1998)
 Francis Galton. "Regression Towards Mediocrity in Hereditary Stature," Journal of the Anthropological Institute, 15:246263 (1886). (Facsimile at: [1])
 Robert S. Pindyck and Daniel L. Rubinfeld (1998, 4h ed.). Econometric Models and Economic Forecasts,, ch. 1 (Intro, incl. appendices on Σ operators & derivation of parameter est.) & Appendix 4.3 (mult. regression in matrix form).
 Kaw, Autar; Kalu, Egwu (2008). "Numerical Methods with Applications". [2]., Chapter 6 deals with linear and nonlinear regression.
External links
Least squares and regression analysis Computational statistics Least squares · Linear least squares · Nonlinear least squares · Iteratively reweighted least squaresCorrelation and dependence Regression analysis Regression as a
statistical modelLinear regressionSimple linear regression · Ordinary least squares · Generalized least squares · Weighted least squares · General linear modelPredictor structureNonstandardNonnormal errorsDecomposition of variance Model exploration Background Mean and predicted response · Gauss–Markov theorem · Errors and residuals · Goodness of fit · Studentized residual · Minimum meansquare errorDesign of experiments Numerical approximation Applications Regression analysis category  Statistics category · Statistics portal · Statistics outline · Statistics topics
Wikimedia Foundation. 2010.