Bedingter erwarteter Wert des Produkts der Normal- und Log-Normalverteilung


7

Könnte jemand bitte die Antwort und die Schritte zur Lösung dieses Ausdrucks geben?

E[(eXY+k)|(eXY+k)>0]
E ist der Erwartungsoperator.
XN(μX,σX2);YN(μY,σY2);Xand Y are independent. Also, k<0

KEY MISSING LINK

Der obige Ausdruck hängt davon ab, ob der folgende allgemeine Ausdruck bewiesen wird oder zumindest gezeigt wird, dass er für den Sonderfall der Normal- und Log-Normalverteilung gilt.

E[UVUV>c]=E[UE[VUV>c]]

Hier sind unabhängige Zufallsvariablen, die diskret oder kontinuierlich sein können und jeder Wahrscheinlichkeitsverteilung folgen. ist der Erwartungsoperator. ist eine Konstante.U,VEc

Natürlich benötigen wir für unsere Hauptfrage nur den Fall, dass einer von ihnen normal verteilt ist und der andere logarithmisch normal ist? Angenommen, ist logarithmisch normal und ist normal. Trifft die obige Identität für diesen Sonderfall zu?UV

Dies wird aufgrund seiner Bedeutung als separate Frage gestellt: Beweis der Vereinfachung der bedingten Erwartung des Produkts zufälliger Variablen

SCHRITTE VERSUCHT

0) GEMEINSAME ZUSTANDSDICHTE

Ich habe Schwierigkeiten, die bedingte Gelenkdichte zu finden, die in der obigen Erwartung verwendet werden kann. Die Gelenkdichtefunktion von nur und ist einfach und folgt aus der Standarddichtefunktion für den bivariaten Normalfall. Wie würden wir den bedingten Aspekt in die Gelenkdichtefunktion einbeziehen?XY

1) NORMALES LOG-NORMALES MISCHPAPIER VON YANG

(Link: http://repec.org/esAUSM04/up.21034.1077779387.pdf )

Dieses Papier hat die ersten vier zentralen Momente ohne Beweis (Gleichung 5 im obigen Papier). Wenn jemand diese Beweise liefern könnte, könnte dies mehr Licht auf das oben genannte Problem werfen.

Die Variablen im Yang-Papier sind korreliert, was oben leicht anzuwenden ist. Sie haben aber auch einen Mittelwert von Null in der Arbeit, was in unserem Fall nicht direkt gilt, da wir einen Mittelwert ungleich Null haben.

2) ANDERE VERWANDTE LINKS

a) Interessante Frage zu einer Erwartung, die eine leicht modifizierte Form des normalen logarithmisch-normalen Gemisches beinhaltet. Diesem fehlt jedoch der bedingte Aspekt, und daher müssen einige Änderungen vorgenommen werden, bevor er für das obige Problem verwendet werden kann.

/math/1142841/covariance-in-normal-lognormal-nln-mixture

b) Eine weitere Frage zur normalen logarithmisch-normalen Mischung, obwohl dies keine tiefere Diskussion enthält.

/math/159818/combination-of-a-normal-rv-with-a-log-normal-one

c) Frage zur bedingten Erwartung des Produkts unabhängiger Zufallsvariablen. Es wäre gut zu wissen, welche Aspekte davon in unserem Fall zutreffen.

/math/544410/result-and-proof-on-the-conditional-expectation-of-the-product-of-two-random-var

d) Weitere interessante Fragen zur bedingten Erwartung unabhängiger Zufallsvariablen.

/math/380866/conditional-expectations-for-independent-random-variables?rq=1

/math/55524/rule-with-independent-random-variables-and-conditional-expectations?rq=1

3) Annäherungen an die Taylor-Serie

Wäre es möglich, hier ernsthafte Näherungswerte zu verwenden? Ich bin wenig verwirrt aufgrund der bedingten Erwartung und der normalen logarithmischen Normalmischung? Hinweise darauf, ob dies möglich ist und wie es weitergehen soll oder ob dies hier nicht anwendbar ist, wären großartig.

