母平均の信頼区間～母分散が既知の場合

/ tau / 5件のコメント

概要

母集団の分散がわかっている場合の、母平均の信頼区間の推定について。

信頼区間の推定の考え方は以下の通り。

  1. サンプルを抽出し、標本平均\overline{x}を求める
  2. 既知の分散σ2から標本平均は正規分布N(μ, σ2/n)に従う
  3. 標本平均をμ, σ2/nで標準化し、標準正規分布の信頼確率αに対する信頼区間を設定
  4. 母平均μの信頼区間を計算

手順

まず、母集団からn個のサンプルx1, …, xnを抽出し、その平均を求める。

(1)   \begin{equation*} \overline{x} = \frac{1}{n}\sum_{i=1}^n x_i \end{equation*}

次に平均と分散で標準化した変数に対して、意図する確率値αに対する標準正規分布の確率変数値zを使って信頼区間を設定。両側の境界を持つ信頼区間の場合は以下のようになる。

(2)   \begin{equation*} z\left( p \le \frac{1 - \alpha}{2} \right) \le \frac{\overline{x} - \mu}{\sqrt{\sigma^2 / n}} \le z\left( p \le \frac{1+ \alpha}{2} \right) \end{equation*}

これを移項してμの信頼区間として表示。

(3)   \begin{align*} \overline{x} - z\left( p \le \frac{1 - \alpha}{2} \right) \sqrt{\frac{\sigma^2}{n}} \le \mu \le \overline{x} + z\left( p \le \frac{1+ \alpha}{2} \right) \sqrt{\frac{\sigma^2}{n}} \end{align*}

信頼確率αに対応する標準正規分布のzを設定してμの信頼区間を算出する。たとえば両側95%信頼区間なら、片側2.5%確率に対応する1.96など、標準正規分布のzの値はこちらを参照。

(4)   \begin{align*} \overline{x} - 1.96 \sqrt{\frac{\sigma^2}{n}} \le \mu \le \overline{x} + 1.96 \sqrt{\frac{\sigma^2}{n}} \end{align*}

例題

e-statの身長・体重に関する国民健康・栄養調査2017年のデータから、40歳代の日本国民の身長の平均171.2cm及び標準偏差6.0cmを母集団のパラメーターとして用いる（データ数は374人）。

このパラメーターから、正規分布に従う10個の乱数を発生させた結果が以下の通り。

これらのデータの平均は170.6となり、これとσ2 = 36、サンプル数10、両側95%に対する1.96を用いて、信頼区間は以下のように計算される。

(5)   \begin{gather*} 170.6 - 1.96 \sqrt{\frac{36}{10}} \le \mu \le 170.6 + 1.96 \sqrt{\frac{36}{10}} \\ 166.9 \le \mu \le 174.3 \end{gather*}

【注】上記のデータはPythonでseed(1)として発生させた。

当初seed(0)で発生させた際には以下のようになり、95%信頼区間が母集団の平均を含まなくなった。

(6)   \begin{gather*} 175.6 - 1.96 \sqrt{\frac{36}{10}} \le \mu \le 175.6 + 1.96 \sqrt{\frac{36}{10}} \\ 171.9 \le \mu \le 179.3 \end{gather*}

seed(0)はよく使う系列だが、このようなこともあるので乱数系列を複数変えて試すのが望ましい。

サンプルサイズに対する信頼区間の傾向

サンプルサイズを大きくしていったときの平均身長の95%信頼区間は以下の通りで、かなりばらつきながら徐々に区間幅は小さくなるが、ある程度サンプルサイズを大きくしてもあまり顕著な区間幅の減少はみられない。

これは信頼区間に現れる1/\sqrt{n}のグラフを描いてみると分かるが、n=20程度まで急激に小さくなり、その後の減少スピードはかなり遅いことがわかる。したがって、信頼区間を狭めようとしても、効果があるのはせいぜいデータ数50程度までということになる。

【補足】

本記事にいただいたコメントの通り、これの考え方は適切ではない。正しくは、1.96 \sqrt{\sigma^2 / 2}などのグラフを描くべき。ご指摘に感謝します。

なお、1つ目のグラフの計算手順は以下の通り。

  1. 母集団の平均・標準偏差から、サンプルサイズを変えながら正規乱数を発生させる
  2. サンプルごとにサンプル平均を計算する
  3. サンプル平均と母分散から母平均推定の信頼区間の上限値と下限値を計算してリストに追加する
  4. 結果をグラフに表示する

 

5件のコメント

  1. データの解析上とても有益な記事を拝見しています。
    当方の浅学のためですが、記事中にある「信頼区間に現れる1/sqrt(n)のグラフ」の意味がよく理解できません。この値はｎのみで決まっていますので、母集団の分散が異なっている場合でも共通ということですか。また、「信頼区間を挟めようと」という意味についてもご教示いただければ幸いです。

    返信

    1. コメントありがとうございます。
      ご指摘の通りです。
      上の考え方は正しくなくて、本来ならσ^2/nの平方根で考えなければいけないですね。

      信頼区間を狭めることについては、平均±平方根の項のレンジが信頼区間の幅なので、平方根の項が小さいほど信頼区間の幅が小さくなる＝より確度の高い推定になる、と考えました。
      ですので、平方根の項の減少が緩やかになって、以後あまりデータ数を増やす効果がなさそうな50程度を（これはご指摘の通り正しくないですが）を目安と考えました。
      俗世の統計もので「データ数が少ないのであてにならない」として幾千ものデータを集めたりしていますが、それが過剰なのでは？との問題意識から考えてみたものです。

      返信

      1. ご丁寧にありがとうございました。
        データ数は多ければいいというものではなく、試験コストにも跳ね返ってきます。まったく同感の問題意識です。

        返信

  2. (5)のμの値が166.9ぐらいになるように思います。
    勘違いであったら申し訳ありませ。

    返信

    1. おっしゃる通りです。
      168.7としていましたが、正しくは166.88…なので166.9とするのが正しいですね。
      この式はPCの電卓で計算した記憶があり、その際のミスです。

      また、右側の方も値が違っていました。
      172.5としていましたが、正しくは174.31…なので174.3が正でした。

      ご指摘いただいてありがとうございました。

      返信

コメントを残す

メールアドレスが公開されることはありません。 * が付いている欄は必須項目です