ぎょーむ日誌 2007-03-(11-20)

2007 年 03 月 11 日 (日)

• 0830 起床． 朝飯． コーヒー． 洗濯 & 怠業． 昼飯． 1250 自宅発． 晴． 1305 研究室着．
• 屋久島照葉樹 22 樹種の葉寿命階層ベイズモデル WinBUGS 化， おおづめのシュート伸長の休眠・二度伸びの事前分布 (これは別データから推定したもの) の導入． なんでコレが必要かといえば照葉樹のシュートは 「見かけの齢」 とホントの齢が異なるから． 先端についてる葉っぱだからといって 「これは今年でた葉っぱだろう」 と見なすのはまちがい． そのシュートが休眠してたらもう何年も前からついてる葉っぱかもしれない， と (下の図は 昨年 1/17 のぎょーむ日誌から)．

• シュート伸長の休眠・二度伸びの分布を生成する 関数は昨日のうちに作っておいた． あるシュート (の節) において齢がどれだけススむか， はゼロ項が「二度伸び」確率となる条件つき幾何分布となる． このようなけったいな確率分布にしたがう乱数 (単なる乱数ではなく条件つき乱数つまり Gibbs sampler) は BUGS 言語の中では dcat() を使うことになる．
• さらに本来なら 10 年とか 100 年連続の休眠， ってのもあつかうべきなんだろうけど， そんなに葉つきシュートは長生きしないだろう， ということで 7 年連続あたりでうちきる． つまり事前分布の規格化による条件つき確率分布化， ということだ． これは R2WinBUGS な R コードの中で計算している．
• 1445 ばぐ取りひととーり終了． 試験運転開始． 2000 step で 838 秒か． DIC = 2123． 葉齢まわりをぜんぶ exp() で囲むと DIC = 2105． あまりかわらんのかなあ． なぜか「あの」窒素濃度の係数がプラスの偏る． こんな「ふつー」な結果はおもしろくないのう．
• 2000 step では短すぎるような気がしたので 5000 step にしてみる． 1528 計算開始． えーと 予想計算時間は 35 分間だから 1600 すぎ終了， か ……
• しかしいったん止めてもひとつ 2000 step 計算やることに． 昨年の牛原さん修論で使ったのはもうちょい別のモデル (齢に依存する部分を「個別に」 exp() した やや奇妙なモデル) で念のため試してみることに．
• お， DIC = 2086 …… しかし収束はがたがただし解釈はめんどうだから， このモデルはヤメときますか． どうせ窒素はプラスのままだし …… これはおそらく窒素濃度を「中央化」して Gibbs sampler にわたしているのが原因か? 一年前はこの手法知らなくて 「窒素まわりの収束わるい」 とかいう理解だったんだよねえ ……
• ということで， 齢の係数は exp(定数項 + 明るさ依存項) モデルでいこう． 明るさが加齢を促進するかいなか， というやつね． 窒素依存項もここにいれるべきか? いや， 「窒素」と呼んでるモノの実体はシュート先端の葉 (年齢不詳ですな) だから， 現状のままでいいか．
• それより， じつは今さら基本的なことなんだけど， 発生から (少なくとも) ひと夏を経過した葉っぱの齢は ゼロとすべきか 1 とすべきか? めんどくさいから試しに 0.5 としてみる?
• 窓の外は雪． しかし次の次の週の月曜日には松山に移動せんといかんのだよねえ ……
• 葉齢 0.5 補正モデル DIC = 2098． やはり改善されたか． 残念ながらというか窒素依存性はいままでと変わらず． 長 step 計算やる前に「窒素ぬき」モデルを試してみるか．
• 「窒素ぬき」 DIC = 2020 …… うーむ． しかも他のパラメーターの収束も改善された (じつは 昨年の計算もこういうかんぢだった)． なぜ「うーむ」なのかといえば DIC で「モデル選択」していいのかどうか， というあたりなんだよね． うーむ， ここはひとまず窒素あり・なしの両モデルの計算やっとくか．
• ということで， 窒素なしでまず 4000 step． 1515 秒． DIC = 2008． 収束だいたいよろしい． leaf.bugs.NoN.RData に保存．
• 窒素あり 4000 step． 1950 秒． DIC = 2050． 収束はさすがにマシにはなったが …… leaf.bugs.RData に保存．
• 作図プログラムをまた手直し． これが窒素なしモデルの推定結果．

• が窒素ありモデルの推定結果．

