ぎょーむ日誌 2007-09-10

2007 年 09 月 10 日 (月)

• 0830 起床． ねぼうしてしまった． 朝飯． コーヒー． 0945 自宅発． 曇． ちょい雨． 1000 研究室着．
• 昨日のゐんばぐす試行錯誤の報告 (つまりぎょーむ日誌) 書いたり， 放置メイルに返信してみたり． そして， 回避不可能感のある査読をひきうけてしまったり …… うーむ， イヤなかんぢだ． たまには楽しい査読とかやりたいものですなぁ ……
• ECOAS プランクトンデータの階層ベイズモデルの改良あれこれ．
  • 観測値はそんなに信頼できない， ということで分散パラメーターの事前分布を修正
  • 標高と湖沼面積の交絡のせいで混交が悪くなってる …… しょうがないので湖沼面積も統計モデル化， つまりベイズモデルの中で標高補正
  湖沼面積の標高補正， こりゃーいいアイデアだとやってみると …… なぜかしらかえって劣悪になった． えーい， 湖沼面積だけでいいや． どうせ標高が高くなると湖沼面積は小さくなるわけだし．
• で， 標高ぬいたら収束するまでの時間が激減， 気分もすっきり．
• 12000 MCMC step (最初の 2000 だけ burn-in) 240 秒ほどで計算終了． たぶんもっともっと短い計算時間でも十分だろう (事後分布表)．



• 昼飯．
• ということで， ここまでの ECOAS プランクトンデータの統計モデリングについて， 占部さんたちに報告かき． 現時点では， ごくごく簡単な報告ですませるつもりなんだけど， 説明にめんどう感つのる …… まあ， どういう説明でもそうなんだけど．
• 二時間ほど費やして ECOAS PukiWiki 上に統計モデル解説メモ書き． けっこう長くなってしまった． うう． 東北方面にメイルかき． ひと休み． 時刻は 1635．



• revised な 母子里 Gibbs 林冠モデルはどうなったのかしらん， と思ったらまだ with Editor どまり， か． Agricultural and Forest Meteorology はこのあたりのんびりしてるのかしらん?
• 論文といえば「アリ」がまだ残存してたな． 早く脱アリしたいんだけど ……
• じつは ECOAS プランクトンデータ， まだ動物プランクトンの組成解析をやってないんだよね ……
• とゆーことで， よろよろと動物プランクトンのバイオマス組成データなどながめてみる．
• データ変換して， また「全データが見える作図」 プログラミング． これは 64 種のプランクトンのうち一種．



• プランクトンデータながめてるうちに気づいたんだけど …… サカナって pH が 6 より小さい (酸性の) 湖沼にはすんでないんだな (動物プランクトンなんかはそういうところにもいる) …… 理由は知らないけど． R の glm() で調べてみたら pH は log(湖沼面積) より良い説明変数になってた (そして log(湖沼面積) は標高より良い)．

[高山? 低山? の酸性湖?]

えーと， このあちこち交絡しまくりデータなので注記しておくと， なぜかしら湖沼の pH は標高が高いほど低くなっている (酸性になってる)． ということで， このデータの場合， 標高 1600 m 以上にも魚はいないわけだが． とはいえ， pH = 3 というレモン汁みたいな湖は標高 1000 m 以下にあるな．

• プランクトン 1 種につき 1 ペイジの図セットを 64 ペイジ PDF 出力する． R の pdf() 関数つかえば簡単に …… しかし， どうも種ごとに見てもよくわからんなあ． プランクトンの知識がないからでもあるけれど．
• 空腹になってきたので， 1920 研究室発． 1935 帰宅． 体重 68.4 kg． 晩飯．
• 冷凍庫なき冷蔵庫の製氷棚の氷はがし作業． 3-4 週間に一度ぐらいやってる．
• 今後のぎょーむ予定． 冷静に考えて
  • 9 月: ECOAS 下うけデータ解析 (プランクトンと八甲田湿原群) でオワりそう， 月末に仙台出張??
  • 10 月: 日光会議の準備と 11 月統計学授業の準備
  • 11 月: 統計学授業やりつつ，脱アリ作文?
  ホントはもっと早く脱アリしたいんだけどね …… そしてフィールド調査の季節がオワると， いろいろとめんどうな問題がもちこまれるんだろうね．
• [今日の運動]
  • うんどう休養日
• [今日の食卓]
  • 朝 (0920): 米麦 0.7 合． キャベツ・ニンジン・キュウリ・タマネギのサラダ． ネギ・ブナシメジ・卵の炒めもの．
  • 昼 (1330): 研究室お茶部屋． 米麦 0.7 合． キャベツ・ニンジン・キュウリ・タマネギのサラダ．
  • 晩 (2030): 米麦 0.7 合． ピーマン・タマネギ・卵炒飯． キャベツ・ニンジン・キュウリ・タマネギのサラダ．