Die anderen Antworten gefallen mir, aber die folgenden hat noch niemand erwähnt. Das Ereignis {U≤t, V≤t} tritt genau dann auf, wenn {max(U,V)≤t} , wenn also U und V unabhängig sind und W=max(U,V) , dann FW(t)=FU(t)∗FV(t) so fürα eine positive ganze Zahl (beispielsweiseα=n ) nehmenX=max(Z1,...Zn) , wo dieZ ‚s sind iid
Für wir auf F Z = F n X umschalten , also wäre X diese Zufallsvariable, so dass das Maximum von n unabhängigen Kopien die gleiche Verteilung wie Z hat (und dies wäre keiner unserer bekannten Freunde) , allgemein). α=1/nFZ=FnXXnZ
Der Fall von einer positiven rationalen Zahl (sagen wir α = m / n ) folgt aus dem vorhergehenden, da
( F Z ) m / n = ( F 1 / n Z ) m .αα=m/n
(FZ)m/n=(F1/nZ)m.
Wählen Sie für irrational eine Folge positiver Rationalen a k, die gegen α konvergieren ; dann konvergiert die Sequenz X k (wo wir unsere obigen Tricks für jedes k verwenden können ) in der Verteilung zu dem gewünschten X.αakαXkkX
Dies ist möglicherweise nicht die gesuchte Charakterisierung, gibt aber zumindest eine Vorstellung davon, wie man über für α angemessen schön denkt . Andererseits bin ich mir nicht sicher, wie viel schöner es wirklich werden kann: Sie haben bereits die CDF, also gibt Ihnen die Kettenregel das PDF und Sie können Momente berechnen, bis die Sonne untergeht ...? Es ist wahr, dass die meisten Z kein X haben , das für α = √ bekannt istFαZαZX , aber wenn ich mit einem Beispiel herumspielen wollte, um nach etwas Interessantem zu suchen, könnte ich versuchen,Zgleichmäßig verteilt auf das Einheitsintervall mitF(z)=z,0<z<1.α=2–√ZF(z)=z0<z<1
EDIT: Ich schrieb einige Kommentare in @JMS answer, und es gab eine Frage zu meiner Arithmetik, also schreibe ich auf, was ich damit gemeint habe, in der Hoffnung, dass es klarer ist.
@ Kardinal richtig im Kommentar zu @JMS Antwort schrieb, dass das Problem zu vereinfacht
,
oder allgemeiner, wenn Z nicht unbedingt N ( 0 , 1 ) ist , wir hat
x = g - 1 ( y ) = F - 1 ( F α ( y ) ) .
g−1(y)=Φ−1(Φα(y)),
ZN(0,1)x=g−1(y)=F−1(Fα(y)).
Mein Punkt war, dass wir , wenn
eine schöne Umkehrfunktion hat, mit der Basisalgebra nur nach der Funktion
y = g ( x ) auflösen können. Ich schrieb in der Bemerkung , dass
g sein sollte
y = g ( x ) = F - 1 ( F 1 / α ( x ) ) .Fy=g(x)gy=g(x)=F−1(F1/α(x)).
Nehmen wir einen speziellen Fall, schließen Sie die Dinge an und sehen Sie, wie es funktioniert. Lasse haben eine Exp (1) -Verteilung, mit CDF
F ( x ) = ( 1 - e - x ) , x > 0 ,
und inversen CDF
F - 1 ( y ) = - ln ( 1 - y ) .
Es ist einfach, alles anzuschließen, um g zu finden . nachdem wir fertig sind, erhalten wir
y = g ( x ) = -X
F(x)=(1−e−x), x>0,
F−1(y)=−ln(1−y).
g
Zusammenfassend lautet meine Behauptung also, dass wenn
X ∼ E x p ( 1 ) und wenn wir
Y = - ln ( 1 - ( 1 - e - X) definieren ) 1 / α ) ,
dann hat
Y eine CDF, die wie folgt aussieht:
F Y ( y ) = (y=g(x)=−ln(1−(1−e−x)1/α)
X∼Exp(1)Y=−ln(1−(1−e−X)1/α),
Y
Wir können dies direkt beweisen (schauen Sie sich
P(Y≤y) anund verwenden Sie die Algebra, um den Ausdruck zu erhalten, im vorletzten Schritt benötigen wir die Wahrscheinlichkeitsintegraltransformation). Nur in dem (oft wiederholten) Fall, dass ich verrückt bin, habe ich einige Simulationen durchgeführt, um zu überprüfen, ob es funktioniert, ... und es funktioniert. Siehe unten. Um den Code zu vereinfachen, habe ich zwei Fakten verwendet:
Wenn X ∼ F, dann ist U = F ( X ) ∼ U n i f ( 0 , 1 )FY(y)=(1−e−y)α.
P(Y≤y)If X∼F then U=F(X)∼Unif(0,1).
If U∼Unif(0,1) then U1/α∼Beta(α,1).
Die Darstellung der Simulationsergebnisse folgt.
Der R-Code, der zum Generieren des Plots verwendet wird (Minuszeichen), lautet
n <- 10000; alpha <- 0.7
z <- rbeta(n, shape1 = alpha, shape2 = 1)
y <- -log(1 - z)
plot(ecdf(y))
f <- function(x) (pexp(x, rate = 1))^alpha
curve(f, add = TRUE, lty = 2, lwd = 2)
Die Passform sieht ziemlich gut aus, finde ich? Vielleicht bin ich nicht verrückt (diesmal)?