[Projectμ] プロジェクトμ ブレーキパッド NS-C フロント用 C-HR NGX10 18/5~ 1.2L 本州は送料無料 北海道は送料500円(税別) 沖縄・離島は送料1000円(税別)

ポアソン分布 (poisson distribution) は平均値 \(\lambda>0\) を唯一のパラメータとする離散確率分布ですが、これをガンマ分布 (gamma distribution) で混合してみることを考えます。つまりポアソン分布を \(p(x|\lambda)\)、

[Projectμ] プロジェクトμ ブレーキパッド NS-C フロント用 C-HR NGX10 18/5~ 1.2L 本州は送料無料 北海道は送料500円(税別) 沖縄・離島は送料1000円(税別) 激安直営店

送料無料![品番SG13] [MF33S MRワゴン] [H23/2-H24/5] [MF33S] [定員4] Regaliaレガリア シートカバー ダイアモンドステッチ;モーターステージ スリップオンマフラー ワイルドブラス ナロースラッシュ 00年-06年 FLSTN 黒 20050141 HD店;【送料無料】[SHOEI] 【NEOTEC II/ネオテック2】 ライトシルバー XL  システムヘルメット(フリップアップ) ショーエイ[4512048475314]

[Projectμ] プロジェクトμ ブレーキパッド NS-C フロント用 C-HR NGX10 18/5~ 1.2L 本州は送料無料 北海道は送料500円(税別) 沖縄・離島は送料1000円(税別);(クーポン配布中)X.A.M A5507-42 525-42T スプロケット A5507-42;【送料無料 ハイエース200系】 195/80R15 15インチ LA STRADA ラ・ストラーダ クロスエイト 6J 6.00-15 YOKOHAMA ヨコハマ RY52 107/105L サマータイヤ ホイール4本セット【YOsum18】;ヨコハマ スタッドレス アイスガード シックス アイスガード6 IG60 215/45R18 & ジャパンsport01 18×7.5J 114.3/5H + 53 アクセラ BYEFP;RSタイチ 2018年春夏モデル ステルス レザーグローブ 黒 Lサイズ (スマホ対応) RST445BK01L JP店

[Projectμ] プロジェクトμ ブレーキパッド レーシングN+ フロント用 ルーチェ HCSS 87/8~ 3.0L 本州・北海道は送料無料 沖縄・離島は送料1000円(税別), Project μ プロジェクト・ミュー ブレーキパッド TYPE PS フロント WRX STI GRB, 【USA在庫あり】 ITP タイヤ サンドスター 21x7-10 2PR フロント 370941 HD店, 175/60R16 82H NEXEN ネクセン N blue HD Plus エヌ ブルー エイチディー プラス WORK EMOTION CR kiwami ワーク エモーション CR 極 サマータイヤホイール4本セット, アルティナ ラバーフロアマット トヨタ アルファード ANH・GGH2# 後期 8人乗/G/システムコンソール付 2011年11月~ ★送料無料★(一部離島等除く) 【キノクニ】ニッサン ホンダ マツダ ミツビシ スバル etc(K100-2058BM、Run・Max ブロック)

、ガンマ分布を \(q(\lambda|a,b)\) として、\begin{align*}r(x|a,b):=\int_0^\infty p(x|\lambda) q(\lambda|a,b)d\lambda\end{align*}を計算すると何が起こるのか、という話です。

確率変数を「サイコロをふる」ことに例えれば、上の積分で書かれている離散確率分布 \(r(x|a,b)\) は、まずはじめにガンマ分布のサイコロを振って平均値を求め、その平均値を使ったポアソン分布のサイコロをふる、という二段階でランダムな値が決定される Artina アルティナ ラバーフロアマット 1台分 レガシィ/レガシィツーリングワゴン BM/BR H21.05~ 全車種適用、ということです。これが確率分布を混合するということのイメージでしょうか[1]

