数理統計の問題
ある電子部品の寿命(単位: 年)は、 が与えられたとき、確率密度関数
の指数分布に従うとする。 は未知の母数で、その事前分布として確率密度関数
のガンマ分布(形状母数 、尺度の逆数 )を考える。
が与えられたとき互いに独立なこの部品 個の寿命を調べたところ、
であった。このデータのもとでの の事後平均 を求めよ。
この問題で使うもの
- 共役事前分布とベイズ推定(説明と入口の問題)
- ガンマ関数・ベータ関数とガンマ分布(説明と入口の問題)
この問題を解けるようになるまで(68 技能)
読み込み中…
右端(太枠)がこの問題で使う技能です。左へドラッグすると、さらに前の前提をさかのぼれます。技能を押すと、そこまでの道筋が光ります。
一覧で見る(前提が先)
- 樹形図で数え上げる
- 積の法則と和の法則
- 和の法則・積の法則と樹形図
- 順列・円順列・重複順列
- 組合せ
- 確率の定義と和事象・余事象
- 期待値
- 確率変数と確率分布
- 平均と合計
- 平均値・中央値・最頻値
- 代表値(平均値・中央値・最頻値)
- 分散と標準偏差
- 変量の変換と平均・分散
- 確率変数の期待値・分散と
- 分数のたし算・ひき算
- 小数の四則計算
- 正負の数の加法・減法
- 分数のかけ算・わり算
- 四則の混じった計算の順序
- 正負の数の乗除・累乗と四則混合
- 分配法則を使った計算のくふう
- 文字式の計算(同類項・分配法則)
- 乗法公式による展開
- 因数分解
- 3 次式までの展開と因数分解
- 微分係数と導関数の定義
- 多項式の微分
- 多項式の不定積分
- 定積分の計算
- 連続型確率変数と確率密度関数
- 等間隔にならんだ数の和
- 等差・等比数列の和と Σ
- 数列の極限
- 関数の極限
- 積と商の微分
- 合成関数の微分
- 円の周の長さと面積
- おうぎ形の弧の長さと面積
- 弧の長さからおうぎ形の中心角と面積を求める
- 三角形・四角形の面積
- 約数の見つけ方
- 素因数分解と約数の個数
- 平方根の意味と根号の簡単化
- 三平方の定理
- 特別な直角三角形の辺の比
- 弧度法と単位円による三角関数
- の極限
- 加法定理と 2 倍角・半角・積和の公式
- 三角関数の導関数
- 根号の乗除と加減
- 指数の拡張と指数関数
- 対数の定義と性質
- 自然対数の底
- 指数関数・対数関数の導関数
- 基本的な関数の不定積分
- 置換積分
- 部分積分
- 定積分の置換・部分積分と偶関数・奇関数
- 広義積分
- 一様分布・指数分布
- ガンマ関数・ベータ関数とガンマ分布
- 確率変数の和と積(独立な場合)
- 重積分と累次積分
- 同時分布・周辺分布と独立性
- 条件付き確率と乗法定理
- 条件付き分布と全期待値の公式
- 原因の確率(ベイズの定理)
- 共役事前分布とベイズ推定
ヒント
行き詰まったら、上から順に 1 つずつ開いてください。
ヒント 1
共役事前分布の考え方を使います。指数分布の尤度にガンマ分布の事前密度を掛けると、事後分布もガンマ分布になります。
ヒント 2
尤度は です。これに を掛けて、 の指数と の肩の係数がどう変わるかを見ましょう。ガンマ分布(形状 、尺度の逆数 )の平均は です。
ヒント 3
データの和は なので、事後密度は に比例し、形状母数 、尺度の逆数 のガンマ分布になります。
答えと解説
答えを見る
どこでつまずきましたか?(選ぶと、戻る場所を案内します)
解説を読む
手順 1: 尤度を書く。
観測値を とすると、 が与えられたとき独立なので
ここで なので 。
手順 2: 事後分布を求める(ベイズの定理)。
事後密度は「事前密度 × 尤度」に比例するから
これは形状母数 、尺度の逆数 のガンマ分布の密度の形である。事前分布がガンマ分布で事後分布もガンマ分布になる(ガンマ分布は指数分布の共役事前分布)。形状母数は (データ数を足す)、尺度の逆数は (データの和を足す)と更新される。
正規化定数は、 と置換して
なので、。
手順 3: 事後平均を計算する。
( を使った。)一般に形状母数 、尺度の逆数 のガンマ分布の平均は であり、 と一致する。
よって事後平均は 。
別解
事前分布は形状母数 、率(尺度の逆数) のガンマ分布である。確かめると、密度は となり、問題文の と一致する。
尤度は である。ここで である。
事後密度は尤度と事前密度の積に比例するので、
これは形状母数 、率 のガンマ分布である。
ガンマ分布の平均は(形状母数)/(率)なので、。
検算: 一般式 で、上の結果と一致する。