ぎょーむ日誌 2010-07-(11-20)
2010 年 07 月 11 日 (日)
-
0840 起床．
首のうしろがカユくてよく眠れなかった
……
朝飯の準備．
朝飯．
コーヒー．
ちょっと怠業．
-
1100 自宅発北大構内走．
今日も曇だけどときどき直射日光．
1145 帰宅．
体重 71.2 kg．
-
13 時すぎにまた自転車二台で自宅発．
大通公園ちかくで昼飯．
-
その後，
またまゆみさんの円山道楽 & うんどう不足解消のおともとして，
円山周辺をうろうろ．
[円山道楽]
こういう円山な古道具屋で椅子やテイブルを吟味すること．
私は店の前の路上に出されていた木製椅子にすわって，
道ゆく人々をながめておりました．
-
その後，
円山北側に位置する道内最大神社である北海道神宮を見物して，
豆腐など買って円山から撤収．
北大構内は高速南北移動に便利．
-
16 時すぎに北 18 西 2 の北辰中学校で参議院投票をすませて，
帰宅．
ちょっと休憩．
-
買いもの．
晩飯の準備．
晩飯．
-
今までナゾに包まれていた
伊庭さんの mixi 日記，
その一端 (数学・科学エッセイもの) が
はてなダイアリー
にて公開進行中!!
伊庭さんタイピングの秘密，
中国暦法の超絶 fitting，
赤池さんのゐすきい選択規準，
算数難題中の保存測，
高次元わーぷアニーリング
……
どれもおもしろいハナシです．
-
JR 札幌駅すぐ北の小映画館
蠍座
で
マイレージ マイライフ
をみた．
アメリカあちこちとびまわる
「りすとら宣告人」
コメディ．
これってホントにそんな職業あるのか?
公儀介錯人
みたいな架空の職業?
といった虚構性をホネぐみにして世界を構築してるのがよいですね．
北米各地の空撮も楽しい．
-
[今日の運動]
-
北大構内走 45 分間．
自転車 1.5 時間ぐらい．
-
[今日の食卓]
- 朝 (0920):
昨日，
南 1 西 11 のパン屋「ピノ」で買ったバタール．
トマト・モッツァレラチーズのサラダ．
ソーセイジ・卵の炒めもの．
- 昼 (1340):
[さえらのサンドイッチ]
大通公園すぐ南の地下三階，
サンドイッチ専門店．
タラバガニとパストラミビーフのくみあわせ．
- 晩 (1850):
米 0.8 合．
トマト．
北海道神宮横の豆腐屋「めめ」で買ったゴマ豆腐．
コマツナ・ニンジン・ショウガ・豚肉の炒めもの．
2010 年 07 月 12 日 (月)
-
0800 起床．
朝飯の準備．
朝飯．
コーヒー．
0900 自宅発．
雨．
0915 研究室着．
-
論文原稿読み作業．
14 時すぎまで．
途中，
昼飯．
-
また別の論文よみ，
すこしばかり．
-
なんといいますか，
当然ながらといいますか，
自分で重要だと思える研究成果は，
自分で論文を書かないといけませんな．
論文かくのがおそいとか，
こんさるに時間をとられたりするから
……
といった状況に流されて論文かきをあとまわしにすると，
次第に「良い研究をしよう」という意欲が失われていく．
このあたり，
ぎょーむ日誌をつけつづけることで，
自分の時間の流れかたがよーやくにしてわかってきた，
ともうしますか．
-
とはいえ，
現時点では岩波本作業をがんばらないといけない．
土日は何も進捗しなかったしね．
すでに 1530 だが．
窓の外はあいかわらず雨．
-
第 10 章「空間構造のある階層ベイズモデル」章の書きなおし，
いったりきたりでオワってしまった．
1920 研究室発．
1940 帰宅．
晩飯．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0820):
米 0.6 合．
ネギ・油揚の味噌汁．
ゴマ豆腐．
- 昼 (1230):
研究室．
食パン．
ナッツ類．
- 晩 (2000):
なぜか焼き鳥．
2010 年 07 月 13 日 (火)
-
1850 研究室発．
1900 帰宅．
晩飯．
-
まゆみさんが極楽湯に行ってるあいだに，
ちょっととびこみデータ解析をどうこしてみる努力．
まずは図示．
library(lattice)
の
histogram()
関数による作図．
-
ゐんばぐすによる解析．
空間スケイルに依存してますな
……
-
[今日の運動]
-
うんどう休養日つづく
……
やはり，
22 時ごろにうんどうできるような時間やりくりが妥当かも
……
-
[今日の食卓]
- 朝 (0810):
ライ麦パン．
ソーセイジ．
トマト．
- 昼 (1230):
研究室．
食パン．
ナッツ類．
- 晩 (2000):
いろいろなスケジュールの都合で，
「ごはんや でんすけ」
で玄米・サバ・かきあげといった晩ごはんいただくことに．
2010 年 07 月 14 日 (水)
-
0740 起床．
朝飯の準備．
朝飯．
コーヒー．
0900 自宅発．
曇．
0910 研究室着．
-
メイルかきあれこれ．
-
そして昨晩の計算ほうこくかき．
-
……
などとうだうだしていて，
岩波本作業にとりつけたのは 14 時すぎ
……
-
第 10 章の空間ベイズモデル章のかきなおしがなかなかおわらない
……
1840 研究室発．
自転車ぱんくを発見．
家までころころ押して帰る．
1900 帰宅．
北 9 東 1 の米屋をめざしたのだけど
……
19 時閉店だった．
そのまま東にいって Homac で自転車の空気いれとか購入
北 12 西 1 の北 12 生協で晩飯の買いもの．
米も 2 kg 買う．
1950 再帰宅．
晩飯の準備．
晩飯．
-
自転車の後輪，
やはり空気ぬけではなくパンクのようで．
小径輪の弱点なのか
……?
-
[今日の運動]
-
[今日の食卓]
- 朝 (0810):
米 0.6 合．
ネギ・油揚の味噌汁．
- 昼 (1310):
研究室．
北大生協サンドイッチ．
- 晩 (2200):
米 0.8 合．
ジャガイモ・ニンジン・タマネギ・豚肉の味噌汁．
コマイ．
コンブのポン酢あえ．
2010 年 07 月 15 日 (木)
-
0740 起床．
朝飯の準備．
朝飯．
コーヒー．
0900 自宅発．
晴．
後輪パンクした自転車をころがして研究室まで輸送．
5 月のひっこし時に工具類をこちらにもってきたので．
0920 研究室着．
-
統計関連学会連合大会
の参加もうしこみをしてみた．
参加費はクレジットカード払い．
-
しかしながら，
確認メイルがとどかず．
「24 時間まて」
……
か．
あとで確認してみよう
……
と思ったら一時間ほどでとどいた．
参加申込者 2 号だったよーで．
-
昨日，
こういう図を
library(lattice)
の bwplot()
で作ってみていきづまったので，
発注者にどうしたもんでしょうかね，
とといあわせてみる．
-
とゆーことで，
また他のことは放置して，
岩波本作業のつづき．
-
しかし，
午前中は進捗よろしくなく，
昼飯．
[ちょっと分解してみる]
Carry me のパンク後輪は取りはずさないとチューブの点検とかできそうにないので，
ナットをいろいろゆるめてみたのだが
……
どうもクランク
(写真左上の歯車)
を外してチェインをゆるめないと後輪が外れないようになっている，
ように見える．
で，
クランクを取りはずすためには直径 10 mm の六角レンチ (棒スパナー)
が必要なんだけど，
そんな太いものは持ってないので，
また Homac に行って調達するしかないな
……
と終了．
-
岩波本作業のつづき．
10 章からよけーな記述をがしがし削除していく
……
-
で，
とりあえず
17 時前にひととーり終わらせてみた．
現時点の本文総ペイジ数は 278 か
……
これは多すぎるので 40-50 ペイジぐらい削りたい．
削らないと
(岩波書店の価格設定を憶測するに)
本の値段が 4000 円突破してしまうからね．
-
まあ，
削るといっても，
有用そうな一部は PDF ファイル化してネットから自由に
ダウンロードできるようにする．
とにかく，
ネットにまわせるもの
(あるいは，まわしたほうがむしろ良いもの，
おもに R や WinBUGS のコマかいハナシ)
は印刷書面には書かない，
とゆーことで
……
-
ここ二ヶ月ぐらい，
アホみたいな効率の悪さで第 7-10 章 (100 ペイジぐらいか)
を書きなおしていた
……
わけか．
-
じつは書きなおし (おそらく，そのついでに内容削除)
が必要とされている章がまだあって
……
それは 4, 5, 6 章だ．
現状はどうなっているかというと，
-
第 4 章のモデル選択・検定章:
講義のーとからあまり改訂していない．
とはいえ，
モデル選択とかについて，
これ以上あまり詳しくは書けそうにない．
一方で，
よけーな記述があるのでそのあたりは削る．
-
第 5 章の GLM 雑章．
これは GLM 説明章のまわりの何度かのリストラで
ハミだしてしまった部分のよせ集め的な
……
がんばって削ってみる．
内容は (少なくとも印刷版では)
あまり増やさないようにする．
-
第 6 章の GLMM 章．
これも講義のーとからあまり変わってないんだけど，
(GLMM については何度も書いてきたことだし)
まあいいでしょう
……
と思っていたのだけど，
階層ベイズ章を書きながらこの第 6 章をちらちらと
参照してるうちに，
どーにもいいかげんなコトが書かれているよーな気がしてきて，
ですね．
なかなかオワリが見えてこない
……
-
1820 研究室発．
買いもの．
1840 帰宅．
晩飯の準備．
晩飯．
-
あ，
北 9 東 5 の
Homac は 21 時閉店だった．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0810):
米 0.6 合．
ジャガイモ・ニンジン・タマネギ・豚肉の味噌汁．
コマイ．
- 昼 (1200):
研究室．
食パン．
クルミ．
- 晩 (2040):
米 0.8 合．
ジャガイモ・ニンジン・タマネギ・豚肉の味噌汁．
麻婆豆腐．
トマト．
2010 年 07 月 16 日 (金)
-
0720 起床．
朝飯の準備．
朝飯．
コーヒー．
0900 自宅発．
晴．
0920 研究室着．
-
あれこれののちに，
岩波本作業のつづき．
さてさて，
今日から第 4 章のモデル選択・検定章の書きなおしなのだが
……
どーにもこーにもモデル選択まわりの説明がひどすぎる．
というか，
説明にもなってないような．
-
いやべつに AIC の導出なんぞを岩波本で示す必要はないんですよ．
そんなのは，
小西・北川「情報量規準」本とか読めばいいんだから．
-
とはいえ，
逆に，
ありがちな「2 かけるパラメーター数はぺなるちーで」
といった
「なんちゃって美しきぱーしもにー実現する厳罰主義」
みたいな安易な説明なんかも通俗宗教的なかんぢで
(職業的うそつきとしては)
おもしろくないな，
と．
今までさんざん使ってきましたが．
-
そこで，
とりあえず，
「とにかく何でも具体的な数値例・図解」
な岩波本のぽりしー
(いつのまにやら?)
の旗をかかげつつ
……
-
構成変更ぷらんを検討してみる．
これまでの第 4 章を分断して，
-
新第 4 章:
旧第 4 章アタマのポアソン回帰の AIC モデル選択はそのまま，
図解・平均対数尤度 (これは従来の教科書ではなかったもの)，
旧第 4 章末尾につけてた
FAQ モデル選択
のちょー圧縮版 (つづきは web 版みろ的な)
-
新第 5 章:
旧第 4 章の尤度比検定パラメトリックブートストラップ版そのまま，
よけーな記述は削除する
といった構成を検討中．
新第 4, 5 章どちらも短い．
現在の構成は章の長さのばらつき大きいのだけど，
主題によって章の長さがかわるのはとうぜん，
などとひらきなおっている今日このごろ．
-
さてさて，
平均対数尤度
&
最大化対数尤度のバイアス補正
なる概念は下のよーな図で説明されることが多く，
これがわかりやすいような気もするのだが
……
-
えーと上の図は，
川崎能典さん
の AIC 導出の説明
スライド
の一部を利用させてもらいました．
-
R をひねくって数値例をいろいろ生成してみると
……
実際のところは，
たとえばこんな描像だ．
-
自分のど忘れ用にメモしとくと，
青曲線が平均対数尤度，
赤点が最尤推定値・最大化対数尤度，
赤破線がその標本平均，
グレイ曲線がある推定における対数尤度関数．
-
1925 研究室発．
買いもの．
2000 帰宅．
晩飯の準備．
晩飯．
[100 均の工具]
いつものごとく北 12 生協で晩飯の買いもの後に帰宅中，
「いまや北大周辺で (自転車修理に使える) 工具うってる店はないな」
(昼飯後にちょっとヨドバシにいったんだけど，
工具はまったく置いてなかった)
……
「そうだ北 12 地下鉄駅そばの 100 円均一店はけっこう工具あった!」
と寄り道してみた．
たしかにいろいろあるんだけど大口径棒レンチはなかった．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0740):
米 0.6 合．
油揚・ワカメの味噌汁．
コマイ．
- 昼 (1300):
研究室．
シリアル．
- 晩 (2200):
米 0.8 合．
ホウレンソウ．
タマネギ・ズッキーニ・ニンジン・ピーマン・トマトの煮物．
ニラ・ショウガ・卵・豚肉の炒めもの
……
らでぃっしゅぼーやで野菜がたくさんとどいたので，
野菜な晩飯に．
2010 年 07 月 17 日 (土)
-
0820 起床．
朝飯の準備．
朝飯．
コーヒー．
ちょっと怠業．
-
1110 自宅発．
晴．
まゆみさん自転車 Patagonia 号を借用して，
北 9 東 5 の Homac へ．
直径 8 mm の六角棒レンチ購入
……
結論からいうと購入する必要なかったんだけど，
まあ一本 200 円ぐらいだったし．
私が子供のころとかに比べると，
工具ってすごく安くなってるような気がする．
ついでにタオルケットも購入．
-
そのまま北大へ．
散髪しようとクラーク会館に行ったんだけど，
こんでたので撤退．
1150 研究室着．
-
昼飯．
ちょっと怠業．
-
岩波本のモデル選択説明の改訂の検討のつづき
……
-
1530 ごろ，
またクラーク会館へ．
理髪店
で散髪．
前回の散髪は 6/5
か．
まあ，
夏は散髪いんたーばるを短縮したい．
-
このあとは自転車の修理にかまけてしまった．
[クランク，外れず]
六角棒レンチで反時計まわりにぐっとまわすと，
ネジはわりとあっさり外れた．
しかしクランクは四角形断面の軸にがっちりささっていて，
ぬけそうにない．
検索してみると専用の工具が必要なようで．
目的はチェインを外すことなんだが，
後述のようにクランクぬきとりは必要でなかった．
左の写真はすでにチェイン外れたあとの図．
-
で，
Carry Me の後輪をとりはずすためには以下のようにすればよい，
とわかった．
また同じような作業するだろうから，
そのメモ．
-
後輪を固定しているネジをはずす
-
後輪ブレイキユニット (モジュール)
ごとぼこっと外す (前方のナットひとつで外れる)
-
車輪が自由に動くようになるので，
前に押しこんで
チェインをゆるめる
……
すると後輪からチェインがはずれる
-
これではまだ後輪は外れないので，
(ゆるんだチェインを利用して)
クランクギアのチェインをはずす
-
これで後輪外れる
[タイヤ分解]
これはけっこうチカラのいる作業だった
……
ホィール横に置いてる 2 本の「車輪はずし」
(といってもプラスチック棒+へらなんだが)
をテコに外していく．
チューブに一ヶ所にアナあいてたので，
紙やすりかけてゴムのりぬってパッチをはる
……
なかなかうまくいかない．
さらに，
その後またタイヤをはめなおすのもなかなか苦闘
……
時間かかる．
[なんとかぱんく修理完了]
けっきょく，
本日は路上運転試験はできなかったけど
……
じつはパンク修理すんのは初めてだったり．
-
2030 研究室発．
Carry me は直立形態のまま，
文字どおり「こわきにかかえて」
まゆみさん自転車 Patagonia 号にのる．
自転車による自転車輸送．
しかし片手運転で重量 8 kg はけっこうクタびれる
……
昔の槍騎兵とか，
けっこうたいへんだったのかもね．
2040 帰宅．
ひとりぶんなので，
いいかげんに晩飯の準備して晩飯．
-
[今日の運動]
-
土曜日なのにうんどうタイミングを逸してしまった
……
-
[今日の食卓]
- 朝 (0930):
イソップベイカリーのクロワッサンとか．
タマネギ・ズッキーニ・ニンジン・ピーマン・トマトの煮物．
- 昼 (1230):
研究室．
イソップベイカリーのライ麦パンとか．
- 晩 (2100):
米 0.8 合．
レタス・牛肉のスープ．
ホウレンソウ．
2010 年 07 月 18 日 (日)
-
0740 起床．
朝飯の準備．
朝飯．
コーヒー．
1135 自宅発．
晴．
暑い．
1150 研究室着．
-
また自転車のパンクの修理．
今朝，
朝飯調達のためパン屋に行くときに，
昨日修理した (つもりになってた)
自転車の後輪の空気がまたもれたので
……
-
すでに手順はわかっているので，
後輪解体は 10 分間ぐらいでできるようになった．
このあとは時間かかったのだけど
……
[破孔パッチ圧着中]
昨日の修理箇所は破損部分を補填するパッチの接着がアマかった，
とわかった
……
そこで，
写真のように椅子の重量で圧着させつつ，
接着剤が乾くまでの時間として 30-40 分まってみることにした．
さらにこのはりつけ作業を一度しくじったので，
またやりなおし．
-
待ち時間は研究室内怠業したり，
モデル選択本をながめたり．
[今度こそ修理完了?]
後輪を組みたてて，
はめこむ作業も 10 分ちょっとでできるようになった．
後輪ささえ後部についてるネジ，
何のためにあるのかしらんといぶかしく思っていたんだけど，
今回よーやくチェインにテンションあたえるものだと理解できた．
これで無事終了かな?
……
結果を書いておくと，
とりあえず帰宅時に乗ってみたけど問題なかった．
-
なぜかしら，
大通公園ちかくで昼飯．
-
「Robert Doweney Jr. の出演作にはお布施をする」
教のまゆみさん「アイアンマン 2」みにいくので，
そのおとも．
すでに札幌周辺ではほとんど上映終了で，
白石区にある
ディノスシネマズ札幌白石
に行ってみる．
本日は
ドニチカ切符
で札幌市地下を移動中．
[ディノスシネマズ札幌白石]
上映室 3 室なのでしねこんとよんでよいのかどうか
……
ここは 48 座席で，
蠍座より小さい．
それはいいんだけど，
アイアンマン 2 は観客数 5 でした．
だいじょうぶなのかな
……
-
アイアンマン 2，
お布施まゆみさんをして「ハナシに何の主張もなかった」
と言わしむ，
まあそういったたぐいのものなのだが
(こういう画像技術誇示的な娯楽動画に
主張を期待するほうもどうかと思うが)
……
前作 (これもお布施のおともで観た)
にも共通するのだが，
製作側の
「ほれほれ，どうせキミたちメカおたくはこういうのが見たいんだろ?
ボクらの超絶技術で可視化してやるよ」
といった「めっせーじ」を受信してしまう自分の感性に困惑させられる．
実際のところ，
こういう映画が大ひっとするのは，
かつて (中学生ぐらい?)
の私のよーにそういった演出に興奮するヒトたちの動員に
成功したためなのだろーが
……
まあ何となくおもしろく思えるのは，
こういった種類のものを「すげー」と無心に楽しめる素直なヒトたちは，
じつは理系おたくにはなれないんぢゃないかな，
いったあたり．
-
ところで，
あとで聞いたところによると，
「RDJ お布施」の理由は
「Robert Doweney がまたクスリに手を出さず，
ちゃんと映画の仕事をしてるな，よしよし」
と確認するため．
この映画もその観点からすると十分だったそうで
……
映画というのはいろいろな見かたがあるものですね．
-
晩飯．
研究室で自転車回収して，
2000 帰宅．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0900):
イソップベイカリーのクロワッサンとか．
タマネギ・ズッキーニ・ニンジン・ピーマン・トマトの煮物．
トマト・卵の炒めもの．
- 昼 (1420):
大通公園駅すぐ横の地下 3F サンドイッチ店「さえら」．
- 晩 (1820):
なぜかしら，
大通公園ちかくでハンバーグ屋．
2010 年 07 月 19 日 (月)
-
0740 起床．
朝飯の準備．
朝飯．
コーヒー．
洗濯のてつだい，
少しばかり．
-
1045 自宅発北大構内走．
完全曇天は良いのだが
……
湿度たかくしかもなかなか暑い．
昨日はフェーン現象のせいで日本海側の札幌市はまた
30 度ごえの真夏日だった．
本日もなかなか
(しかし，あとで調べたら本日の最高気温は 28 度弱).
1130 帰宅．
体重 71.6 kg．
-
昼飯．
-
近所のアオバト珈琲豆店でコーヒー豆を購入．
オリジナルブレンドと黒くローストしたのを半々に混ぜたのを作ってくれました．
どんな味になるんだろう?
そういえば，
店員さんに私のちび自転車 Carry me についていろいろと尋ねられました．
[ぐろびゅーる跡地・続]
珈琲豆調達の帰りに撮影．
スープカレー・グロビュールが 1F に入ってたヘンな建物
(三階だて? 断面みてもよくわからん)，
まるごと派手派手にとりこわし進行中．
昨年の火災以来ずっと閉店してたわけだが
……
どうなることやら?
-
1340 自宅発．
曇．
1350 研究室着．
-
雑用とかしたうけは休日にやったほうが，
なンとなくとりつきやすくて，
ですね
……
だらだらと効率わるくとりくんでみる．
[コース紹介の図]
8/1 のオープンキャンパスで使うので提出しなさい，
と私のところまで命令がおりてきたので．
研究室紹介
の図をてきとーにはりあわせる．
来年度からコース構成が少しはマトモになり，
工藤さんが植物生態に復帰される．
-
まあ，
こんなものかな．
-
1830 研究室発．
買いもの．
1900 帰宅．
晩飯の準備．
晩飯．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0810):
イソップベイカリーのバタール．
レタス・卵の炒めもの．
- 昼 (1220):
素麺．
トマト．
ナメコのダイコンおろしあえ．
- 晩 (2000):
米 0.8 合．
アサリ・スケソウダラのワイン蒸し．
レタス．
鶏レバ・ハツのショウガ煮．
2010 年 07 月 20 日 (火)
-
0755 起床．
ちょっとねぼうした．
朝飯の準備．
朝飯．
コーヒー．
黒くローストした豆の苦味はすごい．
-
昨日の計算結果を図にしてみる．
-
1000 自宅発．
曇．
今日も気温は高くならないよーで，
助かる．
1010 研究室着．
-
午前中はなにやらうだうだと失われてしまった
……
-
昼飯．
-
岩波本のモデル選択章の「図解 AIC 導出」で使う予定の図の一部．
平均対数尤度とバイアス補正
……
バイアス補正ってけっこうコワいもんだ，
と図示されてしまっていますね．
-
1830 研究室発．
買いもの．
1855 帰宅．
晩飯の準備．
晩飯．
-
2245 自宅発北大構内深夜走．
2320 帰宅．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0820):
米 0.6 合．
タマネギ・油揚の味噌汁．
鶏レバ・ハツのショウガ煮．
- 昼 (1250):
研究室．
北大生協サンドイッチ．
- 晩 (2020):
米 0.8 合．
ジャガイモ・タマネギ・ニンジン・鶏肉のスープ．
コマツナ・タマネギ・イカの炒めもの．