• 収束がいまいちなので， 6000 step MCMC 計算を命じるためのしくみを作る． これは窒素なし・ありそれぞれ 2270 秒 + 2913 秒 で終了した．
• 計算まち時間に「葉データぜんぶ見える table」 作ってみる． R による「作図的作表」なんだけど， めんどうなことにと言うか こういう図を作るには library(grid) (grid graphics) を駆使しなければならない．

• また laser printer CANON LBP-2810 の再生トナーカートリッジ (cyan) が不良品だった． 残量 28% でトナーでなくなる． 交換．
• 2210 研究室発． また寒波がきている． 寒い． 2225 帰宅． 明日も真冬日， という予報か． 晩飯の準備． 晩飯．
• [今日の運動]
  • 腹筋運動 30 × 3 回． 腕立ふせ 10 × 3 回． スクワット 100 回．
• [今日の食卓]
  • 朝 (0900): 米麦 0.6 合． 納豆卵焼き．
  • 昼 (1140): うどん．
  • 晩 (2440): 米麦 0.7 合． ネギ・ナメコ・豆腐の味噌汁． コマツナあえもの．

2007 年 03 月 12 日 (月)

• 0910 起床． 朝飯． コーヒー． 1000 自宅発． 雪． 昨晩また 10 cm ぐらい積雪した． 1015 研究室着．
• いろいろメイルやりとりとか． Annika さん PC が生協から配達されるのは水曜日， か …… 午後になって電話連絡あり， これは明日の午前に変更された．
• 来週の松山大会用の屋久島ポスター準備， 計算のたぐいはひととーり終わったので …… そうだ， 昨日の 6000 MCMC step 計算をここに記録しとこう．
• 窒素なしモデルの推定結果 (DIC = 2155) …… なンか DIC が悪くなって， ですね．

• 窒素ありモデルの推定結果 (DIC = 2053) …… こちらのほうは DIC 変化なし， でこちらのほうが「DIC 的に良い」 モデルになってしまった． やはり DIC はばらつきが大きくなりがちなのかもしれない (モデルやデータが不完全なためだと思うんだけど)．

• ポスター準備してるあいだに A801 室の Dell desktop 機には予備的な計算を命じておくか． このデータセットは惜しいことに 根もと側末端の節の葉痕数がいわば「欠測」状態にある． Bayesian なら個体ごとの葉痕数の事後分布だせよってコトになるんだけど …… まあそれは計算時間倍増ぐらいになるので， とりあえず末端における葉数ゼロ効果をもっとも「ひかえめに評価」 した場合に該当する「葉痕数 1」モデルで試しに計算させてみる． 6000 MCMC step．
• ポスター準備つづく (要旨)． いつものごとく LaTeX の a0poster.cls (CTAN) で． まずは 昨年のポスター をひながたにして， ですね …… タイトル両脇の写真の樹木， 母子里のダケカンバ のままですな ……

• なぜかしら dvips した Postscript file が gv (ghostscript viewer) で表示されない． ナマの gs (ghostscript) だったら表示できるんだけどな …… しょうがないので PDF に変換して xpdf で表示．
• 午後になってあれこれ調べてわかったんだけど …… 上の問題解決には psa0toa4 なる bash script で A4 化して gv 表示するのが正解， とわかった． 詳しくは 検索 してみてね …… と (次にポスター作るときにこのへん検索するであろう) 未来の自分に連絡を．
• 1400 すぎ昼飯． 昼飯後もぽすたー作業．
• 根もと側末端「欠測」に架空データいれた場合の 6000 MCMC step (試験的・予備的な) 推定計算が終了した． 窒素なし・ありモデルそれぞれ 3330 秒 + 4300 秒 …… ほーらね， 各個体一ヶ所「データ」を増やしただけでこれだけの計算時間増加． で， 推定結果はこちらのネラいどおり 「明るいところで葉っぱが一斉に脱落 (明るいところで葉の「加齢」は速い)」 というパターンを表現するモノとなった． よしよし．
• 窒素なしモデルの推定結果 (DIC = 2263)．

• 窒素ありモデルの推定結果 (DIC = 2172) …… ということで， またもやこちらのモデルのほうがよい． Ms こと タイミンタチバナ ががんばってるせいか?

