一つの母分散に関する検定と推定

以前の記事で、母分散が未知の場合の一つの母平均の検定と推定について、基本的な考え方と進め方を解説しました。

一つの母平均に関する検定と推定(母分散が未知の場合)

母平均の検定では、t分布の性質を使って検定と推定を行いました。
今回は、母分散に注目し、一つの母分散に関する検定と推定について解説します。

この記事で分かること
  • χ2分布とは
  • 一つの母分散の検定の進め方
  • 一つの母分散の推定の進め方


1. 適用できる場面

以下の事例を使って、母分散の検定と推定を解説します。

【事例1】

ある薬品の1kg中の有効成分の量は、母分散が(0.50g)2の正規分布に従っています。最近、最新鋭の新ラインを立ち上げたことから、新ラインで製造した薬品の1kg中の有効成分の母分散が、従来のラインより小さくなったかどうかを検討することになりました。
新ラインからランダムに10個のサンプルを抜き取り、有効成分の量を測定した結果は、以下の通りでした。(単位:g)

10.2, 9.6, 9.9, 10.3, 10.1, 9.8, 10.2, 10.4, 10.1, 10.3

この事例では、新ライン製造品の有効成分の量が正規分布N(μ,σ2) (母平均μと母分散σ2はともに未知)に従うと考えて、この母集団から抜き取られた10個のサンプルのデータに基づいて、一つの母分散σ2の検定や推定を行うことが目的です。

2. χ2分布について

正規分布に従う母集団からn個のサンプルx1,x2,,xnを抽出すると、母分散に関する情報はサンプルの偏差平方和S=i=1n(xix¯)2やサンプルの分散V=Sn1に含まれています。

ここで、以下の基本事項が知られています。

基本事項(χ2(カイ2乗)分布)

n個のデータx1,x2,,xnが互いに独立に正規分布N(μ,σ2)に従うとき、χ2=Sσ2は自由度ϕ=n1χ2分布に従う。

χ2分布という新しい分布が出てきました。
Zが標準正規分布N(0,12)に従うとき、以下の式で求まるχ2が従う分布をχ2分布と言います。

χ2=Z12+Z22++Zϕ2


χ2分布は自由度のみによって形状が決まる確率分布で、自由度により形状が異なります
χ2分布は、これまで出てきた正規分布やt分布と違って、ゼロ以上の値で定義されていることと、左右非対称の分布形状であることが特徴です。

検定統計量がχ2分布に従うことが分かれば、あとは棄却域を設定して検定統計量が棄却域に入るかどうかを調べることで検定できます。

3. 一つの母分散に関する検定

それでは、母分散の検定手順を見ていきます。
基本的な流れは、以前の記事で解説した母平均の検定の流れと同じですので、異なる点を特に詳しく説明します。

3-1. 一つの母分散の検定手順

手順1. 帰無仮説H0と対立仮説H1を設定する。

検定の目的に応じて、(1)~(3)のいずれかを選択します。

(1) H0:σ2=σ02 (σ02は指定された値)

H1:σ2σ02 (両側検定)

(2) H0:σ2=σ02 (σ02は指定された値)

H1:σ2>σ02 (右片側検定)

(3) H0:σ2=σ02 (σ02は指定された値)

H1:σ2<σ02 (左片側検定)

手順2. 有意水準αを決める。

通常は、α=0.05とします。

手順3. 手順1(仮説)と手順2(有意水準)に対応した棄却域を決める。

(1)棄却域:χ02χ2(ϕ,1α/2)またはχ02χ2(ϕ,α/2)  (両側検定)

(2)棄却域:χ02χ2(ϕ,α)  (右片側検定)

(3)棄却域:χ02χ2(ϕ,α)  (左片側検定)

手順4. 採取したデータx1,x2.,xnから、検定統計量χ02の値を計算する。

   χ02=Sσ02

   ϕ=n1

手順5. 判定する。

χ02が棄却域に入れば、有意水準αで有意と判定し、帰無仮説H0を棄却して対立仮説H1を採択します。
χ02が棄却域に入らなければ、有意水準αで有意でないと判定し、帰無仮説H0を棄却しません。

3-2. 一つの母分散の検定における棄却域の求め方

