2つの独立したポアソン確率変数の加重和

ウィキペディアを使用して、2つのポアソン確率変数の合計から生じる確率質量関数を計算する方法を見つけました。しかし、私が持っているアプローチは間違っていると思います。

LET 平均値を有する2つの独立したポアソン確率変数である、及び、及び定数、次いで確率発生関数でありますによって与えられる $X_1, X_2$ $\lambda_1, \lambda_2$ $S_2 = a_1 X_1+a_2 X_2$ $a_1$ $a_2$ $S_2$ 今、ポアソン確率変数の確率発生関数であるという事実使用、我々は、2つの独立したポアソンランダムの和の確率生成関数を書くことができるが変数として

G_{S_{2}} (z) = E (z^{S_{2}}) = E (z^{a_{1} X_{1} + a_{2} X_{2}}) G_{X_{1}} (z^{a_{1}}) G_{X_{2}} (z^{a_{2}}) .

$G_{S_2}(z) = \operatorname{E}(z^{S_2})= \operatorname{E}(z^{a_1 X_1+a_2 X_2}) G_{X_1}(z^{a_1})G_{X_2}(z^{a_2}).$

G_{X_{i}} (z) = e^{λ_{i} (z - 1)}

$G_{X_i}(z) = \textrm{e}^{\lambda_i(z - 1)}$

確率質量関数は、

導関数を取ることによって回復されるようです

\begin{aligned} G_{S_{2}} (z) & = e^{λ_{1} (z^{a_{1}} - 1)} e^{λ_{2} (z^{a_{2}} - 1)} \\ = e^{λ_{1} (z^{a_{1}} - 1) + λ_{2} (z^{a_{2}} - 1)} . \end{aligned}

$\begin{aligned} G_{S_2}(z) &= \textrm{e}^{\lambda_1(z^{a_1} - 1)}\textrm{e}^{\lambda_2(z^{a_2} - 1)} \\ &= \textrm{e}^{\lambda_1(z^{a_1} - 1)+\lambda_2(z^{a_2} - 1)}. \end{aligned}$

S_{2}

$S_2$

G_{S_{2}} (z)

$G_{S_2}(z)$

、ここで

\Pr (S_{2} = k) = \frac{G_{S_{2}}^{(k)} (0)}{k!}

$\operatorname{Pr}(S_2 = k) = \frac{G_{S_2}^{(k)}(0)}{k!}$

。

G_{S_{2}}^{(k)} = \frac{d^{k} G_{S_{2}} (z)}{d z^{k}}

$G_{S_2}^{(k)} = \frac{d^k G_{S_2}(z)}{ d z^k}$

これは正しいですか？定数とため、確率質量関数を取得するための導関数だけを使用することはできないと感じています。これは正解？別のアプローチはありますか？ $a_1$ $a_2$

これが正しい場合、すべてのkについて無限合計を切り捨てることにより、累積分布の近似を取得できますか？

distributions poisson-distribution

— ミシェル
ソース

なぜと加数をスケーリングしている

と

？合計は、これのない別のポアソン分布です。変数は正の整数の値を取るので、最初の

倍プラス

a_{1}

$a_1$

a_{2}

$a_2$

1

$1$

通常、1秒あたり

回は非常に不自然で、両方の変数の値を復元できます。

\sqrt{2}

$\sqrt{2}$

— ダグラスザレ2013

ここでの難しさがない限り、両方のことである

と

整数であり、1は、ことを確認することはできません

整数値のみを取ります。したがって、あなただけでなく、見つける必要が

の整数値に対して

も

毎に

として表すことができる

非負の整数のための

と

。

a_{1}

$a_1$

a_{2}

$a_2$

S_{2}

$S_2$

P (S_{2} = k)

$P(S_2 = k)$

k

$k$

P (S_{2} = α)

$P(S_2 = \alpha)$

α

$\alpha$

a_{1} m + a_{2} n

$a_1m + a_2n$

m

$m$

n

$n$

— ディリップサルワテ

@DilipSarwateそれは可能ですか？これを行う他のアプローチはありますか？

— ミシェル

@DouglasZare私はこれをしなければなりません...多分私はある種のブートストラップ法に頼らなければなりません。

— ミシェル

S_{2}

$S_2$

α

$\alpha$

P {S_{2} = α} = \sum_{a_{1} m + a_{2} n = α} P {X_{1} = m} P {X_{2} = n} = \sum_{a_{1} m + a_{2} n = α} \exp (- λ_{1} m) \frac{λ_{1}^{m}}{m!} \exp (- λ_{2} n) \frac{λ_{2}^{n}}{n!} .

$P\{S_2 = \alpha\} = \sum_{a_1m + a_2n = \alpha}P\{X_1=m\}P\{X_2=n\} = \sum_{a_1m + a_2n = \alpha} \exp(-\lambda_1m)\frac{\lambda_1^m}{m!}\exp(-\lambda_2n)\frac{\lambda_2^n}{n!}.$

a_{1}

$a_1$

a_{2}

$a_2$

a_{1} = a_{2} = 1

$a_1=a_2=1$

S_{2}

$S_2$

λ_{1} + λ_{2}