• 昨日つくった葉データ 全部みえる表 みると， タイミンタチバナとかあからさまに「明るさ」-「窒素」 が交絡 (confounding) しちゃってるように見えるんだけどねえ? 葉寿命まにあっくなヒトたちは明るさ条件が同一でも 窒素が葉の加齢を促進する (この場合は正確には「加齢」というより脱落時期)， と考えているのだろうか?
• 蛇足ながら La こと バリバリノキ の葉が「明るいところほど長命」なのは「事実」とでも言いますか …… そういうふうに見えるデータだよね．
• このあたりも含めて， また末端側「欠測」データのあつかいに関しても 論文化のときに再検討せんといかんな． 今回のポスター発表では簡単なモデルの推定結果を示すことにして， 昨日の計算のうち 6000 MCMC step 「窒素あり」モデルで得られたあれこれを示して 説明することに．
• 松山移動一週間前にして (発表は来週の火曜日)， よーやくポスター準備用の MCMC 計算ぜーんぶ終了 (のハズ)． Dell 機はいつものことながらごくろーさまでした． 昨年 は自作 Metropolis-Hastings 法 MCMC 計算プログラム (in Perl) で苦闘してたわけだが (いま思うとこれはこれでスゴい計算プログラムだ)， 今回は とりっきーな BUGS プログラミングが必要ではあったけれど， やはり R2WinBUGS 使った計算はぜんぜんラクラクでした …… すでにモデリング試行錯誤は昨年にうちにすませていたから， でもあるんだけど．
• ポスター作業， ひたすら続く．
• 2000 すぎになぜか雑用事務書類提出指示メイルがくる． 2006 年度の兼業状況報告書? 私は 昨年 8 月 に「兼業先: (独) 水産総合研究センター遠洋水産研究所」 でシンポジウム講師として 8 時間働いて 10000 円いただいたことになっている …… よくわからぬまま， このゑくせるファイルを印刷出力してハンコついとく． 明日提出すれば終了．
• ばててきたので撤退． 2130 研究室発． さむー …… -4°C ぐらい． 2145 帰宅． 体重 68.2kg． 晩飯．
• [今日の運動]
  • 腹筋運動 30 × 3 回． 腕立ふせ 10 × 3 回． スクワット 100 回．
• [今日の食卓]
  • 朝 (0940): 米麦 0.5 合． ネギ・ナメコ・豆腐の味噌汁．
  • 昼 (1410): 研究室お茶部屋． 米麦 0.6 合． キムチスープ． リンゴ．
  • 晩 (2250): 米麦 0.8 合． 卵焼き． もずく． ネギ・ナメコ・豆腐の味噌汁．

2007 年 03 月 13 日 (火)

• 0800 起床． 朝飯． コーヒー． 0915 自宅発． 晴． 積雪ある冬景色 …… なんだけど来週には (当地から見ると) 季節が二ヶ月ぐらい進んでいる松山に行かねばならぬことを忘れてはいかん …… 0930 研究室着．
• 雑用メイルいろいろ． Annika さん PC も北大生協から届く．
• ポスター作りの続き． ポスターくだらなワザひとつ発見． 下の図は今つくっているポスターの一部である いつものごとき事後分布図なんだけど， これは R の density() 関数で描かせているんだけど， これをそのまま postscript() 出力すると「重い」図ファイルになる． これは density() の default が n = 512 (the number of equally spaced points at which the density is to be estimated) となっているためだ． そこで下の図では n = 30 にしてみたところファイルサイズは 1/12 (約 550KB → 45KB) になった …… そして「見た目」はほとんど変化してない． default の n = 512 がでかすぎるのである …… しょせん density() 図は「見た目」 だけのモノだし．

• まあ， いまどきそれぐらいのファイルサイズの大小でどうこういうのも ヘンなんだけど， 何となく気になるのである． ついでにいうと上の web 画像ファイルも optipng で (元ファイルサイズの) およそ半分に圧縮されているのである! …… ちょっと病的な気もするな．
• さらにどうでもよいが上の図の「基礎死亡確率」 うんぬんは 2004 年の釧路雑用でつちかわれた (そしてつい先ほどまですっかり忘れていた) LaTeX の picture 0,0 わざ (in ぎょーむ日誌記録) でスーパーインポーズされたものである． いやー， ぎょーむ日誌って役にたちますねえ ……
• 窓の外は雪． ときどき吹雪．
• ポスターの紙面がぢりぢりと図表でうまっていく． 今回は文章の少ないポスターになるかも． 昼飯．
• 夜までひたすらポスター紙面をうめる作業．

• 最後の最後で R の「日本語文字ポストスクリプト出力」でコケた． 苦闘したけど解決せず． R 内の日本語のあつかいには問題ない …… しかし postscript() 出力には問題がある． うーむ， これは R ではなく ghostscript 7.07 in Vine 3.2 の問題かしらん?
• てなことでぢたばたしてるうちにばてきってしまった． ポスター完成度は 80% ぐらい? か． へろへろと 2330 研究室発． 雪． 晩飯． 2440 帰宅．
• [今日の運動]
  • うんどう休養日
