Jumlah dari dua variabel acak gamma independen

13

Menurut artikel Wikipedia tentang distribusi Gamma :

Jika XGamma(a,θ) dan YGamma(b,θ) , di mana X dan Y adalah variabel acak independen, maka X+YGamma(a+b,θ) .

Tapi saya tidak melihat bukti. Adakah yang bisa mengarahkan saya ke buktinya?

Sunting: Terima kasih banyak kepada Zen, dan saya juga menemukan jawabannya sebagai contoh di halaman Wikipedia tentang fungsi karakteristik .

Dexter12
sumber
3
Intuisi: Gamma distribusi timbul sebagai jumlah dari n distribusi eksponensial independen, mana itu segera dalam konteks ini bahwa X + Y akan memiliki Gamma ( a + b , θ ) distribusi disediakan sebuah dan b keduanya bilangan bulat positif. (n)nX+Y(a+b,θ)ab
whuber

Jawaban:

15

Buktinya adalah sebagai berikut: (1) Ingatlah bahwa fungsi karakteristik dari jumlah variabel acak independen adalah produk dari fungsi karakteristik individu mereka; (2) Dapatkan fungsi karakteristik dari variabel acak gamma di sini ; (3) Lakukan aljabar sederhana.

Untuk mendapatkan intuisi di luar argumen aljabar ini, periksa komentar whuber.

Catatan: OP bertanya bagaimana cara menghitung fungsi karakteristik dari variabel acak gamma. Jika , maka (Anda dapat memperlakukan saya sebagai konstanta biasa, dalam hal ini)XExp(λ)i

ψX(t)=E[eitX]=0eitxλeλxdx=11it/λ.

Sekarang gunakan tip Huber: Jika , maka Y = X 1 + + X k , di mana X i adalah independen E x p ( λ = 1 / θ ) . Oleh karena itu, menggunakan properti (1), kita memiliki ψ Y ( t ) = ( 1YGamma(k,θ)Y=X1++XkXiExp(λ=1/θ)

ψY(t)=(11itθ)k.

Kiat: Anda tidak akan mempelajari hal-hal ini menatap hasil dan bukti: tetap lapar, hitung semuanya, cobalah untuk menemukan bukti Anda sendiri. Bahkan jika Anda gagal, penghargaan Anda terhadap jawaban orang lain akan berada pada level yang jauh lebih tinggi. Dan, ya, gagal itu OK: tidak ada yang melihat! Satu-satunya cara untuk belajar matematika adalah dengan bertarung tinju untuk setiap konsep dan hasil.

Zen
sumber
Pernyataan yang direferensikan secara eksplisit menyatakan "asalkan semua Xi independen."
whuber
Satu hal yang saya tidak mengerti, adalah bagaimana kita sampai pada fungsi karakteristik?
Dexter12
Saya akan menambahkannya ke jawabannya. Lihatlah.
Zen
Mungkin Anda bisa menyertakan referensi untuk fungsi karakteristik a untuk nilai non-integer dari a ? Γ(a,θ)a
Dilip Sarwate
14

Berikut adalah jawaban yang tidak perlu menggunakan fungsi karakteristik, tetapi sebaliknya memperkuat beberapa ide yang memiliki kegunaan lain dalam statistik. Kepadatan jumlah variabel acak independen adalah konvolusi kepadatan. Jadi, dengan mengambil untuk kemudahan eksposisi, kita memiliki untuk z > 0 , f X + Y ( z )θ=1z>0

fX+Y(z)=0zfX(x)fY(zx)dx=0zxa1exΓ(a)(zx)b1e(zx)Γ(b)dx=ez0zxa1(zx)b1Γ(a)Γ(b)dxnow substitute x=zt and think=ezza+b101ta1(1t)b1Γ(a)Γ(b)dtof Beta(a,b) random variables=ezza+b1Γ(a+b)
Dilip Sarwate
sumber
3
(+1) It is ideal to have more than one way to prove everything. Maybe someone will post an answer considering the transformation (X,Y)(U,V)=(X+Y,X).
Zen
Can we similarly find the density of XY in a closed form expression? I'm unable to simplify the integrals in that case.
pikachuchameleon
@pikachuchameleon See this answer of mine.
Dilip Sarwate
3

On a more heuristic level: If a and b are integers, the Gamma distribution is an Erlang distribution, and so X and Y describe the waiting times for respectively a and b occurrences in a Poisson process with rate θ. The two waiting times X and Y are

  1. independent
  2. sum up to a waiting time for a+b occurrences

and the waiting time for a+b occurrences is distributed Gamma(a+b,θ).

None of this is a mathematical proof, but it puts some flesh on the bones of the connection, and can be used if you want to flesh it out in a mathematical proof.

Svein Olav Nyberg
sumber