Paradox des Poisson-Prozesses mit mindestens einem Ereignis im Intervall

7

Lassen eine Anzahl von Ereignissen in Poisson Prozess der unitären Rate ( ) innerhalb Intervall der Länge T . Es ist bekannt, dass mindestens ein Ereignis in dem Intervall beobachtet wurde. Ich möchte die Wahrscheinlichkeit ermitteln, dass sich mehr Ereignisse in dem Intervall befinden.XTλ=1T

Meine Intuition ist, dass Pr(XT>1XT>0)=Pr(XT>0) .

Das Grundprinzip dahinter ist das

  1. Wenn das beobachtete Ereignis zum Zeitpunkt t vom Beginn des Intervalls war, reicht es aus, die Wahrscheinlichkeit zu berechnen, dass weder in (0,t) noch in (t,T) offenen Intervallen ein Ereignis aufgetreten ist : Pr(XT=1XT>0)=Pr(Xt=0)Pr(XTt=0)=etetT=eT=Pr(XT=0) ,

  2. Pr(XT>1XT>0)=1Pr(XT=1XT>0)Pr(XT>0)=1Pr(XT=0).

jedoch

Pr(XT>1XT>0)=Pr(XT>1,XT>0)Pr(XT>0)=Pr(XT>1)Pr(XT>0)=1Pr(XT{1,0})1Pr(XT=0)=1TeTeT1eT,

was weder ich noch WolframAlpha als gleich \ Pr (X_T> 0) = 1 - e ^ {- T} beweisen können Pr(XT>0)=1eT.

Da beide Ergebnisse nicht wahr sein können - wo ist mein Fehler?

Ich kann sehen, dass XT>1 und XT>0 stark abhängig sind. Ist das wichtig? Meine Intuition ist, dass XT>0 nur den Abtastraum einschränkt ...

