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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

812 0 Statistical Mathematics<br />

The relationship (0.3.68) allows us to prove properties 1 and 4.<br />

The one property <strong>of</strong> square positive matrices that does not carry over to<br />

square irreducible nonnegative matrices is that r = ρ(A) is the only eigenvalue<br />

on the spectral circle <strong>of</strong> A. For example, the small irreducible nonnegative<br />

matrix<br />

<br />

0 1<br />

A =<br />

1 0<br />

has eigenvalues 1 and −1, and so both are on the spectral circle.<br />

It turns out, however, that square irreducible nonnegative matrices that<br />

have only one eigenvalue on the spectral circle also have other interesting<br />

properties that are important, for example, in Markov chains. We therefore<br />

give a name to the property:<br />

A square irreducible nonnegative matrix is said to be primitive if it<br />

has only one eigenvalue on the spectral circle.<br />

In modeling with Markov chains and other applications, the limiting behavior<br />

<strong>of</strong> A k is an important property.<br />

If A is a primitive matrix, then we have the useful result<br />

k A<br />

lim = vw<br />

k→∞ ρ(A)<br />

T , (0.3.70)<br />

where v is an eigenvector <strong>of</strong> A associated with ρ(A) and w is an eigenvector<br />

<strong>of</strong> A T associated with ρ(A), and w and v are scaled so that w T v = 1. (Such<br />

eigenvectors exist because ρ(A) is a simple eigenvalue. They also exist because<br />

they are both positive. Note that A is not necessarily symmetric, and so its<br />

eigenvectors may include imaginary components; however, the eigenvectors<br />

associated with ρ(A) are real, and so we can write w T instead <strong>of</strong> w H .)<br />

To see equation (0.3.70), we consider A − ρ(A)vw T . First, if (ci, vi) is<br />

an eigenpair <strong>of</strong> A − ρ(A)vw T and ci = 0, then (ci, vi) is an eigenpair <strong>of</strong> A.<br />

We can see this by multiplying both sides <strong>of</strong> the eigen-equation by vw T :<br />

hence,<br />

Next, we show that<br />

civw T vi = vw T A − ρ(A)vw T vi<br />

= vw T A − ρ(A)vw T vw T vi<br />

= ρ(A)vw T − ρ(A)vw T vi<br />

= 0;<br />

Avi = A − ρ(A)vw T vi<br />

= civi.<br />

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

ρ A − ρ(A)vw T < ρ(A). (0.3.71)

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

Saved successfully!

Ooh no, something went wrong!