さて、混合の結果、\(r(x|a,b)\) は負の二項分布 (negative binomial distribution) になります。負の二項分布は(いろいろなパラメータ化があるけども)\[r(x|a,p)=\begin{pmatrix}x+a-1\\a\end{pmatrix}\,\,p^x(1-p)^a\]という感じで定義される離散分布です。\(a\) は一般には正の整数で \(p\) は確率値です。

この証明は以下の通り。ポアソン分布は\begin{align*}p(x|\lambda)=\frac{\lambda^x}{x!}e^{-\lambda}\end{align*}で、ガンマ分布は\begin{align*}q(\lambda|a,b)=\frac{b^a}{\Gamma(a)}\lambda^{a-1}e^{-bx}\end{align*}です(\(\Gamma(\cdot)\) はガンマ関数)。これを直接計算すると\begin{align*}
r(x|a,b)
&=\int \frac{\lambda^x}{x!}e^{-\lambda}\frac{b^a}{\Gamma(a)}\lambda^{a-1}e^{-b\lambda}d\lambda\\
&=\frac{b^a}{x!\Gamma(a)}\int\lambda^{x+a-1}e^{-\lambda(1+b)}d\lambda\\
&=\frac{b^a}{x!\Gamma(a)}\int\bigg(\frac{\mu}{1+b}\bigg)^{x+a-1}e^{-\mu}\frac{d\mu}{1+b}\\
&=\frac{b^a}{x!\Gamma(a)(1+b)^{x+a}}\Gamma(x+a)\\
&=\frac{\Gamma(x+a)}{\Gamma(x+1)\Gamma(a)}b^a(1+b)^{-x-a}\\
&=\frac{\Gamma(x+a)}{\Gamma(x+1)\Gamma(a)}p^x(1-p)^a
\end{align*}
となります。ただし途中で \(\mu=(1+b)\lambda\) 、最後に \(b=(1-p)/p\) と変数変換しています。また \(x!=\Gamma(x+1)\) を使っています。

もし \(a\) が正の整数ならば二項係数を使って 【代引不可】【受注生産品】MRA(エムアールエー):バイク スクリーン ヴァリオ グラデーション付/スライド調整可能 XRV750 AFRICA TWIN(96-00)用 スモーク MV032、
\[r(x|a,p)=\begin{pmatrix}x+a-1\\a\end{pmatrix}\,\,p^x(1-p)^a\]
となります。その前のガンマ関数を使ったバージョンは 【エスペリア】NISSAN モコ ドルチェ MOCO DOLCE MG33S ダウンサス MG33S R06A H25/10~ 2WD NA / ドルチェX フロント 、(一般には正の整数とされている)パラメータ \(a\) を実数に拡張したものといえます。

パラメータ \(a\) が実数だと通常の意味では解釈が困難になりますが 2輪 アグラス リアスライダー 5点セット 品番:P052-3321 ホワイト スズキ GSX-R600 2011年~2012年 JAN:4548664832262、GLM で誤差がポアソン分布としたとき(=ポアソン回帰)に overdispersion になってしまうときは 【店舗塗装サービス】【20 セルシオ ケイブレイク】セルシオ 20系 プラチナステンレスピラー 4P バイザー有り K-BREAK ロゴ入れ 2ヶ所、負の二項分布を誤差分布として \(a\) に相当するパラメータ[2]を調整することで overdispersion を回避することができます(このへんの話はいつか書きたいですね)。こういうケースでは実数で定義されていると便利なんですね。

以上でおしまいです。こういう~分布と~分布を混合すると~分布になる、みたいな話を時間のあるときにまとめていこうと思います。


  1. [1] 英語では mixture distribution とか compound distribution とか言ったりするらしい。
  2. [2] 実際には GLM のときはパラメトリゼーションが違うからこの言い方は厳密ではないのだけど。

北海道は送料500円(税別) 1.2L 沖縄・離島は送料1000円(税別) 18/5~ プロジェクトμ 1.2L [Projectμ] C-HR NGX10 プロジェクトμ NS-C ブレーキパッド 本州は送料無料 フロント用