Verteilung einer quadratischen Form, nicht zentrale Chi-Quadrat-Verteilung

7

Definition . Angenommen, yN(μ,In×n) .

Dann ist

w=yTy=y2χn2(θ=μ2/2=μTμ/2),
dh die χ2 Verteilung mit n Graden von Freiheits- und Nichtzentralitätsparameter θ .

Ich versuche, (am allerwenigsten) findet einen Beweis des folgenden Satzes:

Satz . Angenommen, Σ ist eine n×n positive definitive und symmetrische Matrix, A ist eine n×n symmetrische Matrix mit Rang m und (AΣ)2=AΣ .

Dann yN(μ,Σ)yTAyχm2(μTAμ/2).

Ich sehe nicht, wie diese Definition funktioniert. In meinen Augen müssten wir A in Potenzen von 1/2 "aufteilen" 1/2(was auch immer dies bedeutet), um die Definition, die ich oben gegeben habe, mit dem Satz zu arbeiten , dh (A1/2y)TA1/2y , es sei denn, ich vermisse hier wirklich etwas.

Ich lese nur selbst, daher wären alle Quellen oder Personen, die mir mit einem Beweis dafür helfen können, dankbar.

Klarinettist
quelle
Sind Sie sicher, dass Sie auch durch 2 teilen müssen?
JohnK
@JohnK Nach diesen Notizen, die ich lese, stellen sie fest, dass einige Autoren die Division nicht durch . 2
Klarinettist
1
Letzte Frage: Die Matrix muss symmetrisch sein, oder? A
JohnK

Antworten:

13

Um den Satz zu beweisen, benötigen wir das folgende Zwischenergebnis.

Satz . Sei und sei für eine symmetrische Matrix mit Rang . Wenn dann idempotent ist, hat Q eine -Verteilung.yNn(0,σ2In)Q=σ2yAyArAχ2(r)

Der Satz erstreckt sich auch in die andere Richtung, aber wir brauchen nur die Hinlänglichkeit, also werden wir dies nur beweisen und dies unter Verwendung der Eigenwert-Eigenvektor-Zerlegung einer idempotenten Matrix tun.

Es ist wichtig anzumerken, dass wir zuerst den Fall des Mittelwerts Null betrachten. Wir werden diese Annahme später lockern. Aber denken Sie vorerst daran, dass für eine quadratische Matrix mit Rang rA

A=i=1rλicici

Dabei sind die Lambdas die Eigenwerte und die die entsprechenden Eigenvektoren. Alles ziemlich normal bis jetzt. Wenn wir nun auf symmetrisch beschränken, passieren zwei Dinge:ciA

  1. Die Eigenwerte sind reelle Werte

  2. Eigenvektoren, die unterschiedlichen Eigenwerten entsprechen, sind orthogonal

Dies sind Konsequenzen des sogenannten Spektralsatzes der linearen Algebra, und Sie können jedes gute Lehrbuch für einen Beweis konsultieren. Was bedeutet das für uns? Sie werden gleich sehen, warum Symmetrie erforderlich ist. Schreiben wir die Zerlegung unseres in unserer quadratischen Form auf und sehen, was passiert.A

(1)σ2yAy=σ2y(i=1rλicici)y=i=1rλi(σ1ciy)2

Wir haben unsere quadratische Form als gewichtete quadratische Projektionen auf orthogonale Achsen geschrieben. Untersuchen wir nun die Verteilung von und , . Nach Grundregelnciycjyij

σ1ciyN(0,cici=1)

da die Eigenvektoren nicht eindeutig sind und daher ohne Verlust der Allgemeinheit neu skaliert werden können, um die Länge eins zu haben. Nächster,

Cov(σ1cix,σ1cjx)=σ2ciIncj=0,  ij

