Angenommen, . Show


18

Wie lässt sich am einfachsten feststellen, dass die folgende Aussage zutrifft?

Angenommen, . Zeige .Y1,,YniidExp(1)i=1n(YiY(1))Gamma(n1,1)

Beachten Sie, dass .Y(1)=min1inYi

Mit XExp(β) bedeutet dies, dass fX(x)=1βex/β1{x>0} .

Es ist leicht zu erkennen, dass Y(1)Exponential(1/n) . Darüber hinaus haben wir auch das i=1nYiGamma(α=n,β=1) unter der Parametrisierung

fY(y)=1Γ(α)βαxα1ex/β1{x>0}α,β>0.

Lösung in Xi'ans Antwort : Unter Verwendung der Notation in der ursprünglichen Frage:

i=1n[YiY(1)]=i=1n[Y(i)Y(1)]=i=1nY(i)nY(1)=i=1n{Y(i)Y(i1)+Y(i1)Y(1)+Y(1)}nY(1)=i=1nj=1i{Y(j)Y(j1)}nY(1) where Y(0)=0=j=1ni=jn{Y(j)Y(j1)}nY(1)=j=1n(nj+1)[Y(j)Y(j1)]nY(1)=i=1n(ni+1)[Y(i)Y(i1)]nY(1)=i=2n(ni+1)[Y(i)Y(i1)]+nY(1)nY(1)=i=2n(ni+1)[Y(i)Y(i1)].
Daraus erhalten wir i=2n(ni+1)[Y(i)Y(i1)]Gamma(n1,1)1) .

1
@MichaelChernick 1) Ich bin mir nicht sicher, ob mein Unabhängigkeitsnachweis korrekt ist, und 2) ich bin mir nicht sicher, ob das Ergebnis, das ich oben in Bezug auf den Unterschied der Gamma-Verteilungen vermutet habe, überhaupt korrekt ist. Dies scheint dem zu widersprechen, was hier angegeben ist , aber vielleicht ist diese Situation anders, da dieser Unterschied eine der Ordnungsstatistiken enthält? Ich bin mir nicht sicher.
Klarinettist

1
@Clarinetist, ich bin nicht sicher. Versuchen Sie vielleicht, mit , was eindeutig der Summe entspricht, mit der Sie arbeiten. Die Antwort hier könnte hilfreich sein: math.stackexchange.com/questions/80475/…i=2n(Y(i)Y(1))
Gammer

3
Haben Sie versucht zu beweisen, dass jedes - mit Ausnahme eines , für das und dann unter Verwendung der Tatsache, dass die Summe von iid Exponentialvariablen Gamma-verteilt sein wird? i Y i - Y ( 1 ) = 0 ( n - 1 )(YiY(1))Expon(1)iYiY(1)=0(n1)
Marcelo Ventura

1
Wir haben und bedingt durch , wir teilen dies durch Geben wir , so haben wir . Was mich an diesem Beweis stört: Ich betrachtete als Konstante. Aber ist keine Konstante. Warum sollte das funktionieren? Y ia e - a e - z i ( Z i | Y iA ) ~ Exp ( 1 ) A Y ( 1 )
fZi(zi)=fYi(zi+a)=e(zi+a)
Yiaeaezi(ZiYiA)Exp(1)AY(1)
Klarinettist

1
Der Punkt ist, dass es egal ist, was ist! Die Verteilung ist immer ! Bemerkenswert, nicht wahr? Und daraus können Sie schließen, dass die Verteilung von immer für , unabhängig vom tatsächlichen Wert von . Exp ( 1 ) Y i - Y [ 1 ] Exp ( 1 ) i > 1 Y [ 1 ]aExp(1)YiY[1]Exp(1)i>1Y[1]
Bogenschütze

Antworten:


15

Der Beweis ist in den Büchern der Mutter aller zufälligen Generationen, Devroye's Non-Uniform Random Variate Generation , auf S.211 (und es ist sehr elegant!) Zu finden:

Theorem 2.3 (Sukhatme, 1937) Wenn wir dann leiten sich die normierten Exponentialabstände von ab Ordnungsstatistik einer iid-Exponentialstichprobe der Größe sind selbst iid-Exponentialvariablen( n - i + 1 ) ( E ( i ) - E ( i - 1 ) ) E ( 1 )E ( n )E(0)=0

(ni+1)(E(i)E(i1))
E(1)E(n)n

