Asymptotische Normalität einer quadratischen Form

11

Lassen Sie x ein Zufallsvektor aus gezeichnet sein P . Betrachten wir eine Probe . Definieren Sie und . Sei \ boldsymbol {\ mu}: = \ mathbb {E} _ {\ mathbf {x} \ sim P} [\ mathbf {x}] und C: = \ mathrm {cov} _ {\ mathbf {x} \ sim P} [\ mathbf {x}, \ mathbf {x}] .{xi}i=1ni.i.d.P C :=1x¯n:=1ni=1nxiC^:=1ni=1n(xix¯n)(xix¯n)μ:=ExP[x]C:=covxP[x,x]

Nehmen wir nach dem zentralen Grenzwertsatz an, dass

n(x¯nμ)dN(0,C),

wobei C eine Kovarianzmatrix mit vollem Rang ist.

Frage : Wie kann ich das beweisen (oder widerlegen)?

n(x¯n(C^+γnI)1x¯nμC1μ)dN(0,v2),

für einige v>0 und für einige γn0 so dass limnγn=0 ? Das sieht einfach aus. Aber ich konnte nicht genau herausfinden, wie ich das zeigen sollte. Dies ist keine Hausaufgabenfrage.

Mein Verständnis ist, dass die Delta-Methode es uns ermöglichen würde, leicht zu schließen

n(x¯nC1x¯nμC1μ)dN(0,v2),

oder

n(x¯n(C^+γnI)1x¯nμ(C^+γnI)1μ)dN(0,v2).

Diese sind ein bisschen anders als ich will. Beachten Sie die Kovarianzmatrizen in den beiden Begriffen. Ich habe das Gefühl, dass ich hier etwas sehr Triviales vermisse. Wenn es die Dinge einfacher macht, können wir alternativ auch ignorieren, dh und annehmen, dass invertierbar ist. Vielen Dank.γnγn=0C^

wij
quelle
2
Wir müssen etwas darüber wissen, wie auf 0 geht. Ist es eine Folge von Konstanten? Ich denke, dass Sie zuerst was meiner Meinung nach ein Ergebnis von Slutskys ist. Dann würde ich als schreiben . hat eine einschränkende Verteilung, die mit der Methode gefunden werden kann. Zuletzt können Sie versuchen zu zeigen, dass mit einer Wahrscheinlichkeit auf 0 geht. Obwohl ich nicht sicher bin, ob das gilt ...γnx¯nTγnIx¯np0C^C+bias(C^)x¯nTCx¯nδx¯nTbias(C^)x¯n
AdamO
γn ist eine Folge von Konstanten (nicht zufällig). Die Sequenz kann auf alles eingestellt werden, was die Konvergenz zum Funktionieren bringt (falls eine solche Sequenz vorhanden ist). Ich denke, ist wahr. Ich bin nicht gefolgt, warum wir das zuerst brauchen. Aber lassen Sie mich darüber nachdenken und den Rest mehr. :)x¯nIx¯np0
wij
2
Ich habe nicht erwähnt: Ihr Zögern, die Methode direkt anzuwenden und als erledigt zu bezeichnen, ist berechtigt. Ich denke, Sie können dies sorgfältig aufschreiben. Nützliche Theoreme für diese Art von Beweisen sind Slutskys, das Mann-Wald Continuous Mapping Theorem und das Cramer-Wold Theorem. δ
AdamO
Ich bin damit einverstanden, dass die von Ihnen genannten Ergebnisse nützlich sein könnten. Ich sehe aber immer noch nicht wie. Eigentlich denke ich auch, dass die asymptotische Verteilung möglicherweise keine Normalverteilung ist.
wij
Es scheint, dass dies komplizierter ist, als es scheint. Das arXiv-Papier beschreibt hier , was in hohen Dimensionen passiert. Ich kann kein Analogon mit fester Dimension finden, aber sie haben ein endlichdimensionales Argument in Abschnitt 3.
Greenparker

Antworten:

1

Bei der Verwendung der Delta-Methode treten einige Schwierigkeiten auf. Es ist bequemer, es von Hand abzuleiten.

Mit dem Gesetz der großen Zahl, . Daher . Wenden Sie den Satz von Slutsky an, wir haben Nach dem kontinuierlichen Mapping-Theorem haben wir Daher Nach dem Satz von Slutsky haben wir Die Kombination der beiden oben genannten Gleichheitsausbeuten C^PCC^+γnIPC

n(C^+γnI)1/2(X¯μ)dN(0,C1).
n(X¯μ)T(C^+γnI)1(X¯μ)di=1pλi1(C)χ12.
n(X¯μ)T(C^+γnI)1(X¯μ)P0.
nμT(C^+γnI)1(X¯μ)dN(0,μTC2μ).
n(X¯T(C^+γnI)1X¯μT(C^+γnI)1μ)=n((X¯μ)T(C^+γnI)1(X¯μ)2μT(C^+γnI)1(X¯μ))=2nμT(C^+γnI)1(X¯μ)+oP(1)dN(0,4μTC2μ).
Die verbleibende Aufgabe besteht darin, sich mit zu befassen. Leider konvergiert diese Termdosis NICHT gegen . Das Verhalten wird kompliziert und hängt vom dritten und vierten Moment ab.
n(μT(C^+γnI)1μμT(C)1μ).
0

Um einfach zu sein, nehmen wir unten an, dass normalverteilt sind und . Es ist ein Standardergebnis, dass wobei eine symmetrische Zufallsmatrix mit diagonalen Elementen als ist und nicht diagonale Elemente als . Somit ist durch Matrix-Taylor-Erweiterung haben wir Xiγn=o(n1/2)

n(C^C)dC1/2WC1/2,
WN(0,2)N(0,1)
n(C^+γnIC)dC1/2WC1/2,
(I+A)1IA+A2
n((C^+γnI)1C1)=nC1/2((C1/2(C^+γnI)C1/2)1I)C1/2=nC1(C^+γnIC)C1+OP(n1/2)dC1/2WC1/2.
Also
n(μT(C^+γnI)1μμT(C)1μ)dμTC1/2WC1/2μN(0,(μTC1μ)2).

Also

n(X¯T(C^+γnI)1X¯μTC1μ)dN(0,4μTC2μ+(μTC1μ)2).
kfeng123
quelle
1
Danke für deine Antwort. Es ist genau dieser Begriff, der nicht gegen 0 konvergiert, der das Ganze schwierig macht. Leider kann ich nicht davon ausgehen, dass normal verteilt ist. Aber ich schätze die Antwort immer noch. Wenn Sie kommentieren könnten, wie es vom dritten und vierten Moment abhängt (vielleicht mit Referenzen), wäre das hilfreich. Auch kann ich im Moment nicht erklären. Aber ich habe das Gefühl, dass es langsamer als zerfallen muss . Ich muss mir den Grund genauer überlegen. Xigammano(n1/2)
wij
Ich habe vergessen hinzuzufügen, dass in meinem Fall davon ausgegangen werden kann , dass in einem kompakten Satz lebt (falls erforderlich). Dies könnte bei den momentanen Bedingungen hilfreich sein. Xi
wij