質問: sd / √n =√(2/3) / √3 =√(2/9) 3は標準偏差の公式に当てはめました。 4は標準誤差の公式に当てはめました。 計算間違いあるいは値の代入間違いでなければ、 3の標準誤差と4の標準誤差は同じになるはずだと思うのですが、 なっておりません。 何故なのでしょうか? 3の計算方法が間違っているのでしょうか? お教えていただければ幸いです。 以上、宜しくお願いいたします。
回答: 標準誤差は3のようにして求めるものではないです。 標準誤差は4にあるような s.e. = s.d. / sqrt(n) という計算式によって"推定"されるものです。 確かに、母集団から何回もサンプリングをしてそれぞれの平均値の分布の標準偏差を求めれば、それが標本誤差であるといえます。しかし、どんなにがんばっても、3のような方法(何回もサンプリングして平均値の分布の標準偏差を求める)と4の方法(1つの標本から標本誤差の公式から求める)とで値が一致することはありません。 これはコンピュータで実際にやってみれば実感できます。 3の方法: 1) 平均が50,標準偏差が15の正規分布に従っているデータから25個のサンプルを抽出する。 2) これを100回繰り返して、その度に平均値を求める。 3) 100個の平均値の標準偏差を計算(=標本誤差と呼ばれるもの) rという無料の統計ソフトで、1)と2)をいっぺんにやってしまう関数を定義します(別に意味は分からなくても良いですが)。 my.func <- function(){ dat.mean <- c() for(i in 1:100){ dat.mean[i] <- mean(rnorm(25, mean=50, sd=15)) } return(dat.mean) } #3)を実行 > sd(dat) [1] 2.988013 4の方法: 1) 単純に平均50,標準偏差15の正規分布から25個サンプリングする。 2) その25個のデータから標準偏差を求める。 3) 公式で標準誤差を求める。 #1回目のチャレンジ > x1 <- rnorm(25, mean=50, sd=15) > sd(x1)/sqrt(25) [1] 2.955733 #2回目のチャレンジ > x2 <- rnorm(25, mean=50, sd=15) > sd(x2)/sqrt(25) [1] 2.145187 #3回目のチャレンジ > x3 <- rnorm(25, mean=50, sd=15) > sd(x3)/sqrt(25) [1] 3.807814 ■3の方法ではs.e.=2.988013となりましたが、■4の方法では近い値にはなっても完全には一致しませんよね。そういうものです。