ぎょーむ日誌 2003-01-17

2003 年 01 月 17 日 (金)

• 0820 起床． もう金曜日か． 朝飯． コーヒー． 0900 自宅発． 雪． 0910 研究室着．
• 今日は空間不均質推定問題， 推定部にやっとのことで手をつけられる． とはいえ， しばらくは均質環境での推定試験が続くわけだが．
• 1100 よーやく環境が均質な場合の推定部ができた． このように一様な場合は， とうぜんうまく推定できる． 尤度最大化は Newton-Raphson 法．
• 1220 なんとか最初の sector 分割ができたようだ． とりあえず北大構内雪上走． 1230 研究室発． 晴． 大脳を物理的にゆすってるうちにバグに気づいた． 1305 研究室もどる．
• バグをなおしてみる． まあ， まずは一回のセクター分割． 計算が遅いんで標本数を 6000 から 2000 点に 減らしてみたんだけど， なかなか時間がかかる ……

[ホントは均質環境だけど] 標本数が多いのでふたつに分割してしまった． これはまあしょうがないのかな． ホントの発生確率は 0.12 (全体) なんだけど， 左を 0.12 右を 0.09 としている． モデル選択基準は AIC．

• で， こんどは複数回のセクター分割やらせようとすると， 消したはずのセクターが幽霊のごとく復活しちゃうという 呪われたバグ． いやはや． 昼飯．
• バグとりにちょっとてこずったのちに， 複数回分割の試験． うーむ．

[これもホントは均質] ともかく「切れるところまで切ってみろ」 としてやるとこういう分割をやってくれた． 二次元一様分布， とはこういう不均質性を生成するものなんだろうか? いや，やはり推定のやりかたを改善すべきだろうな

• セクター分割に関してまた手直しして 1720 不均質環境推定の問題を解かせてみる． セクター分割だけ， でどこまで推定できるかどうかを見てみたわけだが ……

[こんなものなのかなぁ] 上が試験問題で下がその回答． なんとなーく， それらしい推定しようとしてるらしい， というのはわかるんだが． 点の個数は 6000． これぐらい多くないと推定できん．

• ということでセクター分割だけでなく， セクター結合という操作も必要のようで．
• 結合は分割よりさらにめんどうで， なかなかプログラミングすすまず．
• M1 大澤君がややこしい作図問題を持ってきて， 今日はそれを片づけておわり． 2100 研究室発． 2120 帰宅． 体重 73.2kg． 晩飯．
• すすまん．
• [今日の運動]
  • 北大構内走 1230-1305． ストレッチング．
• [今日の食卓]
  • 朝 (0830): 昨晩の残り．
  • 昼 (1420): 弁当． 研究室お茶部屋． 米麦 0.8 合． ホウレンソウ・ブナシメジのあえもの．
  • 晩 (2150): 米麦 0.8 合． ハクサイ・ニラ・ワカメ・豆腐の雑炊．