ぎょーむ日誌 2000-02-28

2000 年 02 月 28 日 (月)

• 0930 起床． 曇り． 燃えるゴミ出す．
• 久しぶりに体重測ってみる． 70.6Kg． うーむ， 運動不足のわりには よく食べてるってことかな．
• 最尤推定法プログラムの改良． はい， まだやってます． 尤度関数の計算をモデル依存な部分と 非依存な部分に分離した． これで 「あたし Perl しか教えてもらってません」 と抗議するヒトにもプログラムの改良が容易になった． 利用者が勝手に決めるモデルに依存した部分って これだけだもんな．

  //-- FILE: function.cc //-- CODER: KUBO Takuya (kubo@ees.hokudai.ac.jp) //-- //-- (C)opyright 2000 KUBO Takuya (kubo@ees.hokudai.ac.jp). //-- All rights are reserved. //-- Do not modify without amending copyright, //-- distribute freely but retain this paragraph. //-- //-- CHANGE: Mon Feb 28 10:43:20 2000 // MODIFY: Sat Feb 5 17:58:44 2000 #include "ml.h" double function( Plist* p, ItrSample itr_sample ) { double y = ( *itr_sample )[ 0 ]; // the first column double x = ( *itr_sample )[ 1 ]; // the second column double h = 1.0 / ( 1.0 / p->Val( 2 ) + 1.0 / p->Val( 0 ) * pow( x, - p->Val( 1 ) ) ); // height of tree given by an allometric equation return - ( y - h ) * ( y - h ); // + ( term without p ) }

• 全体 660 行のうちこの 23 行だけ． ついこないだまで二階微分を手で計算してたのに， ずいぶん進歩したもんだ． と， 一人でちょっと感慨にふける．
• 1155 自宅発． 雪． また積もりそう． 1220 JR 琴似駅発． 北七条郵便局に立ち寄ってみるも ThinkPad560E は未だ届かず． 1235 研究室着．
• 浦口さんに尤度比検定プログラムについて 説明させてもらったんだけど， 昨日・今朝直したところで新しいバグが 生じていたのであせる． とりあえず一通り説明．
• 甲山さんが以前とりまとめて投稿してしまった 死亡率の最尤推定法に関する短報の査読結果が 返って来たので対応を考える． 基本的には問題ない． 私がやらねばならぬ修正作業は 推定値のばらつきの偏りに関する二項検定と， 最尤推定プログラムコードをネット上で公開する準備 である．
• 昼飯の弁当食ってから， 浦口さんの前で明らかになってしまった バグ直しに着手し始めたところで， A 棟 3 階に君臨する女王・ 牧さんから電話 「ひまそーな男を集めてください． 私の車が雪で動かなくなりました」
• 清野さんと私の二人でスコップをかついで出動． 目的地まであと少し． ああ， 向うに 歩道をさえぎり車道に頭を出したまま停止している いかにもハタ迷惑な車が． 近づくと牧さんがその周りで雪をひっかいては 途方にくれていたのであった．
• 駐車場から歩道に出ようとして， 除雪で生じた 畝条の段差を強行突破しようとしたものの， 乗り越えた瞬間に車体がその段差によって 持ち上げられて駆動前輪が用をなさなくなった， ということらしい． 清野さん発案によって持参してきた 二本の直線的スコップにて 車体浮かせている雪こぶを破壊し， 前輪下に商品名「おた助くん」なる ヘルパー (ゴム製の板で車輪の下に敷いて摩擦力を増す器材) 噛ませ 全力後進ならびに人力推進の合力をもって ようよう窮地を脱したのであった．
• 聞けば牧さんは， 動揺のあまり 何時間もずっと車の周りで うろうろしておられたとのこと． 低温と疲労で すっかり消耗しておられる． 低体温症による遭難の一歩手前， というところか．
• 時刻はすでに 1630 すぎなので， 郵便局にもう一度行ってみる． ついに英語キーボード ThinkPad560E が届いている． さっそく帰って動作確認． よし問題なし． しかしこいつに Linux を入れるのはかなり難しい， と事前にわかっているので， すぐにはインストールに着手しないことに． とはいえ， ついついあれこれと遊んでしまう．
• 気がつけばすでに空腹になって来たので帰宅することに． ああ仕事が片付いていない． 2210 研究室発． 2229 JR 札幌駅発． なーんか異様に混んでる車内． 2310 帰宅．
• 腹が減ったので飯が炊けるのを待つ間， モヤシ炒めを作って空腹をしのぐ．
• 晩飯食って一息ついたので， 死亡率の偏りに関する二項検定をやってみる． 以前作った自作の検定プログラムを少し修正して 確率を計算してみる． 実現確率 10 のマイナス28乗未満の帰無仮説を棄却する． なぜこういう偏りが生じるのか うまい説明は思い付かないのだが．
• どんどん朝が近付いてくるので， 寝たほうが良さそうな気がしてくる．
• 今日の食卓
  • 朝 (0950): 昨夜の味噌汁の残り． ホウレンソウとエノキダケの 炒めもの．
  • 昼 (1430): 朝の残りの ホウレンソウとエノキダケの 炒めものをおかずにした弁当．
  • 晩 (2340): 味噌汁の残りと ネギ・モヤシの炒めもの．