ぎょーむ日誌 2008-04-10

2008 年 04 月 10 日 (木)

• 0830 起床． うーむ． 朝飯． コーヒー． 0930 自宅発． 曇． 0945 研究室着．
• 以前に送信したデータ解析こんさるメイルでロジスティック回帰 glm(cbind(応答した数, しなかった数) ~ ..., family = binomial) と書くべきところをまちがってたことに気づいたので， あわてて訂正メイルかき．
• そして， アリ論文のつづき …… にはとりかかれず， 天塩川わーるどにもどされる． 河川水質学に関して私を指導する大学院生 Nina さんから reading assignment (論文よみ宿題) を出されたからだ． モデル屋稼業は客商売であり日々勉強でもある．
• …… いかん， 関係ないメイルかきに時間をつかってしまった．
• 論文よみ宿題，

  King RS, Baker ME, Whigham DF, Weller DE, Jordan TE, et al. (2005) SPATIAL CONSIDERATIONS FOR LINKING WATERSHED LAND COVER TO ECOLOGICAL INDICATORS IN STREAMS. Ecological Applications: Vol. 15, No. 1 pp. 137-153 (url)

  なる課題なのだが， まあつまりは
  • 土地利用には空間相関がある (スケイルにもよるけど， これはとうぜんそうでしょう， 相関なければランダムぱたーんだよね)
  • 空間相関あるから各要素の影響を独立に推定するのは難しい
  • 土地利用の割算値だけ見ててはダメで， 地上・地下水が「地図上のどこをとおるか」 (つまり土地利用の空間配置) が重要
  といったことで．
• 天塩川モデリングに関しては私のてもとには集水域の地図はないし， あっても上述のごとき精密モデルからパラメーター推定できるほど観測データがない …… 上のような推定をやるためには川にそってぽつぽつとサンプリングしていてはだめで， もっと二次元的にあちこちの水をくんでくる必要がある．
• まあ， ともかく天塩川モデルでも GIS がはきだす割算値を素朴に使ってるわけではない， ということで． いちおう上流の water & ion flow が下流のそれらに影響あたえる， という点では空間構造明示的なモデルではあるな． ただし， 上流の観測点から下流の観測点に水が流れつく時間おくれを考慮してないけど． 流速データとか無いもんで．
• ともあれ， GIS 普及はたいへんめでたいことだけど， 同時に脱「GIS な割算値」という 統計モデリング問題はいろいろ検討しないといけませんね， ということで．
• といったメイル書き． 昼飯．
• 天塩川 water & ion flow モデル (BUGS code) の図をつくってみる． またメイル連絡．

• これでしばらく天塩川モデリングの世界から脱出できるだろうか， と自分でもあまり信じてもいないコトバを口にしつつ， アリ論文原稿にとりくむ． 半年も放置してた原稿をみるのはツラい …… と書いてみて， 「ホントに半年か?!」 と気になってみたので調べてみた．
• …… なンと 昨年 7 月 23 日 以来の放置だとわかってしまった! R 起動

  > as.Date("2008-4-10") - as.Date("2007-7-23")
  Time difference of 262 days
  

  あらら， いやはや， 262 日も放置していたのか． がっくし．
• ちょっと追跡してみたい． 11 月以降から現在までは 修論をはじめとする下うけデータ解析に圧殺されていた， ということにしてもだ． 7 月末から 10 月末までいったい何をやっていたのか?
• ぎょーむ日誌をたどってみると， 中断前日である 7 月 22 日 の記述に

  ありていに言って， この孤立無援なるアリ作文を放棄してしまえば， ただちにアリ体表面の炭化水素組成データ解析に着手できるんだけど ……

  でその翌日にアリ論文中断後， 体表面炭化水素組成データ解析のための アリ触角ニューラルネットワークモデリング (これは当時は必要とされていて， それなりの成果を得たのだけど， 現時点からの結果論でいえば， この触角モデリングにテをださずに 敵対性行動実験論文を完成させるべきであった)， 沖縄のアリ研究会出席， 数学セミナー原稿執筆． これで 8 月はオワって 9 月はひたすら ECOAS 下うけモデリング． 10 月は日光会議の準備 (つまり昔の研究の発掘)． その後は統計学講義のーと 書きで 11 月になっちゃいました， という次第だ．
• 今年度の下うけ減圧期 (すなわち他の皆さんが野外をうろうろしている季節) はこのようであってはなるまい． 教訓:
  • 一度中断した論文かきはなかなか再開できない
  • 研究発表会に行ってもお座敷ごとに ぜんぜん別のハナシをしなければならず， 準備がたいへんなのでできるだけ回数をへらす
  • 安易にしたうけモデリングをひきうけない， とくにぷろぢぇくとのたぐいにまきこまれない
  • 書きもの安うけあいもほどほどに
  といったところか ……
• なぜかしらたいへん恐縮しておられる事務のかたが来て， なんでも 2 月の試験ぎょーむの手当て +8000 円ふりこまれるべきところを まちがって -12000 円ふりこみになった (マイナスふりこみが可能?!)， 今月は修復できないので来月には …… とのこと． あまりにも恐縮しておられるので， 当方も恐縮もーどにひきこまれつついえいえお気になさらないでください， とお引きとりいただく．
• あり? ネット上で確認できる 給与明細だと「その他手当 (入試)」 +20000 円になっているぞ? まちがって余分にふりこんでしまったのかな? まあ， 当方で解決できる問題ではないので， 気にしないでおこう．
• …… などといつまでもうだうだやってないで， アリ論文アリ論文．
• と思ったら斎藤君のデータ解析こんさる． 「ある調査地で A 種と B 種が同時に存在する確率」 といった相関をデータから推定する問題は階層ベイズモデルむけなハナシ (多変量正規分布 な事前分布をつかったモデリング) にほかならないわけだが …… あとでふと考えてみると， 「fixed effects の影響が似てる」 vs 「random effects の影響が似てる」 はどう考えればよいのやら?
• まあ， 後日にデータもってくるそうなので， そのときにまた考えることにしよう．
• アリ論文原稿よみ …… なかなか進捗せづ， 1950 研究室発． 2010 帰宅． 晩飯の準備． ちょっとうんどう． 晩飯．
• [今日の運動]
  • 腹筋運動 30 × 3 回． 腕立ふせ 10 × 3 回． スクワット 100 回．
• [今日の食卓]
  • 朝 (0910): 米麦 0.6 合． ダイコンの葉の味噌汁．
  • 昼 (1420): 研究室お茶部屋． 米麦 0.6 合． 海藻スープ．
  • 晩 (2150): 米麦 1.0 合． ニラ卵． ダイコンサラダ．