間隔内に少なくとも1つのイベントがあるPoissonプロセスのParadox


7

LET(単一レートのポアソン過程におけるイベントの数であるの長さの間隔内で)。間隔内で少なくとも1つのイベントが観測されたことがわかっているので、間隔内にさらにイベントがある確率を見つけたいと思います。XTλ=1T

私の直感は、です。Pr(XT>1XT>0)=Pr(XT>0)

背後にある理論的根拠は

  1. 観測されたイベントが間隔の最初から時間tにあった場合、(0、t)または(t、T)の開いた間隔でtイベントが発生しなかった確率を計算するだけで十分です:\ Pr(X_T = 1 \ mid X_T> 0)= \ Pr(X_t = 0)\ Pr(X_ {Tt} = 0)= e ^ {-t} e ^ {t-T} = e ^ {-T} = \ Pr(X_T = 0 )(0,t)(t,T)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).

しかしながら

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,

IもWolframAlpha\ Pr(X_T> 0)= 1-e ^ {-T}に等しいことを証明できませんPr(XT>0)=1eT

両方の結果は真実ではあり得ないので-私の間違いはどこですか?

私がいることがわかりますXT>1XT>0大きく依存しています。それは重要ですか?私の直感は、XT>0がサンプリングスペースを狭めるだけであるということです...

[編集#1]

私はサポートするもう1つの方法を見つけました...両方の結果。

が間隔の最初のイベントの時間である場合(間隔の最初から)、その分布密度はとして与えられます。次に、tpdf(t)=et1eT.

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

ただし、同様の手順を均一に分布した()に対して繰り返し、前のイベントも考慮に入れると、まだ pdf(t)=1Tt

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

[編集#2]

@comboのコメントによるフォローアップ(最初のアプローチでの条件付けの喪失について)。

なぜコンディショニングが失われるのかわかりません。

長さ間隔を作成し、その中に単一のポアソンプロセスのイベントが少なくとも1つある状況を想像してみてください。ましょう単一ポアソン過程のランダムイベントであり、均一に分布したランダム変数である。次に、は、間隔の先頭から(均一に分散)で、少なくとも1つのイベントを含む長さの間隔です。イベントの独立性から、間隔内にイベントがなくなる確率はですよね。そして、インターバル中に少なくとも1つのイベントがあったと仮定します。TYt(0, T)(Y  t, Y  t + T)TtPr(Xt==0)Pr(XTt==0)

少なくとも1つのイベントを含む長さ間隔がある場合、状況が異なるのはなぜですか?ランダムに選択されたイベントの時間( ;間隔の最初から)は均一に分布しているため、違いはありません。Tt


なぜですか?なぜそうなのかはわかりませんが、これがトラブルの原因です。私自身で解決すると、あなたのアプローチ2と同じ答えが得られます。P(XT=1XT>0)=P(Xt=0)P(XTt=0)
2017

@StefanJorgensen私はあなたの意見を理解しています。しかし、それはであるべきではありませんか?P(Xt=1)P(XTt=0)=tete(Tt)=teT
abukaj 2017

@StefanJorgensenが2番目のコメントに答えました。これは、イベントを、間隔を未知数のイベントの2つのサブ間隔に分割するポイントと見なしたためです。
abukaj 2017

@StefanJorgensenこれは質問に答えるようなものです。時間にイベントがあった場合、またはどちらにもイベントがなかった確率はどれくらいですか。t(0,t)(t,T)
abukaj 2017

方程式は、またはいずれにもイベントがない確率であり、条件付けをまったく表さないことを除いて。これが、実際には条件付けが結果を変える場合に、少なくとも1つのイベントが存在する可能性(条件付けが最初に失われたため)になるという理由です。P(Xt=0)P(XTt=0)(0,t)(t,T)
2017

回答:


2

ようやくそれがわかりました!

@comboのアドバイスによると、「オカレンス」という用語を使用します。

驚いたことに、私の直感の根拠の最初の部分はほぼ正しかった。

観測された発生が間隔の最初から時間あった場合、またはオープン間隔で発生が発生しなかった確率を計算するだけで十分です:。t(0,t)(t,T)Pr(XT=1t)=Pr(Xt=0)Pr(XTt=0)=etetT=eT=Pr(XT=0)

違いは、をで、と見なすことができます。ここで、は無限小間隔の長さです。以降、我々が見ることができる点における密度としてのみの間隔内での発生。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)

の値は不明なので、ます。これまでのところ良好です。区間内に1回だけ発生する無条件の確率はわかっています。ただし、ですることにより、サンプリングスペースの部分は破棄されました。したがって、残りのすべての確率/密度を係数で正規化する必要があります。。t0Tpdf(t)dt=TeT=Pr(XT=1)XT>0eT11eT

したがって、直感は確率と密度の誤解でした。


1

混乱を避けるための用語の簡単な注記:特定の時間に発生する何かを、サンプルの要素としてのイベントのより厳密な定義との混同を避けるために、「イベント」ではなく「発生」と呼びます。スペース。

ポアソンカウンティングプロセスの定義から始めましょう。してみましょう、時間によって起こる発生回数も。次のプロパティがあります。N(t)T

  1. N(0)=0
  2. N(T1),(N(T2)N(T1))はに対して独立しています(独立したインクリメントプロパティ)T1<T2
  3. 長さ任意の区間での出現回数は、パラメーター(ここでは、)のポアソン確率変数です。tλtλ=1

独立したインクリメントプロパティは、あなたをつまずかせるものです-特に厳密な不等式の影響。

時間で発生する場合、ある確率はどれくらいですか?アプローチに従って、プロセスを3つのセグメントに分割しましょう。一部のでは、次のようになります: そして今、確率 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)
これはあなたが書いたものとまったく同じではないことに注意してください。主な違いは2つあります。
  1. 間隔はばらばらではないため、はから独立していますが、どちらも独立しています。。N(tϵ)N(0)N(T)N(t+ϵ)N(t+ϵ)N(tϵ)
  2. 間隔はすべて正の尺度です。あなたのアプローチでは、間隔を3つの部分分割します。問題は、間隔で発生したイベントに条件を付けていることです、測定値はゼロです(これが問題である理由の詳細については、この投稿を参照してください)。[0,T][0,t),[t],(t,T][t]