• [今日の食卓]
  • 朝 (0840): 米麦 0.7 合． キムチ．
  • 昼 (1350): 研究室お茶部屋． 米麦 0.6 合． 海藻スープ． お茶部屋キャベツつけもの．
  • 晩 (2400): 北 14 西 3 の炭火焼ばー串ごん はこの時間でも営業しており， しかもおいしいパスタをだしてくれる． オオバいりミートソーススパゲッティー．

2007 年 03 月 14 日 (水)

• 0800 起床． 朝飯． コーヒー． 0900 自宅発． 曇． 真冬かと錯覚させるような積雪． 0920 研究室着．

• ポスター作成のつづき． いくらなんでも今日中には終了するはず ……
• 中間さん R postscript() GS_OPTIONS=-DNOKANJI gv hoge.eps 中間さんご指摘によると， Vine Linux の日本語処理はふう変わりというか 「硬派」 ということだそーで ……
• 昨日のぎょーむ日誌の記述があいまいだったので， 現状を整理してみると，
  • R から日本語文字まじりの postscript (EPS) file が出力できる (file1)
  • file1 を gv のたぐいで X11 で見ることができる
  • file1 を LaTeX に \includegraphics{} して postscript 出力できる (file2)
  • file2 も gv のたぐいで見える
  • file2 を lpr すると (ここで gs filter を通過している) エラーはでないけれど，印刷出力には file1 部分だけが印刷されない
  という状況で lpr でつかう gs filter をいじるべきなのかも ……
• ところで file2 をゐんどーづ機に移して Adobe Distiller にほうりこむと file1 部分も含めて正しく PDF ファイル化された． 最終的な A0 版印刷出力のためにはこの PDF file さえあれば良いので， 現状ではそれほど問題ないと考えるべきなのかも．
• ポスター完成にむけての作業つづく． 窓の外はときどき雪．
• あとは葉っぱ脱落モデルの説明だけか …… いったん休憩して昼飯．
• 1545 よーやくのことで， 屋久島葉寿命ポスターひととーりできた …… A4 縮小印刷出力ながめてちょっと修正してから， お茶部屋でひと休み． 院生の皆さんから 「ごちゃごちゃしたポスター」 との講評をうける．
• また内容など手なおし． 共著者の牛原さん・甲山さんに見てもらうべく， 今回の ポスター発表ペイジ 作る．

• 事務からまた何かみょーなことを …… いや， 今までどーりか? まあ， ここにはっておこう． ついでに研究室の雑用 Wiki にも．

【航空機利用時の搭乗半券等徴取に係る取扱い】 2007-03-14

１．搭乗半券等の徴取について
　　　航空機を利用した場合の搭乗半券を徴取することは，『搭乗の証明』『航空券の運賃』の
　　事実が確認でき，旅費の不正請求及び受給等の防止につながります。
　　　また，法人として社会に対する説明責任を果たすうえで必要なことでもあることから，
　　航空券等の領収書の他に搭乗半券の徴取を義務付けるものです。
　　　なお，航空券の運賃は，搭乗半券に印字されている記号により航空運賃種別を判別でき，
　　航空運賃が分かる仕組みになっております。

２．実施開始時期　　平成１９年４月１日出発分の旅行から適用します。

３．提出方法について
　　　本学が認める領収書（原則，氏名・金額・搭乗年月日・搭乗便名が記載されたもの）と
　　合わせて，下記４の搭乗半券等の原本を提出して下さい。なお，概算払の旅行については，
　　旅行完了報告と同時に，別途提出して下さい。これが不可能な場合は，当該旅行に係る航空賃
　　相当額は支給されません。
　　　なお，精算払の領収書等に搭乗半券等が添付されていない場合は，回収いたしませんので，
　　必ず領収書と搭乗半券を一緒に提出してください。

４．搭乗半券について
　　　空港の搭乗口で受け取る搭乗券の半券で，氏名・行先・搭乗日・搭乗便名の記載があります。
　　　旅費システムを利用（電話・ファックスを含む）し，航空券等を手配した場合は，搭乗確認が
　　できるため搭乗半券の提出は不要です。
　　　なお，以下の場合は搭乗半券と同様に扱うこととします。
　　　・ＪＡＬの「ＩＣチェックインサービス」を利用した場合の「ＩＣチェックインお客様控え」
　　　・ＡＮＡの「スキップサービス」を利用した場合の「搭乗券」
　　  ・搭乗証明書

事務部会計係