Beweis. Da the gemeinsame Dichte der Ordnungsstatistik schreibt als Setzen von , die Änderung der Variablen von bis haben einen konstanten Jacobian [übrigens gleichdies muss aber nicht berechnet werden] und damit die Dichte von (E(1),,E(n

i=1nei=i=1ne(i)=i=1nj=1i(e(j)e(j1))=j=1ni=jn(e(j)e(j1))=j=1n(nj+1)(e(j)e(j1))
f( e )=n!(E(1),,E(n)) Y i = ( E ( i ) - E ( i - 1 ) ) ( E ( 1 ) , , E ( n ) ) ( Y 1
f(e)=n!exp{i=1ne(i)}=n!exp{i=1n(ni+1)(e(i)e(i1))}
Yi=(E(i)E(i1))(E(1),,E(n))1 / n ! ( Y 1 , , Y n ) exp { - n i = 1 y i }(Y1,,Yn)1/n!(Y1,,Yn) ist proportional zu was das Ergebnis festlegt. QED
exp{i=1nyi}

Eine Alternative, die mir Gérard Letac vorgeschlagen hat, besteht darin, zu überprüfen, ob dieselbe Verteilung hat wie (aufgrund der memorylosen Eigenschaft), wodurch die Herleitung von einfach.( E 1

(E(1),,E(n))
(E1n,E1n+E2n1,,E1n+E2n1++En1)
k=1n(EkE(1))k=1n1Ek

1
Danke für diese Antwort! Ich möchte für jeden, der dies in Zukunft liest, einige Details : sind beobachtete Werte von und der einfachste Weg, dies zu sehen. schreibt Begriff für Begriff aus. Da die Dichte von proportional zu , trennen Sie das zu sehen, dass die Dichte proportional zu , daher .eiEii=1ne(i)=i=1n(ni+1)(e(i)e(i1))=i=1ne(i)i=1n(ni+1)(e(i)e(i1))(Y1,,Yn)exp(i=1nyi)yii=1neyiY1,,YniidExp(1)
Klarinettist

5

Ich lege hier dar, was in den Kommentaren von @jbowman vorgeschlagen wurde.

Sei eine Konstante . Lassen folgen einem und betrachten . Danna0YiExp(1)Zi=Yia

Pr(ZiziYia)=Pr(YiaziYia)

Pr(Yizi+aYia)=Pr(Yizi+a,Yia)1Pr(Yia)

Pr(aYizi+a)1Pr(Yia)=1ezia1+eaea=1ezi

Dies ist die Verteilungsfunktion von .Exp(1)

Beschreiben wir dies: Die Wahrscheinlichkeit, dass ein rv in ein bestimmtes Intervall fällt (der Zähler in der letzten Zeile), hängt nur von Länge des Intervalls und nicht an der Stelle, an der sich dieses Intervall in der realen Linie befindet. Exp(1)Dies ist eine Inkarnation der Eigenschaft " memorylessness " der Exponentialverteilung, hier in einer allgemeineren Umgebung, frei von Zeitinterpretationen (und gilt für die Exponentialverteilung im Allgemeinen).

Nun wird durch Konditionierung auf wir zwingen nicht-negativ zu sein, und das ist entscheidend, das erhaltene Ergebnis hält . Wir können also folgendes sagen: {Yia}ZiaR+

Wenn , dann . YiExp(1)Q0:Zi=YiQ0 ZiExp(1)

Können wir ein , das frei ist, alle nicht negativen reellen Werte anzunehmen und für das die erforderliche Ungleichung immer gilt (fast sicher)? Wenn wir können, können wir auf das konditionierende Argument verzichten. Q0

Und in der Tat können wir. ist die Statistik minimaler Ordnung , , . So haben wir erhaltenQ=Y(1)Pr(YiY(1))=1

YiExp(1)YiY(1)Exp(1)

Das bedeutet, dass

Pr(YiY(1)yiy(1))=Pr(Yiyi)

Wenn also die Wahrscheinlichkeitsstruktur von unverändert bleibt, wenn wir die Statistik minimaler Ordnung subtrahieren, folgt, dass die Zufallsvariablen und wobei unabhängig sind, sind auch unabhängig, da die mögliche Verbindung zwischen ihnen, keinen Einfluss auf die Wahrscheinlichkeitsstruktur hat.YiZi=YiY(1)Zj=YjY(1)Yi,YjY(1)

Dann wird die Summe enthält uiv Zufallsvariablen (und eine Null), und soi=1n(YiY(1))n1 Exp(1)

i=1n(YiY(1))Gamma(n1,1)
Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.