104
115

Delete article

Deleted articles cannot be recovered.

Draft of this article would be also deleted.

Are you sure you want to delete this article?

More than 3 years have passed since last update.

これでわかった統蚈的怜定Z怜定、怜定、分散のカむ二乗怜定、怜定

104
Last updated at Posted at 2022-12-16

統蚈の勉匷をしおいるず、いろいろな怜定手法が出おきたす。今回は怜定に぀いお代衚的なものをたずめおみたした。怜定は、区間掚定ずよく䌌た考えをベヌスずしおいたす。区間掚定っお䜕ずいう方がいたら、別の蚘事に区間掚定に぀いおたずめおみたしたので、先にこちらを確認しおみおください。

1. 怜定っおそもそも䜕

今回扱う怜定は、区間掚定ずずもに掚枬統蚈孊の最も基本的な手法の぀です。事前に仮説を立お、その仮説が正しいかを怜蚌するための手法です。独特な甚語や蚀い回しが䞊びたすが、これらを䜿いこなせるようになるず、䟋えば、新薬に効果があるかなどに぀いお統蚈的に語れるようになりたす。

s.png

新薬の怜定の䟋をめちゃくちゃざっくりで蚀うず、新薬ず停薬の集団の平均が同じ垰無仮説は䞻匵したい仮説の反察にするので $\mu_1=\mu_2$だずするず、実隓で埗た倀はめちゃめちゃレア、ずいうこずで、新薬の効果があった察立仮説$\mu_1 \neq \mu_2$ず刀断したす。はじめはしっくりこないかもしれたせんが、統蚈的怜定では䞻匵したい仮説ではなく、棄华した垰無仮説を怜定するずいうずころがポむントですいわゆる背理法。

1.1 怜定にはどんなものがあるの

怜定ず䞀蚀で蚀っおも、倚くの方法がありたす以䞋代衚的なものの䟋。今回はその䞭でも正芏分垃、分垃、カむ二乗分垃、分垃を䜿った最も基本的な怜定平均、母比率、分散の怜定に぀いお扱っおいきたいず思いたす。

皮類         ざっくりしたむメヌゞ                  
怜定 統蚈量が暙準正芏分垃に埓うこずを甚いる統蚈孊的怜定法
怜定 統蚈量が分垃に埓うこずを甚いる統蚈孊的怜定法
分散の $\chi^2$ 怜定 統蚈量がカむ二乗分垃に埓うこずを甚いる統蚈孊的怜定法
怜定 統蚈量が分垃に埓うこずを甚いる統蚈孊的怜定法
適合床怜定 床数が適合するか吊かを怜蚌する方法
順䜍和怜定 デヌタの順䜍を甚いお怜定する方法
尀床比怜定 尀床比を怜定量ずしお甚いる統蚈孊的怜定法

このように今回扱う぀の怜定は、確率分垃がすでに分かっおいる統蚈量を䜿っお仮説を評䟡するずいうのが基本の考え方です。やりたいこずに応じお利甚する統蚈量を䜿い分けるこずで統蚈的怜定を行うこずができたす。埌半に怜定したいこずに察応した統蚈量を蚘茉しおいたすので確認しおください。

b.png

2. 怜定の基本

さっそく、基本の甚語ず流れを抑えおいきたいず思いたす。

2.1 怜定で䜿う甚語

怜定で出おくる甚語に぀いお、たずはざっくりしたむメヌゞです。
・垰無仮説無に垰する仮説ずいうこずで、吊定したい仮説。 䟋平均は倉わらない
・察立仮説これが自分の䞻匵に盞圓するもので、本来採択したい仮説。 䟋平均は増えた
・有意氎準垰無仮説を棄华するか受容するかの刀断ラむン。
・怜定統蚈量怜定を目的ずしお、暙本から求める倀。
・p倀垰無仮説が正しいず仮定したずきに、統蚈量が実珟倀より極端な倀になる確率。
・棄华域倀が䞋回ったずきに垰無仮説を吊定する領域。

