研究室のページ/食資源経済統計学/中谷ゼミのページ/2011年度のバックアップ(No.11) - PukiWiki

[ トップ ] [ 新規 | 一覧 | 単語検索 | 最終更新 | ヘルプ ]

バックアップ一覧
ソースを表示
研究室のページ/食資源経済統計学/中谷ゼミのページ/2011年度は削除されています。

中谷ゼミ　トップページへ戻る?

中谷ゼミ　2011年度のページ †

2011年度　参加者
- 修士2年　　2名
- 学部4年目　6名
- 学部3年目　3名

時間　毎週火曜日　午後2時45分から
場所　農学部　共同実験棟　共308

2011年度の教科書 †

前期
- 前園宜彦(2010)詳解演習　確率統計，サイエンス社，ISBN 978-4-7819-1265-3，定価1,890円
  をやる可能性が濃厚です．
  - 上記書籍には正誤表があります．ここからダウンロード（PDF形式）してください．

本年度前期は，問題演習を中心にします．

「事前に問題を解いておき，ゼミ時間中に（私を含め）無作為抽出で選ばれた人がその解説をする」という方式でいきます．

無作為抽出について †

解答者決定の無作為抽出が本当に無作為に行われているのか疑問があるようですので，以下のRプログラムでシミュレーション実験をしてみました．
```
nsim <- 10000000
number <- numeric(nsim)
a <- 1:12
for(i in 1:nsim){
   number[i] <- sample(a, size = 1, replace = T)
}
hist(number, breaks = seq(0, 12, by = 1), freq = F)
```
- sample(...)は，抽出数が一つなので，"replace = T"があってもなくても結果は変わりません．
- 1千万回のシミュレーション実験では，1から12までの番号の出現割合は，上の図のように8%とほぼ均一になりました．ちなみに，1/12 = 0.0833です．

というか，上のプログラムは異常に非効率的でした．forループを回さずに，
```
nsim <- 10000000
a <- 1:12
number <- sample(a, size = nsim, replace = T)
```
とすれば，実行時間が劇的に改善されます．

ちなみに，シミュレーションの回数を50回 (nsim <- 50) にすると
となって，かなりのバラツキがあります．
少ない試行回数（抽選回数）で，出る目の回数を均等化することは，無作為抽出ではありません．つまり，「それぞれの目が出る回数を均等化する」という作為を働かせることになるので，無作為ではなくなってしまいます．