• ばててきたので「今日は洗濯しなくちゃ」という口実で撤退． 1830 研究室発． 吹雪． お， 「ねこや」の横に新しい中華料理屋ができつつある． このあたりはスープカレー屋の高密度地帯として有名であるが， じつはなぜか中華料理屋もけっこう高密度だ． 1845 帰宅． 体重 68.6kg． 洗濯． 晩飯の準備． 晩飯．
• [今日の運動]
  • ばてばてうんどう休養日
• [今日の食卓]
  • 朝 (0830): 米麦 0.5 合． ワカメスープ．
  • 昼 (1340): 研究室お茶部屋． 米麦 0.6 合． こんびにで買ったフリーズドライの ユッケジャン．
  • 晩 (2120): 米麦 0.8 合． ニラ・シイタケの卵とじ．

2007 年 03 月 15 日 (木)

• 0900 起床． うう， 9 時間も寝てしまった …… しかしまだアタマがややへろへろしている． コーヒー． 朝飯． 1020 自宅発． 晴． 1035 研究室着．
• アタマがへろへろしてる理由: (昨日の作業である) LaTeX ポスター作りはある意味ラクでいいんだけど， これってアタマの中でレイアウトの座標計算やりながら 作成することになるんで， このあたりが意外と (自分でも気づかぬうちに) 負担になっているんだろうな ……
• 窓の外はときどき吹雪． 三月になってもこんなに降るとはなかなかのモノだな．
• 午前中は何だかメイル書いたり， 屋久島 web page や自由集会サイトまわりでふらふらしていたり …… いまやありがたいことにネット上をさまようだけで 愛媛大学→「大街道」周辺の地理がアタマに入ったりするんだけど， 自由集会後のうちあげ会場はいつものことながら何とも決めがたい． 人数が把握できないから， という理由もあるんだけど．
• 気分転換にちょっと昼飯調達にでる． 研究室にもどって昼飯．
• 1400 から 2 時間弱ほど石沢さんと豊平川アミノ酸の 階層ベイズモデルのハナシ． じつはこれは (来年度あれこれやらねばならぬ) 群集生態学ハナシのモデリングとも関連ありそうだ …… といった統計モデリングおたくな話はともかく， このアミノ酸うんぬんとサケ母川回帰のハナシがどう関連するか いまいちよくわからない． おそらく石狩川・豊平川溯上してるサケが 分岐点あたり で「えーと，豊平川は右か」とか判断するのに使うんだろう (??)けど， 札幌市豊平川さけ科学館 の豊平川のサケ標識放流調査ってのを見ても 「どれぐらい母川を識別してるか」 はわからない． 「まちがって左の石狩川とかに行ってしまいました」 データはどこにもないんで．
• 粕谷さんと 自由集会 メイルやりとり …… 粕谷さんの一般講演 (ポスター) は「主成分分析と検定」なるスゴそうなタイトルとわかったので， ネット上の要旨 を読んでしまった …… PCA と統計学的検定を連関させんのは意外と難しいような気もする． そりゃ， bootstrap とか使えばどういう場合でも何らかの P 値みたいなものは計算できるでしょうけどね ……
• 経験ベイズ法の例題としては， 縮小推定量 (Stein 推定量) のハナシとかわかりやすいかも …… たとえば， 「100 人のプロ野球選手それぞれ 5 打席ぶんのデータから 各打者の『真の打率』 (?) を推定しなさい」， みたいな． 打者ごとに割りざん推定してしまうと 「0 割打者」とか「10 割打者」とか ぞろぞろ出るんだけど， 階層ベイズモデル化 (打者 100 人の「個体差」考慮しつつ， 全員のデータから「打率の事前分布」を推定しつつ) → 経験ベイズ法で推定計算 (これは最尤推定) するともうちょっとマシな推定結果になる (5 打数 5 安打でも「10 割打者」とはならない)， といったハナシですね．
• 中間さん からのメイルを拝読してるうちに， 今さらながら気づいたんだけど …… 2 年ちょい前 のぎょーむ日誌で R の expression() で日本語文字があつかえない， とかぶちぶち文句を言ってたんだけど， なかまさんたちの 尽力 によってその後の version の R では ちゃんとあつかえるようになっている， と今ごろになって気づいた! device は pdf() や postscript() でも問題ナシ． R 開発者の皆さん， ありがとう ……

  x <- rnorm(100)
  y <- rnorm(100)
  plot(x, y, xlab = expression(Delta * beta["control"]^italic(Z) * " の値"))
  

• 自由集会準備は具体的にはあまりススまぬまま撤退． 1910 研究室発． 1930 帰宅． 体重 68.6kg． 晩飯の準備． 晩飯．
• [今日の運動]
  • 腹筋運動 30 × 3 回． 腕立ふせ 10 × 3 回． スクワット 100 回．