4) VERWENDUNG VON STANDARD NORMAL (SCHEINT WIE EIN TOTES ENDE)

Ich weiß, dass es eine Lösung gibt, wenn wir diese Summe mit der folgenden Standardnormalen ausdrücken können. Bitte beraten Sie auch, wie dies oder andere Alternativen zur Lösung des oben genannten Problems zu tun sind. Dies scheint eine Sackgasse zu sein, wie Experten in diesem Forum bestätigt haben. Aber immer noch hier bleiben, wenn jemand einen Weg findet, diesen Ansatz weiter zu verwenden.

W=(eXY+k)<=>μ+σZ where, ZN(0,1)

[WN(μ,σ2)W=μ+σZ;W>0Z>μ/σ]
Wir müssen dann bestimmen .μ and σ

Wir haben für jede Standardnormalverteilung und für jedes . Hier sind und das normale Standard-PDF bzw. CDF. Daher haben wir, Zu, Pr[Z>-u]=Pr[Z<u]=Φ(u)ϕΦ

E[Z|Z>u]=1Φ(u)[utϕ(t)dt]=1Φ(u)[ϕ(t)|u]=ϕ(u)Φ(u)
E[Y|Y>0]=μ+σE[Z|Z>(μσ)]=μ+σϕ(μ/σ)Φ(μ/σ)
Einstellung, , ψ(u)=u+ϕ(u)/Φ(u)
E[Y|Y>0]=σψ(μ/σ)

@gung, Bitte beachten Sie über den Schritten, die ich versucht habe. Das oben Gesagte als Standardnormalvariable mit geeignetem Mittelwert und Varianz auszudrücken, hat eine Lösung. Der Trick wäre also, das Obige in ein geeignetes Standardnormal umzuwandeln. Wenn es Alternativen gibt, schlagen Sie bitte auch vor. Danke für Ihre Hilfe.
Texmex

Ich habe gerade auch die Mischung Paper on Normal und Log Normal von Yang gefunden, aber sie behandelt korrelierte Variablen und bin sich noch nicht sicher, wie ich die bedingte Erwartung erhalten soll. Bitte beraten Sie über geschlossene Lösungen, die für die oben genannten existieren könnten.
Texmex

Was ist der Index in Ihrer ersten Formel? T
user765195

Entschuldigung, ich hätte das entweder entfernen oder früher klarstellen sollen, dass es Zeit ist. Ich nehme die Erwartung zum Zeitpunkt T an. Bitte lassen Sie mich wissen, wenn Sie weitere Details benötigen. Anregungen zu Antworten und Schritten wären willkommen?
Texmex

2
Jede Referenz, die den Fehler analysiert, würde dies tun. Zunächst empfehle ich Ihnen jedoch, die Frage zu vereinfachen, indem Sie redundante Parameter eliminieren. Sie kann um auf eins reduziert werden, wobei eine Mittelwert-Null-Normalverteilung und eine beliebige Normalverteilung hat. Diese Vereinfachung wird die Art des Problems verdeutlichen und Ihnen möglicherweise dabei helfen, festzustellen, ob eine Taylor-Reihen-Näherung angemessen wäre. eXYXY
whuber

Antworten:


2

Was ist die beabsichtigte Verwendung des Ergebnisses? Dies hängt davon ab, welche Art von Antwort erforderlich ist, um zu berücksichtigen, ob ein stochastischer (Monte-Carlo-) Simulationsansatz angemessen sein könnte. Und selbst das Gesamtbild ist, dass dieses Problem gelöst werden muss, und hat sich jemand dieses Problem ausgedacht ein Problem auf höherer Ebene zu lösen, und es könnte einen besseren Ansatz für das Problem auf höherer Ebene geben, der dies nicht erfordert.

