ぎょーむ日誌 2000-02-12

2000 年 02 月 12 日 (土)

• 0730 起床． 晴れ．
• 久しぶりに風呂に入る． 10 日ぶりに洗濯する． 伸びてしまった不精鬚を剃る． 部屋を少し片付け掃除機をかける． いやはやプログラミングにハマると 何もかも手抜きになってしまいますな …… って私だけか．
• 1035 自宅発． 1059 JR 琴似駅発． 1115 研究室着．
• ここしばらく C++ なんかで 最尤法のコード書いてんだけど …… 琉球大情報工学の河野 真治さん (あの「入門 Perl」 (ASCII) の著者) が fj.comp.lang.c++ というニュースグループで あまりにも面白いこと 書いておられたんで引用 (Message-ID: <18002.950282804@rananim.ie.u-ryukyu.ac.jp>)． C++ を評して曰く，
  • C よりも危険で、
  • ADA よりも複雑で、
  • ML よりも複雑な型理論を持ち、
  • Prolog よりコンパイルが遅く、
  • Java より Version 依存が激しく、
  • Smalltalk より OOP に対する理解を要求し、
  • Fortran よりも、出力コードが予想しづらく、
  • COBOL よりも出力コードがでかい
  …… うーん， めった撃ちというカンジですね．
• プログラミングにはまっていた間に とどこおっていたメイル連絡などを処理していく． TeX を作った Don Knuth は あるプログラム (えーと， たしか「文芸的プログラミング」を 実現するための WEB システム， つまり TeX と Pascal の複合体で インターネットとは関係ない) 開発していた 一年間郵便物をまるで読まなかったので， 終ったときには山のような未処理連絡を まとめて片付けなければならなかったそうだ．
• Gamete シミュレイターを共同開発してる 富樫さん (室蘭の北大理学部附属海藻研究施設) から いろいろと提案のメイルをいただいたんだが …… メス配偶子 410 万個， オス配偶子 910 万個ぐらいの系でどうでしょう， という提案読んでちょっと気絶しそうになる． <蒼龍> でも それだけいっぺんに動かすのは ちょっと苦しい． まぁ， 実際の実験系でも「希釈」という操作が行われるので， 1/1000 とかに「希釈」した系から始めたほうが よいでしょう． それぐらいなら何とかなりそう．
• お茶部屋で 昼飯を食いながら宮本さんとシミュレイターの話などする．
• かんけーないけど， ちかごろ Yahoo などをダウンさせてる 分散型協調攻撃って なかなか凄いと思う． <鳳翔> は大丈夫だろうな． 悪事に加胆してないはずだが …… 攻撃の手口などをネット上で調べてみる．
• 牧さんの当年生シュート動態モデルが思わしくない， ということで サイズ依存な死亡率の関数型をかえてみる． p0 * exp( - p1 * x * x )． 10 反復ほどで収束． これでどーだ．
• けっきょく p0 * exp( - p1 * x ) というシンプルなやつでいいじゃないかということに なりそう． 牧さんも私も C++ と Perl が頭の中でごちゃごちゃに 混乱して， そのたびに意味不明なコードを書いたりする． ということで， 夕方はそういうモデリングにはまる．
• いろいろと議論してるうちに， 奇妙なやりかたを考えついた． ある確率分布をパラメーター推定に使ったとする． このとき， いくつかの条件さえ満たしていれば 「実は推定に用いた時点で すでに両端を切断 (truncate)してました」 というようなことを言っても， すぐには (あるいは永久に) 反論できない場合がある． ま， これもよく考えるとごまかしだと 判断できる場合もある．
• 1945 研究室発． けっこうな降雪． 2022 JR 札幌駅発． まだ早い時間なので， 久しぶりに古本屋をはしごしてみる． 2230 帰宅．
• 雪の琴似をふらふらと歩いているうちに， 最尤推定法がニガテとする 「端っこのある確率分布」 を使ったときでも， Newton 法みたいなやり方で うまく計算する方法を思いついた． データーのばらつきが誤差だけでなく， 確率論的な過程によって生じたという場合に 有効だ． まず， ばらつきをそういう 確率論的モデルとノイズに分離してしまう． 前者は 「端っこのある確率分布」 を使ったとしても， 後者は十分に裾野の長い確率分布である， として 最尤法を適用すればよいだろう． ようするに データー点に長い裾をもったカーネルを かけてやるのである． データーの種類によっては この手法は使えないんだろうけど， インチキなごまかしと断ずるほど ひどいやり方でもないとは思う．
• 今日の食卓
  • 朝 (0830): ダイコンの味噌汁． ダイコン以外なにもない． ジャガイモあるけど， これを混ぜるのはなぁ ……
  • 昼 (1400): お茶部屋にて マーボーナス雑炊． 缶詰もののマーボーナス (100円，北大生協) を温め， 弁当として持って来た飯 (0.5 合ぶん) を投入． じゃんく． だけどおいしい．
  • 晩 (2300): 朝の味噌汁の残りに 油揚げとネギを追加したもの．