ぎょーむ日誌 2009-03-(11-20)
2009 年 03 月 11 日 (水)
-
0720 起床．
朝飯．
コーヒー．
-
朝から不幸のメイル．
Behavioral Ecology and Sociobiology
に出してたアリ論文が reject になった．
うーむ，
また「脱」アリがとおのいてしまった．
-
査読者二名が言ってることはだいたい同じで，
要約すると
-
よけーな discussion を大幅に削除して，
もっと統計モデルにハナシを限定しろ
-
よけーな図を削れ
-
作文がへた
これって二回目の査読なんだから，
このあたりは最初の査読のときに指摘してくれー
……
と言ってもムダでしょうな．
-
査読者の一人は，
この論文が二重の目的を追及しようとして中途半端になってる，
とも指摘している．
つまり階層ベイズモデルを使って動物の行動解析が改善される，
というハナシと，
エゾアカヤマアリの集団構造みたいなハナシ．
ありていに言って，
たしかに「すーぱーころにー」うんぬん議論を全部削除して
行動解析に集中してよいのであれば，
当方だって作文苦労はずいぶん軽減されましたよ．
-
ともあれ，
あちこち手なおしして別のぢゃーなるに再投稿，
となりそうだ．
その際には削除できる図・議論は全部削除する，
と．
-
行動解析のベイズモデルも一点だけ改善したほうがよさそうだ．
今回の reviewer と「蟻」解説を書いたときに恒岡さんに指摘されたのだが，
antennation する → 攻撃成功確率が下がる，
といった現象の説明として，
-
antennation なんて行動に時間をかけているから，
そのすきに侵入アリが逃げだして攻撃失敗におわる
(論文で書いてること)
-
そもそも「敵っぽくない侵入アリ」がいて，
そいつは念入りにチェックされるけど，
「やっぱり敵でなかったか」
ということで無事に釈放されるので攻撃されない
(指摘された別の説明)
のふたつがありうる．
私としてはべつにどっちの説明でもさしつかえないんだけど，
統計モデルにちょっとテを加えることで，
ハナシがしやすくなるかも，
という気になってきた．
-
これは「個体差」のいれかたに関係する．
今まで，
どうせ一個体あたり 5 反復なんだし「個体差」
なんかわかるわけないぢゃん
……
などと考えて antennation と攻撃における「個体差」
は同じ (つまりアヤしい個体は antennation・攻撃されやすい)
といったモデルになっていた．
-
改良はちょっとめんどくさいんだけど，
「個体差」を antennation に関するものと攻撃に関するものに分割，
そしてぽいんとはココなのだが，
これらの事前分布を二変量正規分布とする．
で，
この二変量正規分布の分散共分散行列の事後分布から
相関係数の事後分布を計算するわけだが，
-
相関係数がマイナス:
「敵っぽい侵入アリ」と
「敵っぽくない侵入アリ」がいる説がマシかも，
となる
……
つまり「個体の差」に依存して
「antennation せずに攻撃」
「antennation して釈放」
していますよ，
と
-
それ以外:
攻撃確率の説明変数である「事前の antennation」
という説明がマシでしょう
上で書いたように一個体 5 反復なので，
相関係数はゼロからずれない可能性が高いのだが
(上のハナシが支持されにくい)，
まあいちおうこういう比較ができるようにはなっている，
と．
-
しかし論文がこれ以上長くなるのも問題だな
(すでに 8000 word 超)．
次の投稿前の作業としては，
-
まず可能なかぎり論文を短くする
(投稿可能そうなぢゃーなるの投稿規定を見てると，
これが必要に思われる)
-
モデル改良による再長文化が可能か検討する
といった手順でやるべきだろう．
モデル改良できなかった場合には，
「説明いろいろあります」
式の議論でごまかす展開になるだろう
……
つまりこのあたりは「どーでもいーや (単なる bias 補正でーす)」
といったあつかいにして，
ハナシを実験アリ拉致サイトに依存して攻撃確率が段階的に変化します，
といったあたりに集中する，
と．
これはこれで安易な改訂だな．
-
0920 自宅発．
雪．
まだまだ冬っぽい日にもなる．
0935 研究室着．
-
ちょっとアリ論文など検索したりして，
アリ問題を検討してしまった．
「納得」が必要とされているんでしょうな．
-
小林君の予備審の時間だけど，
まあ何ともよゆーのない状況なので，
ぱす．
まあ，
議論は公聴会でやればいいでしょう．
-
専攻事務から連絡．
あらら，
甲山さんが次期 (生物圏) 専攻長になってしまったか．
どうなることやら．
これ以上，
雑用が増えてもだいじょうぶなんだろうか?
-
とりあえずアリ論文共著者にメイル連絡．
再投稿は早くても 4 月以降．
上のようにいろいろ思惑があるわけだが，
すでに上に書いてしまってアタマの中で整理されたので，
メイルはかなり簡潔になった
……
「ぎょーむ日誌」の効能というやつですか．
-
またべつのアリ問題に関するべつのメイル書き．
どうなることやら
……
-
そしてデータ解析こんさるメイルかき．
昼飯．
-
なんだか進捗ののろくさいウミガメポスター作業のつづき．
-
とりあえず主要部はできた．
紙面の 7 割ぐらいがうまったかんじ．
時刻は 1530．
-
残り 3 割の部分は
……
これから計算が必要だ．
どうなることやら．
まあ，
すでに MCMC 計算で得られた
事後分布標本をひねくるだけですけどね．
-
1997-2006 年の平均対数増殖速度の予測分布，
あっさり計算・作図できた．
こりゃーおもしろいね
……
しかしこれにもとづく「将来予測」なるものは，
とくにおもしろいわけでもない，
と気づいた (動態モデルが単純だから)
……
-
1700 から農学部 (というかフィールド科学センター 2F 会議室)
で ベイズカエル本
輪読会，
最終回．
FSC の森君が第 11 章 Subjective prior で，
著者 McCarthy がやたらとススめる主観事前分布のススめ，
といった内容．
うーむ．
ようするにアンケートとかの結果から事前分布をでっちあげる
その作法に検討の余地いろいろありだと思うんですけどね．
-
そのあと私が一時間ほど，
ベイズカエル本ではすごくカルい
あつかいだった階層ベイズモデルの補講
……
といいますか，
先日の統数研の
公開講座
の流用でありますが．
-
そのあとはその場で輪読会うちあげパーティー．
[うちあげ]
豊田さん・宮田さんがもちこんでくれた料理いろいろ．
ごちそうさまでした．
飯島君ほか院生のみなさんの尽力で，
楽しい輪読会 & うちあげとなりました．
この勉強会，
いろいろ勉強になりましたよ．
著者 McCarthy みたいな安易な
「何でも主観事前分布」
派が増えていきそうな予感もするのだが，
まあ階層な事前分布を普及したいところです．
-
2240 研究室発．
真冬のような吹雪．
2300 帰宅．
-
明日は朝から試験監督だ．
二次試験後期日程．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0750):
米麦 0.8 合．
ネギ卵炒飯．
- 昼 (1300):
研究室．
セイコーマートサンドイッチ．
- 晩 (1920):
うちあげのごちそう．
2009 年 03 月 12 日 (木)
-
0700 起床．
朝飯．
コーヒー．
NHK ラジオのニュース解説によると，
来年度卒業生の就職求人状況はいわゆる「就職氷河期」
ほど悪くはならないはず，
とのことで
……
このあたりが大学院受験者数 (≅
入学者数?)
に影響してる気がする．
0750 自宅発．
雪．
0805 薬学部着．
本日午前中はここで試験監督 (薬学部二次試験後期日程)．
-
まあ，
室外監督員なのでお気楽ともうしますか．
試験前は玄関で受験生の受験票チェック．
試験中は廊下の椅子に座ってるだけ．
一度，受験生をトイレまで案内したけど．
しかも途中で交代するので，
あとは監督控室でのんびり．
[薬学部のセミナー室]
監督控室として使われていた．
棚には漢方薬のビンがずらりと．
-
ウミガメモデル，
一ヶ所修正すべき点が見つかった．
分散パラメーターのいれすぎ．
一個けずる必要がある．
昨日も書いたように，
この random walk 成長なウミガメモデルによる
「将来予測」
ってにはじつにおもしろくないものなんだけど
……
やはり統計モデルは「それを使ってどのように予測できるか」
をよく考える必要があるな．
-
で，
監督控室で再計算などしてみたのだが
……
これってけっこう抜本的な改訂になっている，
といまさらながら気づいた．
-
この改訂版モデルはダメ，
というのも観測データにひきずられすぎてめちゃくちゃになる，
しかも収束が悪い
……
といったことがすでに昨年 1 月のぎょーむ日誌
(2008-01-09)
に書かれており，
その翌日
(2008-01-10)
になって
「水棲動物個体群の時系列データ解析がこんなのでホントによいのか，
まだまだよく検討すべき点はあるけど，
修士論文対策としてはとりあえずこんなモンでしょう」
などと言いわけしつつ今まで使ってきたモデル
(ポスターの準備に使ってたやつ)
を採用している．
-
仮に，
この試験監督控室で試してみた「観測データにひきずられすぎ」
なモデルを「差分 (勾配) 平滑化モデル」
(統数研ふうモデル，とでも申しますか)，
ここまでポスター準備に使ってきた
(「分散パラメーター多すぎ」と言ってた)
動態モデルを
「確率論的密度変化モデル」
とよぶことにしよう．
-
推定モデルとしては，
今までポスター準備に使ってきた
確率論的密度変化モデルのほうが断然に収束がよい．
問題はこのモデルの解釈と動態予測のときの計算のやりかた，
だ．
-
解釈というか他人に説明するときの注意としては，
今までばくぜんと
(対数スケイルでの)
「成長率」とよんでいたものを，
(対数スケイルでの)
「平均成長率」と呼称するべきである，
と．
-
おっと試験終了時刻からまた全員出動か．
答案用紙の回収 → 受験生かえってヨシ，
となるまで 30 分ぐらいかかる．
-
その間もいろいろ検討．
統数研ふうな差分平滑化モデルを
(観測データ追従にならぬよう)
うまくやるためには，
データ数より
「点」の数を多くする
必要がある，
ということだろうな．
たぶん．
そしてウミガメの場合，
「点」を増やしてよい，
と説明する理屈がめんどうだ．
-
ということで，
やはり確率論的密度変化モデルでやってみるか．
ただし用語に関してはいろいろ修正する必要がある．
-
受験会場の忘れものチェックして，
解散．
1220 研究室着．
昼飯．
-
研究室雑用．
そしてしつこく継続される The Lichenologist
編集長殿からの質問メイルへの回答かき．
これは短文なのだけど，
かなり消耗感あり．
-
ウミガメポスター作業のつづき．
いつもながら作業が遅い．
-
日没後にひととーりできた，
ように思える．
あとはウミガメ写真を「安全」なモノにとりかえたいのだが
……
-
ウミガメ写真，
Wikimedia commons 上で
``loggerhead''
などで検索してみる
……
こちらの 水中写真
と宮崎県の 上陸砂浜
を使わせてもらうことにしよう．
-
で，
このふたつは Creative Commons ライセンスなので
……
といったあたりを調べてみたんだけど，
最低限
``(c) 名前, 発表年''
とすればよいのかな?
あとはポスター上の
(少なくとも)
その写真の部分は CC ライセンスにしておく，
と．
うーむ，
出典を Wikimedia Commons と書いておけば，
もとライセンスが調べられるから OK，
とはならないかな?
-
まあ，
とにかくできたことにしよう．
こうやって見ると，
また文字の多いポスターになってしまった．
昨年よりは全体の文字サイズをひとまわり大きくしているんだが
……
-
ポスター PDF ファイルをアップロード，
共著者のお二人にメイル送信して，
2120 研究室発．
買いもの．
2150 帰宅
……
「ぽすたー，ひとまづ終了」
と脱力 & アタマがぽけーっとしていて
ThinkPad を学校におきわすれてしまった．
-
晩飯の準備．
晩飯．
ネットから途絶された夜，
というのも平和でよいかも．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0720):
シリアル．
- 昼 (1300):
研究室お茶部屋．
[試験監督弁当]
試験は午前中で終わったのだけど，
なぜかしら左のような弁当が配給されたので，
研究室にもちかえって食べた．
- 晩 (2240):
米麦 0.6 合．
キャベツ・ニンジン・ソーセイジの蒸し煮．
ヨーグルト．
2009 年 03 月 13 日 (金)
-
昼飯．
その後，
ひさびさにお茶部屋雑談にとらっぷされた．
まあ，
ポスターとかひととーりできてしまうと，
のんびりしてしまうもので
……
生態学会大会 (盛岡)
に出発するのは来週火曜日の朝だしね．
-
ウミガメポスターの細かい修正あれこれ
……
まあ，
こんなもんかな．
一年以上放置していたウミガメ上陸数の統計モデル，
そのポスター準備にはやはり
一週間以上
をついやしてしまった．
-
お茶部屋．
院生たち数名が年度末がらみの雑用にふりまわされている．
院生を雑用戦線に投入する，
というのはまあ学徒動員の一種で，
その場はしのげても後遺症があるように思える．
-
年度末あまり予算でインドからの PD ギータさん用の
英語版ソフトウェアが手配できそうなのは助かるわけだが
……
私などはあまりカネのかからない研究者なんだけど，
それにしても，
ともかく今年度はてもとにそのわずかのカネさえない一年だった．
-
東南アジアから帰ってきた藤沼君のイノシシ &
外来植物のポスターをお茶部屋で一時間ばかり拝見する．
-
GLM 自由集会
(なぜか GLMM 自由集会だとかんちがいしているヒトが多い)
の準備，
よくよく見るとまだまだ手ぬきだと気づいた．
これでは 20 分 + 40 分とやや変則的に分割されたハナシを
のりきることはできん．
-
しかしまだ明日・明後日・しあさってとよゆーあるから，
本日は撤退 (なぜかしらバテぎみなことだし)
……
と思ったのだけど，
この時刻になってなぜか作図質問メイルがくる．
ちょっとややこしい．
一時間ばかり試行錯誤．
返信．
-
1910 研究室発．
1930 帰宅．
晩飯．
-
[今日の運動]
-
晩飯後いきなりばたっと寝てしまった
……
夜中にめざめて
「うん? 今日はなん曜日? 昼? 夜?」
とあたりを見回す
……
-
[今日の食卓]
- 朝 (0730):
米麦 0.6 合．
キャベツ・ニンジン・ソーセイジの蒸し煮．
- 昼 (1300):
研究室お茶部屋．
米麦 0.6 合．
キャベツ・ニンジン・ソーセイジの蒸し煮．
- 晩 (2000):
米麦 0.8 合．
キャベツ・ニンジン・ソーセイジの蒸し煮．
ヨーグルト．
2009 年 03 月 14 日 (土)
-
0640 起床．
先日の試験監督いらい，
ちょっと早めに目がさめる．
次の火曜日の盛岡行きはかなりはやい出発になるので，
これで問題ないんだけど．
えーと昨晩は夜中に目がさめて
……
まあ，
それでも合計 6 時間ちょいは寝ているはず．
朝飯．
コーヒー．
とはいえ，
やはりだるいのでごろごろしながら怠業．
ほんのすこしばかり，
だけどちょっと風邪ぎみなのかも．
-
今日はあわてて仕事にとりかからないことにして
……
すごくひさしぶりに部屋のそうじ，
とか．
-
うんどう．
30 分でへろへろ．
というかぐらぐらする．
-
昼飯．
やはり体調が「ヘン」という気がしたので，
体温はかってみたら
37.3°C
だった．
うーむ．
とりあえず
先月のヒドい風邪
のときに処方してもらった，
解熱剤など飲んでみる．
-
今日は学校にいくのはあきらめて
(少なくともポスター準備が早めに終わっててよかった!)，
横になってみる．
やはりバテぎみなのか，
2 時間ばかりくーくー寝てしまった．
-
で，
体温はかってみたら
36.7°C
までおちてた．
よかった．
ちょっと買いものにでてみる．
といっても隣接する北 12 生協だけど．
外は雪．
-
今日はもうぢたばたしないことにする．
歴史本とか読む．
-
生態学会大会直前なのに平和な一日でした．
-
[今日の運動]
-
エアロバイク 30 分間．
風邪ぎみなのにうんどうしてしまった．
-
[今日の食卓]
- 朝 (0730):
シリアル．
ヨーグルト．
リンゴ．
- 昼 (1300):
スパゲッティー．
キノコソース．
- 晩 (2000):
米麦 0.6 合．
麻婆豆腐．
キャベツ・キュウリ・ネギ・ハムのサラダ．
2009 年 03 月 15 日 (日)
-
0830 起床．
生活周期，
またずれぎみ．
朝飯．
コーヒー．
体温
36.7°C，
とくに問題ナシ，
かな．
まだアタマがちょっとぽーっとしてる
……
といった口実をおもいついたので，
すこしばかり怠業．
-
1030 自宅発．
曇．
気温は高い．
札幌駅の本屋をふらふらしたり，
お茶を買ったり．
1110 研究室着．
-
ちょっとまたデータ解析こんさるメイルかき，
とか．
-
ベイズカエル本の第 10 章の
感想
とか書いてて，
ふと気づいたんだが
……
まあ，
この章で使われているタイプの主観事前分布は
できるだけ使わずにデータ解析したいところなんだけど，
毒をもって毒を制する，
ともうしますか，
甲山さん的な「説明変数の濫用」に対抗するには，
こういった主観事前分布をあえて使ってみせる，
もっと言えば「主観事前分布を階層的な事前分布で生成する」
といった方式が使えるかも，
と気づいた．
-
ハナシの経緯はこうなる．
(今年度ではなく)
昨年度の甲山さん M2 院生 3 名のデータ解析は
階層ベイズモデルを初めて本格導入して，
(これまた「甲山さん的な」としかいいようのない)
多樹種系のめんどうにとりくんだ．
で，
苦闘をカサねて
ひととーりなんとか解析できた
……
と安心した矢先，
3 つのモデルのうち 2 つまでが甲山さんの指示による
「汚染」
というかですね，
ふつーに考えて因果関係が言えそうにない
「説明変数」を組みこまれてしまって，
ですね．
まあ，
統計モデルも推定計算も推定結果も
なんだかヘンてこなものになってしまったわけですよ．
-
このあたりは (一部の) 熱帯生態学データ解析にみられる
悪しき習慣，
つまり相関関係をきちんと定量化すればいいのに，
めんどうだから因果関係にする (説明変数にしちまう)
といった作法なんだよね．
-
昨年度のうちになんらかのテをうてば良かったのだが
……
相関関係を推定しようとすると，
これらの問題の場合には，
二変量正規分布を事前分布をやたらと推定してみせなくてはならないし
(これの推定計算はなかなかたいへん)，
そもそもこれらの説明変数の正体ってのが
「人間がてきとーにカテゴリーわけした樹種ごとの耐陰性」
「めまいがするほどいいかげんな高次モーメント統計量」
「ある基準に従ってずらずらと並べた樹種の順位 (ある時点での)」
といったもので，
こんなものをどうやって
二変量正規分布なんかと対応づけりゃいいのか，
と．
-
で，
McCarthy のカエルベイズ本の第 10 章「主観事前分布」
にもどるわけだが
……
この章であつかってるような
「みんなにアンケート
(例: 森が伐採されると鳥 A, B, C 種それぞれの存続に
どのぐらいいんぱくとがありますか?)
して事前分布を決めました」
といった態度はどうよ，
と思うのだけど
……
この
「アンケート」→「主観的事前分布をでっちあげる」
といった定式化というか数量的操作の手口について考えてみると，
ここで階層ベイズモデル的な発想が使えるかも，
という気もしてきた．
つまりその「あんけーと」とやらの結果が観測データの説明に使えるのではあれば
事前分布が自動的に「狭く」なり，
うさんくさければ
(超パラメーターである)
分散パラメーターを大きくして無情報事前分布にしてしまう，
と．
-
つまり，
うさんくさいものをこねあげて事前分布たち
(その中には無情報事前分布も含まれる)
をでっちあげる，
という操作に超パラメーター・超事前分布を援用できるのであれば，
上述の多樹種系の解析における
「人間がてきとーにカテゴリーわけした樹種ごとの耐陰性」
「めまいがするほどいいかげんな高次モーメント統計量」
「ある基準に従ってずらずらと並べた樹種の順位 (ある時点での)」
といった因子・数量をこねこねして事前分布を作るといったことができ，
そこに階層ベイズモデル的な考えかたをいれることで
(モデル選択というほどきちんとしたものではないにせよ)
何らかの意味で多少はマシなものになるかもね，
と．
-
えーと，
じつのところ，
モデル中の「ケガれた説明変数」
を因果関係として組みこもうが，
主観的事前分布をでっちあげるところに使おうが，
「統計モデルとしての悪さ」
はあまりマシにはならないかもしれない．
とくに，
上述の多樹種系解析は線形予測子にすべての説明変数をくみこんでいるので，
「説明変数としてほうりこむ」
と
「(「切片」とかの?) 事前分布を作るときに使う」
の本質的な差位はほとんどない場合もある．
-
しかしながら，
「統計モデルの説明」
としては多少はマシになりそうだ．
たとえば，
「人間がてきとーにカテゴリーわけした樹種ごとの耐陰性」
をそのまま説明変数としてほうりこんだりするのではなく，
このてきとーなる因子どもから「主観的事前分布を作り」
(階層ベイズ的な発想の助けを借りつつ)
その主観的事前分布にもとづいて事後分布を算出しました，
と書いて説明したほうが「どのあたりがアヤしいのか」
がわかりやすい．
因果関係のアヤしい因子・数量データを使って観測データを
説明したいので
「残念ながら主観的事前分布をもちださねばなりません」
といった態度を明確にできている．
ケガれた部分とそうでない部分が分離されている，
ように見える．
-
むろん，
査読者が「主観的事前分布なんて自然科学の研究ではありえん」
といった理由で論文を却下するかもしれないけど，
それはまた別の問題だ．
-
ここ半年ほどで気づいた「主観事前分布」のつかいみち:
-
アヤしい「事前情報」とやらを使わざるをえないとき:
上述でくだくだしくうだうだとのべてるとーりですが
(うん? 誰もふぉろーできてない?)，
ベイズカエル本の McCarthy と異なるのは，
主観的事前分布をでっちあげるときに
階層ベイズモデル
(階層的な事前分布)
といった考えかたなんかを何か利用できないだろうか，
と検討しているあたり
-
測定誤差のモデリング:
測定誤差と「個体差」の区別がつかなくて困るときは，
測定誤差の大きさを主観的に「えいやー」と決めてしまう
(あるいは測定機器のマニュアルなどにかいてある「精度」
を信じる)，
といった方法をとらざるをえない
-
これがイヤなら (少なくとも一部の)
標本を二度以上独立に測定
してください
-
蛇足なンだけど
……
測定精度がどーのこーの
といったハナシが好きなヒトいますけど，
拝聴していると
「なんだか精密に測定できてるみたいでウレしい」
といった気分的なコトが重要であるらしく，
データ解析への定量的な影響とか統計モデルの中でのとりあつかい，
といったことはあまり念頭にないような
……
-
McCarthy のベイズカエル本
(さきごろ野間口さんによる訳本
「生態学のためのベイズ法」
がでた)，
この本に紹介されてるいろいろな「事前分布のきめかた」
なんかで
ベイズ統計モデリングをやりたくはないんだけど
……
しかし，
ここまで自分の方法とずれていると，
期せずしていろいろ考えさせる本になっているなぁ，
と．
実際のところ，
必ずしも皮肉というわけでもなく，
やくにたちました．
今後は急増すると予想される
McCarthy 亜流なべいぢあん対策についてもいろいろ検討できた．
-
昼飯．
-
GLM 自由集会の準備のつづき．
「いつもの図」
を整理 & 簡素化してみる．
-
いや，
ここからは
「推定方法」のハナシはいっさい削除すべきか?
-
で，
推定方法その他はこういうふうに図示してみる，
と．
-
まだまだ改訂が必要だろうなぁ，
と危惧しつつも
……
1955 研究室発．
2020 帰宅．
晩飯の準備．
晩飯．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0900):
米麦 0.6 合．
キャベツ・キュウリ・ネギ・ハムのサラダ．
- 昼 (1300):
研究室お茶部屋．
米麦 0.8 合．
麻婆豆腐．
- 晩 (2130):
米麦 0.8 合．
ネギ・シュンギク・シイタケ・カジカの鍋もの．
2009 年 03 月 16 日 (月)
-
0750 起床．
朝飯．
コーヒー．
0900 自宅発．
曇．
0915 研究室着．
-
地環研すぐ北の情報基盤センターの大判プリンターで印刷．
[かめぽすたー]
今回は半透過な図を含む PDF ファイルだったせいか
……
Adobe Acrobat からの印刷ができん．
いろいろ試行錯誤してみたんだけど，
Mac OS X の Preview を使うとなぜかしらうまく印刷できた．
また印刷作業に一時間以上を費やしてしまった
……
-
その後，
昼飯前も昼飯後も GLM 自由集会の準備．
-
これは投影資料の終盤あたり．
-
えーと，
LaTeX で「打ち消し線」
いろいろあるみたいだけど
……
お手軽には
\usepackage{cancel}
(検索)
して \xcancel{}
とかすればいいのか．
[消された deviance]
……
とゆーことで，
すみませんが今回は deviance のハナシはやめときます．
deviance のハナシをしようと思ったら，
ポアソン分布な統計モデルの
最尤推定からはじまって，
モデル選択あたりまで紹介しないと
「意味がよくわからん」
となってしまいそうなので
……
-
17 時ごろ，
ひととーりできたので，
アップロードとか．
-
「出発前日の夕方にできた」
というのはふつーのヒトにとっては
「ぎりぎり」なんでしょーけど
……
私としては「こんなに早く
準備がおわってしまっていいのか?」
とただひたすらに呆然としております．
-
とにかく今回の盛岡大会は
「ばてないぐったりしない生態学会大会」
「学会もえつきしない大会」
をめざして蕎麦・豆腐をゆったり楽しむ
のんびり盛岡滞在にしたいと思います．
目標は「したうけうっかり引き受け」数をゼロ件に!
-
いやー，
大会から帰ってもまた〆切・督促 (被督促?)・雑用とかあったりするので，
盛岡にいるあいだは静かに休養します．
-
1835 研究室発．
1850 帰宅．
晩飯の準備．
晩飯．
-
明日は早おきして飛行機にのらんといかん．
-
[今日の運動]
-
いかん，
時間的によゆーあったのにうんどうできなかった
……
-
[今日の食卓]
- 朝 (0820):
シリアル．
ヨーグルト．
- 昼 (1220):
研究室．
北大生協にぎり飯．
- 晩 (2000):
シリアル．
ヨーグルト．
ネギ・卵の炒めもの．
ニンジン
……
冷蔵庫内ののこりもの処分．
2009 年 03 月 17 日 (火)
-
0600 起床．
ねむい．
0630 自宅発．
曇．
今日から
生態学会盛岡大会．
-
0640 JR 札幌駅着．
新千歳空港 1040 円．
0648 同発．
車内こんでる．
ねむい．
0725 新千歳空港着．
ぽーっとしてたので JAL 登場なのに
ANA な保安ゲイトをくぐってしまった．
とくに問題なかったけど．
朝飯．
-
0800 JAL 2900 便 MD81 搭乗．
0810 定刻どーり taxing 開始．
6 割ぐらいの搭乗率．
私のまわりはすいてて助かる．
機内でぐったりしてるうちに
0910 仙台空港に着陸．
えーと 1130 から始まる
(必ずしも出席する必要のない)
委員会にでるべしと命じられたので，
札幌→仙台→盛岡という経路を選んだ．
これが一番早い．
-
しかしながら
……
0938 仙台空港線にのって空港発．
0957 仙台着
……
ここから悪夢のごとき 40 分間となる．
-
手ばやく自販機で新幹線の切符かう．
盛岡まで乗車券 + 自由席特急券 5550 円．
1010 はやてにのって仙台をでようとしたところ
……
はやては全車指定席?!
2-3 年に一度しか
新幹線
に乗らない私はかくのごとき東北新幹線の常識を知らないのであった．
-
とりあえず改札口にもどって
「誤入場」のハンコ押してもらって切符窓口にならぶ．
510 円追加料金はらって指定席にしてもらう．
再入場してベンチでぐったり．
プラットホームにすべりこんできた
1038 発の次のはやてに乗ろうとしたところ
……
ウミガメポスターがない!!
さっきぐったりしてたベンチだ!
-
で，
読者たちはここで私がまた乗りおくれたと考えるだろうけど，
ここは読者たちの想像をはるかに超越するすばらしいダッシュで
エスカレイターを雪崩のごとき勢いでかけおり，
ポスターをひっつかむや上りエスカレイターの上を跳躍しつつ
プラットホームにまいもどりドアのしまりそうなはやてにすべりこんだのであった
……
ますますぐったり．
仙台駅，
悪夢の 40 分間であった．
-
またぐったりしてるうちに 1122 盛岡着．
もう完全に遅刻になっているわけだが，
これほどの波瀾万丈の大冒険を突破してきた私を誰が非難できようか．
1129 いわて銀河鉄道線にのって盛岡発．
1144 滝沢着．
15 分ほどあるくと
岩手県立大学．
[岩手県立大学]
今回の生態学会大会会場，
-
会場についていきなり
「久保さん，やつれましたね」
と指摘されました
……
うーむ，
この冬季繁忙期はいろいろあったから，
なのかな?
-
で，
大会企画委員会というのをやってる
部屋の片隅でちぢこまっていたのだが
……
大会どーするこーするといった議論を聞かされていると，
2004 年
の釧路大会の悪夢がよみがえるばかりだ
……
たとえば講演要旨集の経費に関するハナシで
「1 ペイジの編集料は 1400 円」
とか聞かされると自分はあのとき数十万円ぶんの
タダ働きをしていたのか，
とか．
-
来年度から 3 年なので
2010 年の東京大会，
2011 年の東京大会，
2012 年の近畿大会，
さらに同時に開催される EAFES にまきこまれることも確定．
いやはや．
-
全体会議がおわってがっくりしてると，
大会企画委員会のプログラム部会のヒトたちに，
プログラム部会は単なる季節労働者だからラクですよ，
となぐさめてもらった．
たしかに特に議論すべきこともなく部会は終了．
-
そういえば，
全体会議が終わった瞬間に，
生態学その他で某有名なるヒトが
「大会会長はいるかー?! メイルよんでくれましたか?!」
とどなりこんできて，
「まだ委員会は終わってません」
と手ぎわよく放りだされていた．
皆さん，
慣れておられる
……
-
E 会場で自由集会の準備．
-
1500 より
GLM 自由集会．
今日は初日 (委員会メインの日)
なのでそんなに集まらないかと思っていたら，
今回もほぼ満席 200 名ぐらい集まってくださった．
皆さん，
勉強熱心だなぁ
……
どうもありがとうございます．
-
最初に 20 分ばかり GLM の紹介など．
GLM 自由ペイジ
から投影資料をダウンロードできます．
-
それから 60 分ほど粕谷さんによる交互作用項の解説
……
交互作用は
``much confusion (ひどく混乱)''
``mine field (地雷源)''
とよばれていて，
といった刺激的な導入からはじまる．
とりあげられたおもな点は，
「交互作用の意味」は link 関数依存，
R の交互作用のあつかいは「単純な積」，
交互作用の coding
(あるいは標準化や中央化)
で「主効果の意味」がかわる，
そして主効果を「どこ」で見るかは
統計学だけでは解決しない問題
……
だけど「差」をみたら，
というのが粕谷さんの提案．
近日中に粕谷さん投影資料をアップロードするので，
おたのしみに．
-
私などとは比較にならぬほど生態学会雑用をおしつけられてる
粕谷さんはこのあとすぐに全国委員会に．
-
最後の 40 分間はふたたび私で，
offset 項ワザの説明．
いろいろ質問がでて無事に終わりました．
-
いろいろ質問の例:
ポアソン回帰 (GLM の一部)
の推定結果を利用してシミュレイションを
やるには?
-
R なら
rpois()
関数でポアソン乱数を生成してください
変数変換 & 正規分布な線形モデル (LM) で何でダメなの?
-
(いまだにこういう質問あるわけだが，
「人に教えるときにどう説明するのがよいのか」
といった理由もある)
LM でできることはことごとく GLM ででき，
しかも変数変換 LM には「それでよし」
とする統計学的な根拠がないのに対して
GLM は「カウントデータなのでポアソン分布」
といった明確な説明がある
「割算値な解析するな」というけれど，
○○という割残値は生物学的な解釈が可能であり，
無意味ではない
-
○○という概念が生物学的に意味があることと，
それを観測値 / 観測値というふうに誤差つきな量どうしの
割残値として推定してみせるのが良いかどうかは別の問題です
……
そして生理生態学で使われている呪われ割残値をあつかうには，
ベイズモデルが便利です
offset 項を使ったポアソン回帰をやりたいのだけど，
応答変数にゼロがたくさんふくまれている
-
offset 項つきの
GLMM でうまく対処できるかもしれないし，
zero inflated な統計モデルを導入するといった
工夫も必要かもしれません
-
GLM 集会後，
会場横の休憩コーナーで「臨時診療所」
の開設．
相談希望者が 5-6 人ならび，
「○○で解決すると思います．
はい，つぎのかたー」
-
その後，
1915 まで樹木の水通導の自由集会．
シラカンバやウラジロエノキは「つまった」道管を復旧できるとか，
小笠原のキバンジという外来植物は乾燥しても葉っぱを維持して
次のスコールをまつ
……
といったまにあっくなハナシを楽しむ．
-
1943 滝沢発．
1958 盛岡着．
2005 本日の宿
(東横イン盛岡駅南口開運橋)
着．
部屋に荷物おく．
こういう地方都市は 20 時すぎたら店がしまるんじゃないかと心配して，
ネットで検索したんだけど，
駅地下は 22 時までやっているもよう．
さっそく盛岡の蕎麦を楽しんできました．
-
[今日の運動]
-
まあ，
うんどう休養日だけど，
仙台駅の猛烈なるダッシュとか，
滝沢村の駅-大学往復 30 分とか
……
-
[今日の食卓]
- 朝 (0750):
搭乗ゲイトふきんでサンドイッチ．
- 昼 (1120):
東北新幹線の中でサンドイッチ．
- 晩 (2130):
JR 盛岡駅の地下の
「そば de 十割」．
[天ぷら二色もり]
1200 円．
田舎蕎麦・更級の二種類が食べられてうれしい．
どちらもおいしかった．
揚げたての天ぷらなんて何年ぶりだろー，
と感銘うけた．
2009 年 03 月 18 日 (水)
-
生態学会盛岡大会 2 日目．
0650 起床．
うだうだとぎょーむ日誌かきとか．
[開運橋]
北上川にかかる橋のひとつ．
宿の窓からの風景．
-
階下におりようとエレヴェイターをまつのだが，
上からおりてくるハコの中は生態学研究者過密スシづめ状態．
2000 人ばかりの学会参加者で飽和する盛岡市の宿泊施設，
そのあちこちで今朝はこういう現象が発生しているんだろう．
かくのごとき非日常というほかないハコが
目の前を落下していく様子を三度ばかり見送ってから，
よーやく下におりることができた．
宿の無料給食．
-
0830 宿発．
JR 盛岡駅前のアヤしい行列にならんでしまった．
これまた会場いき臨時バスをまつ生態学研究者の長蛇の列．
0850 駅前発．
0920 会場の岩手県理大学着．
今日は曇．
[長い廊下でポスター発表]
ポスター会場 B と C は長さ 200 m ぐらいの廊下．
私は C に
ウミガメ
ポスターをはった．
まあ，
ポスター間がゆったりしてていいんですけどね．
30 枚準備しておいた縮小版はすぐに消滅してしまった．
[体育館でもポスター]
会場 A にもふらふらと遊びに行った．
ここ数年，
口頭発表やシンポジウムの会場で座っていることはなく，
ポスター会場をひたすらふらふらしていることが多い．
-
1215 から 14 時すぎまでコアタイム，
ということでポスター横でウミガメ上陸数の階層ベイズモデルの
解説につとめてみる．
ポスター発表のよいところは，
聴き手の反応をみながら説明方法を変えていくうちに，
全体的に多少はマシな方法に改善されていく，
といったあたりかな
……
-
ということで，
ウミガメについては「個体群動態の推定」
といった説明ではなく
(ウミガメの個体群は何とも定義しがたいので)，
むしろ
「ウミガメたちはある砂浜を産卵場所として
見捨てつつあるのか?」
といった問題のたてかたのほうがマシなように思えてきた．
これはハナシを聞いてくれたヒトたちが
そのように指摘したわけでは
……
しかし皆さん相手に何度も説明しているうちに，
自分一人では思いつきそうにない改訂案がでました，
というかんぢで．
-
またポスター会場ふらふら放浪を．
おお，
WinBUGS で個体群動態を推定している野心的な研究が二件も!
そして
Levy flight
のハナシはいつも楽しいなあ．
-
本日もまた
「久保さん，やつれてますね」
「目がつかれてますね」
「くまができてますね」
指摘が
……
やはりこの冬季繁忙期におけるうんどう不足がよくなかったのかな?
安楽椅子生態学こそ，
体力しょーぶですな．
-
ときどきウミガメポスターふきんに回帰して，
またウミガメ説明を試みてみたり．
-
GLM は普及しつつある
……
と言いたいところですが，
まだまだですな．
-
17 時すぎにポスターを撤収しようとしたら，
なぜか急にお客さんがぞろぞろと．
撤収期限ぎりぎりに撤収．
-
などとぐずぐずしてたら，
すでに企画集会「樹形変化」の F 会場は
満パイ状態で入りこめなかった．
かかる呪われあろめとりー世界において，
我らが矢澤さんが立派な発表 (私もなぜか共著者の一人)
をしたことを疑いません．
ええ，
まったく疑いませんとも
……
と自分を納得させつつ，
岩手県立大学からすばやく撤退．
-
1750 会場発．
1805 いわて銀河鉄道の滝沢駅着．
1821 同発．
盛岡まで 360 円．
1836 JR 盛岡駅着．
駅前で
A 棟 8F の斎藤君と晩飯．
えー，
何でも斎藤君「念願の」ポスター賞受賞が実現したとのこと
& えーいまたしても本日も某所に置き忘れてしまった
ウミガメポスターを回収してきてくれたので．
2005 宿にもどる．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0810):
宿の無料給食．
- 昼 (1200):
ポスター会場のかたすみでカロリーメイト，
とか．
- 晩 (1900):
盛岡駅前の
盛楼閣
で斎藤君と焼肉 &
盛岡冷麺．
なかなかおいしかった．
お値段もなかなかでしたが．
2009 年 03 月 19 日 (木)
-
生態学会盛岡大会 3 日目．
0650 起床．
宿の無料給食．
雨．
出ようとしたら傘を部屋に置き忘れてきたことに気づいた．
また，
悪夢のごとき生態学会参加者すしづめエレヴェイターにとりくむ．
-
0805 宿発．
今日も盛岡駅前からバスにのる．
0818 駅前発．
0850 会場である岩手県立大学着．
-
今日もひたすらポスター会場をあるきまわっていた．
今までいろいろな下うけ計算やってきたおかげで，
いろいろな分野のポスターたのしめる．
オスの求愛行動の資源分配・血縁度，
なわばり侵入の情報戦，
リター分解速度の MCMC 計算，
植物の分布拡大と集団遺伝学，
開花遺伝子突然変移と多年性
……
-
そういえば，
粕谷さんのポスターも展示されてて，
X vs Y といった回帰をするときに，
X に「誤差」があるのを無視してセンをひっぱると
推定値セットにバイアスが入る
……
といった問題の指摘と，
ベイズで解決するのがよいのでは，
といった提案でした．
たしかにベイズ化がよいでしょう．
-
一日中あるきまわって全体的な展望をえたような気分になったので，
1635 会場発．
あいかわらずの雨．
あるいて滝沢駅へ．
1648 滝沢駅．
このあたりは豆腐が名物らしいんだけど，
今まで食べてなかったので駅の「豆腐おかず屋」
みたいなところで豆腐コロッケをひとつ買って食べた．
-
1658 滝沢発．
1715 盛岡着．
駅の本屋をぷらぷらしたり，
みやげもの屋をのぞいたり，
また蕎麦屋で晩飯くったりして，
1840 宿にもどる．
いきなり寝てしまって，
夜中におきた．
うーむ，
生活周期はどうなることやら
……
-
そして明日は岩手県から撤退するばかりである．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0750):
宿の無料給食．
- 昼 (1200):
ポスター会場周辺でカロリーメイト，
とか．
- 晩 (1820):
[鴨南そば]
二日前と同じく，
JR 盛岡駅の地下の
「そば de 十割」．
鴨南そば (980 円)
の default は更級蕎麦なんだけど，
田舎蕎麦にも変更可能だったのでそちらにしてみた．
2009 年 03 月 20 日 (金)
-
生態学会盛岡大会 4 日目
……
であるが，
私は朝から札幌にもどる．
0630 起床．
生態学会参加者すしづめエレヴェイターに乗るのは
もはやイヤなので早々に撤退．
0650 チェックアウト．
宿の
(塩分過剰なる)
無料給食もイヤになったので盛岡駅のコーヒー屋で朝飯．
-
朝飯くいつつ本日の計画を考える．
花巻空港いきのバスは 0750 発と 0900 発．
遅いほうで問題ナシ．
-
ということで，
けっこう重い荷物を背負い，
くるくると丸めた (だけの) ウミガメポスター片手に盛岡市内を
散歩してみることに．
[07:24 JR 盛岡駅発]
14 年前の 1995 年にも盛岡 (岩手大学)
で生態学会大会があった．
今よりもずっと規模が小さく
……
しかし駅とか何も記憶してない．
今回はあまりよく見えなかった岩手山を少しおぼえているぐらいだ．
[07:26 開運橋]
以前にも宿からの写真をのせたけど，
北上川にかかる橋のひとつ．
これは歩道橋の上から見おろしたもの．
駅から歩いて 2-3 分．
ときどきぱらぱらと雨がふる．
[07:36 盛岡城]
10 分ほどあるくとついた．
建物は何もなく，
石垣のみ．
公園になっている．
別名はは不来方 (こずかた)
城とのことだけど，
Wikipedia
によると厳密には別モノだそうで．
昨年
行った福岡城よりかなり規模は小さい
(福岡城がでかすぎるわけだが)．
一昨年
行った松山城よりやや小さい城が平地にあるかんじ．
[07:44 銅像のない台座]
本丸跡におかれている．
銅像は藩政時代とは何の関係もない，
日露戦争で戦死した騎兵中尉
(手前の写真)
が叙勲
(つまり朝敵だった南部藩の汚名雪辱)
された記念だそーで．
騎兵将校銅像は
1944 年に資源不足な戦争のタシに鋳つぶされてしまった．
それから 2/3 世紀にわたって「本丸あとに空虚なる台座」
をたてまつる盛岡市民たちの，
何やら哲学的なのか仏教的なのか警世的なのか外部のものにはよくわからぬ
ぽりしーについて考えさせられた．
[07:45 盛岡市街]
城の良いところは市街地を俯瞰できるところですね．
殿様の視点で．
ところで城内には新渡戸稲造 (生家は城の近く)
や石川啄木の石碑・歌碑があったり．
後者は例の
「不来方（こずかた）のお城の草に寝ころびて空に吸はれし十五の心」
ですな．
金田一京助の書だった．
その他，
よくわからぬ石碑も多数．
宮沢賢治のもあった．
[07:57 城の堀]
城の東側の堀が公園化されてかくのごとき池になっている．
この堀の池にそって南にむかうと中津川にでる．
[08:06 中津川]
これも盛岡城の防御構造の一部．
ハクチョウやカモなど水鳥が子連れでのんびりしている．
まあ，
休日のこの時間はふつーのヒトが出歩く時間でもないからな．
[08:20 中津川・北上川合流点]
中津川にそって下っていくと北上川との合流点につく．
この写真ではどこで合流してるのかもわからんけど．
ちなみに築城時代の合流点はもうちょい上流にあり，
つまり北上川が城の西側の防御構造になってそーで．
手前の広場は合流点の後背地に作られた公園．
フットサル競技場とか．
[08:25 三つの川の合流点だったか]
地図で説明するとこんなかんじ．
北上川・中津川だけでなく雫石川もここで合流している，
とのこと．
ふーむ，
「強く美しい河川整備」
とは，
ちと自己陶酔的ですな．
[08:25 北上川]
合流点のちょっと上流にある
不来方 (こずかた) 橋からの下流方向のながめ．
さきほどの開運橋のすぐ下流がわにある橋．
[08:31 JR 盛岡駅にもどる]
JR 盛岡駅前にある地図．
中央やや右の緑が盛岡城．
下方やや左の黄色が JR 盛岡駅．
一時間ちょっと歩いた．
荷物がやけに重かったなぁ
……
と調べてみたら，
重さ 1.5 kg ぐらいありそうな生態学会講演要旨集がでてきた．
一度もひらかなかった百数十枚の紙束．
リサイクルされることを期待しつつ，
駅でぽいっと捨てる．
[08:47 じゃじゃ麺]
バス出発まで時間があったので，
駅の立ちぐい蕎麦屋で
じゃじゃ麺
たべてみた．
440 円．
これで盛岡の蕎麦・冷麺・じゃじゃ麺の三つを制覇したわけだが
……
じゃじゃ麺だけがぢゃんくっぽいというか，
学会開催地の「そうるふーど」研究家の田代さんあたりに言わせれば，
3 年前の
新潟
における
``
イタリアン''
(洋風ソースかけ中華麺)
のごときものなのか?
まあ，
立食いで 440 円しはらうのはおススめではなさそう．
究めたいヒトは専門店で検証してください．
[08:53 JR 盛岡駅前]
こんな巨大な看板がたっておりました．
2000 人の学会参加者は盛岡市の経済にどういった
いんぱくとがあったんでしょうなぁ
……
たしかに盛岡駅周辺ではアウトドアなジャケット
(ゴアテックスとか)
を着用し，
デイパックなどかついでうろうろしてる人たちみかけました．
私もデイパックというか，
3-4 日出張用のややでかいバックパックかついでいますけど．
空港いきバス乗り場ふきん，
北大農学部のぷろふぇっさーたちがうろうろしておられたので，
目立たぬ位置にてバスをまつ．
いや，
べつに仕事を押しつけられるわけでもないだろうけど
……
0900 空港いき高速バスにのって盛岡駅発．
1260 円．
[09:46 いわて花巻空港]
どこかの駅みたいですが，
空港です．
-
上の写真など整理してるうちに出発の時間．
1055 JAL 2830 の MD81 に搭乗．
理由はよくわからないんだけど満席．
東北のヒトたちは三連休に北海道に遊びにいくのか?
-
くるくるてきとーに巻いただけ，
そのままのウミガメポスターを機内にもちこんでいたのだが
(往路も，そのときは問題なかった)
……
客室乗務員のボスねーさんがつかつかとやってきて説明されるところによれば，
この長さ 90 cm の丸めポスターは非常不時着時の脱出
に重大なる障害となる危険物であり，
かくのごとき剣呑なるシロモノは隔離したうえで
客室乗務員の監視下におくことに同意されたい，
さもなくば
……
もちろんおびえきった私が
かくかくとアタマを上下するやいなや，
あわれウミガメポスターがいづこにか連行されていったのはのべるまでもない．
-
どうやらこの飛行に「万が一」の不安を感じているのは
客室乗務員の皆様だけでなく，
機内後方にいるとおぼしき一男児も
かくのごとき危険きわまりない運航が予期される機内に
おしこめられているのがいかに不本意でありかつ理不尽なる所業なのか，
といった内容を保護者だけでなくかなり遠方の周囲に対しても主張している．
論点を単純化 & 限定し，
それを執拗に執拗に何度もくりかえし，
さらに大音声
(なかなかよくとおる声質でもあり，
機内アナウンスをうち負かしている)
を発しつづけるすばらしい持久力を発揮している
……
ふーむ，
将来はある種の研究者として成功できるのではなかろうか．
-
人々のさまざまな不安を機内にかかえたまま定刻より 10 分遅れの
1100 taxing 開始．
離陸．
飛行時間は短いので飲物はでない．
-
機内では粕谷さんの論文原稿を拝読していた．
統計学もので，
その統計学的結論は明快なのだが，
「さて，
このような問題のある状況が生じるのは具体的にはどういう場合なのか?
問題が生じているかどうかデータを見ただけでわかるのか?」
といったあたりを検討してみると，
なかなかの難問のように思えてきた．
-
この粕谷さん問題を検討してるうちに，
多くの人々の予感をうらぎり，
なにごともなく
1150 新千歳空港着．
ウミガメポスター，
かえしてもらえた．
1204 JR にのって同発．
車内，
かなりこんでる．
立ったまま，
また粕谷さん問題を検討してるうちに
1240 札幌着．
気温ゼロ度ふきん．
かなり強い吹雪．
また冬の世界にもどった．
-
昼飯かって 1255 研究室着．
いきなりお茶部屋でお留守番院生たちから
30 分間にわたって盛岡大会尋問．
釈放されて昼飯．
-
「盛岡ばてしない」
「学会もえつきしない」
とゆーのが今次大会で実現すべき最大の目標であったわけだが
……
とりあえずこの日は仕事がススまず．
-
1940 研究室発．
なかなか寒い．
買いもの．
2010 帰宅．
晩飯の準備．
晩飯．
-
[今日の運動]
-
うーん，
まあ，
重荷 (12-13 kg?)
をせおって一時間強の盛岡散歩，
とか?
歩く速度はそれなりに速かったんで．
-
[今日の食卓]
- 朝 (0710):
JR 盛岡駅のターリーズでエッグサンドセット．
480 円．
- 昼 (1350):
研究室．
幕の内弁当．
500 円．
- 晩 (2130):
米麦 0.8 合．
キャベツ・ミズナ・ネギ・ハムのサラダ．
麻婆豆腐．