Hier ist eine stochastische (Monte Carlo) Simulationslösung in MATLAB.

a = 1; b = 2; c = 3; d = 4; k = -1; % Made up values for illustrative purpose
n = 1e8; % Number of replications
mux = 10; sigmax = 4; sigmay = 7; % Made up values for illustrative purposes
X = mux + sigmax * randn(n,1); Y = sigmay * randn(n,1); Y1 = a + b + c + d * Y;
success_index = exp(X).*Y1 > 0; % replications in which condition is true
num_success = sum(success_index);
Cond_Sample = exp(X(success_index)) .* Y1(success_index) + k;
disp([num_success mean(Cond_Sample) std(Cond_Sample)/sqrt(num_success)])
1.0e+09 *
0.058475265000000   1.502775087443930   0.057342191058931

Ich bin der Meinung, dass dies eine so grundlegende Kombination aus Normal und Log-Normal ist, dass es eine Lösung geben muss, wenn keine Annäherung verfügbar ist. Wenn dies nicht der Fall ist, wären angesichts der eher einfachen Natur des Problems weitere Fortschritte oder Durchbrüche für viele hilfreich, wie es scheint. Ich würde mich freuen, wenn ich so schnell wie möglich ein Kopfgeld dafür bekomme, nachdem ich fünf weitere Punkte erhalten habe.
Texmex

Zu diesem Zeitpunkt habe ich wohl genug Werbegeschenke für einen Milliardärs-Hedgefonds-Tycoon gegeben. Sie können Pony spielen, wenn sie spezialisiertere Hilfe und Beratung wünschen. Oder ist das ein akademisches Projekt?
Mark L. Stone

Dies ist ein akademisches Projekt. Ich bin kein Milliardär Hedge Fund Tycoon :-) Wenn dies nicht akademisch wäre, hätte eine Simulation es längst getötet. Eine bessere Antwort zu finden, wäre hauptsächlich im Interesse, Wissen zu erlangen und es zu verbreiten… Richtig?
Texmex

Vielleicht sollten wir zusammenarbeiten. Sie können eine Einladung an mein LinkedIn-Konto senden, das in meinem Profil verlinkt ist. Und ich dachte nicht, dass Sie der Milliardärs-Hedgefonds-Tycoon sind - eher, dass Sie ein Angestellter sind, der für den Milliardärs-Hedgefonds-Tycoon arbeitet :).
Mark L. Stone

Da ich das große Ganze immer noch nicht wirklich kenne, kann ich nicht beurteilen, welche Art von Lösung, wenn überhaupt, wirklich dafür benötigt wird. Und warum ist eine dieser Normalen oder Lognormalen und nicht einige dickere Schwanzverteilungen? Der Simulationsansatz bleibt bestehen, wenn sich Verteilungen ändern, Abhängigkeiten eingeführt werden usw.
Mark L. Stone

2

Bemerkungen:

  1. Die Verbindungsdichte wird durch Multiplizieren der Dichten angegeben, da sie idp sind. Eine Variable ist nur ein Parameter für die andere.
  2. YeX ist nicht normal verteilt, daher funktioniert Ansatz (4) nicht.

  3. Mit den folgenden Ausdrücken können Sie möglicherweise eine Annäherung finden. Wenn nicht, sind sie mit einem Computer relativ einfach zu bewerten.

Sei , bezeichne und . Sei , dann und wenn und andernfalls. Da also nur ist, wenn , ist XN(μ,σ2)c=k/yd(y)=log(k/y)Z=eX

