ぎょーむ日誌 2002-08-22

2002 年 08 月 22 日 (木)

• 0940 起床． いやはや． 0950 自宅発． 曇． 1000 研究室着．
• 昨日調べた統計学 freeware R の挙動に関して， ネット上の文書をいろいろと眺めてみる． 同じ S 言語で動く S-Plus 用の文書も有効である．
• なんとなく Math::MatrixReal をインストールする …… なんとなんと， Perl で「演算子上書き」という悪夢のごときワザが使えるとは知らなかった ……!! 標準ライブラリ， use overload．
• 世の中には CPAN-Mac などというものもあるらしい． これまたすごいワザだ．
• いやはや， 一昨日 の「ぎょーむ日誌」がトんでたことに気づく． 生成機構のバグというか操作ミスだった． fool-proof になってないなぁ． なんともはらだたしい． いらいらさせられる． とりあえず「間違いにくくする」応急処置． そのうち機会をみて このあたりのシステムの安全性を改善しないとだめだな． バックアップから失われたデータを完全復旧する． ついでに hosho バックアップスクリプトの小さなバグも修復．
• 空腹になったけど昼飯の弁当がないので， 北大生協に買いだしにいく． 1440 かとーオフィス発． ついでに図書館に本を返却する． 北大生協書籍部をふらふら． 「樹に咲く花 (合弁花・単子葉・裸子植物の巻)」 という図鑑を買う． 表紙中央に我が宿敵たる Abies の枝分かれ構造がはっきり写ってる写真があったんで …… Abies にかぎらず樹冠全体像のカタチがわかる写真が多く掲載されている (特に針葉樹に関して)．

[その図鑑から勝手に引用] 石槌山のシラビソだそーで．

それから例の 「空間データモデリング」 も買う． 1520 かとーオフィスにもどる．
• お茶部屋雑談などやっててなかなか本業にかかれない．
• 苫小牧の田中健太さんからメイルいただき， 来月に富良野の 東大演習林 で LTER の研究会があるんで参加しませんか， とのご案内いただく． 参加してみることにする．
• 牧さん (すでに道北に移動)・下野さんとロジスティック回帰メイルあれこれ．
• 夕方すぎになって， にせアカマツどらふとの続きを再開する． あいかわらずぢりぢりたる進捗なれど 構成転換によってこの対象に関する最善説明法は確立しているんで， 気分的にはひどく楽だ． 言うなれば， 地図なき山岳地帯の踏破を終えて， あとはひたすらに林道をくだり歩いてふもとの町をめざすばかり．
• 以前に書いた結果説明を再利用できるところまで合流して本日は撤退． 2100 研究室発． 2115 帰宅． 体重 72.0kg．
• 晩飯前に運動． 晩飯．
• 一般化してみた logistic な二項分布最尤推定って， ぢつは案外かんたんに数値計算できるんじゃないの? という昼間ふと思いついたアイデアを検討 …… Ek = exp(-Σibixi,k) として数式をヒネってみる．
  • 最大化する対数尤度： log(L) = Σk=1N-M log(Ek) - Σk=1N log(1 + Ek) + (定数項)
  • その偏微分： Fi = ∂log(L)/∂bi = -Σk=1N-M xi,k + Σk=1N xi,kEk /(1 + Ek)
  • そのまた偏微分： Gi,j = ∂Fi/∂bj = -Σk=1N xi,kxj,kEk /(1 + Ek)2
  おお， なんかイケてますよ， これは．
• [今日の素読]
  • NHK ラジオ「ビジネス英会話」 Lesson 21 ``Pet Care'' (8/19-8/23)．
• [今日の運動]
  • 腹筋運動 30 × 3 回．
  • 自宅空自転車． 2145 より 28 分間 300kcal． ストレッチング．
• [今日の食卓]
  • 朝 (1020): 研究室お茶部屋． コーヒーだけ， か．
  • 昼 (1540): 研究室お茶部屋． パン 2 枚． ブルーベリージャム． コーヒー．
  • 晩 (2240): 米麦 0.7 合． モロヘイヤ・ナス・マイタケの炒めもの．