統計の基礎

平均値

$n$ 個の数値 $x_1$，$x_2$，...，$x_n$ が与えられたとき，これらの値を代表する値（代表値）として最もよく使われるのが平均値（mean，average）

です。$\overline{x}$（エックス・バー）は $x$ の平均値を意味します。

平均値を意味する英語は mean と average の二通りがあります。平均値を求めるのにExcelでは average() という関数を使いますが，Rでは短いほうの mean() が関数名になっています。

> x = c(0, 7, 8, 9, 100)
> mean(x)
[1] 24.8

運動競技などの採点で，極端な点数を付ける審査員の影響を少なくするために，点数を大きさの順に並べて，両側から同数ずつ削除してから平均を求めることがあります。このような平均を，トリム平均（トリムド平均），調整平均，刈り込み平均（trimmed mean）などと呼びます。Rでは mean() に trim=... というオプションを与えて計算します。次の例では20%ずつ両側から（合わせて約40%）外します（正確には，個数×0.2 を切り捨てて整数にした個数ずつ両側から外します）。

> mean(x, trim=0.2)
[1] 8

このトリム平均を推し進めて，両側からどんどん外していき，残りが1個または2個になったときのトリム平均を，中央値またはメジアン（メディアン，median）といいます（学習指導要領では中学校1年の数学に「中央値」という用語が現れます）。つまり，$x_1\leqq x_2\leqq \cdots \leqq x_n$ となるように並べ替えて，$n$ が奇数なら $x_{(n+1)/2}$，$n$ が偶数なら $(x_{n/2}+x_{n/2+1})/2$ がメジアンです。Rでは median() という関数で求めます。

> median(x)
[1] 8

これら三つの代表値の使い分けとしては，データの分布が後で説明する正規分布に近いなら平均値，後で説明するコーシー分布のように極端な値が非常に多いならメジアンのほうが安定した結果が得られます。トリム平均は，分布にかかわらず，最善ではないけれども妥当な代表値となるものとして，近年よく使われています。

分布の形が左右非対称の場合，例えば所得（収入）の分布のような場合に，平均値を使うかメジアンを使うかで結果が大きく異なります。どちらが正しいということはなく，全員の収入をプールして均等に再分配したらどうなるかという意味では平均値ですし，「平均的な」人の収入という意味ではメジアンのほうが妥当です。

中央値は全体の50%の点です。同様に，0%の点（最小値），25%の点（第一四分位数），75%の点（第三四分位数），100%の点（最大値）もよく使われます。これらは quantile() という関数で求められます。例：

> quantile(1:10)
   0%   25%   50%   75%  100% 
 1.00  3.25  5.50  7.75 10.00 

確率変数，乱数，母集団，標本

ある変数 $X$ の値を調べようとすると，調べるたびに値が変わってしまうような，変な変数のことを，数学屋さんは確率変数（random variable）と呼び，コンピュータ屋さんは乱数（random numbers）と呼びます。確率変数は $X$ や $Y$ など大文字で表すことが多いようです。

確率変数 $X$ から $n$ 個の値 $X_1$, $X_2$, ..., $X_n$ を引き出してくることは，統計学で，非常に大きな集団（母集団，population）からランダムに $n$ 個のもの（標本，sample）を選び出してくるのと同じことだと考えられます。

確率変数 $X$ の平均値（または期待値，expectation value または expectation）とは $\overline{X}=\frac{1}{n}(X_1+X_2+\cdots +X_n)$ の $n\to \mathrm{\infty }$ での極限です。これを記号で $E(X)$ または $⟨X⟩$ と書きます。$E(X)$ の値を1文字で表したいときは，よく $\mu$（ミュー）という文字を使います。$\mu$ は英語の m に相当するギリシア文字で，mean（平均）の頭文字です。

統計学のことばでは，$\mu =E(X)$ を母平均（population mean），$n$ が有限なときの $\overline{X}=\frac{1}{n}(X_1+X_2+\cdots +X_n)$ を標本平均（sample mean）といいます。$E(X)$ は一定の値ですが，$X$ から取り出した $X_1$ や $X_2$ などの個々の値，およびそれから計算した $\overline{X}$ は，標本ごとに違う値になります。つまり，$X_1$ などや $\overline{X}$ も確率変数といえます。それらの期待値はいずれも同じ値です：