1. 最後の段落とについて少し混乱しています。サイズ間隔のオーバーラップがあるということですか?2.混乱の影響である可能性がありますが、およびでの発生が与えられた場合:。ϵ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 2017

それについて申し訳ありませんが、私はそれらが同じであることを意味しましたが、編集中に私のコンピュータが死にました。違いがまさに混乱の原因です。イベントはイベントに依存しません、。(制限ではなく等価)を作成しようとすると、間隔はシングルトンであり、メジャーがゼロのセットを条件としています。{N(tϵ)N(0)=0,N(T)N(t+ϵ)=0}{N(t+ϵ)N(tϵ)=1}ϵ>0ϵ=0
コンボ

1
さて、依存関係はわかりません。しかし、私はなんとか間違いを見つけました(私の答えを参照してください)-そして、間隔が重要だったため、+ 1となりました。[tϵ,t+ϵ]
abukaj 2017

0

問題へのアプローチの2つの方法のうち、2番目の方法は適切であるようです。それははるかに厳格で、私には完全に理にかなっています。しかし、私は最初のアプローチを理解するのにもう少し困難があり、これがあなたの間違いの原因であると信じる理由を与えました。

まず、好奇心から、なぜを計算するのですか?以下で使用したのとまったく同じアプローチを考えると、この確率(関連する限り)はと等しくなるはずですこれは、定義上とは等しくありません。。お役に立てれば!P(XT=1|XT>0)TeT1eTeT


私はまだ最初のアプローチで私の間違いを見ることができません。私は、以来、それを計算しています -編集を参照してください。P(XT>1|XT>0)=1P(XT=1|XT>0)
abukaj 2017
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.