• [今日の食卓]
  • 朝 (0930): 米麦 0.5 合． ワカメスープ．
  • 昼 (1320): 研究室． 北大生協にぎり飯．
  • 晩 (2220): 米麦 0.7 合． タマネギ・ブナシメジ・ホタテのカレー． ダイコンのゴマサラダ．

2007 年 03 月 16 日 (金)

• 0820 起床． なンかあまりよく眠れなかった． コーヒー． 朝飯． 0930 自宅発． 晴． 0945 研究室着．
• 松山大会用のポスターもうちょい修正して， 「これで印刷してみます」版を 屋久島葉ペイジ にアップロード． つくば (!) にいる牛原さんと ブラジル (!!) にいる甲山さんにメイルで連絡． 地環研から北に 100 m ほどいったところにある 情報基盤センター北館 に印刷にいく．

[基盤センター北館]

玄関． 玄関右側 (東側) が最後に利用報告書を提出する 受付 (事務室) である． プリンター利用のためには アカウント取得が必要． われわれは甲山さんのアカウント使わしてもらってるわけだが ……

[2F 西側のディスプレイ端末室]

というよくわからぬところに A0 プリンター二台が設置されている． 端末はゐんどーづ機と iMac． インストールされてる ソフトウェアは使いそうなものひととーり． Adobe Illustrator なんかも．

[HP の大判インクジェットプリンター]

大判プリンタ ペイジも参照． この時期は混雑している． 私も 20 分ほど順番まち． いざ印刷にかかると， ここで LaTeX の a0poster.cls によるサイズぎめを私がよく理解していない弊害がでて， ですね …… 私のポスターの場合「ぬり」が少なかったせいか， 15 分ぐらいで印刷終了．

[印刷終了]

終了すると「インク乾燥中」とやらでその状態で 5 分間ほど作業中止． 5 分間経過するとカッターがロール紙を切断する (ただしこの段階でもまだインクは完全には乾いていない)． 左右などの余白は同じ部屋にある大型カッターを使って切断できる．

[利用者受付]

1F にある事務室に利用報告書 (これはプリンター横の端末のそばにある) に使用枚数など書いて提出． A0 光沢紙だと一枚 900 円． お金はあとからアカウント保持者に請求される．

[屋久島ポスター]

用紙サイズ設定その他にいろいろ課題を残してしまった …… いつもは年度末予算消化もかねて印刷屋さんにおねがいしてたので， むこうでかなり適切なサイズに調節してくれてたわけだが． すんづまりポスター， ということで「見るヒトの腰にやさしい」 「後ろにならんで説明きいてるヒトたちにやさしい」 けれど， 文字が小さく行間がつまっているので 「目にやさしくない」 ポスターですなぁ．

• 用紙設定の教訓としては， dvips などで出力される preview はそれなりに正しく， a0poster.cls の挙動はあまり斟酌せずにタテ・ヨコの range を決めればよい． 上記 A0 プリンターのドライヴァーが最終的な縮小・拡大を かなり適切にやってくれる． また， このドライヴァーまわりをうまく設定できれば， A0 よりも 10cm ほどハバの広い印刷が可能になるだろう．
• ただし …… 私にとってはまあどうでもいいんだけど …… RGB → CMYK 計算はちょっと偏りがあるのか， 青なんかがかなり紫っぽくなる．
• 研究室にもどって昼飯． もう 1400 か．
• メイルかきあれこれ． jeconet に 自由集会 の案内を投稿してみる …… ああ， 準備ができてない．
• しかしかんぢんの自由集会準備そのものはあまり進捗せず． 2015 研究室発． 2030 帰宅． 晩飯．
• 一週間後にはもう松山から札幌に戻ってるのか ……
• [今日の運動]
  • うんどう休養日
• [今日の食卓]
  • 朝 (0840): 米麦 0.6 合． タマネギ・ブナシメジ・ホタテのカレー． ダイコンのゴマサラダ．
  • 昼 (1400): 研究室お茶部屋． 米麦 0.5 合． ワカメスープ．
  • 晩 (2130): 米麦 0.7 合． タマネギ・ブナシメジ・ホタテのカレー．

2007 年 03 月 17 日 (土)