E[YZB]=E[YE[ZB]],
E[ZZ>c]=eμ+σ2/2P(X>dσ2)P(X>d),
c>0E[Z]=eμ+σ2/2c>0Y<0
E[YZB]=12eμ+σ2/2(E[Y1(Y>0]+E[YP(X>d(Y)σ2)P(X>d(Y))1(Y<0)]).

Der erste Teil ist einfach und der zweite

0yfY(y)dy,
0yP(X>log(ky)σ2)P(X>log(ky))fY(y)dy.

Danke für die hilfreiche Antwort. Hier ist , aber ich nehme an, wir ändern nur das Vorzeichen und alles folgt. Könnten Sie bitte auch die Schritte näher erläutern, um die folgenden beiden Ergebnisse zu erhalten? Folgt der erste einfach aus der Unabhängigkeit oder handelt es sich um einen anderen besonderen Zweck? # 1) # 2)k<0
E[YZB]=E[YE[ZB]],
E[ZZ>c]=eμ+σ2/2P(X>dσ2)P(X>d),
texmex

Könnten Sie bitte auch erklären, wie Sie hier die bedingte Dichte einrichten?
E[ZYZYB]=BzyfZYBdzdy=y(k/yzfZdz)fYdy
texmex

Es gab einen Fehler im zweidimensionalen Integral. Wenn Sie es sich jetzt ansehen, sollte klar sein, woher stammt - ansonsten lassen Sie es mich wissen! Die bedingte Erwartung ist ein Standardergebnis, das beispielsweise im Wikipedia-Artikel für lognormal zu finden ist. Sie finden es, indem Sie e ^ X f_X mit einigen Einschränkungen integrieren. E[YZB]
Hunaphu

Sehr dankbar für Ihre Klarstellungen. Ich sehe, wie wir das bekommen würden, aber ich habe diesen Link gefunden, als ich versucht habe, dies besser zu verstehen, und vielleicht verwirre ich etwas, würde aber die Antwort unter diesem Link gilt in unserem Fall nicht. math.stackexchange.com/questions/544410/…
E[YZB]=E[YE[ZB]],
texmex

Wir habenE[YE[ZB]E[Y]E[ZB]
Hunaphu

0

Diese Lösung ist auf die Vorschläge und Korrekturen von @Hunaphu, @whuber und anderen zurückzuführen. Könnte jemand bitte überprüfen, ob alle Schritte sinnvoll sind?

ANTWORT SCHRITTE STARTEN

Verwenden einiger Notationsverknüpfungen,

Betrachten Sie: Hier ist die Wahrscheinlichkeitsdichtefunktion für , hier ist die Wahrscheinlichkeitsdichtefunktion für ,

E[(eXY+k)|(eXY+k)>0]=E[k|(eXY+k)>0]+E[(eXY)|(eXY+k)>0]
=k+E[(YeX)|(YeX+k)>0]
=k+yexf(yex|{yex+k}>0)dxdy
f(w)wf(w)w
=k+yexf(yex;{yex+k}>0)f({yex+k}>0)dxdy
[We note that, yex>k>0y>0]
=k+yexf(y)f(ex;{yex+k}>0)f({yex+k}>0)dxdy
=k+y[exf(ex;{ex>ky})f(ex>ky)dx]f(y)dy
=k+y[exf(ex|{ex>ky})dx]f(y)dy
=k+0(y<k)y[exf(ex|{ex>1})dx]f(y)dy+(y>k)y[exf(ex|{ex<1})dx]f(y)dy
=k+0(k)y[E(W|W>c)]f(y)dy+(k)y[E(W|W<c)]f(y)dy;here, W=eX and c=1
Vereinfachung der inneren Erwartungen,
E(W|W>c)=1P(eX>c)cw1wσX2πe12[ln(w)μXσX]2dw
Setzen , haben wir,t=ln(w)dw=etdt
E(W|W>c)=1P(X>ln(c))ln(c)etσX2πe12(tμXσX)2dt
t12(tμXσX)2=12σX2(t(μX+σX2))2+μX+σX22
E(W|W>c)=e(μX+12σX2)P(μX+σXZ>ln(c))ln(c)1σX2πe12[t(μX+σX2)σX]2dt;ZN(0,1)
Put und wir haben,s=[t(μX+σX2)σX]b=[ln(c)(μX+σX2)σX]ds=dtσX
E(W|W>c)=e(μX+12σX2)P(Z>ln(c)μXσX)b12πe12s2ds
=e(μX+12σX2)P(Z<ln(c)+μXσX)[12πe12s2dsb12πe12s2ds]
=e(μX+12σX2)P(Z<ln(c)+μXσX)[1Φ(b)];Φ is the standard normal CDF
=e(μX+12σX2)Φ(ln(c)+μXσX)[Φ(b)]
In ähnlicher Weise gilt für den anderen Fall Setzen , wir haben,
E(W|W<c)=1P(eX<c)0cw1wσX2πe12[ln(w)μXσX]2dw
t=ln(w)dw=etdt
E(W|W<c)=1P(X<ln(c))ln(c)etσX2πe12(tμXσX)2dt
t12(tμXσX)2=12σX2(t(μX+σX2))2+μX+σX22
E(W|W<c)=e(μX+12σX2)P(μX+σXZ<ln(c))ln(c)1σX2πe12[t(μX+σX2)σX]2dt;ZN(0,1)
Put und wir,s=[t(μX+σX2)σX]b=[ln(c)(μX+σX2)σX]ds=dtσX
E(W|W<c)=e(μX+12σX2)P(Z<ln(c)μXσX)b12πe12s2ds
=e(μX+12σX2)P(Z<ln(c)μXσX)[Φ(b)];Φ is the standard normal CDF
=e(μX+12σX2)Φ(ln(c)μXσX)[Φ(b)]
Verwenden Sie die Ergebnisse für die inneren Erwartungen.
E[(eXY+k)|(eXY+k)>0]=k+0(k)y[e(μX+12σX2)Φ(ln(c)+μXσX)[Φ(b)]]f(y)dy+(k)y[e(μX+12σX2)Φ(ln(c)μXσX)[Φ(b)]]f(y)dy
=k+e(μX+12σX2)[0(k)y{Φ(μX+σX2σX)Φ(μXσX)}f(y)dy+(k)y{Φ([μX+σX2σX])Φ([μXσX])}f(y)dy]
=k+e(μX+12σX2)[0(k)y{Φ(μX+σX2σX)Φ(μXσX)}f(y)dy+(k)y{1Φ(μX+σX2σX)1Φ(μXσX)}f(y)dy]
=k+e(μX+12σX2)[{Φ(μX+σX2σX)Φ(μXσX)}μYσY(k+μYσY)(μY+σYz)12πe12z2dz+{1Φ(μX+σX2σX)1Φ(μXσX)}(k+μYσY)(μY+σYz)12πe12z2dz];ZN(0,1)
=k+e(μX+12σX2)[{Φ(μX+σX2σX)Φ(μXσX)}{μY[Φ([k+μYσY])Φ(μYσY)]σY2π[e12(k+μYσY)2e12(μYσY)2]}+{1Φ(μX+σX2σX)1Φ(μXσX)}{μY[1Φ([k+μYσY])]+σY2π[e12(k+μYσY)2]}]

Ja, das oder einfach eine stochastische Simulation. Wer sind diese Leute, die im akademischen Urteil sitzen und keine simulationsbasierten Lösungen zulassen? Es gibt Professoren an den renommiertesten Universitäten der Welt, Zeitschriften und Bücher, die sich nur mit stochastischer Simulation befassen. Hinweis: Nennen wir es stochastische Simulation, nicht Monte-Carlo-Simulation - es klingt raffinierter und weniger "vulgär".
Mark L. Stone

@ MarkL.Stone Du bist lustig und hast Recht mit der stochastischen Simulation. Wie ich bereits erwähnt habe, besteht die Frage zu diesem Zeitpunkt darin, keinen Professor oder Hedgefonds-Manager glücklich zu machen. Es ist Wissen um seiner selbst willen ... :-)
Texmex
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.