g.png

2.2 怜定の流れ

䞀般的には以䞋の流れで怜定を行いたす。
n_k.png

3. 怜定の抂芁

Z怜定、怜定、分散の $\chi^2$怜定、怜定のいずれの怜定もパラメトリックず呌ばれる、䜕かしらの確率分垃を前提ずした怜定手法です。それぞれ利甚する確率分垃をもずに以䞋のような考え方で怜定を行いたす。

3.1 怜定の考え方

怜定では、垰無仮説が正しいず仮定した䞖界で、今手元にある珟実デヌタから求めた倀が、その統蚈量が埓う確率分垃の䞭のどの䜍眮にあるかで評䟡を行いたす。これは、珟実から埗られた倀が、垰無仮説が正しいずいう仮説の䞖界で、どの皋床起きやすいのかを確率で考え評䟡するこずを意味したす。

垰無仮説が正しいずいう仮説の䞖界で、珟実から埗られた倀以䞊に極端な倀がでる確率を p倀ず呌び、垰無仮説を棄华するかどうかの刀断の境目に甚いる倀を有意氎準 $\alpha$ ず呌びたす。

p倀が有意氎準 $\alpha$ よりも小さい堎合、レアなこずが起きおいる偶然ではなく意味のあるこずが起きおいるず考え、垰無仮説を棄华し、察立仮説を採択したす統蚈的に有意な差がみられた状態。䞀方、p倀が有意氎準 $\alpha$ よりも倧きい堎合、この差は単なるばら぀きから生じたものであり、垰無仮説を棄华するこずはできない統蚈的に有意な差はみられなかった。ずいうこずになりたす。

r.png

ここで泚意しなくおはいけないのは、統蚈的に有意な差がみられた堎合でも、察立仮説が絶察に正しいずいうこずではなく、察立仮説を支持する぀の蚌拠が埗られた。ずいう意味を瀺しおいるにすぎないずいうこずです。

3.2 䞡偎怜定ず片偎怜定

怜定においお、䞀般には䞡偎怜定分垃の䞡方の裟の和を䜿う怜定を甚いたすが、片偎の確率だけで有意氎準 $\alpha$ ずする片偎怜定もありたす。

k_.png

片偎怜定を行うずきは、怜定を行う察象が 偎か ヌ偎か論理的に決たっおいる堎合に䜿いたす。その堎合、仮説の䞍等号ず棄华域に぀いおは、以䞋のように蚭定したす。
 $\mu  \mu_0$  $\mu$ は基準の倀よりも倧きいず思っおいるのだから、右裟に入っおきたら有意。
 $\mu  \mu_0$  $\mu$ は基準の倀よりも小さいず思っおいるのだから、巊裟に入っおきたら有意。


䞡偎で怜定すべきずころで安易に片偎で怜定をしおしたうず、本来反察偎で棄华できた垰無仮説が棄华されないため差があるにも関わらず差がないず刀定をくだしおしたうケヌスや、以䞋のような同じ怜定統蚈量でも片偎怜定にするこずで察立仮説を採択できおしたうケヌスもありたす。結果に差が出るため、片偎怜定を䜿う際は泚意が必芁です。
k.png

4. いろいろな怜定ず怜定量

ここでは各皮怜定に぀いお䞡偎怜定のみを蚘茉しおありたす。片偎怜定を䜿いたい堎合は、棄华域を前章を参考に倉えお䜿っおください。

4.1 平均の怜定暙本の怜定

たずは、最も基本的な母集団の平均の怜定からはじめたいず思いたす。やりたいこずは、母平均が倉わったずいえるのかを評䟡するこずです。母集団の平均の怜定は、区間掚定ず同様、暙準化倉量を利甚したす。母分散がわかっおいるずきわからないずき、暙本サむズが倧きいずき小さいずきの条件によっお蚈算方法が倉わるのも区間掚定ず同様です。もし区間掚定に぀いお埩習したい方は、こちらを確認しおみおください。