• 0810 起床． 今朝も冷えこんでて雪． しかし明後日にはサクラとか咲きそうな ところ に行くわけだよなぁ …… コーヒー． 朝飯． 洗濯 …… で， 怠業するわけにもいかず， 自由集会 準備にとりくもうとする努力のつづき． 1050 自宅発． 曇． 1110 クラーク会館地下 理髪店 着． 1130 散髪．
• なンで自由集会の準備もできてないのに髮なんか切ってるんだよう， と指摘されそうだけど …… まあ， こういう状況ではいったん計算機の前から離れてあれこれ考えるのが有効， ということで． じっさい， ここですごして一時間ちょいのあいだに 「ベイズ推定と最尤推定のちがい」 に関する説明方法いろいろと整理できたし．
• そして「学会後」になれば時間によゆーがあるのかと言えば， 必ずしもそうではない． たとえば
  「久保さん，○○の計算もうできましたか?」
  「いやー，松山大会の準備とかがいそがしくて，ですね」
  「いつごろできるんですか?」
  「松山から帰ったらすぐにとりかかりますよ」
  というやりとりが， 最近ひんぱんになってしまって …… つまり生態学会大会を「時間かせぎ」の口実に使っている罰として， 「学会後」と呼ばれる期間にその報いを受けねばならなくなりつつある， ということで．
• 1220 クラーク会館地下の散髪屋でる． 2625 円． 雪．

[雪の中のクラーク先生]

たまには北大的な写真でも．

• 研究室お茶部屋で昼飯たべてから， また自由集会準備のつづき．

• うう， まだ終わらん． しかしばててきたので 2120 研究室発． 2140 帰宅． 体重 68.0kg． 準備やせ? 晩飯．
• [今日の運動]
  • 腹筋運動 30 × 3 回． 腕立ふせ 10 × 3 回． スクワット 100 回．
• [今日の食卓]
  • 朝 (0910): 米麦 0.6 合． 納豆． ネギ・豆腐・卵の炒めもの．
  • 昼 (1330): 研究室お茶部屋． 米麦 0.6 合． ワカメスープ．
  • 晩 (2220): 米麦 0.7 合． 納豆． ネギ・豆腐・卵の炒めもの．

2007 年 03 月 18 日 (日)

• 0750 起床． コーヒー． 朝飯． 自宅で ベイズ自由集会 の準備作業のつづき．

• 1755 自宅発． 雪． 1810 研究室着． また作業のつづき．

• ああ， 粕谷さんはすでに松山で ぐるめ な報告を …… うう， 当方は松山における自由集会うちあげ の店をまだ探せていないのに ……
• 自由集会準備， さすがにほぼでそろってきたような気もする …… けど， まだまだみたいな気も． 2330 研究室発． 晩飯． 2450 帰宅．
• [今日の運動]
  • うんどう休養日
• [今日の食卓]
  • 朝 (0850): 米麦 0.5 合． ネギ・卵の炒めもの．
  • 昼 (1330): もりもとのサンドイッチ．
  • 晩 (2410): この時間， 北大周辺で飯がくえる場所は少なく …… また北 14 の串ごんでパスタ．

2007 年 03 月 19 日 (月)

• 0810 起床． コーヒー． 朝飯． 今日から松山への出張ということになる． しかし飛行機の出発は 1420 なんで， 時間的によゆーある …… ということで洗濯しつつ旅じたく．
• 0950 自宅発． 曇． 1005 研究室着． しつこく ベイズ自由集会 準備のつづき．
• 正午すぎてしまったので 自由集会 ペイジに「とりあえず版」をアップロード． あわただしく昼飯．
• 1250 研究室発． 1255 JR 札幌駅着． 自動チェックイン機で航空券うけとり． 1310 快速エアポートにのって同発． 1040 円． 1346 新千歳空港着． 荷物あずける． ターミナルからバスで移動． 1400 ANA866 に搭乗． 搭乗率 7 割ぐらい? 道内の生態学関係者多数．

[Airbus A320-200]

北海道発松山直行便． 私はかなり後方の席えらんだ． 北に離陸 → 時計まわりの旋回で進路南に． 珍しく窓際えらんだら， 新千歳空港 がよく見えたので 撮影してしまった． 奧の滑走路が千歳基地．

[瀬戸内上空]

二時間以上とんでよーやく瀬戸内に． この多島海を飛行機からながめるのはいつも楽しい (といっても 10 年ぶりぐらいなんだが)． 飛行場が見えたので 「松山空港ってあんな山奥だっけ?」 と思ったんだけど， これは広島空港だった． 手前の 河口都市は三原市． ここから南方に左旋回して 海ぞいの松山空港へ．

