ぎょーむ日誌 2001-01-22
2001 年 01 月 22 日 (月)
- 0730 起床．
昨晩，
寝る前に魔がさしてメイルを取ったところ，
パラメーター推定に関するなかなか難しい質問来て，
目がさめてしまったのである．
やはり 2400 すぎたらメイル読んではいかん．
- 朝飯・弁当の準備．
朝飯．
コーヒー．
眠い．
- 0805 自宅発．
曇．
気温高い．
0822 東京モノレイル流通センター発．
浜松町の気温 5 ℃．
0840 研究所着．
- 気になってた推定問題に関して，
浜松町に置いてる統計学本調べてみる
……
まあ，
一般的なコトはこれでわかったような気分に
なれるんだが．
- 調べているうちに，
しだいに不安になってきたので，
検定 Perl スクリプトを一個書いてみることにする．
- 質問されたのは，
パラメーター推定のときの外れ値の取り扱いである．
調べてみると，
確率分布として正規分布を仮定している場合には，
Grubbs-Smirnov 検定
(Smirnov-Grubbs 検定?)
すれば良い，
とわかった．
- 北大・宮本さんの
統計学 Perl モジュール
を利用させていただいて
(t-値の計算が必要なのである)，
その Grubbs 検定コードを作っていく．
- どうにか動いて，
試験運転も終わったのは 1250．
ふう．
試験運転では Box-Muller 法で発生させた
標準正規乱数のサンプルによけーな外れ値を
つけたデータ読ませて調べてみた．
gs_test.pl は stepwise 的に
このよけーな値を取り除いていくか．
調べてみると，
正規乱数標本に付け加えたよけーな値が
そんなに外れていない場合は，
当然ながら「外れ値」とは認定されない．
- ふう．
つかれた．
昼飯にしよう．
- 飯を食いつつ細かい手直し．
飯食ったら歯磨き．
今日も歯医者に行かねばならない．
ああ．
- メイル連絡して，
外れ値話はいったんおしまい．
- 今週は領域セミナーでしゃべらないとダメなんで
その準備と，
あと研究計画書なるものを提出しなければならないんで，
ちょっとあせりぎみである
(どちらも何の準備もしてないんで)．
- しかしながら，
先週末にやりかけていた，
開空度計算の評価試験のとっかかり部分は
もうちょっとだけ進めておこう．
今日，
何とかそれに一区切りつけておかないと，
もはや月末までそれにかかわってる
時間がなさそうである．
- 小川シミュレイターを観測点開空度評価モードで
呼び出せるような Perl スクリプトを書いていると，
北大から「外れ値」追求メイル．
なけなしの勇気をふるって
時間稼ぎメイルを送ってその場をしのぎつつ，
午前中に作った gs_test.pl の
入出力まわりをイッキに改良．
すごくぷれっしゃーのかかる作業でした．
<鳳翔>の FTP サイトに置いちゃって，
解説メイル送信．
gs_test.pl で使わせてもらった
Pstatlib 作者の宮本さんにも Cc で送る．
どうにかこうにか終了したのが 1605．
- まぁ，
今までデータ解析する上で「外れ値」ってのを
放置しておいたむくいでこのような
ドロ縄作業になってしまったわけだ．
とはいえ，
Grubbs-Smirnov 検定って
正規分布の母集団から標本した場合しか使えんからなぁ．
考え方を応用すれば，
他の確率分布に関しても検定を構成できるんだろうか．
- 小川「とりあえずここまでのとりまとめ」
作業に復帰．
- さて，
結果出そろったんが
……
なんかうまく描画できんな．
あ，
つまらんミス見つけた．
三つのテイブルを結合するところだ．
しかしもはや歯医者に行かねばならぬ時刻．
1740 研究所発．
1746 JR 浜松町発．
1756 JR 大森着．
ここまで 160 円．
1800 駅の近くの歯科医院着．
- 待つことしばし．
今日もすごい勢いで患者が回転している．
1825 診療．
今日は麻酔なしで処理されたんで，
自分の歯がどうされてるか，
よくわかった．
つまりドリルででかい穴開けられて，
そこに薬剤つめられて，
何か可塑的な素材で一時的にフタされてる，
という状況だったのである．
今日もそれが繰り返された．
しばらくは歯科医通いになりそうである．
- 1835 診療終了．
今日は 410 円．
次は木曜日の 1815 か．
1855 帰宅．
- (一言近況)
虫歯で歯医者通いやってます．
痛いのを長期放置しておいたら，
なかなかひどいことになりました．
もっと早めに行けばよかった．
(20010122)
- う．
Grubbs-Smirnov 検定づかれか
歯医者づかれか，
晩飯食ったら 3 時間ほど寝てしまった．
- 私が寝てる間に，
北大では gs_test.pl の改良がすすみ，
直線回帰対応版までが開発されてしまったようである．
うーむ，
おそるべし
……
- 真夜中ちかくなってから，
夕方の作業の続き．
いかんな，
これは．
ちゃんと比べてみると，
三次元モデルを使った開空度の計算値と
小川の観測値があまり一致してない，
とわかった．
すなわち三次元モデルがあまりよろしくないのである．
- そこから 0400 ごろまでパラメーター調節など．
全然ダメ，
という領域からはぬけだす．
が，
やはりよくない．
図なんかのアップロードは明日，
ということで．
- やはり計算方法変えるかな．
しかし変えても
良くなるのか，
判然としないトコロが．
- う，
「今日の朗読」
できなかった．
蝋山本のおそるべき第 1 章も
あとは短い 1.8 と 1.9 で終わるのに
……
- もう 0500．
さっさと寝るか．
- 今日の食卓
- 朝 (0740):
米 0.5 合．
豆腐・キャベツ・エノキダケを炊いたもの．
麺つゆ味．
- 昼 (1300):
弁当．
米 0.8 合．
朝と同じ．
- 晩 (1920):
米 0.7 合．
朝と同じ．