ぎょーむ日誌 2004-11-30

2004 年 11 月 30 日 (火)

• 0620 起床． 朝飯． コーヒー． 0750 自宅発． 晴． 0800 研究室着．

[いきなり積雪]

昨日の朝はまったく雪なんか積もっていなかったんだけど．

• 全国的に修論・D 論しーづんということで， 週に何通か統計学相談メイルいただく． 今朝も一通返信．
• なーんか北大でも学術リポジトリを作ろうとしてるみたいね …… とりあえず， アンケート に答えてみる． 図書館 関係は， 北大の学術的な基盤の中ではかなりマトモに機能してる部門だと思う．
• また北大ネット不調． 学外につながったりつながらなかったり via http． 北大の学術的な基盤の中ではかなりマトモに機能してないへっぽこ部門．
• それだけじゃなく， ちかごろは接続先とはほぼ無関係にこういう ``Access denied'' しやがるし． これは北大の proxy server (強制通過型) の設定・挙動がおかしいのでこういう動作になる．

• うーむ …… さっさと次の仕事にかかるべきなんだが， 「よーやく授業おわった」 脱力からのりはびり状態で， 何もてきぱきとススまん．
• 苫小牧直径成長モデル …… 次の一手をどうするべきか判断が難しいので， お尋ねメイル送ってから， べつに急ぎでもない推定プログラム改善作業なんぞやってみたり． 名目的には， 推定プログラムの点検とでもなるのか …… とサボタージュ労働者は考えるのであった．
• そういや， アカマツ作図宿題もあったなぁ ……
• 苫小牧モデル推定プログラムの手直しだけで昼飯前の仕事終了． これは次の高速化の布石． 雪がふってるんで， しばらくは北大構内走は様子見． 昼飯． 昼飯くってるといきなり煙まみれな大学院生があらわれたり (注: 本人が燃焼していたわけではない)．
• うーむ． 苫小牧モデル作文の私が担当すべき箇所に関しては， 時間的よゆーが生じているような …… DEFCON 一段さげてよいのだろうか．
• とりあえず， 午前の続きで， 推定プログラムのアヤしげ改善のつづき． ともあれ， アヤしげな部分ごとに関数に隔離されてるんで， 困ったことにはならないんだけど．
• しかし， 気分的にはかなりラクなはずの修正作業もススまん． ということで， 先日つくった書籍リストを送って本の発注．
• 苫小牧モデル推定高速化問題， 予想があたった． この複雑なモデルのパラメーター推定には glm() のような単純なやりかたは使えない． しかしながら， 直径成長年変動の説明に使う代謝気象値・光合成気象値を 数千・数万の候補の中からしぼりこむには (特定するわけではない)， glm() みたいな計算で十分， とわかった． これでうまくいけば 7-8 時間かかっていた計算を 10-20 分ぐらいに 短縮できる．
• で， 一回試運転やってみると， ほとんどうまくいった． ただしミズナラでは AIC が 0.2 だけ大きくなってしまった． ということで， 絞りこむ候補数を 10 → 20 にしてみるか．
• 20 にしてもミズナラだけだめだ …… ちくしょう， 食植昆虫に食われまくりの変態樹木め．
• アルゴリズムを変更する． ノミネイト数固定方式はヤメ． 良さそうなやつから順番に試していく． そして， ある気象値セットが N 回にわたって「勝ちぬく」 ことができたら， その気象値を採用するという「勝ちぬき戦」型を採用する． とりあえず N = 10 ぐらいでやってみるか …… ミズナラ以外はこんなことやる必要もないんだけどね．
• ミズナラ， N = 10 でもダメか． ならば N = 20 なら?
• ダメだ． N = 30 もダメ …… そうだな， 粗い選抜であっても glm(..., family = poisson) はまずかったか? glm(..., family = quasipoisson) にしてみる． 平均 vs 分散の関係はこんなふうだし …… あれ? ミズナラでは分散がだいたい平均に比例してるぞ?

• quasipoisson 推定で使うのは擬似尤度だから aic は計算できん． そこで deviance (これも擬似っぽいな) 見て選抜かけてみる …… ダメだ． ミズナラだけはどうしてもうまくいかん． 他の 5 樹種は完ペキなのに．
• ということで， 人間に (少なくともこの私に) 悪意をもっているとしか思えんミズナラ一樹種にわずらわされるのはヤメ． AIC だって 574.5 → 574.7 と 0.2 悪化するだけだし． これ以外は以前とまったく同じ． 計算時間は 7.5 時間 → 5 分， と短縮． 90 倍の高速化 (60-70 時間を必要としてたころと比べると 800 倍以上の高速化)． すばらしいことですね．
• といった結果どもを wiki にアップロード．
• 1830 研究室発． 2000 帰宅． 体重 73.4kg． 晩飯．
• 苫小牧直径成長モデル， PDF 出力式の作図プログラム作りなおし． うう． われながら健忘症とゆーか， ちょうど 一ヶ月前 にやらかしたのと同じ間違いを連発． いやはやと思いつつバグとり続く．
• [今日の運動]
  • 腹筋運動 30 × 3 回． 腕立ふせ 5 × 3 回．
• [今日の食卓]
  • 朝 (0720): 米麦 0.7 合． ネギ卵炒飯． キャベツ・ダイコン・サバぶしの味噌汁．
  • 昼 (1240): 研究室お茶部屋． 米麦 0.7 合． キャベツ・ダイコン・サバぶしの味噌汁．
  • 晩 (2050): 米麦 1.0 合． チンゲンサイ・タマネギ・豆腐の炒めもの． キャベツ・ダイコン・サバぶしの味噌汁．