• 1640 松山に着陸． 荷物回収． 伊予鉄空港リムジンバスにのって 1658 同発． 駅まで 300 円． お， 松山城の天守閣は小さいけど山城だから遠くからよく見えるもんだな …… 1715 JR 松山駅着． うーむ， 意外と暖かくない?
• 1720 市電の循環線 (これって単線だから時計まわりにしかまわれない?) に乗って駅前発． 1735 清水町で下車． 運賃は 150 円で固定されてる． 1740 松山城北側にある宿， ホテル泰平 着．
• チェックインして部屋へ． すぐにネットにつないで， 明日のベイズ自由集会後の「うちあげ」会場をさがす． 大学からの近さと大人数対処能力たかそうなところを選んで， 電話 …… わりとあっさりと居酒屋「庄屋」 に決定．
• 自由集会うちあげは人数不確定なんでいつも店に嫌がられるんだけど， 今回は伊予気質とでもいうか「ええよ，ええよ」 でハナシが進んでしまった …… 自由集会うちあげペイジ を更新 & アップロード． すでに松山 ぐるめ 生活を楽しんでおられる粕谷さんにメイルで連絡． 一件落着．
• この宿は温泉地から離れてるんだけど， よくわからん方法で温泉を「引きこんで」いるそうで …… ともかく私の部屋の階に大浴場がある． ということで， お手軽にとりあえず「温泉」． 温泉はともかく， びぢねすホテルに大浴場ってのはなかなかおトク感ありますなぁ ……
• 晩飯． 晩飯はすぐ下のレストランで． ほてるのれすとらん， と言ってもこれまたお手軽なところで， 1050 円しはらってセルフサーヴィス (ただし時間は 2000 まで)． 料理は多種多様 (なぜか肉料理が多かったような)． 本日は「ビールなど酒類が飲み放題の日 (月に 8 回)」 らしいんだけど， エタノール服用しない私にはご利益ないよな ……
• しつこく， 自由集会の準備のつづき， とか …… ちょっと変えては 自由集会 サイトにアップロード …… あ， もう 2600 すぎか ……

• [今日の運動]
  • 腹筋運動 30 × 3 回． 腕立ふせ 10 × 3 回． スクワット 100 回．
• [今日の食卓]
  • 朝 (0830): イソップパンの巨大クロワッサン． ネギ・卵の炒めもの …… 冷蔵庫内， ほぼくりあー．
  • 昼 (1240): 研究室お茶部屋． 米麦 0.5 合． ワカメスープ．
  • 晩 (1930): ホテルのセルフサーヴィスレストラン． 1050 円． けっこう良いかんぢの料理でした．

2007 年 03 月 20 日 (火)

• 0715 起床． 松山滞在 2 日め． 4.5 時間ぐらいねた． なンかあまり食欲ないけど朝飯． コーヒー． しつこく 自由集会 の準備．
• 0920 宿発． 徒歩． 平和通を東へ． 0930 大会会場の愛媛大学着． うけつけ． すでに混雑時間帯はとっくにすぎている．
• ポスター会場でポスターはりだし …… そのまま 1610 までポスター説明が続いた．
• ようやく客足がとだえたので， ほかの部屋のポスター見物． 植物生理の部屋 (今回のポスター会場は細かく分断されている) で種子田さんから実生の資源分配問題のモデル & 実測のハナシなど拝聴する．
• 自由集会準備のつづきを， と休憩室にいったらなぜか九大の佐々木さんに Bayes な説明を． あいかわらずスゴい記憶力に驚く (ずっと以前に読んだはずの James-Stein 推定量 の解説を覚えておられた)．
• 1700 自由集会会場である A 会場へ．

[A 会場]

でかい． しかし最終的にはほぼ満席になってしまった ……

• 冒頭のあいさつの準備をしてなかったので， あわてて PDF file を生成．
• 1800 データ解析で出会う統計的問題 -- ベイズ統計学の考えかた・使いかた 開始．

[べいづ自由集会]

粕谷さんの発表． 参加者である森本さん撮影のもっと良い写真は 自由集会ペイジ に．

• 2000 ぶじに終了． 参加者の皆さんありがとうございました． 2030 から 自由集会うちあげ．

[うちあげ]

居酒屋「庄屋」． 参加者 14 名． 楽しく 2 時間ほど．

• 2320 宿着．
• [今日の運動]
  • うんどう休養日
• [今日の食卓]
  • 朝 (0740): 宿のセルフサーヴィス朝飯．
  • 昼 : 食ってない …… 夕方， 飲食禁止の自由集会会場でクラッカーかじる．
  • 晩 (2030): 自由集会うちあげ． すきやきとかすごくひさしぶりに食べた． 「坊っちゃん」 の地で「牛鍋」というのも良いかも (べつに当方は芋ぜめ・豆腐ぜめだったわけではないが)．