ぎょーむ日誌 2002-12-16

2002 年 12 月 16 日 (月)

• 0820 起床． 朝飯． 0850 自宅発． 曇． 0900 研究室着．
• なにかばたばたと計算してたんだけど思いだせん …… また quasi-likelihood 計算だったか?
• ひと区切りついたのでお茶部屋で休憩． M2 井田君と一般化線形モデルのギロン． なんか近ごろあちこちでこの計算方法 すすめてまわってるような気がする．
• かとー先生に MacOS X の挙動など教えてもらう．
• M1 大澤君と prunusdb のススめかたの相談． ともあれ現時点では Buzzaz 本にたたられていて， いそがしそうだ．
• D3 浦口さんのカエデ成長 quasi-likelihood モデルの ここまでの報告など． さーて， このハナシの統計モデリングは これで最後であってほしいところであるが ……
• 1250 北大構内走発． 積もった雪が踏み固められてしかもとけかけている． 気をつけないとすべりそうだ． 1325 研究室帰着． 昼飯． 昼飯後に大澤君のアメフト概論を拝聴する．
• 土曜日のぎょーむ日誌に TWINSPAN がどーのこーのと書いたら， もっとしっかり勉強せよとの御鞭撻と 関連情報を教えていただく …… ふーむ． この名義変数としか思えぬ優占度を 連続変数であるかのよーに扱うお作法， というのは ……
• 2000 年 3 月にメキシコのバハカリフォルニアで遭難した 中野繁さんの論文集「川と森の生態学」 (6000 円，ISBN4-8329-8021-1) がとどく． 中野論文を苫小牧の村上さんたち共著者たちが和訳してまとめたもの． 私もこれのお金だしてたんだなあ (すでに忘却)． 巻末用語集に関してコメントしたのは覚えてる． 出版社は 北海道大学図書刊行会 なんだけど， 同サイトにはまだ情報があがっていない．
• 1530 研究室出て理学部・数学科 (理 3 号館) 5F の図書室へ． かなり充実している． Quasi-likelihood 本 あった …… 書いた本人としてはすごく具体的に書いてるつもりなんだろうけど， 私からするとかなり数学よりな内容． ともあれ借りだしていく．
• ついでに農学部の図書室にも立ち寄る． ここの中央図書室はなかなか立派なんだが， 探してる Forest Science 古雑誌は無い． なんでも森林科学の図書室に収蔵とのことだが …… えーと北館 2F か． うすぐらく迷宮のごとき農学部内をさまよう． 森林科学図書室． はりがみ． 利用可能なのは火曜日・金曜日の 1200-1300 …… 週に 2 時間しか開いてないのか!! さすがは北大・農．
• 帰りがけに図書館本館にもたちよりごそごそと．
• もどってくると 1730 すぎ． 浦口さんから推定宿題メイル (本人は所用ですでに下校)． イタヤカエデはうまくいったけどオオモミジはうまくいかない? ふーむ． 一見するとオオモミジのほうがラクそうだけどな ……
• 推定計算やってみるとたしかにうまく収束しない． 謎． いろいろ試行錯誤してみる．
• その結果わかったことは， パラメーター推定の反復試行の途中において， 平均値がゼロ近くでぢたばたするとまずいことが生じうる， とゆー可能性だ． Quasi-score のマイナス側の極値が消滅して …… Newton-Raphson 法が破綻してるんだろうな．
• これを回避するには …… 数値計算上の都合として quasi-score 関数を某所でへし折ってしまうか． 一見すると乱暴なようだけど 今までやってることとは無矛盾であり， 疑似尤度を尤度そのものとみなして これを最大化するという方針とも合致している．
• さーて， 結果は …… この簡潔な処置で どうやらうまく推定できるようになったみたいだ． データのばらつきがでかい場合においても， いいかげんな初期値からスタートして 不安感なく収束するようになった． 一件落着．
• カエデの推定とは無関係なんだけど， かとー先生から R における ロバスト推定あれこれを教えていただく．
  • MASS ライブラリの の rlm 関数: これは M 統計量を使う方法のよーで
  • lqs ライブラリの の lqs 関数: これはどうも平均値のまわりの， 何と言うか「重みづけカーネル」 みたいなものを使って計算するみたいで …… しかし例をみるとうまく計算しているように見える
• 私はこの種の問題はぶーつすとらっぷ的に解決しよう， ともくろんでいたんだけど …… どちらが統計学的にはマシなんだろうか． 2110 研究室発． 雪どかどか． 2130 帰宅． 体重 73.6kg． 晩飯．
• [今日の運動]
  • 北大構内走 1250-1225． ストレッチング．
• [今日の食卓]
  • 朝 (0830): 米麦 0.6 合． ハクサイ・タマネギ・ピーマン・油揚の炒めもの．
  • 昼 (1335): 弁当． 研究室お茶部屋． 米麦 0.8 合． ハクサイ・タマネギ・ピーマン・油揚の炒めもの．
  • 晩 (2200): 米麦 0.8 合． ネギ・ニンジン・豆腐の味噌汁．