生態学のデータ解析 - 統計学授業 2017

• 統計学の授業やります (2017 年度前期後半, 2017 年 6 月)
  • 「植物生態学特論 I」の一部, 2017 年の 6/5 (月) から開始
    • 誰でも参加できます
    • 事前のもうしこみなどは何も必要ありません (聴講のみで，単位を必要としない場合)
    • 全部ではなく部分的に参加してもらってもかまいません
  • 授業メーリングリスト:
    • 授業ごとの課題 (単位をとる人のため) と「講義資料ダウンロードできます」案内などを通知します
    • 単位をとる人は必ず ML に加入，そうでない人は加入してもしなくてもよい
  • 教科書「統計モデリング入門」 (あるいはやや古いですが: 2008 年の講義のーと)
  • 月曜日・水曜日の 3 講目 (13:00-14:30)
  • 教室: 地球環境科学研究院 A 棟 A809 (エレヴェイターでてすぐ前の部屋)
  • 短縮 URL: http://goo.gl/76c4i
  • 注意: ここの例題に使われている「データ」はすべて架空のものです

---

---

第 1 回: 6/05 (月) 観測されたパターンを説明する統計モデル  

• 投影資料 kubostat2017a.pdf (slide), HOkubostat2017a.pdf (handout)
• 教科書「統計モデリング入門」
  • 参考: 2008 年の講義のーと
• R のインストール R download site

(*** 6/07 (水) は授業なし -- NO class!! ***) 

第 2 回: 6/12 (月) 確率分布と最尤推定 

• 投影資料 kubostat2017b.pdf (slide), HOkubostat2017b.pdf (handout)
• Example data: data.RData
  • how to read the data using R? -- load("data.RData")

第 3 回: 6/14 (水) 一般化線形モデル: ポアソン回帰 

• 投影資料 kubostat2017c.pdf (slide), HOkubostat2017c.pdf (handout)
• 例題データ: data3a.csv
  • うまく読めない場合は R 上で以下の命令を実行してください:

  d <- read.csv("http://hosho.ees.hokudai.ac.jp/~kubo/stat/2017/Fig/poisson/data3a.csv")

第 4 回: 6/19 (月) モデル選択と検定 

• 投影資料 kubostat2017d.pdf (slide), HOkubostat2017d.pdf (handout)
  • 参考: FAQ モデル選択
• 例題データ: data3a.csv
• parametric bootstrap 尤度比検定: pb.R (関数 get.dd() と pb() の定義)
• うまく読めない場合は R 上で以下の命令を実行してください:

•   d <- read.csv("http://hosho.ees.hokudai.ac.jp/~kubo/stat/2017/Fig/poisson/data3a.csv")
• source("http://hosho.ees.hokudai.ac.jp/~kubo/stat/2017/Fig/lrtest/pb.R")

第 5 回: 6/21 (水) 一般化線形モデル: ロジスティック回帰 

• 投影資料 kubostat2017e.pdf (slide), HOkubostat2017e.pdf (handout)
• 例題データ: (logistic 回帰) data4a.csv, (offset 項) d2.RData
• CSV file をうまく読めない場合は R 上で以下の命令を実行してください:
  • read.csv("http://hosho.ees.hokudai.ac.jp/~kubo/stat/2017/Fig/binomial/data4a.csv")
• 参考: 生態学会大会2012 内の kubo2012table.pdf (分割表の作図・GLM・ベイズモデル)

第 6 回: 6/26 (月) 階層ベイズモデル 

• 投影資料 kubostat2017f.pdf (slide), HOkubostat2017f.pdf (handout)
  • 本/データ解析のための統計モデリング入門 の第 8-10 章の例題なども参考にしてください
• MCMC の例題ファイルなど:
  • data.RData (3. と同じ)
    • あるいは data <- c(4,3,4,5,5,2,3,1,4,0,1,5,5,6,5,4,4,5,3,4)
  • R で install.pacakges(c("rjags", "R2WinBUGS")) が必要
  • runjags.R
  • JAGS の入手先: http://sourceforge.net/projects/mcmc-jags/files/
• 階層ベイズの例題ファイルなど:
  • R で install.pacakges(c("rjags", "R2WinBUGS")) が必要
  • data7a.csv
  • runjags.R (上のファイルと同じファイル名だけど別もの)
  • JAGS の mcmc.list を bugs データに変換する関数: mcmc.list2bugs.R

第 7 回: 6/28 (水) 階層ベイズモデルと時間変化モデル 

• 投影資料 kubostat2017g.pdf (slide), HOkubostat2017g.pdf (handout)
• 例題 1 (植木鉢上の植物個体の種子数) のデータ: d1.csv
• 例題 2 (給食と身長ののび) 関連ファイル: 岩波DS01

第 8 回: 7/03 (月) 生態学の時系列データ解析入門 

• 投影資料 kubostat2017h.pdf (slide), HOkubostat2017h.pdf (handout)

---

• 昨年度 2016 統計モデリング入門授業，改善してほしいところの要望・質問
  • R 実習を多くしてほしかった
  • 統計モデルの R コードの説明をもっと詳しく
  • 階層ベイズモデルの説明がわかりにくい（説明の数値例とか？)
  • 個体差 r が分布の期待値を変えるところがわかりにくい…階層ベイズモデルがよくわからない
  • 統計学用語について…ページや統計モデリングの本の該当ページを課題メールに書いてほしい
  • JAGS の使いかたをきちんと説明してほしい