ぎょーむ日誌 2002-10-04

2002 年 10 月 04 日 (金)

• 0830 起床． 自業自得とはいえ眠い． シャワー． 0840 自宅発． 曇． 0850 研究室着． 朝飯． コーヒー．
• にせアカマツどらふと読み直し． あとはちょっと修正して敵性語検閲とおしてから投稿するだけ …… のはず．
• 文献引用の書式を動的に決定する .bst ファイルのカスタマイズ …… なぜか上のようなペイジは参照しないで， 構文から文法を自力で推定する， という方法に没頭してしまった． 読めばだんだんわかってくるわけだが …… このいかれた文法は誰の趣味なんだろうか． Lamport?
• 小林さんに要旨前半の実験に関する部分を書いていただいたんだけど， やはりワタクシ的には少し高邁ではないかとの懸念が …… つまりひじょーに正確ではあるけれどやや難解に思えて， ですね． 当方の趣味でかなりしろーとじみたものに勝手に「翻案」してみる． たいていの読者は題名とせいぜい要旨ぐらいしか読まないんで (ついでに図表とかながめる)， ここは少し努力を傾注してよいところだ． 要旨書いててヤなのは時制をどうすべきかよくわからん， というところ (過去形つかうなという暴論もある)．
• 1350 ひとくぎりついたんでやや遅い昼飯．
• 昼飯後ももうちょっと直してから make pdf． 連絡 & コメントメイル書き． 1520 おわり． お茶部屋でお茶など飲みつつ雑談．
• 粕谷さんから出された Mathematica 宿題に取り組む …… 内容はぎょーかいの標準統計学本 Sokal & Rohlf (1995) Biometry (3rd edition) 418 ペイジにのっているような Box-Cox 変換のパラメーターλと その対数尤度 (正規分布へのあてはまりの良さ，か?) の作図問題にすぎないんだけど …… λちいさくしていくとみょーな虚数解とかが出るんだよね．
• 2 時間ほどあれこれ試行錯誤してみて variance の計算に問題がある， と特定できた． λが小さいときに Box-Cox 変換した量の絶対値は小さくなるんで

  Var(x) = E(x²) - E(x)²

  方式で計算すると (ふつーはこの式で計算する)， この差がうちきり誤差にうもれてしまうことがある (そして Mathematica における実数の取り扱いには かなり気をつけないといけない)．
• 応急の対応策として

  Var(x) = E((x - E(x))²)

  というふうに variance を計算してみるように変更した． これで作図に関しては， とりあえず問題ない．
• お茶部屋に行くと， かとー先生たちが野外作業終了ビール飲んでいて (だれかの差入れかも)， 北大内所在の某研究林は「終っている」とのこと．
• 敵国語メイルを書く． 短いものなのにいつものごとく 30 分以上かかる． 2050 研究室発． 2100 帰宅． 体重 72.2kg．
• 晩飯後いきなり寝てしまった． 昨晩の夜ふかし読書のたたりか．
• [今日の素読]
  • ということで生活習慣がなかなか修復しない．
• [今日の運動]
  • 腹筋運動 30 × 3 回．
  • こちらも，だ．
• [今日の食卓]
  • 朝 (0910): 研究室お茶部屋． パン． コーヒー．
  • 昼 (1350): 弁当． 研究室お茶部屋． 米麦 0.7 合． モロヘイヤ・エノキダケ・豆腐の炒めもの．
  • 晩 (2140): 米麦 0.7 合． タマネギ・油揚・煮干の味噌汁． タマネギ・エノキダケ・卵の炒めもの．