06.06.2013 Views

Theory of Statistics - George Mason University

Theory of Statistics - George Mason University

Theory of Statistics - George Mason University

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Optimal Properties <strong>of</strong> the Moore-Penrose Inverse<br />

5.5 Applications 425<br />

The solution corresponding to the Moore-Penrose inverse is unique because<br />

that generalized inverse is unique. That solution is interesting for another<br />

reason.<br />

Theorem 5.9<br />

Let b ∗ be any solution to the normal equations (5.72), that is,<br />

and let<br />

then<br />

b ∗ = (X T X) − X T Y,<br />

β = (X T X) + X T Y<br />

β2 ≤ b ∗ 2.<br />

Pro<strong>of</strong>.<br />

To see that this solution has minimum norm, first factor Z, as<br />

X = QRU T ,<br />

and form the Moore-Penrose inverse as<br />

X + −1<br />

R<br />

= U 1 0<br />

<br />

Q<br />

0 0<br />

T .<br />

Now let<br />

β = X + Y.<br />

This is a least squares solution (that is, we have chosen a specific least squares<br />

solution).<br />

Now, let<br />

Q T Y =<br />

c1<br />

where c1 has exactly r elements and c2 has n − r elements, and let<br />

U T <br />

t1<br />

b = ,<br />

where b is the variable in the norm Y − Xb2 that we seek to minimize, and<br />

where t1 has r elements.<br />

Because multiplication by an orthogonal matrix does not change the norm,<br />

we have<br />

c2<br />

t2<br />

<br />

,<br />

Y − Xb2 = Q T (Y − XUU T b)2<br />

<br />

<br />

c1 R1<br />

= <br />

0 t1 <br />

<br />

−<br />

<br />

c2 0 0 t2<br />

<br />

<br />

c1<br />

= <br />

− R1t1 <br />

<br />

<br />

.<br />

<strong>Theory</strong> <strong>of</strong> <strong>Statistics</strong> c○2000–2013 James E. Gentle<br />

c2<br />

2<br />

2

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!