[EDIT # 1]

Ich habe einen weiteren Weg gefunden, um ... beide Ergebnisse zu unterstützen.

Wenn die Zeit des ersten Ereignisses im Intervall ist (ab dem Beginn des Intervalls), wird seine Verteilungsdichte alsDann isttpdf(t)=et1eT.

Pr(XT=1XT>0)=0Tpdf(t)Pr(XTt=0)dt=0TetetT1eTdt=TeT1eT.

Wenn ich jedoch ähnliche Schritte für gleichmäßig verteilte ( ) Zufallsereignisse im Intervall wiederhole und auch Ereignisse vor berücksichtige, ich immer noch pdf(t)=1Tt

Pr(XT=1XT>0)=0Tpdf(t)Pr(Xt=0)Pr(XTt=0)dt=0TeTTdt=eT.

[EDIT # 2]

Follow-up aufgrund des Kommentars von @combo (über den Verlust der Konditionierung im ersten Ansatz).

Ich verstehe nicht, warum die Konditionierung verloren geht.

Stellen Sie sich eine Situation vor, in der wir ein Intervall der Länge mit mindestens einem Ereignis eines einheitlichen Poisson-Prozesses erstellen . Sei ein zufälliges Ereignis eines einheitlichen Poisson-Prozesses und ist eine Zufallsvariable, die gleichmäßig in . Dann ist ein Intervall der Länge das mindestens ein Ereignis enthält, bei (gleichmäßig verteilt) vom Beginn des Intervalls. Aufgrund der Unabhängigkeit der Ereignisse ist die Wahrscheinlichkeit, dass sich keine Ereignisse mehr im Intervall befinden, , nicht wahr? Und es wird vorausgesetzt, dass es in dem Intervall mindestens ein Ereignis gab.TYt(0, T)(Y  t, Y  t + T)TtPr(Xt==0)Pr(XTt==0)

Warum ist die Situation anders, wenn ich ein Intervall der Länge das mindestens ein Ereignis enthält? Die Zeit eines zufällig ausgewählten Ereignisses ( ; ab Beginn des Intervalls) ist gleichmäßig verteilt, daher sehe ich keinen Unterschied.Tt

abukaj
quelle
Warum geben Sie an, dass ? Ich verstehe nicht, warum dies der Fall sein sollte, und dies ist die Quelle Ihrer Probleme. Wenn ich es für mich selbst ausarbeite, bekomme ich die gleiche Antwort wie Ihr Ansatz Nr. 2. P(XT=1XT>0)=P(Xt=0)P(XTt=0)
Combo
@StefanJorgensen Ich verstehe deinen Standpunkt. Sollte es jedoch nicht ? P(Xt=1)P(XTt=0)=tete(Tt)=teT
Abukaj
@StefanJorgensen beantwortet Ihren zweiten Kommentar: Ich habe das Ereignis als einen Punkt angesehen, der das Intervall in zwei Teilintervalle mit unbekannter Anzahl von Ereignissen unterteilt.
Abukaj
@StefanJorgensen Es ist wie die Beantwortung einer Frage: Wenn zum Zeitpunkt ein Ereignis aufgetreten ist, wie hoch ist die Wahrscheinlichkeit, dass weder in noch in ein Ereignis aufgetreten ist . t(0,t)(t,T)
Abukaj
Außer dass die Gleichung einfach die Wahrscheinlichkeit ist, dass weder in noch in Ereignis vorliegt, und Ihre Konditionierung überhaupt nicht ausdrückt . Aus diesem Grund ist Ihre Antwort einfach die Wahrscheinlichkeit, dass es mindestens ein Ereignis gibt (weil die Konditionierung zu Beginn verloren gegangen ist), obwohl die Konditionierung tatsächlich das Ergebnis ändert. P(Xt=0)P(XTt=0)(0,t)(t,T)
Combo

Antworten:

2

Ich habe es endlich herausgefunden!

Nach dem Rat von @ combo werde ich den Begriff "Vorkommen" verwenden.

Überraschenderweise war der erste Teil der Begründung für meine Intuition fast richtig.

Wenn das beobachtete Auftreten zum Zeitpunkt vom Beginn des Intervalls war, reicht es aus, die Wahrscheinlichkeit zu berechnen, dass weder in noch in offenen Intervallen ein Auftreten aufgetreten ist : .t(0,t)(t,T)Pr(XT=1t)=Pr(Xt=0)Pr(XTt=0)=etetT=eT=Pr(XT=0)

Der Unterschied besteht darin, dass durch , was als kann. Dabei ist eine Länge des unendlich kleinen Intervalls . Da , können wir als Dichte im Punkt des einzigen sehen Auftreten innerhalb des Intervalls .Pr(XT=1XT>0)Pr(XT=1t)Pr(XT=1Xdt=1)dt[tdt2,t+dt2]Pr(XT=1,t)=eTdtpdf(t)=Pr(XT=1t)t(0,T)

Da der Wert von unbekannt ist, integrieren wir . So weit so gut - wir kennen die bedingungslose Wahrscheinlichkeit, genau ein Vorkommen im Intervall zu haben. Durch Konditionieren auf der Anteil des Abtastraums verworfen - daher ist es notwendig, jede verbleibende Wahrscheinlichkeit / Dichte um einen Faktor von zu normalisieren .t0Tpdf(t)dt=TeT=Pr(XT=1)XT>0eT11eT

Die Intuition war also der Irrtum, Wahrscheinlichkeit mit Dichte zu verwechseln.

abukaj
quelle
1

Ein kurzer Hinweis zur Terminologie, um uns nicht zu verwirren: Ich werde etwas, das zu einem bestimmten Zeitpunkt passiert, eher als "Ereignis" als als "Ereignis" bezeichnen, um Verwechslungen mit der strengeren Definition eines Ereignisses als Element der Stichprobe zu vermeiden Platz.

Beginnen wir mit der Definition des Poisson-Zählprozesses. Sei die Anzahl der Vorkommen, die zum Zeitpunkt . Es hat die folgenden Eigenschaften:N(t)T

  1. N(0)=0
  2. N(T1),(N(T2)N(T1)) sind unabhängig für (unabhängige Inkrementeigenschaft)T1<T2
  3. Die Anzahl der Vorkommen in einem Intervall der Länge ist eine Poisson-Zufallsvariable mit dem Parameter (für unsere Zwecke ).tλtλ=1

Die Eigenschaft der unabhängigen Inkremente ist das, was Sie auslöst - insbesondere die Auswirkungen der strengen Ungleichung.

Wir stellen die Frage, wie hoch die Wahrscheinlichkeit ist, dass ist , wenn ein Ereignis zum Zeitpunkt auftritt. Lassen Sie uns den Prozess nach Ihrem Ansatz in drei Segmente unterteilen. Für einige haben wir: und jetzt betrachten wir die Wahrscheinlichkeit t[0,T]N(T)=1ϵ<min{t,Tt}

N(T)=(N(T)N(t+ϵ))+(N(t+ϵ)N(tϵ))+(N(tϵ)N(0))
limϵ0P(N(T)=1N(t+ϵ)N(tϵ)=1)=limϵ0P(N(tϵ)N(0)=0,N(T)N(t+ϵ)=0N(t+ϵ)N(tϵ)=1)
Beachten Sie, dass dies nicht ganz das ist, was Sie geschrieben haben. Es gibt zwei wesentliche Unterschiede:
  1. Die Intervalle sind nicht disjunkt, so dass unabhängig von , und auch nicht unabhängig von .N(tϵ)N(0)N(T)N(t+ϵ)N(t+ϵ)N(tϵ)
  2. Die Intervalle haben alle ein positives Maß. Bei Ihrem Ansatz teilen Sie das Intervall in drei Teile auf: . Das Problem besteht darin, dass Sie jetzt auf das Ereignis konditionieren, dass ein Ereignis im Intervall auftritt mit dem Maß Null (weitere Informationen dazu, warum dies ein Problem ist, finden Sie in diesem Beitrag ).[0,T][0,t),[t],(t,T][t]
Combo
quelle
1. Ich bin etwas verwirrt über und im letzten Absatz. Meinen Sie damit, dass es eine Intervallüberlappung von size ? 2. Es mag ein Effekt der Verwirrung sein, aber ich sehe keinen Unterschied zwischen und vorausgesetzt, das Auftreten bei ist gegeben: . ϵsϵslims0P(N(ts)N(0)=0,N(T)N(t+s)=0N(t+s)N(ts)=1)lims0P(N(ts)N(0)=0,N(T)N(t+s)=0)tP(N(t+s)N(ts)=1)=1
Abukaj
Tut mir leid, ich wollte, dass sie gleich sind, aber mein Computer ist während der Bearbeitung gestorben. Der Unterschied ist genau die Quelle Ihrer Verwirrung: Das Ereignis ist nicht unabhängig vom Ereignis für . Wenn Sie versuchen, zu machen (mit Gleichheit, nicht mit dem Limit), ist das Intervall ein Singleton und Sie konditionieren eine Menge mit dem Maß Null. {N(tϵ)N(0)=0,N(T)N(t+ϵ)=0}{N(t+ϵ)N(tϵ)=1}ϵ>0ϵ=0
Combo
1
Nun, ich sehe die Abhängigkeit nicht. Es gelang mir jedoch, den Irrtum zu finden (siehe meine Antwort) - und da das Intervall wichtig war, also +1. [tϵ,t+ϵ]
Abukaj
0

Von Ihren beiden Herangehensweisen an das Problem scheint die zweite genau richtig zu sein. Es ist viel strenger und macht für mich vollkommen Sinn. Ich hatte jedoch etwas größere Schwierigkeiten, den ersten Ansatz zu verstehen, was Anlass zu der Annahme gab, dass dies die Quelle Ihres Fehlers ist.

Warum berechnen Sie aus Neugier zunächst ? Angesichts des gleichen Ansatzes, den Sie unten verwendet haben, sollte diese Wahrscheinlichkeit (soweit relevant) gleich was per Definition nicht gleich . Hoffe das hilft!P(XT=1|XT>0)TeT1eTeT

Felix van Doorn
quelle
Ich kann meinen Fehler im ersten Ansatz immer noch nicht sehen. Ich berechne es, da - siehe die Bearbeitung. P(XT>1|XT>0)=1P(XT=1|XT>0)
Abukaj