ぎょーむ日誌 2006-01-20

2006 年 01 月 20 日 (金)

• 0700 起床． 東京大手町の宿． 飛行機出発まで時間あるけど， こんな都心部にいてもやることないので， さっさと撤収することに． 宿 1F で無料給食やってるんでそこで朝飯． このあたりも東横イン的だな． よく知らないんだけど， 最近のびぢねすホテルとやらは全部こうなのか? 東横インだと「型ぬき」握り飯 + みそ汁 or バターロール + コーヒー， ここは安そうなパン三種 + ぬりもの + サラダ + ゆで卵 + コーヒー． 0740 チェックアウト． 一泊 10350 円 (!)， まあこの値段なので朝飯にも原価で 150 円ぶんぐらいの上積みしてますってことだろう．

[大手町]

左から昨日の会議場所の三菱総研， 日経新聞本社， 経団連会館． 東京駅にむかう途中 ここ からの南方の視界なんだけど， 地図みたら JR 神田にいったほうがよっぽど近いぢゃん． いやはや．

• 0801 JR 東京駅発． 0806 浜松町で東京モノレイルに乗り換え． ここまで 150 円． ここから 470 円． 0812 同発の 0837 羽田空港第二ターミナル着． ここでチェックインしてみると， 予定より一時間早い便である 0900 発に乗れそうなんでそちらに変更．
• すでに搭乗がはじまってるので急ぐ． こういうときにかぎって， ですね． 私はカバンの中とかに検問でひっかかるモノはいれないように注意してるんだが …… 近ごろの荷物チェック流行はなんと「ツメ切り」． 二度にわたる X 線照射ではツメ切りを十分に破壊できなかったとのことで， ツメ切りだせと命じられた． 動作緩慢なる係員あれこれとひねくりまわしたあげく， これは航空テロ犯たちが愛用している 大量殺戮型ツメ切りとはいささか異なるタイプらしいとの見解に到達したよーで …… よーやくのことで釈放された．
• 今日の運動: 羽田空港第二ターミナルビル走 300m， 荷物つき．
• 0857 ANA 55 便 B747-400 搭乗． 今日は最後尾の席にしてみたんだけど， このあたりは人口密度が少なくていいね． 0918 羽田空港 C 滑走路を南→北と離陸． 晴． すでにばてぎみ． 1026 新千歳空港着陸． 1051 快速エアポートで同発． 1131 札幌着． ここまで 1040 円． 本日も 5 分おくれぐらいで運行中． 1145 研究室着． ばてた．
• ばてつつ昼飯．
• 鈴木牧画伯が イヌ植物だじゃれ年賀状 公開しておられる． 牧さんわーるどというかんぢで良いですよね．
• とーとつに R のハナシ． lmer() という GLMM solver が Matrix package 内に入ってるな． 挙動を調べてみなくては． しかし， ここに書いててもどわすれしそうだ．
• 900 億回 乱数を発生させる必要のある計算に一時間ぐらいかかるんだけど， もっと速くなる計算アルゴリズムない? との質問くる． そんなのあるわけないぢゃん． 乱数発生だけでも一秒間に 2500 万回になる …… ホントに一時間で計算できてんのか?
• と思ったら， ぢつは私の ThinkPad X31 (Pentium M 1.6GHz) ですら 1 秒未満で 1000 万個の乱数を発生させることができた． R で調べると，

  > system.time(runif(10^7, 0, 1))
  [1] 0.97 0.11 1.08 0.00 0.00
  

  とすると 900 億個だと 150 分間か． ところが， ここで単純に R で一億回乱数発生やらせると，

  > system.time(runif(10^8, 0, 1))
  [1] 20.06  8.43 55.26  0.00  0.00
  

  10 秒ではなく 20 秒かかる． これは計算結果格納 vector のメモリ確保とかに時間くっているのか?
• ところで， そんなに莫大なる個数の乱数発生させるならちゃんとした 擬似乱数 発生法つかってるんでしょうね (R の場合は めるせんぬついすたー) ……
• ふーむ， 多項版の超幾何分布は The Multivariate Hypergeometric Distribution とゆーのか．
• などとうだうだやってて， 来週の授業の準備はススまず． 1855 研究室発． 1910 帰宅． 晩飯．
• [今日の運動]
  • あいかわらずな
• [今日の食卓]
  • 朝 (0720): 宿の無料給食．
  • 昼 (1250): 研究室お茶部屋． こんびに握り飯二個．
  • 晩 (1940): 食パン． ハクサイ・モヤシ・エノキダケ・卵の炒めもの．