ぎょーむ日誌 2009-08-27

2009 年 08 月 27 日 (木)

• 0730 起床． ねむい． 朝飯． チャイ． 0850 自宅発． 曇． 0910 研究室着． コーヒー．
• 今日は大学院入試なのだが試験監督もなければ， 採点ぎょーむもナシ． ただし， 得点のとりまとめ， という雑用を命じられた． すでに残業弁当まで発注されている．
• うーむ， ぎょーむ起動状態がよくない ……
• 関係ないけど， /etc/modprobe.conf で

  alias net-pf-10 off
  alias ipv6 off
  

  としてみる．
• データ解析こんさるメイル． ふーむ，

  Hartigan, JA and Hartigan, PM. 1985. The dip test of unimodality. Annals of Statistics 13:70-84

  …… これですか．
• 例によって， R の CRAN には diptest なる package があるのだが (library(help = diptest)) …… これは上記論文 Hartigan & Hartigan (1985) を読まないかぎり， 何をやってるのかまったく理解できない package．
• で， Hartigan & Hartigan (1985) の Table 1 に 「標本数が X なら sqrt(X) * dip で補間せよ」 ってなコトが書かれてて， 「何だ，こりゃ?」 と思ったのだが …… 上記 library(diptest) をヒネくりまわしているうちに， こういう意味ぢゃないかな， とわかった．

  > library(diptest)
  > qDiptab[,"0.95"]
          4         5         6         7         8         9        10        15        20 
  0.2072698 0.1863918 0.1647696 0.1599034 0.1539783 0.1466038 0.1396114 0.1187608 0.1051302 
         30        50       100       200       500      1000      2000      5000 
  0.0881689 0.0702738 0.0511279 0.0368418 0.0237295 0.0169344 0.0120381 0.0076506 
  > x <- as.numeric(rownames(qDiptab))
  > plot(sqrt(x), 1 / qDiptab[,"0.95"], type = "b")
  > grid()
  

• また別件のデータ解析こんさるメイル． 極値分布 うーむ ……
• そして 15 時からまた別のデータ解析こんさる ……
• これは一時間弱で終了したので， さっさとサル論文原稿だの岩波本作業にもどらなくては …… と思ったのだけど， 例によって「ベイズ化による抜本的な改善」 案がおもいついてしまって， そちら方面にはまりこんでしまった ……
• そうこうしているうちに， 大学院入試の採点結果集計作業がはじまる． といっても， 私は全科目合計 200 枚ばかりの答案用紙に書かれた 受験番号と得点を読み上げる (そして工藤さんが入力していく)， ばらばらと届く採点結果をひたすらまつ， といったぎょーむにすぎなかったんだけど．
• 1930 ごろ， めでたく終了． 残業晩飯．
• またベイズ化の問題にはまる …… うーむ， このゐんばぐすの trap 窓は …… で苦闘したのだけど， 本日は解決しないまま撤退． 1850 研究室発． 1900 帰宅．
• 部屋にもどって 「うーん， MCMC が混交しないよー」 とウナされつつごろごろしてると …… とつぜん解決策がわかってやってみたらうまくいった (こういう場合はうまくいかない事例のほうが圧倒的に多いんだけど)． いやー， すごい!
• しかし ThinkPad はいよいよ不調だ． さっさと system board (メイン基板) を交換したほうがよいのかな ……
• [今日の運動]
  • うんどう休養日．
• [今日の食卓]
  • 朝 (0740): 米麦 0.6 合． ナス・油揚の味噌汁．
  • 昼 (1230): 研究室． 米麦 0.8 合． 梅干．
  • 晩 (1940): 残業晩飯． 北大生協に発注されたもの．