Ist in einem mit einer gewissen Effizienz gemessenen Poisson-Prozess die gemessene Anzahl immer noch Poisson?

8

Situation:

Angenommen, ich habe einen Poisson-Prozess wie den radioaktiven Zerfall, bei dem R- Partikel pro Sekunde erzeugt werden. Ich messe mit einem Detektor. Es besteht eine Wahrscheinlichkeit P, dass ein Partikel vom Detektor erfasst wird.

Dinge, die ich zu wissen glaube:

  1. Die Zwischenankunftszeit der Partikelemission ist exponentiell mit Parametern verteilt, die auf R basieren .
  2. Die Anzahl der vor dem Nachweis emittierten Partikel wird durch ein negatives Binom auf Basis von P angegeben .
  3. Wenn eine Zahl N aus (2) abgetastet wird, kann eine einzelne Probe der Zwischenankunftszeit für detektierte Partikel durch die Summe von N Proben aus (1) gegeben werden. Diese Summe kann durch Abtasten aus einer Gammaverteilung mit Parametern erhalten werden, die auf N und R basieren .

Meine Frage:

Wenn eine einzelne Zwischenankunftszeit durch Abtasten von einem Gamma basierend auf N und R berechnet werden kann , wie wird dann die Anzahl der Detektorzählungen in einem Intervall wieder Poisson? (Um Poisson zu sein, muss die Zwischenankunftszeit für den Detektor exponentiell sein und darf nicht nach einer seltsamen Gammasache verteilt sein.) Natürlich schwankt N , aber ich kann nicht sehen, wie das funktioniert.

Ich bin mir jedoch fast sicher, dass die Detektorzahlen tatsächlich Poisson-verteilt sind. Könnte mir jemand die Mathematik zeigen? Danke für die Hilfe!

BEARBEITEN:

Ich fand dieses Papier: Fried, DL "Rauschen im Photoemissionsstrom." Applied Optics 4.1 (1965): 79 & ndash; 80.

Dies zeigt das Ergebnis, dass eine binomial ausgewählte Poisson-Zufallsvariable auch Poisson mit einer von PR gegebenen Rate ist. Dies bestätigt den Kommentar von jbowman. Trotzdem würde mich interessieren, wie mein Prozess der Erzeugung der Zwischenankunftszeit am Detektor unter Verwendung der negativen Binomial- und Gammaverteilung falsch ist. Dies ist mein größter geistiger Schluckauf. Vielen Dank.

EDIT 2:

Ich habe dieses Matlab-Skript geschrieben, um zu testen, ob das, was ich mit der Gamma-Verteilung versucht habe, funktioniert. Es stellt sich heraus, dass die mit einem geometrisch verteilten N erzeugten Gamma-Zwischenankunftszeiten exponentiell sind und mit den von Poisson (PR) vorgeschlagenen Zwischenankunftszeiten übereinstimmen. (ia2 und ia3 sind identisch verteilt). Irgendeine Idee, wie das analytisch funktioniert? Es war mir nicht intuitiv klar!

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');
user487100
quelle
PRPNNN
Ich verstehe nicht, woher Sie wissen, dass die Rate Poisson (RP) ist. Könntest du mir zeigen? Das ist das Herzstück dieser Frage, denke ich. In # 2 habe ich vermutlich, wenn ich eine P-Chance habe, den Detektor zu treffen, die Anzahl der Partikel, die vor dem Auftreffen auf den Detektor emittiert werden, ist geometrisch mit einem Mittelwert von 1 / P verteilt. Somit kann ich aus dieser geometrischen Verteilung eine Stichprobe berechnen, um N zu erhalten, und dann N Zwischenankunftszeiten summieren, um eine einzelne Zwischenankunftszeit am Detektor zu erhalten. Können Sie den Fehler in dieser Logik erklären? Ich denke, Ihre Aussage über die Rate als Poisson (RP) ist wichtig. Vielen Dank!
user487100
Kennen Sie sich mit Momenterzeugungs- / Charakteristikfunktionen aus? Ich würde es mit diesem Ansatz aufschreiben, da es einfach ist, es sei denn, es ist auch nicht hilfreich.
Jbowman
Nein, ich habe nicht mit momentgenerierenden Funktionen gearbeitet. Haben Sie eine Vorstellung davon, wie Sie zeigen können, dass Poisson + eine feste Akzeptanzwahrscheinlichkeit nur die Poisson-Rate skaliert? Ich bin bereit, den momentan generierenden Funktionsansatz zu lernen, wenn Sie zeigen könnten, wie dies funktioniert.
user487100
1
Ich fürchte, es wird heute viel später sein (Pacific Std Time). Ich kann es auch auf einfache Weise tun, was weniger undurchsichtig sein wird.
Jbowman

Antworten:

4

Rp0iP

λobs=RP

O(t)N(t)PP(r)N(t)O(t)psN(s)=nO(s)n,p) Verteilung.

Dieser Ansatz verwendet wahrscheinlichkeitsgenerierende Funktionen:

E[zO(t)|N(t)=n]=j=0nzj(nj)pj(1p)nj=(1p+pz)n

Letzte Gleichheit nach dem Binomialsatz. Dann bedingungslos, da :N(t)Poisson(rt)

E[zO(t)]=E[E[zO(t)|N(t)=n]]=n=0(1p+pz)nrtnn!ert=ertert(1p+pz)=erpt(z1)

Welches ist die Wahrscheinlichkeit erzeugende Funktion einer Poisson ( ) Zufallsvariablen.rpt

Vermutungen
quelle