期待値については次の式が成り立ちます（$E(\cdots )$ の線形性）：

ここで $a$，$b$ は定数です。これは $X$ と $Y$ が独立でなくても成り立つ便利な式です。さらに $X$ と $Y$ が独立ならば

も成り立ちます。

分散と標準偏差

確率変数 $X$ の期待値（母平均）を $\mu =E(X)$ とするとき，$X$ と母平均との差の2乗の期待値 $E((X-\mu {)}^2)$ を $X$ の分散（variance），統計学のことばでは母分散（population variance）といいます。記号ではよく $V(X)$ と表しますが，分散の値を文字で表すときはよく ${\sigma }^2$（シグマ2乗）と表します：

この平方根

を標準偏差（standard deviation）といいます。$\sigma$（シグマ）は英語の s に相当するギリシア文字で，standard deviation の頭文字です。S.D. あるいは SD と書くこともあります。

確率変数 $X$，$Y$ が独立なら，和の分散は分散の和です。もっと一般に，$a$ と $b$ を定数として，$V(aX+bY)=a^{2}V(X)+b^{2}V(Y)$ が成り立ちます。特に，$V(X\pm Y)=V(X)+V(Y)$ です。標準偏差 $\sigma$ に焼き直せば，${\sigma }_{X\pm Y}^2={\sigma }_X^2+{\sigma }_Y^2$ という「三平方の定理」が成り立ちます。

$X$ から引き出した $n$ 個の値 $X_1$，$X_2$，...，$X_n$ については，

ですが，母集団の平均値 $\mu$ は実際には知ることのできない値ですので，この代わりに標本の平均値 $\overline{X}=\frac{1}{n}(X_1+\cdots +X_n)$ を使うと，$\mu$ を使う場合より小さめの値になります（実際，$\overline{X}$ は $(X_1-\overline{X}{)}^2+\cdots +(X_n-\overline{X}{)}^2$ を最小にします）。後で証明するように，

です。そこで，

と置くと，$s^2$ の期待値は母分散に一致します（これは $X$ の分布によらず成り立ちます。正規分布などは仮定していません）：

期待値が母分散に一致するということは，偏りがないということなので，$s^2$ は母分散 ${\sigma }^2$ の不偏推定量（an unbiased estimator of the population variance ${\sigma }^2$）であるといいます。あるいは，これを略して不偏分散（unbiased variance）ということもあります。

数学的な話が続いてしまいましたが，要するに，$X$ の平均 $\mu$ が最初からわかっている場合は，分散は

で求め，同じデータから求めた平均 $\overline{X}=(X_1+X_2+\cdots )/n$ を使う場合は

で求めるのが正しいのです。しかし，高校の数学の教科書では $n-1$ の説明が難しいのでこれらを区別せずにつねに $n$ で割っています。違いは小さいので，目くじらをたてるほどのことではありませんが，大学ではちゃんと区別しましょう。以下では高校教科書の流儀を「$n$ で割る分散」，正しい流儀を「$n-1$ で割る分散」と呼ぶことにします。

Excelには $n$ で割る分散を求める varp() という関数があります（Excel 2010以降は var.p() も同義）。$n-1$ で割る分散はExcelでもRでも var() という関数で求めます（Excel 2010以降は var.s() も同義）。Excelの varp() に相当するものはRにはありません（だってそんなものを統計学者は使わないので）が，次のように打ち込めば作ることができます：

> varp = function(x) { var(x) * (length(x)-1) / length(x) }

これらを試してみましょう：

> x = 1:3   # x = c(1,2,3) でも同じ
> var(x)
[1] 1
> varp(x)
[1] 0.6666667

