Situasi:
Katakanlah saya memiliki proses Poisson, seperti peluruhan radioaktif, menghasilkan partikel R per detik. Saya mengukur dengan detektor. Ada probabilitas P bahwa partikel akan dideteksi oleh detektor.
Hal-hal yang saya rasa saya tahu:
- Waktu antar kedatangan emisi partikel secara eksponensial didistribusikan dengan parameter berdasarkan R .
- Jumlah partikel yang dipancarkan sebelum deteksi diberikan oleh binomial negatif berdasarkan P .
- Jika sejumlah N sampel dari (2), satu sampel waktu antar kedatangan untuk partikel yang terdeteksi dapat diberikan dengan jumlah sampel N dari (1). Jumlah ini dapat diperoleh dengan sampling dari distribusi gamma dengan parameter berdasarkan N dan R .
Pertanyaan saya:
Jika waktu antar kedatangan tunggal dapat dihitung dengan pengambilan sampel dari gamma berdasarkan N dan R , bagaimana jumlah detektor dalam interval akhirnya menjadi Poisson lagi? (Untuk menjadi Poisson, waktu antar kedatangan untuk detektor harus eksponensial, tidak didistribusikan menurut beberapa hal gamma aneh.) Tentu saja N berfluktuasi, tetapi saya tidak bisa melihat bagaimana ini bekerja.
Namun, saya hampir sepenuhnya yakin jumlah detektor sebenarnya didistribusikan Poisson. Bisakah seseorang menunjukkan matematika kepada saya? Terima kasih untuk bantuannya!
EDIT:
Saya menemukan makalah ini: Fried, DL "Noise in photoemission current." Optik Terapan 4.1 (1965): 79-80.
Yang menunjukkan hasil bahwa variabel acak poisson terpilih secara binomi juga Poisson dengan nilai yang diberikan oleh PR. Ini mengkonfirmasi komentar oleh jbowman. Namun, saya akan tertarik melihat penjelasan tentang bagaimana proses saya menghasilkan waktu antar kedatangan di detektor menggunakan distribusi binomial dan gamma negatif tidak benar. Ini adalah gangguan mental utama saya. Terima kasih.
EDIT 2:
Saya menulis skrip matlab ini untuk menguji apakah apa yang saya coba dengan distribusi gamma berhasil. Ternyata, entah bagaimana, waktu kedatangan antar gamma yang dihasilkan dengan N yang didistribusikan secara geometris adalah eksponensial dan sesuai dengan waktu antar kedatangan yang disarankan oleh Poisson (PR). (ia2 dan ia3 didistribusikan secara identik). Adakah yang tahu bagaimana ini bekerja secara analitik? Bagi saya itu tidak begitu jelas bagi saya!
close all
n = 100000;
ia1 = exprnd(1,n,1); % create exponentially distributed inter-arrival times
t1 = cumsum(ia1); % running sum (the real experiment time)
mask = (rand(n,1) > 0.5); % flip a coin
t2 = t1(mask); % get only the events for which "the coin landed on heads"
ia2 = diff(t2); % calculate the inter-arrival times at the detector.
% plot the distributions
figure; hist(ia1,100); title('exponential inter-arrival times');
figure; hist(ia2,100); title('binomial sampled inter-arrival times');
%%
spacing = geornd(0.5,n,1) + 1; % how many events before we get heads
ia3 = gamrnd(spacing,ones(n,1)); % generate the interarrival times with gamma
figure; hist(ia3,100); title('geom/gamma inter-arrival times');
sumber
Jawaban:
Argumen non-teknis yang cepat mungkin menggunakan jaringan Jackson . Dalam kasus Anda, total kedatangan eksternal adalah tarifR , dan tidak ada transisi internal (partikel yang diamati tidak beralih ke antrian yang tidak teramati). Proporsi pemisahan antara node yang diamati dan tidak teramatihal0 i adalah P , sehingga
Jika Anda mencari prinsip pertama, teleponlahO ( t ) proses penghitungan yang diamati, dan N( t ) ∼ P.P( r ) total proses penghitungan. Di mana setiap kedatanganN( t ) masuk O ( t ) dengan probabilitas hal . Jadi kalau untuk beberapas kita punya N( s ) = n kemudian O ( s ) memiliki binomial (n , p ) distribusi.
Pendekatan ini menggunakan fungsi-fungsi yang menghasilkan probabilitas:
Kesetaraan terakhir oleh teorema binomial. Kemudian, tanpa syarat, sejak ituN( t ) ∼ P.o i s s o n ( r t ) :
Yang merupakan fungsi menghasilkan probabilitas dari Poisson (r p t ) variabel acak.
sumber