durch die zweite Implikation des Spektralsatzes. Somit sind unsere Summanden nicht korreliert. Durch die Normalität sind sie auch unabhängig. Es ist dann leicht zu erkennen, dass die Summe aus gewichteten Zufallsvariablen besteht (gewichtet mit den Eigenwerten). In diesem Thread wurde gefragt, ob dies unabhängig von der Verteilung folgt . Die Antwort ist natürlich nein.χ2χ2

Geben Sie die Idempotenz ein. Aus der Definition einer idempotenten Matrix und dem Eigenwert / Eigenvektor-Problem folgt leicht, dass eine idempotente Matrix Eigenwerte hat, die entweder eins oder null sind. Da die Matrix unter der Annahme den Rang , gibt es Eigenwerte gleich eins. DeshalbArr

σ2yAy=i=1r(σ1ciy)2χ2(r)

Und das vervollständigt den Beweis.


Was passiert jetzt, wenn der Vektor einen Mittelwert ungleich Null hat? Kein Schaden angerichtet, wir werden nur die Definition der nicht zentralen Verteilung verwenden, um zu folgern, dass wennyχ2

yNn(μ,σ2In)

dann

σ2yAy χ2(r,μAμ)

wobei der zweite Term den Nicht-Zentralitätsparameter angibt. (Aus Gewohnheit werde ich die Division um überspringen, aber Sie können es einfach auf Ihre Weise tun).2

Wir sind jetzt bereit, das gewünschte Ergebnis zu beweisen.

Satz . Angenommen, ist eine positive definitive und symmetrische Matrix, ist eine symmetrische Matrix mit Rang und .Σn×nAn×nr(AΣ)2=AΣ

DannyN(μ,Σ)Q=yAyχ2(r,μAμ).

Das ist uns gegeben

AΣAΣ=AΣ

woraus folgt das

AΣA=A

und daher können wir unsere quadratische Form umschreiben als

(2)Q=(Σ1/2y)Σ1/2AΣ1/2(Σ1/2y)

Da eine positiv definierte Matrix ist, ist ihre Quadratwurzel immer gut definiert und Sie können dies anhand der Eigenwert- / Eigenvektor-Zerlegung überprüfen (Sie können auch überprüfen, ob Gleichung ( ) Gleichung ( ) entspricht )!)Σ21

Natürlich würden Sie zustimmen, dass . Wenn wir vorübergehend annehmen, dass wir uns in der Situation des ersten Satzes, richtig? Nun, nicht genau - wir müssen immer noch zeigen, dass die mittlere Matrix idempotent ist, aber das ist unter unseren Annahmen einfach genug.Σ1/2yNn(Σ1/2μ,In)μ=0

Σ1/2AΣ1/2Σ1/2AΣ1/2=Σ1/2AΣ1/2

Wenn also , was impliziert, dass wir jetzt wechseln zurück zur Situation von ungleich Null bedeuten wiryNn(0,Σ)Qχ2(r)

Q=yAyχ2(r,μAμ)

wie erforderlich.

JohnK
quelle
Wenn Sie "Mittelwert" sagen (z. B. "... wir betrachten zuerst den Fall des Mittelwerts Null"), meinen Sie damit den Nichtzentralitätsparameter?
Klarinettist
@Clarinetist Ich meine den Mittelwert der Normalverteilung, , was wiederum bedeutet, dass der Nichtzentralitätsparameter ebenfalls Null ist. μ
JohnK
Ah ja. Okay. Soll das bei der ersten Satzanweisung ? yny
Klarinettist
@Clarinetist Ja, soll die Dimension betonen.
JohnK
Sehr gut geschrieben. Ich verfolge ungefähr 98% davon beim ersten Lesen, was angesichts der Erklärungen, die ich zu diesem und ähnlichen Themen versucht habe, erstaunlich ist. Ich muss ein paar Dinge lernen und überprüfen, aber ich verstehe diese Antwort größtenteils. Vielen Dank und viel Spaß mit den +12.
Klarinettist