# Deriving a Polynomial Regression Model

• May 5th 2012, 09:56 AM
Halty
[SOLVED] Deriving a Polynomial Regression Model
Hi,

I've been learning about Polynomial regression from Internet resources. I came across this video: Polynomial Regression Model: Derivation: Part 1 of 2 - YouTube
which derives the matrix everyone uses. I think I understand most it but I have a couple of questions.

Overview of derivation
Link to a shorter derivation at the bottom or the overview.

If you can find $a_0, a_1$ and $a_2$ such that:
$S_y=\sum^{n}_{i=1}(y_i-(a_0+a_1x_i+a_2x^2_i))^2$
is minimised. You will have found the polynomial equation of the form:
$y=a_0+a_1x+a_2x^2$
with the least squared error.

To do this the partial derivative is taken for $a_0, a_1$ and $a_2$ which gives:
$\frac{\partial S_y}{\partial a_0} =\sum^{n}_{i=1}2(y_i-a_0-a_1x_i-a_2x^2_i)(-1)$
$\frac{\partial S_y}{\partial a_1} =\sum^{n}_{i=1}2(y_i-a_0-a_1x_i-a_2x^2_i)(-x_i)$
$\frac{\partial S_y}{\partial a_2} =\sum^{n}_{i=1}2(y_i-a_0-a_1x_i-a_2x^2_i)(-x^2_i)$

Setting the derivatives to 0 minimises $S_y$.
$\frac{\partial S_y}{\partial a_0} = 0;\frac{\partial S_y}{\partial a_1} = 0;\frac{\partial S_y}{\partial a_2} = 0;$

After expanding and simplifying you're left with:
$a_0\sum^{n}_{i=1}n + a_1\sum^{n}_{i=1}x_i + a_2\sum^{n}_{i=1}x^2_i = \sum^{n}_{i=1}y_i$
$a_0\sum^{n}_{i=1}x_i + a_1\sum^{n}_{i=1}x^2_i + a_2\sum^{n}_{i=1}x^3_i = \sum^{n}_{i=1}x_iy_i$
$a_0\sum^{n}_{i=1}x^2_i + a_1\sum^{n}_{i=1}x^3_i + a_2\sum^{n}_{i=1}x^4_i = \sum^{n}_{i=1}x^2_iy_i$

Which leaves some linear simultaneous equations which can solved.

Shorter alternate explanation:
Least-Squares Parabola

Questions
Q1
When the partial derivative is taken, the square used to find the least square error is changed to a multiply by two, which is then cancelled out later. As far as I can tell you would get the same answer trying to find the least cubed error or the least square root error or any (non zero) error. Does this mean that the resulting function would be the same for any (non zero) powered error?

If it is the case I can't quite wrap my head around how that would work. I would have thought the higher the power the more outliers would pull the function towards themselves.

Q2
It's never explained in the video, but why does setting the partial derivatives to 0 minimise $S_y$?

I think I've got this one.
$\frac{\partial^2 S_y}{\partial a_0^2} = n$
$\frac{\partial^2 S_y}{\partial a_1^2} = \sum^{n}_{i=1}x^2_i$
$\frac{\partial^2 S_y}{\partial a_2^2} = \sum^{n}_{i=1}x^4_i$

All three must always be greater or equal to 0. For it to be equal to 0 the entire data set needs to be populated by 0s or have no data. In either case polynomial regression probably wouldn't be used. Thus it is always at a minima.

Matthew
• May 5th 2012, 10:39 AM
SpringFan25
Re: Deriving a Polynomial Regression Model
Q1
its differenciation, so the power reduces by 1. if you used a cubic error you would be left with 3 * (some quadratic)=0 which will not necessarily end up at the same answer as 2*(some linear function)=0.
Q2
im not sure i follow your reasoning here. Are you trying to prove that the stationary point you just solved for could never be a maximum?
• May 5th 2012, 11:31 AM
Halty
Re: Deriving a Polynomial Regression Model

Q1
Sorry, I don't think I explained this well enough. Using just one of the equations, if we start with:
$S_y=\sum^{n}_{i=1}(y_i-(a_0+a_1x_i+a_2x^2_i))^2$
and take the partial derivative of $a_0$ we get:
$\frac{\partial S_y}{\partial a_0} =\sum^{n}_{i=1}2(y_i-a_0-a_1x_i-a_2x^2_i)(-1)$
From $\frac{\partial S_y}{\partial a_0} = 0$ we get:
$\sum^{n}_{i=1}2(y_i-a_0-a_1x_i-a_2x^2_i)(-1) = 0$
From here we can now divide both sides by 2 or whatever the power was in the first equation to get:
$\sum^{n}_{i=1}(y_i-a_0-a_1x_i-a_2x^2_i)(-1) = 0$
The same can be (and is) done to each of the equations.

Q2
At this point I should probably come clean: I'm not really a mathematician (more of a scientist) so it could be completely wrong. You're nearly right about what I was trying to do though. I was trying to show that it was a local minima (and thus actually minimising the squared error) using the second derivative test. Second derivative test - Wikipedia, the free encyclopedia

If it's wrong then I'd really like to know the actual method.

Thanks,
Matthew
• May 5th 2012, 12:06 PM
SpringFan25
Re: Deriving a Polynomial Regression Model
Q1
for a cubic error the equvalent would be:
$u_i^3 = (y_i - a_0 -a_1x_i - a_2x_i^2)^3$

$\frac{\delta S_y}{\delta a_0} = 0 = -3 \sum (y -a_0 -a_1x -a_2x^2)^2$

This does not have the same solutions as the normal equations in post #1. in fact, it basically says that the sum of squared residuals is 0; which has no solutions unless all data points are exactly on the predicted line.

Q2 your first post was right, i didn't read it properly ;)
• May 5th 2012, 12:18 PM
Halty
Re: Deriving a Polynomial Regression Model
Ah, yes. You're right. I feel a little bit stupid, but at least I got the answer.

Thanks. :)