ぎょーむ日誌 2002-09-24

2002 年 09 月 24 日 (火)

• 0710 起床． シャワー． 朝飯・弁当の準備． 朝飯． 0745 自宅発． 曇． 0755 研究室着．
• 「ぎょーむ日誌」まわりの CSS を少し変更してみる …… がやはり NN4.x では表示がかなりヘンになるのでやめ．
• にせアカマツどらふと， もう少しだけ直してから再アップロード．
• 昼飯． インドネシアから帰ってきた塩寺さんハナシとか， 赤岩青巌峡谷で岩登りデビューした浦口さんハナシとかうかがう．
• 一般化線形ロジスティック回帰がとても簡単にできる logistic.pm の使いかたを浦口さんにご説明． いつも久保ソフトウェアの評価試験者となってくださる， じつに奇特な存在だ．
• その後， このモジュールに汎用性を付加しようとする改造． 分散が一定ではない， という統計モデルの組み込みが難しい …… 分散の関数の reference つまり \&function で 渡せばいいか， と思ってたんだけど， コトはそれほど単純ではない． その一階微分・二階微分も計算せんといかんわけで．
• そもそも分散のパラメーターが， この一般化線形モデルの世界では異端， というか． パラメーター行列のどこにほうりこめばよいのやら．
• 小林さんから返信 …… いつもながら速すぎ …… こちらもなんとか対応してみようと苦慮したり．
• あれこれと交渉の結果として， このような割り算指標図を復帰させる (小林案)， そのかわり怪しげ検定は廃止 (久保案)， というあたりで合意に達する．

[これも] Perl+Gnuplot 作画機構を使ったもの． ラベルなんか自動位置決めをやってる． T#3 は少しずらしたかったので 「もし T#3 なら右に 0.3 ポイント移動」 という命令を追加．

• Perl 最尤推定法プログラムに関して， しょうがないので分散一定の場合の GaussianSample モジュールを作る． これを汎用 Newton-Raphson 計算機構にほうりこんでしまえば， 結果が出てくるというところは LogisticSample と同様である．
• 最小二乗法で計算で記述すればすぐにコード書けるところを， あえて最尤法にこだわり …… よーし， とりあえず計算結果得られた． 2030 研究室発． 2040 帰宅． 体重 72.0kg．
• 夜中の 2600 までひたすら計算ばかり． 今度は分散が変化する場合なんだけど …… さすがに真夜中すぎるとアタマがくらくらしてくる．
• [今日の素読]
  • ということで今日はナシ．
• [今日の運動]
  • こっちも．
• [今日の食卓]
  • 朝 (0720): 米麦 0.7 合． コマツナ・エノキダケあえもの．
  • 昼 (1240): 弁当． 研究室お茶部屋． 米麦 0.7 合． コマツナ・エノキダケあえもの．
  • 晩 (2140): 米麦 0.7 合． ターツァイ・マイタケ・油揚の炒めもの．