WIP: In Arbeit
Nach p. 370 von Cramers mathematischen Methoden der Statistik von 1946 definierenHier ist die kumulative Verteilungsfunktion der Standardnormalverteilung . Infolge seiner Definition ist uns fast sicher, dass .Ξn=n(1−Φ(Zn)).
ΦN(0,1)0≤Ξn≤n
Betrachten Sie eine gegebene Realisierung unseres Probenraums. In diesem Sinne ist sowohl eine Funktion von und als auch eine Funktion von und . Für ein festes können wir als deterministische Funktion von und als deterministische Funktion von und , wodurch das Problem vereinfacht wird. Wir wollen Ergebnisse zeigen, die mit ziemlicher Sicherheit für alleω∈ΩZnnωΞnZn,nωωZnnΞnZnnω∈ΩDies ermöglicht es uns, unsere Ergebnisse von einer nicht deterministischen Analyse auf die nicht deterministische Umgebung zu übertragen.
Nach p. 374 von Cramers mathematischen Methoden der Statistik von 1946 gehen für den Moment davon aus (ich möchte später zurückkommen und einen Beweis liefern), dass wir zeigen können, dass (für jedes gegebene ) die folgende asymptotische Expansion gilt (unter Verwendung von Teilintegration und Definition von ):ω∈ΩΦ
2π−−√nΞn=1Zne−Z2n2(1+O(1Z2n)) as Zn→∞.(~)
Es ist klar, dass für jedes , und ist mit ziemlicher Sicherheit eine zunehmende Funktion von als , daher behaupten wir im Folgenden für (fast sicher alle) feste :Zn+1≥ZnnZnnn→∞ωZn→∞⟺n→∞.
Daraus folgt, dass wir haben (wobei Bezeichnet asymptotische Äquivalenz ):∼
2π−−√nΞn∼1Zne−1Z2n as Zn→∞n→∞.
Wie wir im Folgenden vorgehen, entspricht im Wesentlichen der Methode des dominanten Gleichgewichts , und unsere Manipulationen werden formal durch das folgende Lemma gerechtfertigt:
Lemma: Es sei angenommen , daß als , und (also ). Dann müssen wir bei jeder Funktion die durch Zusammensetzungen, Additionen und Multiplikationen von Logarithmen und Potenzgesetzen gebildet wird (im Wesentlichen jede " Polylog " -Funktion), auch die als :Mit anderen Worten, solche "Polylog" -Funktionen bewahren die asymptotische Äquivalenz .f(n)∼g(n)n→∞f(n)→∞g(n)→∞hn → ∞ h ( f ( n ) ) ∼ h ( g ( n ) )n→∞h(f(n))∼h(g(n)).
Die Wahrheit dieses Lemmas ist eine Konsequenz von Satz 2.1. wie hier geschrieben . Beachten Sie auch, dass das Folgende meistens eine erweiterte (detailliertere) Version der Antwort auf eine ähnliche Frage ist, die hier zu finden ist .
Wenn wir die Logarithmen beider Seiten nehmen, erhalten wir Folgendes:
log(2π−−√Ξn)−logn∼−logZn−Z2n2.(1)
Hier ist Cramer etwas käfig; er sagt nur "vorausgesetzt ist begrenzt", wir können schließen. Aber zu zeigen, dass ziemlicher Sicherheit angemessen begrenzt ist, scheint eigentlich etwas nicht trivial zu sein. Es scheint, dass der Beweis dafür im Wesentlichen Teil dessen ist, was auf den Seiten 265-267 von Galambos besprochen wird, aber ich bin mir nicht sicher, da ich immer noch daran arbeite, den Inhalt dieses Buches zu verstehen.ΞnΞn
Unter der AnnahmelogΞn=o(logn) , dass man zeigen kann, dass , folgt daraus (da der Term dominiert ), dass:−Z2n/2−logZn
−logn∼−Z2n2⟹Zn∼2logn−−−−−√.
Das ist etwas schön, da es bereits das meiste ist, was wir zeigen wollen, obwohl es wieder erwähnenswert ist, dass es im Wesentlichen nur die Dose die Straße hinunter tritt, da wir jetzt eine gewisse fast von . Andererseits hat die gleiche Verteilung für jedes Maximum von iid kontinuierlichen Zufallsvariablen, so dass dies nachvollziehbar sein kann.ΞnΞn
Wie auch immer, wenn as ist, kann man eindeutig auch schließen, dass für jedes Das ist als . Unter Verwendung unseres Lemmas über Polylog-Funktionen, bei denen die asymptotische Äquivalenz oben erhalten bleibt, können wir diesen Ausdruck wieder in , um Folgendes zu erhalten:Zn∼2logn−−−−−√Zn∼2logn−−−−−√(1+α(n))α(n)o(1)n→∞(1)
log(2π−−√Ξn)−logn∼−log(1+α)−12log2−12loglogn−logn−2αlogn−α2logn.
⟹−log(Ξn2π−−√)∼log(1+α)+12log2+12loglogn+2αlogn+α2logn.
Hier müssen wir noch weiter gehenlogΞn=o(loglogn) as n→∞ und davon ausgehen, dass fast sicher ist . Wieder sagt Cramer nur "vorausgesetzt, ist begrenzt". Aber da alles, was man a priori über sagen kann , as ist, scheint es kaum klar zu sein, dass man fast sicher haben sollte, was die Substanz von Cramers Behauptung zu sein scheint.ΞnΞn0≤Xin≤nΞn=O(1)
Aber wie auch immer, ein unter der Annahme , glaubt , dass, dann folgt daraus , dass die dominierende Begriff, der nicht enthält heißt . Da , folgt, dass und eindeutig , so enthält der dominante Term IS . Daher können wir das neu anordnen und (indem wir alles durch oder teilen ) feststellenα12loglognα=o(1)α2=o(α)log(1+α)=o(α)=o(o(αlogn))α2αlogn12loglogn2αlogn
−12loglogn∼2αlogn⟹α∼−loglogn4logn.
Wenn wir dies wieder in das Obige einsetzen, erhalten wir Folgendes:
Zn∼2logn−−−−−√−loglogn22logn−−−−−√,
wieder vorausgesetzt, wir glauben bestimmte Dinge über .Ξn
Wir wiederholen dieselbe Technik erneut. da , folgt auch diesZn∼2logn−−−−−√−loglogn22logn√
Zn∼2logn−−−−−√−loglogn22logn−−−−−√(1+β(n))=2logn−−−−−√(1−loglogn8logn(1+β(n))),
wenn . Lassen Sie uns ein wenig vereinfachen, bevor wir direkt in (1) zurückkehren. wir bekommen das:β(n)=o(1)
logZn∼log(2logn−−−−−√)+log(1−loglogn8logn(1+β(n)))log(O(1))=o(logn)∼log(2logn−−−−−√).
Z2n2∼logn−12loglogn(1+β)+(loglogn)28logn(1β)2o((1+β)loglogn)∼logn−12(1+β)loglogn.
Wenn wir dies wieder in (1) einsetzen, stellen wir fest, dass:
log(2π−−√Ξn)−logn∼−log(2logn−−−−−√)−logn+12(1+β)loglogn⟹β∼log(4πΞ2n)loglogn.
Daraus schließen wir fast sicher
Zn∼2logn−−−−−√−loglogn22logn−−−−−√(1+log(4π)+2log(Ξn)loglogn)=2logn−−−−−√−loglogn+log(4π)22logn−−−−−√−log(Ξn)2logn−−−−−√.
Dies entspricht dem Endergebnis auf S.374 von Cramers Mathematical Methods of Statistics von 1946, außer dass hier die genaue Reihenfolge des Fehlerterms nicht angegeben ist. Anscheinend ergibt die Anwendung dieses weiteren Terms die genaue Reihenfolge des Fehlerterms, aber es scheint ohnehin nicht notwendig zu sein, die Ergebnisse über die Maxima der iid-Standardnormalen zu beweisen, an denen wir interessiert sind.
In Anbetracht des oben genannten Ergebnisses, nämlich mit ziemlicher Sicherheit:
Zn∼2logn−−−−−√−loglogn+log(4π)22logn−−−−−√−log(Ξn)2logn−−−−−√⟹Zn=2logn−−−−−√−loglogn+log(4π)22logn−−−−−√−log(Ξn)2logn−−−−−√+o(1).(†)
2. Dann folgt aus der Linearität der Erwartung:
EZn=2logn−−−−−√−loglogn+log(4π)22logn−−−−−√−E[log(Ξn)]2logn−−−−−√+o(1)⟹EZn2logn−−−−−√=1−E[logΞn]2logn+o(1).
Deshalb haben wir das gezeigt
limn→∞EZn2logn−−−−−√=1,
solange wir das auch zeigen können
E[logΞn]=o(logn).
Dies ist möglicherweise nicht allzu schwierig , da für jede kontinuierliche Zufallsvariable wieder dieselbe Verteilung hat . Somit haben wir das zweite Ergebnis von oben.Ξn
1. In ähnlicher Weise haben wir auch von oben, dass fast sicher:
Zn2logn−−−−−√=1−log(Ξn)2logn+o(1),.
Wenn wir das zeigen können:
log(Ξn)=o(logn) almost surely,(*)
dann haben wir das erste Ergebnis von oben gezeigt. Ergebnis (*) würde auch eindeutig ein Fortiori implizieren, dass , wodurch wir auch das erste Ergebnis von oben erhalten.E[log(Ξn)]=o(logn)
Beachten Sie auch, dass wir im obigen Beweis von ( ) ohnehin davon ausgehen mussten, dass fast sicher (oder zumindest etwas Ähnliches) ist, damit wir dann ( ) zeigen können Wir werden höchstwahrscheinlich auch den Prozess haben, der erforderlich ist, um fast sicher zu zeigen, und wenn wir daher beweisen können wir höchstwahrscheinlich in der Lage sein, sofort alle folgenden Schlussfolgerungen zu ziehen.†Ξn=o(logn)†Ξn=o(logn)(†)
3. Wenn wir jedoch dieses Ergebnis haben, dann verstehe ich nicht, wie man auch das , da . Aber zumindest scheint es wahr zu sein, dassEZn=2logn−−−−−√+Θ(1)o(1)≠Θ(1)EZn=2logn−−−−−√+O(1).
Dann können wir uns also darauf konzentrieren, die Frage zu beantworten, wie manΞn=o(logn) almost surely.
Wir müssen auch die Grunzarbeit machen, einen Beweis für (~) zu liefern, aber nach meinem besten Wissen ist das nur Kalkül und beinhaltet keine Wahrscheinlichkeitstheorie, obwohl ich mich noch nicht hinsetzen und es noch versuchen muss.
Lassen Sie uns zunächst eine Reihe von Trivialitäten durchgehen, um das Problem so umzuformulieren, dass es leichter zu lösen ist (beachten Sie, dass per Definition ):Ξn≥0
Ξn=o(logn)⟺limn→∞Ξnlogn=0⟺∀ε>0,Ξnlogn>ε only finitely many times⟺∀ε>0,Ξn>εlogn only finitely many times.
Man hat auch das:
Ξn>εlogn⟺n(1−F(Zn))>εlogn⟺1−F(Zn)>εlognn⟺F(Zn)<1−εlognn⟺Zn≤inf{y:F(y)≥1−εlognn}.
Definieren Sie entsprechend für alle :n
u(ε)n=inf{y:F(y)≥1−εlognn}.
Daher zeigen uns die obigen Schritte, dass:
Ξn=o(logn) a.s.⟺P(Ξn=o(logn))=1⟺P(∀ε>0,Ξn>εlogn only finitely many times)=1⟺P(∀ε>0,Zn≤u(ε)n only finitely many times)=1⟺P(∀ε>0,Zn≤u(ε)n infinitely often)=0.
Beachten Sie, dass wir schreiben können:
{∀ε>0,Zn≤u(ε)n infinitely often}=⋂ε>0{Zn≤u(ε)n infinitely often}.
Die Sequenzen werden mit abnehmendem gleichmäßig größer , sodass wir schließen können, dass die Ereignisse abnehmen (oder zumindest irgendwie monoton), da auf geht . Das Wahrscheinlichkeitsaxiom bezüglich monotoner Ereignissequenzen lässt daher den Schluss zu:u(ε)nε{Zn≤u(ε)n infinitely often}
ε0
P(∀ε>0,Zn≤u(ε)n infinitely often)=P(⋂ε>0{Zn≤u(ε)n infinitely often})=P(limε↓0{Zn≤u(ε)n infinitely often})=limε↓0P(Zn≤u(ε)n infinitely often).
Daher genügt es zu zeigen, dass für alle ,ε>0
P(Zn≤u(ε)n infinitely often)=0
denn natürlich ist die Grenze jeder konstanten Folge die Konstante.
Hier ist so etwas wie ein Vorschlaghammer-Ergebnis:
Satz 4.3.1., P. 252 von Galambos, The Asymptotic Theory of Extreme Order Statistics , 2. Auflage. Sei iid Variablen mit gemeinsamer nicht entarteter und kontinuierlicher Verteilungsfunktion , und sei eine nicht abnehmende Folge, so dass ebenfalls nicht abnehmend ist. Dann wird für ,
nach wie
X1,X2,…F(x)unn(1−F(un))un<sup{x:F(x)<1}P(Zn≤un infinitely often)=0 or 1
∑j=1+∞[1−F(uj)]exp(−j[1−F(uj)])<+∞ or =+∞.
Der Beweis ist technisch und dauert ungefähr fünf Seiten, aber letztendlich stellt sich heraus, dass er eine Folge einer der Borel-Cantelli-Deckspelzen ist. Ich werde vielleicht versuchen, den Beweis zu verdichten, um nur den für diese Analyse erforderlichen Teil sowie nur die Annahmen zu verwenden, die im Gaußschen Fall gelten, der kürzer sein kann (aber vielleicht nicht), und tippe ihn hier ein. Es wird jedoch nicht empfohlen, den Atem anzuhalten. Man beachte , dass in diesem Fall , so daß die Bedingung vacuous ist, und ist somit eindeutig nicht abnimmt.ω(F)=+∞n(1−F(n))εlogn
Der Punkt ist jedenfalls, dass wir diesen Satz ansprechen, wenn wir das zeigen können:
∑j=1+∞[1−F(u(ε)j)]exp(−j[1−F(u(ε)j)])=∑j=1+∞[εlogjj]exp(−εlogj)=ε∑j=1+∞logjj1+ε<+∞.
Da das logarithmische Wachstum für jeden Exponenten des positiven Potenzgesetzes langsamer ist als jedes Potenzgesetzwachstum (Logarithmen und Exponentiale bleiben monoton), ist und die frühere Ungleichung gelten immer für alle groß genug sind, da und eine Änderung der Variablen) Folgendes haben:loglogn≤αlogn⟺logn≤nαnlogn≤n
∑j=1+∞logjj1+ε≤∑j=1+∞jε/2j1+ε=∑j=1+∞1j1+ε/2<+∞,
da bekannt ist, dass die p-Reihe für alle konvergiert und natürlich impliziert .p>1ε>01+ε/2>1
So unter Verwendung der obigen Satz haben wir gezeigt , dass für alle , , was zu rekapitulieren bedeuten soll , dass fast sicher.ε>0P(Zn≤u(ε)n i.o.)=0Ξn=o(logn)
Wir müssen noch zeigen, dass . Dies folgt nicht aus dem Obigen, da z.logΞn=o(loglogn)
1nlogn=o(logn),−logn+loglogn≠o(logn).
Wenn man jedoch bei einer gegebenen Folge zeigen kann, dass für beliebiges , folgt daraus, dass . Idealerweise möchte ich dies für mit dem obigen Lemma zeigen können (vorausgesetzt, es ist sogar wahr), kann es aber (noch) nicht.xnxn=o((logn)δ)δ>0log(xn)=o(loglogn)Ξn