ぎょーむ日誌 2000-12-05

2000 年 12 月 05 日 (火)

• 0745 起床． 寝たのは 0300 ごろかなぁ． 寝ようとすると地震で建物が ぎしぎしと揺れるし．
• 朝飯・弁当の準備． 朝飯． おしまい． うう．
• 0815 自宅発． 晴れ． 0835 東京モノレイル流通センター発． 浜松町の気温 5 ℃． 0852 研究所着．
• まったくどうでもよいことなのだが …… ちかごろ Yahoo! オークションで 英語キーボードのノートPC が出品されなくなったな， と思ってたら …… オークション検索エンジンの 仕組みが変わったようで， 「英語」ではひっかからないけど 「英語 キーボード」というような指定をすると ぞろぞろと見つかるのである． 一般的すぎるキーワードはダメ， ということなのか．
• 開空度計算ライブラリーの汎用性を改善して， さらには 実測との対応を調節したいんだけど， 小川シミュレイター解説と作図に関して 待機命令が解除されていない状態なんで， この方面では動きまわりにくい．
• これまたどうでもいいのだが， 上のような状況の原因である 某あやしげ編集部はまたしても音信不通である． たぶん何か破綻しているんだろう． ま， 私は困らんハズだけど ……
• ということで， いつまでたっても終わらないが もはや誰も何も期待していない 最尤法初歩の初歩解説の 数値解析篇をてきとーにごまかすべく， 数値計算本などを調べ直してみる．
• 多次元の Newton-Raphson 法ややこしそう …… などと言ってたけど， 実は単なる Taylor 展開， とわかった． ぼけぼけだな．
• 数値計算に関しては， ひろくいろいろな方法を検討するなら ``Numerical Recipes'' (Press, Teukolsky, Vetterling and Flannery) がどの手法でも楽しく紹介されており (驚くべきことに， なんと前記 URL では本文が 完全に公開されてしまってる !!)， ちょっとつっこみたいというか 基本から再検討しませうというときは 「数値計算の常識」 (伊理・藤野) などが読むと「そうか !」と感銘うけたりする．
• …… などとやってるうちに昼飯時間． とりあえず弁当食う． 食いつつ， どーも具合の思わしくない openssh 周りの 情報をネット上で漁ってみるも， やはりよくわからない …… いや， わかった． /etc/ssh/sshd_config で，

	IgnoreRhosts no
	RhostsRSAAuthentication yes

