Steply

Least squares line through (0, 1), (1, 2), (2, 2), (3, 4)

Step by step, with the reason for every step · Least squares line · Linear Algebra

Step by step

The question

Question: find the least squares line through (0, 1), (1, 2), (2, 2) and (3, 4).

(0, 1), (1, 2), (2, 2), (3, 4)

We want a line y = b + m x.

These 4 points are not all on one line, so we find the line that fits best.

line: y = b + m x (b: where it starts, m: slope)

1 One equation per point

Each point gives one equation: put its x and y into y = b + m x.

4 points, 4 equations, but only 2 unknowns, b and m.

(0, 1): b + 0m = 1
(1, 2): b + 1m = 2

As a matrix: A has a column of 1s for b and a column of the x values for m.

The right side holds the y values.

No line goes through all 4 points exactly.

[1, 0; 1, 1; 1, 2; 1, 3] × [b; m] = [1; 2; 2; 4]

2 Write the normal equations

The best b and m solve A transpose A, times x, equals A transpose y.

These are called the normal equations.

With them, the leftovers are at a right angle to the columns of A: as small as they can be.

ATA x = ATy

A transpose: the rows of A become its columns.

AT = [1, 1, 1, 1; 0, 1, 2, 3]

Careful!

A transpose goes in front of A, and in front of y.

Both sides get it.

Careful: Aᵀ in front, on both sides

3 Work out the two sides

A transpose times A.

ATA = [1, 1, 1, 1; 0, 1, 2, 3] × [1, 0; 1, 1; 1, 2; 1, 3]
ATA = [1 × 1 + 1 × 1 + 1 × 1 + 1 × 1, 1 × 0 + 1 × 1 + 1 × 2 + 1 × 3; 0 × 1 + 1 × 1 + 2 × 1 + 3 × 1, 0 × 0 + 1 × 1 + 2 × 2 + 3 × 3]

The top left counts the points: 4.

ATA = [4, 6; 6, 14]

A transpose times y.

ATy = [1, 1, 1, 1; 0, 1, 2, 3] × [1; 2; 2; 4]
ATy = [1 × 1 + 1 × 2 + 1 × 2 + 1 × 4; 0 × 1 + 1 × 2 + 2 × 2 + 3 × 4]
ATy = [9; 18]

4 Solve for b and m

The two equations as an augmented matrix: the unknowns are b and m.

= [4, 6, 9; 6, 14, 18]

We call the rows R 1 and R 2.

→ [1, 3 over 2, 9 over 4; 6, 14, 18]
→ [1, 3 over 2, 9 over 4; 0, 5, 9 over 2]
→ [1, 3 over 2, 9 over 4; 0, 1, 9 over 10]
→ [1, 0, 9 over 10; 0, 1, 9 over 10]

Row 1 says b = 9/10.

b = 910

Row 2 says m = 9/10.

m = 910

So the least squares line is y = 9/10 times x + 9/10.

y = 910x + 910

5 Check

Let's check.

The leftover of a point is its y minus the line's y.

For the best line, the leftovers add to 0, and x times leftover adds to 0.

Each y minus the line's y at that x: 1/10, 1/5, −7/10 and 2/5.

leftovers: 1 − 910 = 110, 2 − 95 = 15, 2 − 2710 = −710, 4 − 185 = 25

Added: 0.

added: 110 + 15 + (−710) + 25 = 0

x times leftover, added: 0.

Both are 0, so the answer is right!

x × leftover: 0 × 110 + 1 × 15 + 2 × (−710) + 3 × 25
= 0

Let's look back.

First: one equation per point.

Then: write the normal equations.

Then: work out the two sides.

Then: solve for b and m.

Last: check.

Answer
y = 910x + 910

Your turn

Try this one yourself, then check:

Least squares line: (0, 1), (1, 2), (2, 2), (3, 5)
Show the answer
y = 65x + 710

More from Least squares line

The whole lesson: Least squares line →