\text{暙本平均の暙準化倉量 }  z  = \frac{\bar{x}-\mu}{\frac{\sigma}{\sqrt{n}}} \\  

䞊の匏では、「暙本平均」の暙準化倉量を考えおいるので、分母の分散郚分は暙本サむズで割ったものになっおいたす。この蟺が「」ず思う方は、こちらの3章あたりに暙本平均が䟿利なこずに母集団の分散や平均を䜿っお正芏分垃で衚珟できる、ずいうこずを蚘茉しおみたしたので、ぜひ芋おみおください。繰り返しにはなりたすが、ここで考えおいるのは、暙本平均の埓う分垃、ずいうこずを考えればしっくりくるのではないかず思いたす。

垰無仮説 $H_0$  $\mu=\mu_0$ 、察立仮説 $H_1$  $\mu≠\mu_0$

条件          怜定統蚈量     怜定統蚈量の分垃         棄华域
$母集団は正芏分垃、\sigma^2は既知$ $z=\frac{\bar{x}-\mu_0}{\frac{\sigma}{\sqrt{n}}}$ 暙準正芏分垃 $ z> z_\frac{\alpha}{2}$
$倧暙本 (n≧30)$ $z=\frac{\bar{x}-\mu_0}{\frac{u}{\sqrt{n}}}$ 暙準正芏分垃 $ z> z_\frac{\alpha}{2}$
$母集団は正芏分垃、\sigma^2は未知小暙本$ $t=\frac{\bar{x}-\mu_0}{\frac{u}{\sqrt{n}}}$ 自由床 $n-1$ の分垃 $ t> t_\frac{\alpha}{2}(n-1)$

※$u$ : 暙本䞍偏分散

4.2 母比率の怜定1暙本の怜定

母比率の怜定も、平均の怜定ず基本的には同じ考えになりたす。こちらもやりたいこずは母比率が倉わったずいえるのかを評䟡です。母比率ずは䜕かに぀いおもくこちらの蚘事に蚘茉しおありたすので、埩習したい方はこちらの3.2章を確認しおみおください。

垰無仮説 $H_0$  $p=p_0$ 、察立仮説 $H_1$  $p \neq p_0$

条件              怜定統蚈量             怜定統蚈量の分垃        棄华域
倧暙本
$n>30 , np>5 , nq>5$
$z=\frac{\hat{p}-p_0}{\sqrt{\frac{p_0(1-p_0)}{n}}} ここで \hat{p}=\frac{x}{n}$ 暙準正芏分垃 $ z> z_\frac{\alpha}{2}$

4.3 母分散の怜定

母分散の怜定は、母分散が倉わったずいえるのかを評䟡したす。

垰無仮説 $H_0$  $\sigma^2 = \sigma_0^2$ 、察立仮説 $H_1$  $\sigma^2 \neq \sigma_0^2$

条件        怜定統蚈量      怜定統蚈量の分垃     棄华域          
母集団は正芏分垃 $\chi^2=\frac{(n-1)u^2}{\sigma^2}=\frac{ns^2}{\sigma^2}$ 自由床$n-1$の $\chi^2$ 分垃 $\chi^2\chi^2_{1-\frac{\alpha}{2}}(n-1)$
たたは
$\chi^2\chi^2_{\frac{\alpha}{2}}(n-1)$

4.4 平均の差の怜定2暙本の怜定

぀の母集団からそれぞれ取っおきた぀の暙本の暙本平均に差があった時、぀の母平均に差があるずいえるのかを評䟡したす。

垰無仮説 $H_0$$\mu_1=\mu_2$ 、察立仮説 $H_1$$\mu_1 \neq \mu_2$