$\lambda_1+\lambda_2$

回答:

この線形結合のすべての確率が単一の値に集中していない場合、コーニッシュフィッシャー展開が（逆）CDFに適切な近似を提供するように見えます。

$S_2$ $\beta_1$

\frac{a_{1}^{3} λ_{1} + a_{2}^{3} λ_{2}}{{(\sqrt{a_{1}^{2} λ_{1} + a_{2}^{2} λ_{2}})}^{3}}

$\frac{a_1^3 \lambda_1 + a_2^3 \lambda_2}{\left(\sqrt{a_1^2 \lambda_1 + a_2^2 \lambda_2}\right)^3}$

$\beta_2$

\frac{a_{1}^{4} λ_{1} + 3 a_{1}^{4} λ_{1}^{2} + a_{2}^{4} λ_{2} + 6 a_{1}^{2} a_{2}^{2} λ_{1} λ_{2} + 3 a_{2}^{4} λ_{2}^{2}}{{(a_{1}^{2} λ_{1} + a_{2}^{2} λ_{2})}^{2}} .

$\frac{a_1^4 \lambda_1 + 3a_1^4 \lambda_1^2 + a_2^4 \lambda_2 + 6 a_1^2 a_2^2 \lambda_1 \lambda_2 + 3 a_2^4 \lambda_2^2}{\left(a_1^2 \lambda_1 + a_2^2 \lambda_2\right)^2}.$

$\alpha$ $S_2$

w_{α} = z + \frac{1}{6} β_{1} (z^{2} - 1) + \frac{1}{24} (β_{2} - 3) (z^{2} - 3) z - \frac{1}{36} β_{1}^{2} z (2 z^{2} - 5 z) - \frac{1}{24} (β_{2} - 3) β_{1} (z^{4} - 5 z^{2} + 2)

$w_\alpha = z +\frac{1}{6} \beta _1 \left(z^2-1\right) +\frac{1}{24} \left(\beta _2-3\right) \left(z^2-3\right) z-\frac{1}{36} \beta _1^2 z \left(2 z^2-5 z\right)-\frac{1}{24} \left(\beta _2-3\right) \beta _1 \left(z^4-5 z^2+2\right)$

$z$ $\alpha$ $S_2$

a_{1} λ_{1} + a_{2} λ_{2} + w_{α} \sqrt{a_{1}^{2} λ_{1} + a_{2}^{2} λ_{2}} .

$a_1 \lambda_1 + a_2 \lambda_2 + w_\alpha \sqrt{a_1^2 \lambda_1 + a_2^2 \lambda_2}.$

$\lambda_1$ $\lambda_2$ $5$ $\lambda_1 = 5,$ $\lambda_2=5\pi/2,$ $a_1=\pi,$ $a_2=-2$

$S_2$

— whuber
ソース

λ_{1}

$\lambda_1$

λ_{2} \leq 5

$\lambda_2 \leq 5$

畳み込みを使用します。

$f_{X_1}(x_1)= \dfrac{\lambda^{x_1}e^{-\lambda}}{x_1!}$ $x_1 \geq 0$ $f_{X_1}(x_1)= 0$ $f_{X_2}(x_2)=\dfrac{\lambda^{x_2}e^{-\lambda}}{x_2!}$ $x_2 \geq 0$ $f_{X_2}(x_2)= 0$

$Z=X_1+X_2\rightarrow X_1=Z-X_2$

f_{Z} (z) = \int_{- \infty}^{\infty} \int_{- \infty}^{\infty} f_{x_{1}, x_{2}} (z - x_{2}, x_{2}) d x_{1} d x_{2}

$f_Z(z)=\int\limits_{-\infty}^{\infty}\int\limits_{-\infty}^{\infty}f_{x_1,x_2}(z-x_2,x_2)dx_1dx_2$

$X_1$ $X_2$

f_{Z} (z) = \int_{- \infty}^{\infty} \int_{- \infty}^{\infty} f_{X_{1}} (z - x_{2}) f_{X_{2}} (x_{2}) d x_{1} d x_{2}

$f_Z(z)=\int\limits_{-\infty}^{\infty}\int\limits_{-\infty}^{\infty}f_{X_1}(z-x_2)f_{X_2}(x_2)dx_1dx_2$

f_{Z} (z) = \sum_{x_{2} = 0}^{z} \frac{λ_{1}^{z - x_{2}} e^{- λ_{1}}}{(z - x_{2})!} \frac{λ_{2}^{x_{2}} e^{- λ_{2}}}{x_{2}!}

$f_Z(z)=\sum\limits_{x_2=0}^{z} \dfrac{\lambda^{z-x_{2}}_1e^{-\lambda_1}}{(z-x_2)!}\dfrac{\lambda^{x_2}_2e^{-\lambda_2}}{x_2!}$

= e^{- (λ_{1} + λ_{2})} \frac{(λ_{1} + λ_{2})^{z}}{z!}

$= e^{-(\lambda_1+\lambda_2)}\dfrac{(\lambda_1+\lambda_2)^z}{z!}$

λ_{1} + λ_{2}

$\lambda_1+\lambda_2$

— QAChip
ソース

