ぎょーむ日誌 2014-12-12

2014 年 12 月 12 日 (金)

• 東京出張 3 日目． 0600 起床． 宿の朝飯．
• ちょっとばてぎみかも． のんびりしてみる．
• 生態学会さーば雑用を一件かたづけて， 0930 宿発． 0955 統数研着．
• 今日も統数研共同研究集会 「野外生態データの観測過程と統計モデリング」． 午前は「動物個体の発見率」の推定が重要になるハナシ …… 端的にいって， 推定はむずかしい．
• 昼飯．
• library(unmarked) について調べたほうが良さそうだな …… 動物生態学でやたらと使われてるが， package がおかしいか， ユーザーが使いかたをまちがっているか， どちらかだという気がする．
• N-mixture モデルというやつなんだけど， ひじょーにかんたんで，
  • 発見された個体数 ~ dbin(総個体数, 発見確率)
  • 総個体数 ~ dpois(平均人口密度×面積)
    (あれ? ちゃんと密度×面積で計算してたヒトいたかな?)
  • (だれもきちんと説明してなかったけど，たぶん)
    log(平均人口密度) ~ 無情報事前分布
    logit(発見確率) ~ 無情報事前分布
  というもの 発見個体数が観測データで， 総個体数はわかっていない．
• 別にこんなんでいいぢゃん …… とゆー気もするけど …… ちょっと気になるのは， 観測データである発見個体数 ≤ 総個体数， という制約があること．
• 直接観察できない， 総個体数という離散値の変数の (事後分布の組み合わせにもとづく) 予測分布みたいなものを知りたいわけだが， それは「切断ポアソン (truncated Poisson) 分布」を積分した確率分布になるのでは? これはけっこうへんてこなカタチの確率分布になるはず．
• 16 時すぎ， 立川を出発して氷点下の北海道に帰った．
• [今日の運動]
  • ラジオ体操: 0 回
• [今日の食卓]
  • 朝 (0630): 宿の朝飯．
  • 昼 (1200): 納豆まき．
  • 晩 (1900): 羽田空港の蕎麦．