> x = 1:10
> var(x)                   # n-1で割る分散
[1] 9.166667
> mean(combn(x, 4, var))   # xから4個選んだ組合せのvar()の平均
[1] 9.166667
> mean(combn(x, 6, var))   # xから6個選んだ組合せのvar()の平均
[1] 9.166667
> varp(x)                  # nで割る分散
[1] 8.25
> mean(combn(x, 4, varp))  # xから4個選んだ組合せのvarp()の平均
[1] 6.875
> mean(combn(x, 6, varp))  # xから6個選んだ組合せのvarp()の平均
[1] 7.638889

標準偏差 $\sigma =\sqrt{{\sigma }^2}$ を標本から求める場合には，一般に次の式を使います。

$E(s^2)={\sigma }^2$ でしたが，一般に $E(s)\ne \sigma$ です。$E(s)$ と $\sigma$ がどれくらい違うかは分布の形によります。

Rで標準偏差 $s$ を求めるには sd() を使います。

> x = 1:10
> sd(x)
[1] 3.027650

確率変数 $X$ の平均値 $\mu =E(X)$，分散 ${\sigma }^2=E((X-\mu {)}^2)$ がわかっているとします。$X$ から引き出した数を $X_1$, $X_2$, ..., $X_n$ とすると，その平均値

の期待値は $\mu$ に一致します。一方，$\overline{X}$ の分散 $E((\overline{X}-\mu {)}^2)$ は，$\mu =0$ の場合だけ考えればいいので，

となり，$\overline{X}$ の分布は平均値が元と同じ $\mu$ で，分散が ${\displaystyle \frac{{\sigma }^2}{n}}$，したがって標準偏差は ${\displaystyle \frac{\sigma }{\sqrt{n}}}$ になります。

一般に，$\overline{X}$ のような統計量の標準偏差を標準誤差といいます。測定回数 $n$ を増やすと，平均値 $\overline{X}$ の標準誤差は $1/\sqrt{n}$ に比例して小さくなります。これが，できるだけ何度も測定して平均をとるべき理由です。測定回数が100倍になれば，標準誤差は 1/10 になります。

中心極限定理と正規分布

一番簡単な乱数は，ある範囲の数がまんべんなく出る一様乱数（uniform random numbers）です。Rの関数 runif(n) は，$0<x<1$ の範囲の一様乱数を $n$ 個取り出します（Rはデフォルトではメルセンヌ・ツイスターというアルゴリズム（日本の松本眞先生らによって考案されました）で乱数を発生します。RNGkind() で変更できます）。

> runif(1)
[1] 0.388267
> runif(10)
 [1] 0.2146394 0.2765450 0.5433135 0.4784538 0.8147103 0.1141375
 [7] 0.6488306 0.7947468 0.1698610 0.2440027

もっとたくさんの数，例えば百万個でやってみましょう。

> x = runif(1000000)  % 百万個の乱数
> hist(x, freq=FALSE) % ヒストグラム（度数分布図）を描く

ヒストグラム（度数分布図）を描く hist() を使いました。オプション freq = FALSE は縦軸の目盛りを，デフォルトの個数（frequency）でなく，密度（density）にすることを意味します。密度とは，棒グラフの棒の面積の和が 1 になるように付けた目盛りです。下で説明する確率密度関数の縦軸と同じものになります。

このヒストグラムを関数で表したものを確率密度関数（probability density function，略してp.d.f.）または単に密度関数といいます。数学的にいえば，任意の $a<b$ について，$a\leqq x\leqq b$ の範囲に乱数が入る確率が ${\int }_a^{b}f(x)dx$ になるように選ばれた関数 $f(x)$ が確率密度関数です。runif() の確率密度関数は

と表すことができます。

この runif() を二つ加えたらどうなるでしょうか。2倍するのではなく，独立な乱数を二つ加えるのです。

> x = runif(1000000) + runif(1000000)
> hist(x, freq=FALSE)

三個ならどうなるでしょうか。

> x = runif(1000000) + runif(1000000) + runif(1000000)
> hist(x, freq=FALSE)

だんだん釣り鐘の形になってきました。

runif() は平均値が0.5ですので，0.5を引いて範囲を $-0.5\leqq x\leqq 0.5$ にして考えましょう。この分散は