母平均の検定では、正規分布またはt分布を使って検定しました。
正規分布とt分布は左右対称なので、両側検定では左右の棄却限界値の絶対値は同じで、単にーが付くかつかないかの違いでした。

一方、χ2分布は左右非対称のため、両側検定の場合、上側と下側の棄却限界値をそれぞれ求める必要があります
χ2分布の棄却限界値をExcelで求める場合、以下のようにセルに入力すれば求められます。(有意水準α=0.05の場合)

両側検定

確率5%を上下に2.5%ずつ割り当てます。
下側:「=CHISQ.INV(0.025,ϕ))」と入力して得られる値を採用します。

(ϕは自由度n1)
上側:「=CHISQ.INV.RT(0.025,ϕ))」と入力して得られる値を採用します。


片側検定

下側:「=CHISQ.INV(0.05,ϕ))」と入力して得られる値を採用します。

上側:「=CHISQ.INV.RT(0.05,ϕ))」と入力して得られる値を採用します。

3-3. 一つの母分散の検定の実施例

事例1について、検定手順に従って検定してみましょう。

手順1. 帰無仮説H0と対立仮説H1を設定する。

母分散が従来より小さくなったかどうかを知りたいので、左片側検定で帰無仮説と対立仮説を設定します。

H0:σ2=σ02  (σ02=0.502)

H1:σ2<σ02

手順2. 有意水準αを決める。

α=0.05

手順3. 棄却域を決める。

棄却域:χ02χ2(9,0.05)=3.33

(ϕ=n1=101=9)

手順4. 検定統計量χ02の値を計算する。

得られたデータより、

S=0.5690 (Excelの「DEVSQ」関数で求められる)

χ02=Sσ02=0.56900.502=2.276

手順5. 判定する。

χ02=2.276χ2(9,0.05)=3.33で検定統計量χ02は棄却域に入るので有意です。

よって帰無仮説H0を棄却して、新ライン製造品の有効成分の量の母分散σ2は、従来ラインの母分散0.502より小さくなったと判断できます。

4. 一つの母分散に関する推定

以前の記事で解説した通り、母分散σ2の推定も点推定区間推定の2種類があります。

4-1. 一つの母分散の推定手順

点推定はデータの分散V=Sn1を使えばよいです。

区間推定については、基本事項からχ2=Sσ2は自由度ϕ=n1χ2分布に従うことから、以下が成立します。

Pr(χ2(ϕ,1α/2)<Sσ2<χ2(ϕ,α/2))=1α

これを変形すると、以下のようになります。

  Pr(Sχ2(ϕ,α/2)<σ2<Sχ2(ϕ,1α/2))

  =1α

左項の括弧内のσ2の範囲が、信頼率1αの信頼区間です。

一つの母分散の推定手順をまとめると、以下のようになります。

点推定

σ^2=V (Excelの「DEVSQ」関数で求められる)


区間推定:信頼率1αの信頼区間

(Sχ2(ϕ,α/2),Sχ2(ϕ,1α/2))

4-2. 一つの母分散の推定の実施例

事例1について、点推定と区間推定を行ってみましょう。

点推定

σ^2=V=0.06322=0.2512

区間推定:信頼率95%の信頼区間を求めます。

(Sχ2(9,0.025),Sχ2(9,0.975))

=(0.569019.02,0.56902.70)

=(0.02992,0.2107)

=(0.1732,0.4592)

5. 実践のためのアドバイス

モノづくりにおいては、特性値のばらつきに着目することはとても重要で、通常は工程能力指数で規格値に対する余裕度を判断します。
工程能力指数の算出には母分散が必要なので、母分散が基準値と比べて変化しているかどうか、あるいは母分散はどの程度かを検定や推定したい場面は多いでしょう。
そのときは、今回のχ2分布を使った検定と推定を使ってください。

6. おわりに

今回は、一つの母分散の検定と推定について解説しました。
一つ母分散の検定と推定では、左右非対称のχ2分布を使うので、棄却限界値を求めるときに注意が必要ですが、基本的な考え方と手順は変わりません。
工程能力を改善するためにはばらつきの低減が必要なので、対策の効果の有無を確認するときに、母分散の検定と推定を使うとよいでしょう。

コメントを残す

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