Wenn und unabhängige Normalvariablen mit jeweils dem Mittelwert Null sind, ist ebenfalls eine Normalvariable

11

Ich versuche die Aussage zu beweisen:

Wenn und unabhängige Zufallsvariablen sind,X N ( 0 , σ 2 1 ) XN(0,σ21)Y N ( 0 , σ 2 2 )YN(0,σ22)

dann ist X Y.X 2 + Y 2XYX2+Y2 auch eine normale Zufallsvariable.

Für den Sonderfall σ 1 = σ 2 = σσ1=σ2=σ (sagen wir) haben wir das bekannte Ergebnis, dass X Y.X 2 + Y 2N(0,σ24 )XYX2+Y2N(0,σ24) wenn X.X und Y.Y unabhängige N ( 0 , σ 2 )N(0,σ2) -Variablen sind. Tatsächlich ist allgemeiner bekannt, dass X Y.X 2 + Y 2 ,X2-Y22 X 2 + Y 2XYX2+Y2,X2Y22X2+Y2 sind unabhängige N ( 0 , σ 24 )N(0,σ24) Variablen.

Ein Beweis des letzten Ergebnisses folgt unter Verwendung der Transformation ( X , Y ) ( R , Θ ) ( U , V )(X,Y)(R,Θ)(U,V) wobei x = r cos θ , y = r sin θx=rcosθ,y=rsinθ und u = r2 sin(2& thgr;),v = r2 cos(2&thgr;)u=r2sin(2θ),v=r2cos(2θ) . In der Tat ist hier U = X Y.X 2 + Y 2U=XYX2+Y2 und V = X 2 - Y. 22X 2 + Y 2V=X2Y22X2+Y2 . Ich habe versucht, diesen Beweis für das vorliegende Problem nachzuahmen, aber es scheint chaotisch zu werden.

Wenn ich keinen Fehler gemacht habe, erhalte ich für ( u , v ) R. 2(u,v)R2 die Gelenkdichte von ( U , V )(U,V) als

f U , V ( u , v ) = 2σ 1 σ 2 π exp[-u 2 + v 2 (u 2 + v 2 +vσ 2 1 +u 2 + v 2 -vσ22)]

fU,V(u,v)=2σ1σ2πexp[u2+v2(u2+v2+vσ21+u2+v2vσ22)]

Ich habe den Multiplikator oben, da die Transformation nicht eins zu eins ist.22

Die Dichte von würde also durch , was nicht ohne weiteres ausgewertet werden kann.UURfU,V(u,v)dvRfU,V(u,v)dv

Jetzt bin ich interessiert zu wissen, ob es einen Beweis gibt, bei dem ich nur mit und kein berücksichtigen muss, um zu zeigen, dass normal ist. Die CDF von sieht für mich im Moment nicht so vielversprechend aus. Dasselbe möchte ich auch für den Fall tun .UUVVUUUUσ1=σ2=σσ1=σ2=σ

Das heißt, wenn und unabhängige möchte ich zeigen, dass ohne Änderung von Variablen. Wenn ich irgendwie argumentieren kann, dass , dann bin ich fertig. Also hier zwei Fragen, der allgemeine Fall und dann der besondere Fall.XXYYN(0,σ2)N(0,σ2)Z=2XYX2+Y2N(0,σ2)Z=2XYX2+Y2N(0,σ2)Zd=XZ=dX

Verwandte Beiträge zu Math.SE:

X2Y2/X2+Y2N(0,1)X2Y2/X2+Y2N(0,1) wenn unabhängig voneinanderX,YN(0,1)X,YN(0,1) .

Wenn iid , zeigen Sie, dass sind iidX,YX,YN(0,1)N(0,1)XYX2+Y2,X2Y22X2+Y2XYX2+Y2,X2Y22X2+Y2N(0,14)N(0,14) .

Bearbeiten.

Dieses Problem ist in der Tat auf L. Shepp zurückzuführen, wie ich in den Übungen von Eine Einführung in die Wahrscheinlichkeitstheorie und ihre Anwendungen (Band II) von Feller zusammen mit einem möglichen Hinweis herausgefunden habe:

enter image description here

Sicherlich ist und ich habe die Dichte von zur Hand.U=XYX2+Y2=11X2+1Y2U=XYX2+Y2=11X2+1Y21X21X2

Mal sehen, was ich jetzt tun kann. Abgesehen davon ist auch eine kleine Hilfe mit dem oben genannten Integral willkommen.

HartnäckigAtom
quelle
1
Der MGF-Ansatz für das Gelenk ist zwar ähnlich, aber etwas einfacher. Siehe die letzte Antwort von: math.stackexchange.com/a/2665178/22064 und: math.stackexchange.com/questions/2664469/…(U,V)(U,V)
Alex R.
@AlexR. Ja, ich hatte den gemeinsamen mgf-Ansatz gesehen, der recht gut funktioniert, wenn ich die gemeinsame Verteilung für den Fall gleicher Varianz finden möchte. Aber ich habe in diesem Fall bereits den Beweis durch Änderung von Variablen, was meiner Meinung nach einfacher ist. Ich versuche, nur mit zu arbeiten , da dies die Verteilung ist, nach der ich strebe. UU
Hartnäckig
1
Der Trick besteht darin, dass die Summe von und , die inverse Chi-Quadrat-Verteilungen skaliert sind, auch eine skalierte inverse Chi-Quadrat-Verteilung ist (d. H. die Eigenschaft stabiler Verteilungen). Die Magie geschieht also in der dritten Gleichung der folgenden: 1X21X21Y21Y2U=XYX2+Y2=11X2+1Y2=11Z2=Z
U=XYX2+Y2=11X2+1Y2=11Z2=Z
Sextus Empiricus
@MartijnWeterings Anscheinend ist das der ursprüngliche Beweis von Shepp.
Hartnäckig am
Ich hätte mir das nicht ausgedacht, wenn Sie den Kommentar von Shepp nicht erwähnt hätten. Aber ich hatte die Idee, dass Sie diesen Beweis nicht bekommen haben. Zumindest war unklar, ob dies der Fall war.
Sextus Empiricus