条件          怜定統蚈量           怜定統蚈量の分垃      棄华域
$母集団は正芏分垃、\sigma_1^2,\sigma_2^2は既知$ $z=\frac{\bar{x_1}-\bar{x_2}}{\sqrt{\frac{\sigma_1^2}{n_1}+\frac{\sigma_2^2}{n_2}}}$ 暙準正芏分垃 $ z> z_\frac{\alpha}{2}$
$倧暙本n_1,n_2≧30$ $z=\frac{\bar{x_1}-\bar{x_2}}{\sqrt{\frac{u_1^2}{n_1}+\frac{u_2^2}{n_2}}}$ 暙準正芏分垃 $ z> z_\frac{\alpha}{2}$
$母集団は正芏分垃$
$\sigma_1^2=\sigma_2^2=\sigma^2等分散で\sigma^2は未知$
$t=\frac{\bar{x_1}-\bar{x_2}}{\sqrt{\frac{u^2}{n_1}+\frac{u^2}{n_2}}}$
$u^2=\frac{(n_1-1)u_1^2+(n_2-1)u_2^2}{n_1+n_2-2}$
自由床 $n_1+n_2-2$ の分垃 $ t> t_\frac{\alpha}{2}(n_1+n_2-2)$
$母集団は正芏分垃察応のある堎合$ $t=\frac{\bar{d}}{\sqrt{\frac{u_d^2}{n}}}$
$\bar{d}=\frac{1}{n}\sum{d_i}$
$u_d^2=\frac{1}{n-1}\sum{(d_i-\bar{d})^2}$
自由床$n-1$の分垃 $ t> t_\frac{\alpha}{2}(n-1)$

※$n_1$、$n_1$ :䞍偏分散

補足察応のある堎合は、いわば暙本の怜定ずしお「平均に差がない」を垰無仮説ずしたものです。たずえば、Aさんが薬を飲む前ず飲んだ埌の差などそれぞれのサンプルごずの差dを利甚しお怜定を行うものです。察応のあるデヌタの堎合、察応のある怜定を甚いた方が、第二皮の過誀が起きにくい、すなわち怜出力が䞊がる傟向があるず蚀われおいたす。

参考母分散が未知で等しくない堎合、怜定をそのたた行うこずはできたせん。怜定を改良したりェルチの怜定で算出するこずもできたすが、぀の矀の分散が異なる堎合は、平均の差を怜定するこずに意味があるのかずいう問題ベヌレンス・フィッシャヌ問題が生じるず蚀われおいたすので衚内の蚘茉は省略したす。ちなみにりェルチの怜定の自由床 $df$ は、以䞋の匏で求められるようです。耇雑 ですね。。。

df=\frac{\Bigl({(\frac{u_1^2}{n_1})+(\frac{u_2^2}{n_2})}\Bigr)^2}
   {\frac{\bigl(\frac{s_1^2}{n_1}\bigr)^2}{n_1-1}+\frac{\bigl(\frac{s_2^2}{n_2}\bigr)^2}{n_2-1}}

4.5 比率の差の怜定2暙本の怜定

぀の母集団の母比率に差があるずいえるのかを評䟡したす。

垰無仮説 $H_0$  $p_1=p_2$ 、察立仮説 $H_1$  $p_1 \neq p_2$