ですので，12個加えるとちょうど分散が1になります。

次の図が，runif() - 0.5 を12個加えた分布のヒストグラムです。

滑らかな曲線は

を描いたものです。ここで $\exp (z)=e^z$ は指数関数です。$e=2.718\dots$ は自然対数の底（てい）です。${\int }_{-\mathrm{\infty }}^{\mathrm{\infty }}\exp (-x^2/2)dx=\sqrt{2\pi }$ ですので，これは確率密度関数です。この分布を標準正規分布（standard normal distribution）といいます。

> integrate(function(x){exp(-x^2/2)}, -Inf, +Inf)
2.506628 with absolute error < 0.00023
> sqrt(2*pi)
[1] 2.506628

乱数を足し合わせる上の実験は，一様分布の乱数 runif() から出発しましたが，どんな分布から出発しても，まったく同じことが成り立ちます。

数学的に言えば，平均 $\mu$，分散 ${\sigma }^2$ の確率変数 $X$ がどんな分布であっても，そこから引き出した数 $X_1$, $X_2$, ..., $X_n$ の平均値

の分布は，平均値が元と同じ $\mu$ で，分散が ${\sigma }^2/n$ になりますので，

の分布は平均 0，分散 1 になります。ここまでは $n$ の値にかかわらず言えることです。ここで，$n$ が十分大きくなると，この分布は標準正規分布に近づきます。このことを中心極限定理といいます。

$X$ が標準正規分布ならば，$\sigma X+\mu$ は平均 $\mu$，分散 ${\sigma }^2$ の正規分布（normal distribution）になります。平均 $\mu$，分散 ${\sigma }^2$ の正規分布を $N(\mu ,{\sigma }^2)$ と表します。標準正規分布は平均 0，分散 1 の正規分布ですので $N(0,1)$ と表されます。

なお，正規分布を研究した数学者ガウス（Carl Friedrich Gauss，1777-1855）に敬意を表して，正規分布のことをガウス分布（Gaussian distribution）ともいいます。

変数がある範囲に入る確率を求めるには，密度関数を積分しなければなりません。そこで，あらかじめ密度関数 $f(x)$ を積分した

を求めておくと便利です。この $F(q)$ を累積分布関数（cumulative distribution function，略してCDF）あるいは単に分布関数と呼びます。$F(q)$ がわかっていれば，確率変数 $X$ が $a\leqq X\leqq b$ の範囲に入る確率は $F(b)-F(a)$ で求められます。

分布関数 $p=F(q)$ の逆関数 $q=F^{-1}(p)$ を分位関数または分位点関数（quantile function）と呼びます。

Rには標準正規分布 $N(0,1)$ の

• 密度関数 dnorm(x)
• 分布関数 pnorm(q) $={\int }_{-\mathrm{\infty }}^q\mathrm{dnorm}(x)dx$
• 分位関数 qnorm(p)
• 乱数を $n$ 個発生する rnorm(n)

があります。下の図の色を付けた部分の面積が dnorm(x) です。

x = (-35:15)/10
y = dnorm(x)
plot(NULL, xlim=c(-3.5,3.5), ylim=c(0,0.4), xlab="", ylab="")
polygon(c(x,x[length(x)],x[1]), c(y,0,0), col="gray")
curve(dnorm, lwd=2, add=TRUE)

これらは，dnorm(x, mean=μ, sd=s) のように平均と標準偏差を与えることもできます。省略すれば dnorm(x, mean=0, sd=1) の意味になります。

正規分布の乱数を正規乱数といいます。

物理現象の測定では，測定誤差がほぼ正規分布することがよくあります。これに対して，他の多くの分野では，観測値そのものが正規分布することはまずありませんが，中心極限定理のおかげで，観測値の平均値はほぼ正規分布します。

正規分布 $N(\mu ,{\sigma }^2)$ に従う変数 $X$ が平均±標準偏差の範囲内に入る確率，つまり $\mu -\sigma <X<\mu +\sigma$ になる確率は，標準正規分布 $N(0,1)$ に従う変数 $Z$ が $-1<Z<1$ になる確率に等しく，Rで次のようにして求められます。

