ぎょーむ日誌 2012-12-19

2012 年 12 月 19 日 (水)

• 清水出張第 3 日目， 0540 起床． またどろなわで今日の準備あれこれ ……
• 宿の朝飯を食って 0830 チェックアウト． 会場が近くて助かる …… 0840 清水テルサ着．
• 統計セミナー最終日， 今日はとりあえず分割表 と統計モデルを対応づけるハナシを紹介して， そのあとは皆さんから出していただいた漁業な例題をあつかったのだが …… まあ， 今朝つくったペイジ に示しているように， なかなか難しかった ……
  • DeLury 法: 捕獲された個体数にもとづいて集団サイズを推定する という問題なのだが …… 集団サイズを潜在変数としてベイズモデル化してみたんだけど， 「捕獲前の個体数」に関しては主観的事前分布とかもちださないと， うまくパラメーター推定できないみたいだ (工夫が足りないだけかもしれないが)
  • CPUE 標準化: CPUE は漁獲高にもとづく「個体群密度」みたいな 量なのだが (標準化とか漁具の違いによる効果などをのぞくこと) …… これまでこの推定は， 「データどうしの割算」 「各年ごとの効果をいちいち最尤推定する」 といった， どうみても「まずいやりかた」の寄せ集めなので， ここは架空データを使ってそのまずさを指摘しようとしたのだが …… やってみると， その「まずいやりかた」でも， そこそこ悪くない結果が得られると判明してしまった． これは将来予測とか興味のない状況， 「過去のある時点での CPUE はどれぐらいなのか」 といった問題にとりくむ場合， かつデータがかなりたくさんある場合ならば， 「まずいやりかた」でもそんなにまずくない …… ということのようで．
  • 体長でランクわけされたマグロいち場のハコ数データから もとのサイズ分布を憶測する問題: これについては「ひきのばし法」という， まあカーネルみたいなものを使ったスムージングで対処していたとのことで， このあたり何かもうちょっとマシな方法は …… と調べてみたのだが， けっきょくのところ BUGS 言語で「ひきのばし法」をやったら， 架空データの「真のサイズ分布」がそれっぽく再構成されてしまったので， まあ特に新しい工夫はいらないかも， でおわってしまった．

• ということであしかけ 3 日間にわたった， 清水の統計モデリングセミナーは午前中で終了．
• 会場ちかくで昼飯． 1309 清水発． 1337 静岡， ひかりに乗り換えて 1432 品川， 京浜東北線に乗り換えて 1440 浜松町． 1448 東京モノレイルの特急にのって同発 …… たしかにこれは速いな． 1502 羽田空港．
• 1520 JAL521 20D 登場． ひさびさの B777-200． この時間なので搭乗率は 6 割ぐらい? 1655 新千歳空港． 空港 2F の出発フロアで AIRDO の道民会員の延長てつづき． 1719 JR の快速エアポートにのって同発． 温暖な清水からもどると …… 札幌は寒い! 1840 帰宅． 晩飯．
• [今日の運動]
  • 帰宅してからとらちゃんえくささいず， 清水にいるあいだは「背筋の疲労感」 がない状態だった， と気づいた．
• [今日の食卓]
  • 朝 (0820): 宿の朝飯．
  • 昼 (1200): 清水魚いち場ちかくの某店で「静岡丼」 1300 円， これはマグロ・サクラエビ・シラスなどのせたばら寿司． 清水は何やら「アナゴ」が名産らしいので， アナゴ天丼を食べたかったのだが …… 「時間がかかる」と言われてあきらめた． 来年 3 月に静岡市で生態学会大会があるので， そのときにチャンスがあるか? この店はけっこうひどくて， 「時間のかからない品がほとんどない」 「会計は全員一括でないとうけつけない」 「両替もしない」 というないないづくし．
  • 晩 (1900): 米 0.7 合． 肉じゃがなど．