Antworten:

6

Die ursprüngliche Lösung des Problems von Shepp verwendet das Konzept des stabilen Rechtseigentums, das mir im Moment etwas fortgeschritten erscheint. Daher konnte ich den Hinweis in der Übung, die ich in meinem Beitrag zitiert habe, nicht verstehen. Ich denke, ein Beweis, der nur die einzelne Variable und keine Änderung von Variablen verwendet, ist schwierig zu finden. Daher teile ich drei Open-Access-Dokumente, die eine alternative Lösung für das Problem darstellen:U=XYX2+Y2U=XYX2+Y2

Der erste hat mich überzeugt, den Integrationspfad, den ich mit dieser Wahl der Variablen , um die Dichte von abzuleiten, nicht zu beschreiten . Es ist das dritte Papier, das wie etwas aussieht, dem ich folgen kann. Ich gebe hier eine kurze Skizze des Beweises:VVUU

Wir nehmen ohne Verlust der Allgemeinheit an und setzen . wir nun feststellen, dass und unabhängig sind, haben wir die gemeinsame Dichte von . Wir bezeichnen es mit .σ21=1σ21=1σ22=σ2σ22=σ2X2χ21X2χ21Y2σ2χ21Y2σ2χ21(X2,Y2)(X2,Y2)fX2,Y2fX2,Y2

Betrachten Sie die Transformation so, dass und . Wir haben also die Gelenkdichte von . Bezeichnen wir es mit . Nach dem Standardverfahren integrieren wir wrt in , um die Randdichte von .(X2,Y2)(W,Z)(X2,Y2)(W,Z)W=X2Y2X2+Y2W=X2Y2X2+Y2Z=X2+Y2Y2Z=X2+Y2Y2(W,Z)(W,Z)fW,ZfW,ZfW,ZfW,ZzzfWfWWW

Wir finden, dass eine Gamma-Variable mit den Parametern und , so dass . Wir stellen fest, dass die Dichte von um symmetrisch ist . Dies impliziert, dass und damit .W=U2W=U212122(1+1σ)22(1+1σ)2(1+1σ)2Wχ21(1+1σ)2Wχ21UU00(1+1σ)UN(0,1)(1+1σ)UN(0,1)UN(0,(σσ+1)2)UN(0,(σσ+1)2)

HartnäckigAtom
quelle
0

demzufolge

Zwei normale Zufallsvariablen transformieren

X=rcos(θ)Y=rsin(θ)X,Ynormal(0,1)θUniform(0,2π)r2chi(2)X=rcos(θ)Y=rsin(θ)X,Ynormal(0,1)θUniform(0,2π)r2chi(2) . und sind unabhängig und sind unabhängig.
XXYY θθrr

auch dass da sin(θ)cos(θ)sin(2θ)2sin(θ)cos(θ)cos(2θ)cos(2θ)fsin(θ)cos(θ)sin(2θ)2sin(θ)cos(θ)cos(2θ)cos(2θ)ff(z)=1π(1z2)I[1,1](z)f(z)=1π(1z2)I[1,1](z)z=sin(θ)f(z)=|ddzsin1(z)|fθ(sin1(z))+|ddz(πsin1(z))|fθ(πsin1(z))=1(1z2)12π+1(1z2)12π=1π(1z2)z=sin(θ)f(z)=|ddzsin1(z)|fθ(sin1(z))+|ddz(πsin1(z))|fθ(πsin1(z))=1(1z2)12π+1(1z2)12π=1π(1z2)

ähnlich für andere.

2XY(X2+Y2)=2r2cos(θ)sin(θ)r=2rcos(θ)sin(θ)=rsin(2θ)rsin(θ)N(0,1)

so können wir zeigen:

X=σrcos(θ) undY=σrsin(θ)

so

2XY(X2+Y2)=2r2σσcos(θ)sin(θ)rσ=2σrcos(θ)sin(θ)=σrsin(2θ)σrsin(θ)σN(0,1)=N(0,σ2)

unabhängig zu zeigen

2XY(X2+Y2)=σrsin(θ)

X2Y22(X2+Y2)=r2σ2(cos2(θ)sin2(θ))2rσ=12rσ(cos2(θ)sin2(θ))12rσcos(2θ)12rσcos(θ) und leicht zu sagen, dass sie unabhängig sind.

masoud
quelle
Was ist, wenn ? σXσY
Sextus Empiricus
Ich habe nicht darüber nachgedacht. Aber einige Berechnungsprobleme sqrt(X2+Y2)
treten