考査１_解説
考査のポイント.
分散の重要性.
標準偏差 sd は 分散の平方根です.
- 平均±２(標準偏差)の範囲にデータの95%が収まり.
- 平均±３(標準偏差)の範囲にデータの99%が収まります.
つまり、 (m,sd)=(平均,標準偏差)
が分かれば、データの全体像が分かります。
このことを理解するには R を使ったシミュレーションが有効です.
Google Colabo でシミュレーションしてみましょう.
Google Colabo で検索すると、「Colabo
へようこそ」という画面になります。
「ノートブックを新規作成」をクリックします。
ランタイムから「ランタイムのタイプを変更」をクリックし、Rに変更します.
＋コードをクリックするとコード入力欄が現れます.
「3章 分散」に関連した内容を練習しましょう.
rnorm(個数,平均,標準偏差)
rnorm(10, 5, 2)=rnorm(個数,平均,標準偏差) は平均5,標準偏差2の正規乱数を10個出力します
rnorm(10, 5, 2) ## [1] 3.841037 7.232155 4.984596 3.384862 7.498565 7.556739 3.208485 3.598070
## [9] 9.461185 6.725545個数を100個に増やし、hist関数でヒストグラムを描きます.
var(a)で分散を計算します
a <- rnorm(100, 5, 2)
hist(a)
var(a)## [1] 4.714295個数が増えれば、分散は理論値４に近くなります。
課題
- 平均や標準偏差を変えてみましょう.
- 個数を1000個に増やしてみましょう.
混合分布.
混合分布について考えましょう。
A：平均３、標準偏差１の正規乱数３０個.
B：平均５、標準偏差１の正規乱数３０個.
Google Colab にコピー・ペーストして実行してみましょう。
A <- rnorm(30,3,1)#平均3,標準偏差1,個数30の乱数
B <- rnorm(30,5,1)#平均5,標準偏差1,個数30の乱数
U <- c(A,B)#A,B合わせた混合分布
par(mfrow=c(3,1))
hist(A,prob=T,breaks=c(seq(-0.35,13.5,by=0.5)),col="#ff00ff40",ylim=c(0,1))
hist(B,prob=T,breaks=c(seq(-0.35,13.5,by=0.5)),col="#0000ff40",ylim=c(0,1))
hist(U,prob=T,breaks=c(seq(-0.35,13.5,by=0.5)),col="grey88",ylim=c(0,1))
課題
- A群の個数を60個にしてみよう.
- A群の標準偏差を２にしてみよう.
Google Colab にコピー・ペーストして実行してみよう。
A <- rnorm(60,3,2)#平均3,標準偏差1,個数30の乱数
B <- rnorm(30,5,1)#平均5,標準偏差1,個数30の乱数
U <- c(A,B)#A,B合わせた混合分布
par(mfrow=c(3,1))
hist(A,prob=T,breaks=c(seq(-5.5,13.5,by=0.5)),col="#ff00ff40",ylim=c(0,1))
hist(B,prob=T,breaks=c(seq(-5.5,13.5,by=0.5)),col="#0000ff40",ylim=c(0,1))
hist(U,prob=T,breaks=c(seq(-5.5,13.5,by=0.5)),col="grey88",ylim=c(0,1))
まとめ.
統計学の第一の目的は、母集団全体の分布を把握することです.
実際の調査は各県ごとに行われ、それらを合わせた混合分布が母集団ということになります。
データを適切に理解するには、
各群の分布と混合分布の分布がどのような関係になるかを理解しておくことが必要になります。
研究「混合分布の平均・分散」.
２群 \(U_A, U_B\) について、
\(U_A\) : \(n_A\)個 平均 \(\mu_A\) 標準偏差 \(\sigma_A\)
\(U_B\) : \(n_B\)個 平均 \(\mu_B\) 標準偏差 \(\sigma_B\)
\(U=U_A\cup U_B\) ：\(n\)個 平均 \(\mu\) 標準偏差 \(\sigma\) について、
混合率 \(p_A=\dfrac{n_A}{n}~,~~p_B=\dfrac{n_B}{n}~~~\)とすると、
\(\mu=p_A~\mu_A+p_B~\mu_B\)
\(\sigma^2=p_A~\sigma_A^2+p_B~\sigma_B^2+p_Ap_B(\mu_A-\mu_B)^2\).
中間考査 問４では、
\(n_A=4,~~n_B=4,~~n=8\) より、\(p_A=\dfrac{1}{2},~~p_B=\dfrac{1}{2}\).
\(p_A~\mu_A+p_B~\mu_B=\dfrac{1}{2}\cdot 2+\dfrac{1}{2}\cdot 3=2.5\)
\(p_A~\sigma_A^2+p_B~\sigma_B^2+p_Ap_B(\mu_A-\mu_B)^2\\=\dfrac{1}{2}\cdot 1+\dfrac{1}{2}\cdot 2+\dfrac{1}{4}\cdot (2-3)^2=1.75\).
課題
中間考査 問５において、この方法で計算しましょう。
混合分布については、2つのグループの差を統計的に判断する方法として、将来、分散分析の中で学ぶことになります。
考査解答.
1. 次の各問について平均と分散を計算しなさい。
(1) x=1,7.
\(m=\dfrac{1+7}{2}=4\) \(V=\dfrac{(1-4)^2+(7-4)^2}{2}=9\)
(2) x=3,4,5.
\(m=\dfrac{3+4+5}{3}=4\) \(V=\dfrac{(3-4)^2+(4-4)^2+(5-4)^2}{3}=\dfrac{2}{3}\)
(3) x=2,3,3,4.
\(m=\dfrac{2+3+3+4}{4}=3\) \(V=\dfrac{(2-3)^2+(3-3)^2+(3-3)^2+(4-3)^2}{4}=\dfrac{1}{2}\)
2.
次の各問について、第１四分位数、第２四分位数、第3四分位数を求めなさい。また、四分位範囲、四分位偏差 を求めなさい。
(1) \(A=\{1,2,2,4,4,5,5,5,6,6,7\}\)
第１四分位数 \(Q_1=\)（ 2 ） 第２四分位数 \(Q_2=\)（ 5 ） 第３四分位数 \(Q_3=\)（ 6 ）
四分位範囲＝（ 4 ） 四分位偏差＝（ 2 ）
(2) \(A=\{1,2,2,4,4,5,5,6,6,7\}\)
第１四分位数 \(Q_1=\)（ 2 ） 第２四分位数 \(Q_2=\)（ 4.5 ） 第３四分位数 \(Q_3=\)（ 6 ）
四分位範囲＝（ 4 ） 四分位偏差＝（ 2 ）
3. 表を利用して、平均と分散を計算しなさい。
(1)
| 階級値 | 度数 | 階級値・度数 | 偏差 | 偏差平方・度数 |
|---|---|---|---|---|
| 10 | 1 | 10 | -20 | 400 |
| 20 | 2 | 40 | -10 | 200 |
| 30 | 3 | 90 | 0 | 0 |
| 40 | 2 | 80 | 10 | 200 |
| 50 | 1 | 50 | 20 | 400 |
| Total | 9 | 270 | 1200 |
平均 \(m=\dfrac{270}{9}=30\) 分散 \(V=\dfrac{1200}{9}=\dfrac{400}{3}\)
(2)
| 階級値 | 度数 | 階級値・度数 | 偏差 | 偏差平方・度数 |
|---|---|---|---|---|
| 30 | 8 | 240 | -15 | 1800 |
| 40 | 12 | 480 | -5 | 300 |
| 50 | 14 | 700 | 5 | 350 |
| 60 | 4 | 240 | 15 | 900 |
| 70 | 2 | 140 | 25 | 1250 |
| Total | 40 | 1800 | 4600 |
平均 \(m=\dfrac{1800}{40}=45\) 分散 \(V=\dfrac{4600}{40}=115\)
4. 次の各問に答えなさい.
(1) \(A=\{1 ,3,1,3\}\)について、平均と分散を計算しなさい。
\(m=\dfrac{1+3+1+3}{4}=2\) \(V=\dfrac{(1-2)^2+(3-2)^2+(1-2)^2+(3-2)^2}{4}=1\)
(2) \(B=\{1,3,3,5\}\)について、平均と分散を計算しなさい。
\(m=\dfrac{1+3+3+5}{4}=3\) \(V=\dfrac{(1-3)^2+(3-3)^2+(3-3)^2+(5-3)^2}{4}=2\)
(3) 研究問題.
A,Bを合わせた8個のデータ\(U=\{1 ,3,1,3,1,3,3,5\}\)について、平均と分散を計算しなさい。
| Xの値 | 度数 | 階級値・度数 | 偏差 | 偏差平方・度数 |
|---|---|---|---|---|
| 1 | 3 | 3 | -1.5 | 6.75 |
| 3 | 4 | 12 | 0.5 | 1 |
| 5 | 1 | 5 | 2.5 | 6.25 |
| Total | 8 | 20 | 14 |
平均 \(m=\dfrac{20}{8}=2.5\)
分散 \(V=\dfrac{14}{8}=\dfrac{7}{4}=1.75\)
５. 次の各問に答えなさい.
(1) \(A=\{1,2,4,7\}\)について、平均と分散を計算しなさい。
\(m=\dfrac{1+2+4+7}{4}=3.5\)
\(V=\dfrac{(1-3.5)^2+(2-3.5)^2+(4-3.5)^2+(7-3.5)^2}{4}=\dfrac{21}{4}=5.25\)
(2) \(B=\{5,8,3,6\}\)について、平均と分散を計算しなさい。
\(m=\dfrac{5+8+3+6}{4}=\dfrac{11}{2}=5.5\)
\(V=\dfrac{(5-5.5)^2+(8-5.5)^2+(3-5.5)^2+(6-5.5)^2}{4}=\dfrac{13}{4}=3.25\)
(3) 研究問題.
A,Bを合わせた12個のデータ\(U=\{1,2,4,7,5,8,3,6\}\)について、平均と分散を計算しなさい。
（同じ値がない場合は、表を使わないほうがよい）
平均 \(m=\dfrac{36}{8}=4.5\)
分散 \(V=\dfrac{21}{4}=\dfrac{21}{4}=5.25
\)
混合分布のヒストグラム.
- 考査 問４、問５のヒストグラムは次の図のようになります。
問４のヒストグラム
a=c(1,3,1,3)
A <- rep(a,1)
b=c(1,3,3,5)
B <- rep(b,1)
U <- c(A,B)
U## [1] 1 3 1 3 1 3 3 5par(mfrow=c(3,1))
hist(A,prob=T,breaks=c(seq(-3.5,5.5,by=1)),col="#ff00ff40",ylim=c(0,1))
hist(B,prob=T,breaks=c(seq(-3.5,5.5,by=1)),col="#0000ff40",ylim=c(0,1))
hist(U,prob=T,breaks=c(seq(-3.5,5.5,by=1)),col="grey88",ylim=c(0,1))
問5のヒストグラム
a=c(1,2,4,7)
A <- rep(a,1)
b=c(5,8,3,6)
B <- rep(b,1)
U <- c(A,B)
U## [1] 1 2 4 7 5 8 3 6par(mfrow=c(3,1))
hist(A,prob=T,breaks=c(seq(-5.5,13.5,by=0.5)),col="#ff00ff40",ylim=c(0,1))
hist(B,prob=T,breaks=c(seq(-5.5,13.5,by=0.5)),col="#0000ff40",ylim=c(0,1))
hist(U,prob=T,breaks=c(seq(-5.5,13.5,by=0.5)),col="grey88",ylim=c(0,1))