> pnorm(1) - pnorm(-1)
[1] 0.6826895

正規分布が左右対称であることを使えば，次のようにするほうが簡単です。

> 1 - 2 * pnorm(-1)
[1] 0.6826895

つまり，ほぼ68%の確率で正規分布は $\mu \pm \sigma$ の範囲に入ります。

同様にして，$\mu \pm 2\sigma$，$\mu \pm 3\sigma$ に入る確率を求めてください。

逆に，$\mu \pm k\sigma$ の範囲に入る確率が与えられたときに $k$ の値を求めるには qnorm() を使います。例えば 95% の場合は

> qnorm(0.025)
[1] -1.959964

ですので，ほぼ $\mu \pm 1.96\sigma$ です。ざっくり $\mu \pm 2\sigma$ で 95% と覚えておくと便利です。

現実のテストの分布は正規分布といえるでしょうか。2015年全国学力テスト中学理科の正答数の度数分布をrika_hist.csv（SJIS）・rika_hist-utf8.csv（UTF-8）として置いておきましたので，ヒストグラムを描き，平均と分散の同じ正規分布の密度関数を重ね書きしてみてください。

Scrapboxのほうに書いた正規分布の何がいいの？もご参照ください。

本当に困る分布

どれだけ正規分布から離れていても，分散が有限であれば，中心極限定理が成り立つので，その平均値は正規分布に近づきます。

しかし，世の中には中心極限定理の前提が満たされない本当に困る分布もあります。その有名な例が，密度関数が $1/(1+a{x}^2)$ に比例するコーシー分布（Cauchy distribution）です。

この分散を計算しようとしても，$x^2/(1+a{x}^2)$ の積分は発散してしまいます。そもそも $x/(1+a{x}^2)$ も積分できません。

Rには密度関数が $1/(1+x^2)$ に比例するコーシー分布の

• 密度関数 dcauchy(x)
• 分布関数 pcauchy(q) $={\int }_{-\mathrm{\infty }}^q\mathrm{dcauchy}(x)dx$
• 分位関数 qcauchy(p)
• 乱数を $n$ 個発生する rcauchy(n)

があります。

実験をしてみましょう。100万個のコーシー分布の乱数を作ります。

x = rcauchy(1000000)

この度数分布を表示するために hist(x) としても，まともに表示されません。x = sort(x) で並べ替えて，最初の数個を head(x) で表示させたり，最後の数個を tail(x) で表示させたりすると，6〜7桁のとんでもない値が含まれることがわかります。そのため，平均 mean(x) を求めると，0からかなり外れてしまいます。このような分布は，分散はもとより，平均値を求めることにも意味がありません。強いて求めるなら，メジアンが最適です。

分散についても，平均を求める操作をメジアンに焼き直したMAD（median absolute deviation）と呼ばれるものを使えば，このような分布でも求めることができます。これは，データからメジアンを引いたものの絶対値のメジアン，つまりRで書くと

median(abs(x - median(x)))

になります。

ところで，標準正規分布の median(abs(x - median(x))) は qnorm(0.75) ≒ 1/1.4826 なので，正規分布の場合に標準偏差と比較しやすくするために，MADの定義を1.4826倍しておくことがあります。Rで mad(x) で求められる値はデフォルトでは

1.4826*median(abs(x - median(x)))

になっています。

MAD以外に，四分位範囲（interquartile range）もよく使われます。これは，順位が全体の1/4と3/4の値の差です。Rでは IQR() で求められます。補正しないMADのほぼ2倍に相当します。

正規分布に比べて外れ値（outliers）が非常に多い場合の分布のモデルにコーシー分布が使われることがあります。

正規分布から導かれる分布

以下の分布は頻繁に現れるので，Rの練習にもなりますので，乱数を発生させて度数分布を描き，理論的な密度関数と比べてみてください。

カイ2乗分布

$\nu$ 個の数値 $X_1$, $X_2$, ..., $X_{\nu }$ が標準正規分布に従うとき，