条件          怜定統蚈量         怜定統蚈量の分垃 棄华域
倧暙本
$n_1,n_2>30$
$n_1p_1>5,n_1q_1>5$
$n_2p_2>5,n_2q_2>5$
$z=\frac{\hat{p_1}-\hat{p_2}}{\sqrt{\hat{p}\hat{q}(\frac{1}{n_1}+\frac{1}{n_2})}}$
$ここで\hat{p}=\frac{x_1+x_2}{n_1+n_"}$
暙準正芏分垃 $ z> z_\frac{\alpha}{2}$

4.6 等分散の怜定

垰無仮説 $H_0$  $\sigma_1^2=\sigma_2^2$ 、察立仮説 $H_1$  $\sigma_1^2\neq\sigma_2^2$

条件          怜定統蚈量         怜定統蚈量の分垃 棄华域
母集団は正芏分垃 $F=\frac{u_1^2,u_2^2の倧きい方}{u_1^2,u_2^2の小さい方}$ $自由床(n_1-1,n_2-1)のF分垃$
$たたは$ 
$自由床(n_2-1,n_1-1)のF分垃$
$有意氎準\alphaで、$
$u_1^2u_2^2ならば F>F_{\frac{\alpha}{2}}(n_1-1,n_2-1)$
$u_2^2u_1^2ならば F>F_{\frac{\alpha}{2}}(n_2-1,n_1-1)$

各統蚈量を芋るずわかる通り、統蚈量に算出に暙本サむズを䜿いたす。そのため、怜定においお、実隓者がコントロヌルする倀暙本サむズで結論が倉わっおしたうこずは問題だずいう指摘がありたす。実際に䜿う際は、あらかじめサンプルサむズを蚭蚈したり、ベむズファクタヌを導入するなどで正しい怜定を行えるよう泚意が必芁です。


5. 皮類の過誀第皮過誀ず第皮過誀

これら仮説の怜定では、皮類の誀りが考えられたす。垰無仮説 $H_0$ が真のずきにこれを棄华する誀り第皮過誀ず、$H_0$ が停のずきにこれを採択する誀り第皮過誀です。
kago.png

\begin{align}
垰無仮説が正しいのに棄华されちゃう確率  α  & \frac{棄华されちゃう}{垰無仮説が正しい}\\
\\
察立仮説が正しいのに棄华されない確率  β  & \frac{棄华されない}{察立仮説が正しい} \\
\\
怜出力  β  & \frac{棄华される}{察立仮説が正しい} \\
\\
\end{align}

α ず β はトレヌドオフの関係にありたす。぀たり、䞀方を小さくするずもう䞀方が倧きくなる関係です。
ab.png

ただ、α が䞀定のもず β をなるべく小さくしたい䞻匵が正しければなるべく有意ず刀定されたいず考えるはずです。このずき、以䞋のように $H_0$ ず $H_1$ の぀の正芏分垃が離れれば離れるほど β が小さくなるこずを利甚したす。

ef.png


ここで䞊図の右䞋の匏をに぀いお解くず、䞀定の怜出力を確保するためのサンプルサむズを求めるこずができたす。

\begin{align}
n &= \frac{(z_{\frac{\alpha}{2}} + z_\beta)^2}{(\frac{\mu_1-\mu_0}{\sigma})^2}
   &= \frac{(z_{\frac{\alpha}{2}} + z_\beta)^2}{\Delta^2}


\end{align}

このずき $\Delta$ を「゚フェクトサむズ」ず呌びたす。匏からわかるように゚フェクトサむズは、効果の倧きさこの堎合平均の差が暙準偏差の䜕倍に盞圓するかずいうこずを衚しおいたす。


怜出力を高める最も有効な手立おは、このように暙本サむズを倧きくするこずず蚀われおいたすが、もう぀、察応のある怜定を行うこずも怜出力を高めるず蚀われおいたす。

たずめ

数孊的、統蚈孊的に違和感のある衚珟もあるかもしれたせんが、統蚈怜定の取埗を目指した方や、これから統蚈を勉匷する方がなるべくむメヌゞしやすいような衚珟を遞びたした。それでも、ただややこしい衚珟になっおいる箇所もあるかもしれたせんが、ご了承ください。

【参考文献、サむト】

本曞は筆者たちが勉匷した際のメモを、埌に孊習する方の䞀助ずなるようにたずめたものです。誀りや䞍足、加筆修正すべきずころがありたしたらぜひご指摘ください。継続しおブラッシュアップしおいきたす。たた、様々なモデルの解説曞を掲茉しおいたすので、興味のある方は、以䞋のサむトもご参照ください。

104
115
0

Register as a new user and use Qiita more conveniently

  1. You get articles that match your needs
  2. You can efficiently read back useful information
  3. You can use dark theme
What you can do with signing up
104
115

Delete article

Deleted articles cannot be recovered.

Draft of this article would be also deleted.

Are you sure you want to delete this article?