a_{1} = a_{2} = 1

$a_1=a_2=1$

a_{1} = a_{2}

$a_1 = a_2$

a_{1} \neq a_{2}

$a_1 \ne a_2$

S = \sum_{i = 1}^{N} X_{i}

$S = \sum_{i=1}^N X_i$

N

$N$

X_{i}

$X_i$

i i d

$iid$

N

$N$

X_{i} = k

$X_i=k$

k N

$k N$

k

$k$

N

$N$

E [s^{k N}] = E [(s^{k})^{N}] = G_{N} (s^{k}) = \exp (λ (s^{k} - 1))

$E[s^{k N}] = E[(s^{k})^N] = G_N(s^{k}) = \exp(\lambda(s^k-1))$

Z = k_{1} N_{1} + k_{2} N_{2}

$Z = k_1 N_1 + k_2 N_2$

G_{Z} (s) = \exp (λ_{1} (s^{k_{1}} - 1) + λ_{2} (s^{k_{2}} - 1)) .

$G_Z(s) = \exp(\lambda_1(s^{k_1}-1) + \lambda_2(s^{k_2}-1)).$

λ = λ_{1} + λ_{2}

$\lambda = \lambda_1 + \lambda_2$

G_{Z} (s) = \exp (λ (\frac{λ_{1}}{λ} (s^{k_{1}} - 1) + \frac{λ_{2}}{λ} (s^{k_{1}} - 1)) = \exp (λ (\frac{λ_{1}}{λ} s^{k_{1}} + \frac{λ_{2}}{λ} s^{k_{1}} - 1)) .

$G_Z(s) = \exp(\lambda ( \frac{\lambda_1}{\lambda}(s^{k_1}-1)+ \frac{\lambda_2}{\lambda}(s^{k_1}-1)) = \exp(\lambda (\frac{\lambda_1}{\lambda}s^{k_1}+ \frac{\lambda_2}{\lambda}s^{k_1}-1)).$

λ = λ_{1} + λ_{2}

$\lambda = \lambda_1 + \lambda_2$

X_{i}

$X_i$

k_{1}

$k_1$

λ_{1} / λ

$\lambda_1/\lambda$

k_{2}

$k_2$

λ_{2} / λ

$\lambda_2/\lambda$

$k_1$ $k_2$ $0$

ディスカッション後に編集：

あなたができる最善はMCだと思います。これは複合ポアソン分布であるという導出を使用できます。

$Pois(\lambda)$
$i=1,\ldots,N$ $X_1$ $X_2$ $\lambda_1/\lambda$ $\lambda_1/\lambda$ $1$ $k_1$ $k_2$

あなたは秒単位で例えば100 000のサンプルを持っているでしょう。

あるいは、最初の表現の2つの加数を別々にサンプリングすることもできます...これも同じくらい簡単です。

定数係数k1とk2が完全に一般的である場合、その他すべて（FFT）は複雑になります。

— リック
ソース

そして、因子が整数の場合、最終的な分布はPanjerアルゴリズムによって見つけることができます。

— Ric

G_{S_{2}} (z) = e^{λ_{1} (z^{a_{1}} - 1)} e^{λ_{2} (z^{a_{2}} - 1)}

$G_{S_2}(z) = \textrm{e}^{\lambda_1(z^{a_1} - 1)}\textrm{e}^{\lambda_2(z^{a_2} - 1)}$

a_{1}, a_{2} \in R^{1}

$a_1,a_2 \in R^1$

P {S_{2} = α} = \sum_{a_{1} m + a_{2} n = α} P {X_{1} = m} P {X_{2} = n} = \sum_{a_{1} m + a_{2} n = α} \exp (- λ_{1} m) \frac{λ_{1}^{m}}{m!} \exp (- λ_{2} n) \frac{λ_{2}^{n}}{n!},

$P\{S_2 = \alpha\} = \sum_{a_1m + a_2n = \alpha}P\{X_1=m\}P\{X_2=n\} = \sum_{a_1m + a_2n = \alpha} \exp(-\lambda_1m)\frac{\lambda_1^m}{m!}\exp(-\lambda_2n)\frac{\lambda_2^n}{n!},$

a_{1}, a_{2}

$a_1,a_2$

こんにちはミシェル、私は私の応答を編集しました。はいPanjerの使用は制限されています。しかし、フーリエ変換アプローチを試すことができます。しかし、非整数ユニットは問題があります...この場合、何をすべきかについてもっと考えなければなりません。どちらの方法でも、結果が複合ポアソン分布になることに注意することが重要です（「単純な」ポアソン分布ではありません）。

— Ric

P r (S_{2} = x) = \frac{1}{2 π} \int_{R} e^{- i t x} G_{S 2} (e^{i t}) d t

$Pr(S_2=x) = \frac{1}{2\pi}\int_{\mathbf{R}} e^{-itx}G_{S2}(\mathrm{e}^{it})dt$

なんらかの方法で...特性関数を計算できる連続分布があった場合（そうすることにより）、これにより迅速で優れた結果が得られます。私たちの場合、それについて考える時間がもっと必要です。何かもっと簡単なはずです。

— Ric