…… こうだ． /etc/rc.d/init.d/sshd restart で daemon を再起動． これで $HOME/.shosts が反映され「お友だち関係の輪」 を広げていくことができる．
• 火星に堆積構造， すなわち大量の水が存在してた痕跡を発見． うーん， 楽しそう． 「火星の暴走の夏」って アイデアありましたなぁ ……
• ファクシミリで某怪しい雑誌カラーペイジの 校正刷が遅れられてくる． うーん， やっぱり あやしげ．
• さらに怪しげな 「技術成果の外部発表等管理票」 なる 無意味としか思えぬ 紙キレになにごとかを書き込んで NASDA に出さねばならない． ついでに原稿も． 小川シミュレイター解説記事の中に 「なぜ NASDA ロケットは失敗続きなのか」 といった悪口を書かれていないか 検閲せずにはいられないんだろーね．
• この web site などは 「技術成果」だけでなく 浜松町独房群のありさまなんかを ずいぶんと 「外部」に「発表」しちゃってるんだけどな ……
• さてさて， 数値最尤法解説の続き …… と Numerical Recipe の あちこちを眺めつつ， 説明方法を考えて， ちょっと LaTeX ファイルを書きかけて， 打鍵の指が固着した …… この mlfitting のやりかたって …… 全然まちがえてるのではなかろーか ?
• そんなバカなと思いつつ， mlfitting のソースコードを当たる． うーん， ちょっと待て． 私は Newton-Rapthon 法で ナニがゼロになるパラメーターヴェクターを 計算してたのか …… ?
• アタマを冷やすために独房から離れ …… このフロアーには 逃げ出す先がないんだよね． しょうがないんで， 体を動かすために 一階まで階段降りて自販機で缶コーヒー買う． その場で飲む． ふう．
• 冷静に考えると数値微分してたのは， 対数尤度関数． 対数尤度関数の微分とは， すなわちスコアー関数． これはゼロにすべき値だ． したがって問題ない， と頭の中が整理された． 間違いではなかったワケだ．
• 先ほどのように ついつい錯覚してしまった原因は …… うん， ふたつほど挙げられるだろう．
• ひとつは mlfitting のソースコードの わかりにくさに由来するものである． f だの df などという名前のヴェクターや行列が 何をあらわしているか わからなくなったんで， 焦ってしまったのである．
• このコードはずっと前に書き始めたもので， 昔のコードの影響を引きずっているのである． なにしろもともとは最尤法でもなければ Newton 法でもなかった． 九大理学部 3 年のときの 生体高分子の授業で配られた Polytope 法 (あるいは Simplex 法) による非線形最小二乗法だったんだよね．
• それはポケットコンピューター (プログラムできる電卓， といいますか) 用の BASIC で書かれていて． それを Fortran77 で書き直して …… そっか， ここで一度「アルゴリズム事典」(奥村) に掲載されてた C 版に全面的に置き換えてしまったな． いま思い出した． で， それを C++ にして． また C++ で書かれたコードは何回か改定されるうちに 分岐して最尤法になって， 最後に (今年の 2 月なんだが) 数値微分・逆行列・Newton-Raphson 法連携型の 数値最尤法コードになったんだよね．
• コードに歴史あり． 静的な配列 (…… のハズなんだが，なぜか実行時にサイズが決まる) が使われているかと思えば， STL だの動的連結関数 (この改造は今年 9 月) だの 使ったりしてるところもある． そろそろ設計を一からやりなおしたほうがよさそう．
• それから 「mlfitting って全然間違ってんじゃないの」と 錯覚したもーひとつは， 計算方法そのものがわかりにくい， というトコロにありそうだ．
• 最尤法の「定石」が スコアー関数→ゼロかつ 対数尤度関数→最大， だからといって， いつもいつも 計算機の中でばか正直にそれを Newton-Raphson 法などで再現してやる 必要などない．
• たとえばいきなり修正 Powell 法か何かで 対数尤度関数→最大 とやったっていいはずだ． というか， たいていの場合， こちらのほうが計算量も危険性も少ないのでは あるまいか．
• えーい， いろいろと改良せんといかんなぁ． どうして最初からスジ道が見えんのだろう．
• うん? なにやら ESTO から封筒． 全員に配られているのか． CD-ROM …… うーん， 私にこんなデーター配布されてもなぁ．

[観測衛星からみた温度分布]
……って何に使えっての．

• 1755 研究所発． モノレイル定期券買う． 10300 円． 浜松町駅の本屋 dan． プログラミング言語 Ruby の本の種類増えてる． 1830 東京モノレイル浜松町発． 1910 帰宅．
• 乳脂肪分のおおい安直な晩飯食ったせいか， 食後に強烈に眠くなる． 少し濃いコーヒー飲む． 眠いけど覚醒しつづけてる状態．
• …… というアタマの状態では， 仕事もはかどらんので， ちょっと大森町商店街のほうにでかける． セヴンイレヴンで電話代の支払い． 6910 円． また額を抑えることに成功した …… しかし， 来月の請求はまた高くなるだろう． 今月はむやみに意味もなく 長時間接続やっちまったからなぁ．
• しかし今週の金曜日には大森にも文明の曙光 …… 東京めたりっく通信 の ADSL 回線網が開局するはず． 申し込んでみよーかな ……
• 買物などして帰宅． 明日以降に食べるための ダイコン・ニンジンの味噌汁を作る． 冬のほうが野菜はむしろ豊かだ． 北海道だと「輸入モノ」ばかりになるので， 冬の野菜は値段が高かったような気がする．
• あ， もう 0120． いかん． コーヒーが効きすぎて すぐには眠れそうにない．
• 今日の食卓
  • 朝 (0755): 米 0.4 合． ダイコン・ニンジン味噌汁の残り． コマツナのおひたしゴマあえ．
  • 昼 (1220): 弁当． 米 0.8 合． コマツナのおひたしゴマあえ．
  • 晩 (1930): 米 0.6 合． 洋風雑炊． タマネギ・ニンニクの薄切りを炒める． 牛乳いれて加熱． 冷や飯・顆粒状クリームシチューの素を投入．