12
$\begingroup$

Let’s say we have the input (predictor) and output (response) data points A, B, C, D, E and we want to fit a line through the points. This is a simple problem to illustrate the question, but can be extended to higher dimensions as well.

Problem Statement

enter image description here

The current best fit or hypothesis is represented by the black line above. The blue arrow ($\color{blue}\rightarrow$) represents the vertical distance between the data point and the current best fit, by drawing a vertical line from the point till it intersects the line.

The green arrow ($\color{green}\rightarrow$) is drawn such that it is perpendicular to the current hypothesis at the point of intersection, and thus represents the least distance between the data point and the current hypothesis. For points A and B, a line drawn such that it is vertical to the current best guess and is similar to a line which is vertical to the x axis. For these two points, the blue and the green lines overlap, but they don’t for points C, D and E.

The least squares principle defines the cost function for linear regression by drawing a vertical line through the data points (A, B, C, D or E) to the estimated hypothesis ($\color{blue}\rightarrow$), at any given training cycle, and is represented by

$Cost Function = \sum_{i=1}^N(y_i-h_\theta(x_i))^2$

Here $(x_i, y_i)$ represents the data points, and $h_\theta(x_i)$ represents the best fit.

The minimum distance between a point (A, B, C, D or E) is represented by a perpendicular line drawn from that point on to the current best guess (green arrows).

The goal of least square function is to the define an objective function which when minimized would give rise to the least distance between the hypothesis and all the points combined, but won’t necessarily minimize the distance between the hypothesis and a single input point.

**Question**

Why don’t we define the Cost Function for linear regression as the least distance between the input data point and the hypothesis (defined by a line perpendicular to the hypothesis) passing through the input datapoin, as given by ($\color{green}\rightarrow$)?

$\endgroup$
  • 5
    $\begingroup$ Simple linear regression assumes that there is no error in the values of the observations' x-coordinates (e.g. because they are experimental manipulations). If there are errors on the x-axis one can account for them by minimizing a cost function similar to the one you propose; this requires to set the ratio between the variance of the errors on the x and y axis. If the ratio $=1$, this amounts to minimizing the perpendicular distance between the points and the line (orthogonal regression). If the ratio $\neq1$ it is called Deeming regression $\endgroup$ – matteo Dec 15 '17 at 18:41
  • $\begingroup$ See this post on PCA: cerebralmastication.com/2010/09/… $\endgroup$ – James Dec 16 '17 at 13:31
12
$\begingroup$

When you have noise in both the dependent variable (vertical errors) and the independent variable (horizontal errors), the least squares objective function can be modified to incorporate these horizontal errors. The problem in how to weight these two types of errors. This weighting usually depends on the ratio of the variances of the two errors:

  1. If the variance of the vertical error is extremely large relative to the variance of the horizontal error, OLS is correct.
  2. If the variance of the horizontal error is extremely large relative to the variance of the vertical error, inverse least squares (in which $x$ is regressed on $y$ and the inverse of the coefficient estimate for $y$ is used as the estimate of $\beta$) is appropriate.
  3. If the ratio of the variance of the vertical error to the variance of the horizontal error is equal to the ratio of the variances of the dependent and independent variables, we have the case of "diagonal" regression, in which a consistent esti­mate turns out to be the geometric mean of the OLS and inverse least squares estimators.
  4. If the ratio of these error variances is one, then we have the case of "orthogonal" regression, in which the sum of squared errors measured along a line perpendicular to the estimating line is minimized. This is what you had in mind.

In practice, the great drawback of this procedure is that the ratio of the error variances is not usually known and cannot usually be estimated, so the path forward is not clear.

$\endgroup$
  • $\begingroup$ I tried to edit to change "dependent" to "independent" in the first sentence but edits must be 6 characters. Maybe update the answer to fix the typo? $\endgroup$ – Ryan Stout Dec 15 '17 at 21:30
  • $\begingroup$ @RyanStout Thanks, and done. I think inserting spaces will get you around that. $\endgroup$ – Dimitriy V. Masterov Dec 15 '17 at 21:33
  • $\begingroup$ Now I'm a bit confused: aren't the vertical errors the errors in the dependent variable (y) and horizontal errors in the independent variable (x) ? $\endgroup$ – Ryan Stout Dec 15 '17 at 21:36
  • $\begingroup$ @RyanStout I messed it up again $\endgroup$ – Dimitriy V. Masterov Dec 15 '17 at 21:48
9
$\begingroup$

One reason is that $$\sum_{i=1}^N(y_i-h_\theta(x_i))^2$$ is relatively easy to compute and optimize, while the proposed cost $$\sum_{i=1}^N \min_{x,y}\big[(y_i-h_\theta(x))^2+(x_i-x)^2\big]$$ has a nested minimization problem which may be quite hard depending on the choice of family for $h_\theta(x)$.

$\endgroup$
  • $\begingroup$ Thats a good point. I was thinking how to go about calculating the Cost Function in general. $\endgroup$ – alpha_989 Dec 15 '17 at 19:02
  • $\begingroup$ I am not necessarily sure how to evaluate the distance between the point and a non-linear plane/surface, but to evaluate the distance between a point and a linear surface/plane, we may not need nested minimization: mathinsight.org/distance_point_plane $\endgroup$ – alpha_989 Dec 15 '17 at 19:02
  • $\begingroup$ Secondly, when we do use regression, our goal is to evaluate the weights to find the best fit. From what I understand, during the actual computation, we rarely evaluate the cost function, but some derivative of the cost function? $\endgroup$ – alpha_989 Dec 15 '17 at 19:02
  • 1
    $\begingroup$ @whuber. I see. Once we establish those meanings for those two terms, I agree that the problems being solved are different (is there or is there not the possibility of there being error in x). I don't think that you will get broad agreement from knowledgeable persons on the meaning of those terms, but that's a side point. $\endgroup$ – stochastic Dec 15 '17 at 23:39
  • 1
    $\begingroup$ @Stochastic I agree there may be fuzziness about the concept of "curve fitting," but the concept of regression I am invoking appears in writings by the best authorities. $\endgroup$ – whuber Dec 15 '17 at 23:41
2
$\begingroup$

The oversimplified version is that X is assumed to have no error. So if you look at point E in your plot for example, it's assumed that its X coordinate is precisely accurate. Typically this is the case when we can control X, in other words when we can set it to a specific value. In that case, the only error that can exist is in the Y direction, and that's why the error / cost function only includes the Y direction.

Whenever that's not the case, whenever we don't control X and X can have error, people do incorporate the X direction in the error function in something called type II or model II regression, and its variants. It can be tricky to do this if X and Y have different scales, so then you have to think about normalizations and such.

$\endgroup$
1
$\begingroup$

At the risk of being prosaic, the reason for the error function is that the standard interpretation is that the x is given and one is trying to best describe (or predict) the y component. So there is no error in the 'x'. For example you might try and understand (or predict) the closing price of a stock tomorrow based on today's closing price . Similarly one could try and understand the average temperature tomorrow in terms of today's average temperature. Obviously these examples are simple minded, but that is the idea. Incidentally something most people don't realize, but I think is clear from your examples, is that if one regresses y against x the regression line doesn't have to have any particular resemblance to the regression of x against y. Orthogonal regression is the term for a regression where one tries to find the line that minimizes the distance of points from a line. For example if one was trying to understand the relationship between the price of IBM stock and the price of AAPL stock, that would be the appropriate method.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.