の分布を「自由度 $\nu$ のカイ2乗分布（${\chi }^2$ 分布）」（chi-square distribution with $\nu$ degrees of freedom）と呼びます（2乗は「にじょう」以外に「じじょう」と読むこともあります。これはもともと「自乗」（自分を乗算する）から来ています）。記号では ${\chi }^2(\nu )$ と書きます。$\chi$ はギリシア文字「カイ」（chi）です。自由度を表すのによく使う文字 $\nu$ はギリシア文字「ニュー」（nu）です。Rのマニュアルでは $\nu$ ではなく df と書かれています。

自由度 $\nu$ の ${\chi }^2$ 分布は，平均 $\nu$，分散 $2\nu$ で，その密度関数 $f(x)$ は $x^{\nu /2-1}{e}^{-x/2}$ に比例します。

${\chi }^2$ は $X_1^2$, $X_2^2$, ... という独立な確率変数の和ですので，中心極限定理のため $\nu$ が大きいと正規分布に近づきます。

Rには自由度 $\nu$ の ${\chi }^2$ 分布の

• 密度関数 dchisq(x, ν)
• 分布関数 pchisq(q, ν) $={\int }_0^q\mathrm{dchisq}(x,\nu )dx$
• 分位関数 qchisq(p, ν)
• 乱数を $n$ 個発生する rchisq(n, ν)

が備わっています。

$X_1$, $X_2$, ..., $X_n$ が正規分布 $N(\mu ,{\sigma }^2)$ に従うとき，標本分散

の分布は（定数倍を除いて）${\chi }^2$ 分布です。正確には，$(n-1)s^2/{\sigma }^2$ が自由度 $n-1$ の ${\chi }^2$ 分布になります。

t 分布

$X$ が標準正規分布 $N(0,1)$ に従い，$Y$ が自由度 $\nu$ の ${\chi }^2$ 分布に従うとき，

の分布を，自由度 $\nu$ のスチューデントの $t$ 分布（Student's $t$ distribution with $\nu$ degrees of freedom）あるいは単に自由度 $\nu$ の $t$ 分布と呼びます（Student はゴセット（W. S. Gosset，1876〜1937年）のペンネームです）。

自由度 $\nu$ の $t$ 分布の密度関数 $f(t)$ は $(1+t^2/\nu {)}^{-(\nu +1)/2}$ に比例します。

特に，$X_1$, $X_2$, ..., $X_n$ が $N(\mu ,{\sigma }^2)$ に従うとき，${\sigma }^2$ の値にかかわらず，

は自由度 $n-1$ の $t$ 分布に従います。

tt = function() {
    x = sample(1:5, 50, replace=TRUE)
    (mean(x) - 3) / sqrt(var(x) / 50)
}
t = replicate(1000000, tt())

> t0 = -qt(0.025, 50-1)
> mean(abs(t) > t0)
[1] 0.050288

Rには自由度 $\nu$ の $t$ 分布の

• 密度関数 dt(x, ν)
• 分布関数 pt(q, ν) $={\int }_{-\mathrm{\infty }}^q\mathrm{dt}(x,\nu )dx$
• 分位関数 qt(p, ν)
• 乱数を n 個発生する rt(n, ν)

が備わっています。

F 分布

$u_1$ と $u_2$ が独立な確率変数で，それぞれ自由度 ${\nu }_1$，${\nu }_2$ の ${\chi }^2$ 分布に従うとき，

の分布を，自由度 $({\nu }_1,{\nu }_2)$ の $F$ 分布（$F$ distribution with $({\nu }_1,{\nu }_2)$ degrees of freedom）といいます。

自由度 $({\nu }_1,{\nu }_2)$ の $F$ 分布の密度関数 $f(F)$ は $F^{{\nu }_1/2-1}(1+{\nu }_{1}F/{\nu }_2{)}^{-({\nu }_1+{\nu }_2)/2}$ に比例します。

特に，$X_1$, ..., $X_n$, $Y_1$, ..., $Y_m$ が $N(\mu ,{\sigma }^2)$ に従うとき，$X_1$, ..., $X_n$ の標本分散を $s_1^2$，$Y_1$, ..., $Y_m$ の標本分散を $s_2^2$ とすれば，

