ความแตกต่างของตัวแปรสุ่มแกมมา


10

รับตัวแปรสุ่มอิสระสองตัวและการกระจายความแตกต่างคือคืออะไร?Y ∼ G a m m a ( α Y , β Y ) D = X - YX∼Gamma(αX,βX)Y∼Gamma(αY,βY)D=X−Y

หากผลลัพธ์ไม่เป็นที่รู้จักฉันจะไปหาผลลัพธ์ได้อย่างไร


ฉันคิดว่าอาจมีความเกี่ยวข้อง: stats.stackexchange.com/q/2035/7071
— Dimitriy V. Masterov

4
น่าเสียดายที่ไม่เกี่ยวข้องกันโพสต์นั้นพิจารณาผลรวมถ่วงน้ำหนักของตัวแปรสุ่มแกมมาซึ่งน้ำหนักเป็นบวกอย่างเคร่งครัด ในกรณีของฉันน้ำหนักจะเป็น +1 และ -1 ตามลำดับ
— FBC

กระดาษ Moschopoulos อ้างว่าวิธีการนั้นสามารถขยายไปยังการผสมแบบเส้นตรง แต่คุณพูดถูกว่าการลดอัตราการรอดชีวิตดูเหมือนจะถูก จำกัด น้ำหนักที่มากกว่า 0 ฉันยืนแก้ไข
— Dimitriy V. Masterov

มีความหวังเล็กน้อยว่าจะได้รับสิ่งใด ๆ ที่เรียบง่ายหรืออยู่ในรูปแบบปิดเว้นแต่ว่าปัจจัยสองระดับจะเหมือนกัน
— whuber

3
เพียงคำพูดเล็ก ๆ : สำหรับกรณีพิเศษของการแจกแจง rvs ด้วยพารามิเตอร์เดียวกันผลลัพธ์คือ Laplace ( en.wikipedia.org/wiki/Laplace_distribution )
— Ric

คำตอบ:


19

ฉันจะสรุปว่าปัญหาสามารถแก้ไขได้อย่างไรและระบุสิ่งที่ฉันคิดว่าผลลัพธ์จะเป็นกรณีพิเศษเมื่อพารามิเตอร์รูปร่างเป็นจำนวนเต็ม แต่ไม่ได้กรอกรายละเอียด

  • ครั้งแรกที่ทราบว่าจะใช้เวลาในค่าใน( - ∞ , ∞ ) และอื่น ๆฉX - Y ( Z )มีการสนับสนุน( - ∞ , ∞ )X−Y(−∞,∞)fX−Y(z)(−∞,∞)

  • ประการที่สองจากผลมาตรฐานที่มีความหนาแน่นของผลรวมของสองตัวแปรสุ่มอิสระอย่างต่อเนื่องเป็นบิดความหนาแน่นของพวกเขา, ที่อยู่, และความหนาแน่นของตัวแปรสุ่ม - Yคือ ฉ- Y ( α ) = F Y ( - α ) , อนุมานว่า ฉX - Y ( Z ) = F X + ( - Y ) ( Z ) = ∫ ∞ - ∞ f X ( x ) f - Y ( z - x )

    fX+Y(z)=∫−∞∞fX(x)fY(z−x)dx
    −Yf−Y(α)=fY(−α)
    fX−Y(z)=fX+(−Y)(z)=∫−∞∞fX(x)f−Y(z−x)dx=∫−∞∞fX(x)fY(x−z)dx.
  • ประการที่สามสำหรับที่ไม่ใช่เชิงลบตัวแปรสุ่มและYทราบว่าการแสดงออกดังกล่าวข้างต้นจะช่วยลดความยุ่งยาก ฉX - Y ( Z ) = { ∫ ∞ 0ฉX ( x ) ฉY ( x - Z )XY

    fX−Y(z)={∫0∞fX(x)fY(x−z)dx,z<0,∫0∞fX(y+z)fY(y)dy,z>0.
  • ในที่สุดใช้ parametrization เพื่อหมายถึงตัวแปรสุ่มที่มีความหนาแน่น λ ( λ x ) s - 1Γ(s,λ)λ(λx)s−1Γ(s)exp⁡(−λx)1x>0(x)X∼Γ(s,λ)Y∼Γ(t,μ)z>0

    fX−Y(z)=∫0∞λ(λ(y+z))s−1Γ(s)exp⁡(−λ(y+z))μ(μy)t−1Γ(t)exp⁡(−μy)dy(1)=exp⁡(−λz)∫0∞p(y,z)exp⁡(−(λ+μ)y)dy.
    z<0
    fX−Y(z)=∫0∞λ(λx)s−1Γ(s)exp⁡(−λx)μ(μ(x−z))t−1Γ(t)exp⁡(−μ(x−z))dx(2)=exp⁡(μz)∫0∞q(x,z)exp⁡(−(λ+μ)x)dx.

s=t

∫0∞xs−1(x+β)s−1exp⁡(−νx)dx
βfX−Y(z)

stp(y,z)yz(s+t−2,s−1)q(x,z)xz(s+t−2,t−1)

  • z>0(1)sy1,z,z2,…zs−1X−YΓ(1,λ),Γ(2,λ),⋯,Γ(s,λ)z>0t

  • z<0X−YΓ(1,μ),Γ(2,μ),⋯,Γ(t,μ)(μ|z|)k−1exp⁡(μz)(μz)k−1exp⁡(−μz)s


2
+1: เมื่อพิจารณาปัญหานี้มาก่อนฉันพบว่าคำตอบนี้น่าสนใจ
— Neil G

ฉันจะยอมรับคำตอบนี้แม้ว่าจะไม่มีโซลูชันแบบปิด มันใกล้เคียงกับที่ได้รับขอบคุณ!
— FBC

f−Y(α)≠fY(−α)

f−Y(α)=fY(−α) P{Y>0}=1−Y01

1
YfY(α)fY(α)α<0fY(α)0α<0f−Y(α)=fY(α)=0αYY−−fYR+

7

จากความรู้ของฉันการกระจายความแตกต่างของสองแกมมาอิสระได้ถูกศึกษาครั้งแรกโดย Mathai ในปี 1993 เขาได้รับโซลูชั่นแบบปิด ฉันจะไม่ทำซ้ำงานของเขาที่นี่ แต่ฉันจะชี้ให้คุณไปที่แหล่งต้นฉบับ วิธีการแก้ปัญหาแบบปิดสามารถพบได้บนหน้า 241 เป็นทฤษฎีบท 2.1 ในกระดาษของเขาที่ไม่กลางทั่วไป Laplacianness รูปแบบสมการกำลังสองตัวแปรปกติ

โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.