Differenz der Gamma-Zufallsvariablen

10

Wie ist die Verteilung der Differenz bei zwei unabhängigen Zufallsvariablen und Y G a m m a ( α Y , β Y ) , dh D = X - Y ?XGamma(αX,βX)YGamma(αY,βY)D=XY

Wie würde ich das Ergebnis ableiten, wenn das Ergebnis nicht bekannt ist?

FBC
quelle
Ich denke, kann relevant sein: stats.stackexchange.com/q/2035/7071
Dimitriy V. Masterov
4
Leider nicht relevant, berücksichtigt dieser Beitrag die gewichtete Summe von Gamma-Zufallsvariablen, bei denen die Gewichte streng positiv sind. In meinem Fall wären die Gewichte +1 bzw. -1.
FBC
Das Moschopoulos-Papier behauptet, dass die Methode auf lineare Kombinationen erweitert werden kann, aber Sie haben Recht, dass die Neuskalierung auf Gewichte größer als 0 beschränkt zu sein scheint. Ich stehe korrigiert da.
Dimitriy V. Masterov
Es gibt wenig Hoffnung, etwas Einfaches oder Geschlossenes abzuleiten, wenn die beiden Skalierungsfaktoren nicht gleich sind.
whuber
3
Nur eine kleine Bemerkung: Für den Sonderfall exponentiell verteilter Wohnmobile mit demselben Parameter lautet das Ergebnis Laplace ( en.wikipedia.org/wiki/Laplace_distribution ).
Ric

Antworten:

19

Ich werde skizzieren, wie das Problem angegangen werden kann, und angeben, was meiner Meinung nach das Endergebnis für den Sonderfall sein wird, wenn die Formparameter Ganzzahlen sind, aber die Details nicht ausfüllen.

  • Beachten Sie zunächst, dass Werte in ( - , ) annimmt und daher f X - Y ( z ) Unterstützung ( - , ) hat .XY(,)fXY(z)(,)

  • Zweitens von der Standard - Ergebnisse , dass die Dichte der Summe zweier unabhängiger kontinuierlichen Zufallsvariablen ist die Faltung ihrer Dichten, das heißt, und dass die Dichte der Zufallsvariablen - Y ist f - Y ( α ) = f Y ( - α ) , ableitendass f X - Y ( z ) = F X + ( - Y ) ( Z ) = - f X ( x ) f - Y ( z - x )

    fX+Y(z)=fX(x)fY(zx)dx
    YfY(α)=fY(α)
    fXY(z)=fX+(Y)(z)=fX(x)fY(zx)dx=fX(x)fY(xz)dx.
  • Drittens ist für nicht negative Zufallsvariablen und Y zu beachten, dass der obige Ausdruck zu f X - Y ( z ) = { 0 f X ( x ) f Y ( x - z ) vereinfacht wird. XY

    fXY(z)={0fX(x)fY(xz)dx,z<0,0fX(y+z)fY(y)dy,z>0.
  • Schließlich bedeutet die Parametrisierung eine Zufallsvariable mit der Dichte λ ( λ x ) s - 1Γ(s,λ)λ(λx)s1Γ(s)exp(λx)1x>0(x)XΓ(s,λ)YΓ(t,μ)z>0

    fXY(z)=0λ(λ(y+z))s1Γ(s)exp(λ(y+z))μ(μy)t1Γ(t)exp(μy)dy(1)=exp(λz)0p(y,z)exp((λ+μ)y)dy.
    z<0
    fXY(z)=0λ(λx)s1Γ(s)exp(λx)μ(μ(xz))t1Γ(t)exp(μ(xz))dx(2)=exp(μz)0q(x,z)exp((λ+μ)x)dx.

s=t

0xs1(x+β)s1exp(νx)dx
βfXY(z)

stp(y,z)yz(s+t2,s1)q(x,z)xz(s+t2,t1)

  • z>0(1)sy1,z,z2,zs1XYΓ(1,λ),Γ(2,λ),,Γ(s,λ)z>0t

  • z<0XYΓ(1,μ),Γ(2,μ),,Γ(t,μ)(μ|z|)k1exp(μz)(μz)k1exp(μz)s

Dilip Sarwate
quelle
2
+1: Nachdem ich mir dieses Problem schon einmal angesehen habe, finde ich diese Antwort faszinierend.
Neil G
Ich werde diese Antwort akzeptieren, obwohl es anscheinend keine geschlossene Lösung gibt. Es ist so nah wie es nur geht, danke!
FBC
fY(α)fY(α)
fY(α)=fY(α) P{Y>0}=1Y01
1
YfY(α)fY(α)α<0fY(α)0α<0fY(α)=fY(α)=0αYYfYR+
7

Meines Wissens wurde die Verteilung der Differenz zweier unabhängiger Gamma-RVs erstmals 1993 von Mathai untersucht. Er leitete eine Lösung in geschlossener Form ab. Ich werde seine Arbeit hier nicht reproduzieren. Stattdessen werde ich Sie auf die Originalquelle verweisen. Die Lösung für geschlossene Formen finden Sie auf Seite 241 als Satz 2.1 in seiner Arbeit über nicht zentrale verallgemeinerte Laplace-Eigenschaften quadratischer Formen in normalen Variablen .

Nathan Crock
quelle