は自由度 $(n-1,m-1)$ の $F$ 分布に従います。

Rには自由度 $({\nu }_1,{\nu }_2)$ の $F$ 分布の

• 密度関数 df(x, ν1, ν2)
• 分布関数 pf(q, ν1, ν2) $={\int }_0^q\mathrm{df}(x,{\nu }_1,{\nu }_2)dx$
• 分位関数 qf(p, ν1, ν2)
• 乱数を $n$ 個発生する rf(n, ν1, ν2)

が備わっています。

仮説検定と区間推定

高校「数学B」では，おもに正規分布に基づいて，仮説検定や区間推定の考え方を学びます。

例えば，歪んでいない硬貨を100回投げて表（おもて）が $X$ 回出たとすると，$X$ は2項分布 $B(100,0.5)$ に従う確率変数です。2項分布 $B(n,\theta )$ の平均は $n\theta$，分散は $n\theta (1-\theta )$ ですので，この2項分布の平均は $\mu =100\times 0.5=50$，分散は ${\sigma }^2=100\times 0.5\times (1-0.5)=25$，標準偏差は $\sigma =\sqrt{25}=5$ です。ここまでは厳密ですが，ここでこの2項分布を，これと同じ平均・分散の正規分布 $N(50,25)$ で近似します。正規分布は平均を中心とする $\pm 1.96\sigma$ の範囲に 95% が入ります。この範囲外の事象が起きる確率は 5% 以下ですので，もしそのような事象が起きたら，「硬貨は歪んでいない」という最初の仮定（帰無仮説）は有意水準 5% で棄却されると考えます。つまり，100回投げて表が40回以下または60回以上出たら，硬貨は歪んでいると結論します。帰無仮説が棄却されない場合は，結論を保留します（歪んでいないと結論するのではありません）。

この「棄却する」という考え方は，誤解されやすいので，注意が必要です。例えば平均から $1.96\sigma$ 以上離れているデータは異常な「外れ値」なので削除するという考え方は間違いです。偏差値 30 以下や 70 以上の人を抹殺してはいけません。何度も繰り返せば，必ず $\mu \pm 1.96\sigma$ の範囲外の事象は起こります。

実際には，仮説検定を適用するのは，たいへん費用がかかるなどの理由で何度も繰り返すことができない場合に限られます。例えば新型コロナのワクチンが効くかどうかを調べるには大規模な治験が必要ですので，あらかじめ「まったく効かないという帰無仮説から（正規分布と仮定した場合）$1.96\sigma$ 以上離れた結果が出たら効果ありと判断する」といったルールを決めて治験を行います。

言い換えれば，データを中心として $\pm 1.96\sigma$ 以内の範囲に帰無仮説があれば，5% 水準で帰無仮説は棄却されません。この範囲を 95% 信頼区間といいます——という説明では雑すぎるので，以下でもうちょっとちゃんと説明します。

高校「数学B」でよく扱われるのは，データが正規分布 $N(\mu ,{\sigma }^2)$ に従うと仮定して，平均 $\mu$ は未知，分散 ${\sigma }^2$ は既知で $\mu$ によらず一定とする場合です。この正規分布から大きさ $n$ の標本 $x_1,x_2,\dots ,x_n$ を取り出してその平均値 $\overline{x}$ を求めます。この $\overline{x}$ は正規分布 $N(\mu ,{\sigma }^2/n)$ に従うはずですから，$\mu -1.96\sigma /\sqrt{n}\le \overline{x}\le \mu +1.96\sigma /\sqrt{n}$ を $\overline{x}$ が満たす確率は 0.95 です。この式を変形して $\overline{x}-1.96\sigma /\sqrt{n}\le \mu \le \overline{x}+1.96\sigma /\sqrt{n}$ を $\overline{x}$ が満たす確率は 0.95 と言っても同じことです。ここで現れる区間 $[\overline{x}-1.96\sigma /\sqrt{n},\overline{x}+1.96\sigma /\sqrt{n}]$ を $\mu$ の 95% 信頼区間といいます。$\mu$ は 95% の確率でこの中に入りますが，確率変数は $\mu$ ではなく区間のほうです。

