CDF zur Macht erhoben?


14

Wenn FZ eine CDF ist, sieht es so aus, als wäre ( ) ebenfalls eine CDF.FZ(z)αα>0

F: Ist das ein Standardergebnis?

F: Gibt es einen guten Weg , um eine Funktion zu finden mit st , wobeiX g ( Z )gXg(Z)FX(x)=FZ(z)αxg(z)

Grundsätzlich habe ich eine andere CDF in der Hand, . In einer reduzierten Form möchte ich die Zufallsvariable charakterisieren, die diese CDF erzeugt.FZ(z)α

EDIT: Ich würde mich freuen, wenn ich ein Analyseergebnis für den Sonderfall . Oder zumindest wissen, dass ein solches Ergebnis nicht umsetzbar ist.ZN(0,1)


2
Ja, das ist ein ziemlich bekanntes Ergebnis und lässt sich leicht verallgemeinern. (Wie?) Sie können g auch zumindest implizit finden. Es handelt sich im Wesentlichen um eine Anwendung der inversen Transformationstechnik, die üblicherweise zum Erzeugen von Zufallsvariablen einer beliebigen Verteilung verwendet wird.
Kardinal

2
@ Kardinal Bitte antworten Sie. Das Team beschwert sich später darüber, dass wir nicht mit einer niedrigen Beantwortungsquote kämpfen.

1
@mbq: Danke für deine Kommentare, die ich sehr verstehe und respektiere. Bitte haben Sie Verständnis dafür, dass ich manchmal aus Zeit- und / oder Ortsgründen keine Antwort veröffentlichen kann, aber einen kurzen Kommentar zulassen kann, der das OP oder andere Teilnehmer an den Start bringt. Seien Sie versichert, dass ich dies in Zukunft tun werde, wenn ich eine Antwort posten kann. Ich hoffe, dass meine weitere Teilnahme durch Kommentare auch in Ordnung ist.
Kardinal

2
@ Kardinal Einige von uns sind aus den gleichen Gründen auch schuldig ...
whuber

2
@brianjd Ja, das ist ein bekanntes Ergebnis , das zur industriellen Herstellung von „verallgemeinerten“ Verteilungen verwendet wurde, sehen . Es gibt viele Transformationen wie diese, und die Leute benutzen sie für diesen Zweck: Sie finden eine parametrische Transformation, wenden sie auf eine Distribution an und schon haben Sie eine Arbeit, indem Sie nur ihre Eigenschaften berechnen. Und natürlich ist der Normalfall das erste "Opfer".

Antworten:


10

Die anderen Antworten gefallen mir, aber die folgenden hat noch niemand erwähnt. Das Ereignis {Ut, Vt} 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=FXnXnZ

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=(FZ1/n)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 istFZαα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.α=2ZF(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 ) ) .

g1(y)=Φ1(Φα(y)),
ZN(0,1)
x=g1(y)=F1(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)g
y=g(x)=F1(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)=(1ex), x>0,
F1(y)=ln(1y).
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(1ex)1/α)
XExp(1)
Y=ln(1(1eX)1/α),
Y Wir können dies direkt beweisen (schauen Sie sichP(Yy) 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)=(1ey)α.
P(Yy)
If XF then U=F(X)Unif(0,1).
If UUnif(0,1) then U1/αBeta(α,1).

Die Darstellung der Simulationsergebnisse folgt.

ECDF and F to the alpha

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)?


See my comment in @JMS answer. For ZN(0,1) it looks like the answer is g(z)=Φ1(Φ1/α(z)) which isn't closed form but can be calculated easily enough. And you can make it easier by recognizing that the input to the inverse CDF is a suitably chosen Beta distribution. The answer will be nice in cases where the inverse CDF is nice, and there are some of those running around.

It would be good to double-check your arithmetic.
cardinal

@cardinal errr... OK, I did,... and it's right? Would you please point out the error?

(+1) Apologies. I'm not sure where my head was when I first looked at this. It's obviously (well, should have been!) correct.
cardinal

@cardinal, no harm, no foul. I admit, though, you really had me sweating for a minute! :-)

14

Proof without words

enter image description here

The lower blue curve is F, the upper red curve is Fα (typifying the case α<1), and the arrows show how to go from z to x=g(z).


Nice pic! Q: What was that drawn in? TikZ?
lowndrul

@brianjd: If I recall, @whuber does many of his plots using Mathematica.
cardinal

3
@cardinal You're right. Actually, I use whatever is handy and seems like it will do a good job quickly. FWIW, here's the code: Module[ {y, w, a = 0.1, z = 3.24, f = ChiDistribution[7.6], xmin=0, xmax=5}, y = CDF[f,z]; w = InverseCDF[f, y^(1/a)]; Show[ Plot[{CDF[f, x],CDF[f,x]^a} , {x, xmin, xmax}, Filling->{1->{2}}], Graphics[{ Dashed, Arrow[{{z,0}, {z,y}}], Arrow[{{z,y}, {w,y}}], Arrow[{{w,y}, {w,0}}] }] ] ]
whuber

6

Q1) Yes. It's also useful for generating variables which are stochastically ordered; you can see this from @whuber's pretty picture :). α>1 swaps the stochastic order.

That it's a valid cdf is just a matter of verifying the requisite conditions: Fz(z)α has to be cadlag, nondecreasing and limit to 1 at infinity and 0 at negative infinity. Fz has these properties so these are all easy to show.

Q2) Seems like it would be pretty difficult analytically, unless FZ is special


@JMS: What about ZN(0,1) ?
lowndrul

2
@brianjd: I don't believe so. Let g be a continuous strictly monotonic function (hence, having a well-defined inverse g1) that satisfies your conditions. Then, it must be that Φα(u)=P(g(Z)u)=P(Zg1(u))=Φ(g1(u)) and so g1(u)=Φ1(Φα(u)). So the inverse is identified fairly explicitly, but not g itself. This is what I meant in my previous comment about g being found implicitly.
cardinal

@brianjd - What @cardinal said :) I couldn't even think of a special case for FZ where you'd get a closed form (not to say there isn't one of course).
JMS

@JMS: ZU[0,1] would be one positive example.
cardinal

@cardinal I never would have thought of such a rare distribution... but now that you mention it a Beta(a,1) should work in general, giving you back a Beta(aα,1).
JMS
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.