ぎょーむ日誌 2003-10-29

2003 年 10 月 29 日 (水)

• 0800 起床． 朝飯． コーヒー． 0850 自宅発． 曇． 0900 研究室着．
• 苫小牧 D3 鍋嶋さんから出されている宿題をこなすべく， とりあえず課題論文にとりくんでみる …… 多変量解析ねえ． 使うにしろ使わんにしろ勉強せんといかん技法なんでしょうなあ． なんでこのへんが好きではないかとゆーと， これらの手法の確率・統計学的側面がかなりいいかげんに思えるからだ． にゅーろとか SVM とか発展してくると， こういうのって衰退していくんぢゃないの?
• とはいえ， 使ってるヒトたちはたくさんいるわけで．
• 岩熊研 (陸水) COE PD の石川さんの奇襲攻撃 …… 問答無用でまるっきり別の統計学的問題の鉄砲水に ざざーと流されていく．
• すでにいろいろ試行錯誤されてるよーで …… 群馬大の 青木さん R ペイジ は大人気ですなぁ． しかし， うまく動かなかった， と
• とりあえず青木関数の換骨奪胎からとりくむ …… む， R における無名関数とかイキなわざが使われてるぞ． しかしこの折れ線回帰， 二直線をまったく独立に推定してるぢゃん (分割点において二直線が接続してない) …… というのが判明したところで午前終了．
• 雨のあいまをぬって北大構内走． 昼飯． D3 谷さんがタニギキョウの発芽実生をお茶部屋にもってくる． なんでも， 東研の院生がオシダ (これも谷さんの研究対象のひとつ) の地下にのみ造巣するアリを研究しており， その巣の中にはこのタニギキョウの種子ばかりがある， とのこと． ふーむ， 不思議な ……
• 午後も R プログラミング． 青木関数の改竄は放棄してゼロから作り直す． ところどころかとー先生に助けていただく．
  • なんであれ階層構造のあるデータ構造は list …… list の中は何でもよい
  • tag つき list をうまく使う
  • 変数のスコープ， ブロックの中のブロックに対してぐろーばるナリ …… ケガれてるけど便利だ
  • 複数行に一個の式を書きたけりゃ (...) でくくるしかない?
  • sort() より order() のほうがよい場合がある
  • 最尤推定はちょい面倒になると optim() する以外に策がない …… っつーか， optim() 初使用 (どきどき)

[hockey regression?]

可能なすべての折れ曲がり点において 最尤推定を行う． 次に全折れ曲がり点の中で最大の最大化対数尤度を 実現している点を探索する． 最後に折れ曲がり点ナシモデルと比較してから 作図する …… という手順をふんでいく R プログラムを書いた．

• 1600 計算とかひとまず終了． アップロードできるようにする． 計算メモは書くことがない． しかしプログラムみなおしてると， 間違いがぼろぼろと見つかっていく．
• 1930 研究室発． 2000 帰宅． 体重 72.0kg． 晩飯．
• 晩飯後も間違い直しがつづく． 幸いにして間違い直すたびに結果はマシになる． しかし， もういいかげんに「新発見」は止まってほしー ……
• [今日の素読]
  • Salsburg, D. 2001. `` The Lady Tasting Tea -- How statistics revolutionized science in the twentieth century''. Owl Book.
  • Chapter 18. Does Smoking Cause Cancer?
  • In 1958, R.A. Fisher published a paper entitled ``Cigarettes, Cancer, and Statistics'' in the Centennial Review, and two papers in Nature entitled ``Lung Cancer and Cigarettes?'' and ``Cancer and Smoking.'' He then pulled these together, along with an extensive preface, in a pamphlet entitled ``Smoking: the Cancer Controversy. Some Attempts to Assess the Evidence.'' In these papers, Fisher (who was ofted phtographed smoking a pipe) insisted that the evidence purported to show that smoking caused lung cancer was badly flawed.
• [今日の運動]
  • 北大構内走 1305-1340． ストレッチング．
• [今日の食卓]
  • 朝 (0800): 飯たくの忘れてた …… コーンフレイクス．
  • 昼 (1350): 研究室お茶部屋． 「北欧」サンドイッチ．
  • 晩 (2130): うどん． ハクサイ・タマネギ・ワカメ・豆腐・煮干の味噌汁． さつまあげ．