上の議論では，データが正規分布に従うと仮定しました。でも，現実のデータは正規分布に従うわけではありません。それなら上の議論は使えないかというと，そうではありません。データが正規分布から外れていても，そこからとった大きさ $n$ の標本の平均値 $\overline{x}$ は $n$ がある程度大きければほぼ正規分布に従います（中心極限定理）。上の議論は $\overline{x}$ が $N(\mu ,{\sigma }^2/n)$ に従えば十分なので，データが正規分布という仮定を外しても近似としては十分正確です。

上の議論のもっと大きな問題点は，${\sigma }^2$ すなわち母分散（母集団の分散）が既知と仮定しているところです。実際には母分散はわからないので，データから計算した分散の推定値で間に合わせることになります。すると，推定した母分散にも誤差が含まれますので，実際の 95% 信頼区間は上で計算したものより若干広くなることが考えられます。

このあたりの議論をもうちょっと詳しくすれば，このような場合に信頼区間を求めるには，正規分布ではなく，それより少し裾（すそ）の広い $t$ 分布を使えばいいことがわかっています。このため，高校数学以外では，データが正規分布に従う場合でも，この $t$ 分布を使って検定します。正規分布を使って検定することはまずありません。

R には $t$ 分布を使った検定や区間推定をするための t.test() という関数が備わっています。例えば大きさ 5 のサンプル $\{1,2,3,4,5\}$ が得られたとき，母集団の平均値の 95% 信頼区間を求めるには，次のようにします。

> t.test(c(1, 2, 3, 4, 5))

	One Sample t-test

data:  c(1, 2, 3, 4, 5)
t = 4.2426, df = 4, p-value = 0.01324
alternative hypothesis: true mean is not equal to 0
95 percent confidence interval:
 1.036757 4.963243
sample estimates:
mean of x 
        3 

データの平均値は $\overline{x}=3$ で，母集団の平均値 $\mu$ の 95% 信頼区間は [1.036757, 4.963243] です。なお，p-value = 0.01324 と出力されますが，これは母集団の平均値が 0 という帰無仮説を設定した場合の $p$ 値ですので，ここでは無視します。

なお，仮説検定で t.test() がよく使われるのは次のような場合です。集団 X から大きさ 5 の標本を取り出したら $\{1,2,3,4,5\}$ でした。集団 Y から大きさ 6 の標本を取り出したら $\{3,4,5,6,7,8\}$ でした。分散については何の仮定もしません。両集団の平均値は等しいという帰無仮説を検定してください。

> x = c(1, 2, 3, 4, 5)
> y = c(3, 4, 5, 6, 7, 8)
> t.test(x, y)

	Welch Two Sample t-test

data:  x and y
t = -2.4019, df = 8.9894, p-value = 0.0398
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
 -4.8549526 -0.1450474
sample estimates:
mean of x mean of y 
      3.0       5.5 

これは $t$ 検定に基づいてウェルチ（Welch）という人が考えた方法を使っています。「両群の母分散が等しいと仮定しない場合の $t$ 検定」ということもあります。$p$ 値は 0.0398 ですので，水準 0.05 で両集団の母平均は有意に異なります。また，両集団の母平均の差の 95% 信頼区間は [-4.8549526, -0.1450474] です。

なお，X，Y 両群の母分散が何らかの理由で等しいと仮定できる場合は，次のようにします。

> t.test(x, y, var.equal=TRUE)

	Two Sample t-test

data:  x and y
t = -2.3619, df = 9, p-value = 0.04247
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
 -4.8944379 -0.1055621
sample estimates:
mean of x mean of y 
      3.0       5.5 

結果がほんの少し変わりました。通常は母分散が等しいと仮定できる理由はないので，先ほどの Welch の検定のほうを使います。

より詳しくは t検定 のページをご覧ください。

奥村 晴彦

Last modified: 2022-05-29 17:50:45 JST