ぎょーむ日誌 2009-02-(21-28)
2009 年 02 月 21 日 (土)
-
0810 起床．
なんだか，
うーむ，
なかんぢ．
朝飯．
コーヒー．
0940 自宅発．
大雪．
0955 研究室着．
-
月曜日の公開講座の準備のつづき．
まにあうのか?
-
窓の外は大雪．
-
うー
……
よーやくにして 90 分 + 60 分の公開講座のうち，
うしろの 60 分のほうができた
……
時刻はすでに 1745．
日没．
いったいどうなることやら．
-
明日は 1640 千歳離陸，
か
……
-
うーむ，
オワリがみえん
……
2045 研究室発．
雪はやんでる．
2115 帰宅．
ちょっとうんどう．
晩飯の準備．
晩飯．
-
晩飯くったら仕事する気がメルトダウン
……
-
[今日の運動]
-
腹筋運動 30 ×
3 回．
背筋運動 30 ×
3 回．
腕立ふせ 10 ×
3 回．
スクワット 100 回．
ストレッチ．
-
[今日の食卓]
- 朝 (0850):
シリアル．
ヨーグルト．
- 昼 (1420):
研究室お茶部屋．
米麦 0.6 合．
レトルトパウチドカレー．
- 晩 (2220):
米麦 0.8 合．
ホウレンソウ・豆腐の味噌汁．
サバの電磁波酒蒸し．
2009 年 02 月 22 日 (日)
-
0750 起床．
朝飯．
コーヒー．
旅じたく
……
するほどのこともないか．
洗濯．
怠業
……
ではなく明日の公開講座の準備のつづき．
-
……
なかなか終わらない．
1400 自宅発．
けっこう雪がふってる．
飛行機，
とぶのだろうか?
1415 研究室着．
-
ちょっと準備をススめて
……
-
1520 研究室発．
あいかわらず雪．
1525 JR 札幌駅着．
空港で時間があまりないと思ったので
(これはじつにおめでたい錯覚であった)
統数研へのおみやげなど買う．
1540 札幌発．
1617 新千歳空港着．
1040 円．
-
空港は今日も混乱のきわみであった．
雪のせいであちこちで「遅延」が．
ヒトもびっしり．
苫小牧の平尾君にあったけど，
その後すぐに 2 時間おくれの便で出発していった．
-
まちあい室の片隅で電力窃盗しつつ，
ぢりぢりと準備をススめる
……
しかし何ともおちつかない．
2 時間以上おくれるらしいんで，
サンドイッチ買って簡単に晩飯．
-
1850 ADO 22 便 B-767-300 に登場．
2 時間 15 分おくれで 1855 taxing 開始．
1910 離陸．
機内ではひたすら休む．
2040 羽田着．
-
2055 京急にのって羽田発．
切符はどうしようかと思ったんだけど，
とりあえず渋谷まで 560 円．
2112 品川着．
JR に乗り換え．
その前にこんびにで朝飯など調達
(あとで宿の近くにもこんびにあるとわかった)．
2119 山手線に乗って同発．
渋谷から歩くのがイヤになったので
2129 恵比寿着．
-
地下鉄日比谷線に乗って
2135 恵比寿発．
2139 広尾着．
2145 宿着．
-
そして準備はまだまだつづく
……
いつになったらオワるのやら
……
-
[今日の運動]
-
[今日の食卓]
- 朝 (0820):
米麦 0.6 合．
ホウレンソウ・豆腐の味噌汁．
- 昼 (1220):
米麦 0.8 合．
卵キムチ炒飯．
- 晩 (1800):
新千歳空港でサンドイッチ．
2009 年 02 月 23 日 (月)
-
あ，
これは昨晩の雪の新千歳空港だ．
-
0500 起床．
東京西麻布．
3 時間弱ねてしまった．
本日 1000 から公開講座の準備．
はい．
まだできてません．
-
ぢりぢりと準備のつづき
……
-
なんとかできたよーに思えたので，
0900 宿発．
雨．
0910 南麻布の 統数研
着．
2F 講堂の公開講座
マルコフ連鎖モンテカルロ法の基礎と実践
会場へ．
公開講座を支えてくださってるスタッフの皆さんにご挨拶．
-
伊庭さん
はまだみえていない
……
開始時間直前にあらわれ，
教室前方で PC の準備．
そして立ったまま朝食のおにぎりとお茶を．
なんだかスゴい．
-
受講者は 80 名ぐらい．
北海道からくるとちょっと異様に思えたのは，
マスク着用なヒトが多いとゆーあたり．
スギ花粉がすでにとびちってるようで．
-
公開講座開始．
1000-1100 伊庭さん「アルゴリズムとしての MCMC」．
内容的には
岩波「統計科学のフロンティア12
計算統計 II
マルコフ連鎖モンテカルロ法の基礎」
の p.27 あたりまで．
Java によるイジングモデルのデモとか．
最後にベイズと MCMC の関係について．
MCMC は「統計学の範囲に限っても階層ベイズへの応用だけではない」
わけだが，
「階層ベイズ + MCMC のセットが当たったのは事実
(+ WinBUGS の 3 点セット)」．
Metropolis-Hastings 法は，
ちょっと Gibbs sampler よりな Metropolis，
といったハナシも．
-
1100-1230 久保「ベイジアンモデリングと MCMC」．
生態学ねたではない架空例題で階層ベイズモデルを説明しようとしたのだが
……
いまいちだったか?
90 分のうち 75 分あたりで終わってしまったのだが，
その後に質問がいろいろと出た．
いきなり「ベイズでモデル選択するには」
といったちょー難問だったので，
いささかの迷いもなく伊庭さんに回答をお願いする．
伊庭さんの回答の結論もやはり「難しいですね」
ということなのだが，
そこにいたるまでにいろいろとありがたいハナシが聞けるのである．
[久保説明風景]
受講者内の知人のかたがとってくださいました．
どうもありがとうございます．
うーむ，
いかにもアヤしい講師だ
……
しゃべってる本人にはちょっと「へろへろ」疲労感あり，
コトバつまりぎみ．
辞任したさきの財務大臣
(この時事ねたは一年後には忘却されて意味不明だろう)
のキモチが少しわかった．
私はカゼ薬も酒ものんでないけど．
-
昼休み．
昼飯くいつつ午後の投影資料のチェック．
-
午後の部．
1330-1430 久保「R と WinBUGS の使いかた」
ということで午前にやった「打率の推定」
階層ベイズモデルを推定計算させるための
ひたすらコマかいハナシ．
午前・午後の久保のハナシの投影資料・ファイルは
こちら
で公開してます
(配布資料はかなり修正する必要があるのでまだアップロードしてない)．
-
これまた 60 分もち時間のうち 45 分ぐらいで終わったのだけど，
質疑応答があって，
伊庭さんからのまとめコメントで無事に終了．
質問時など，
スタッフのかたがマイクの支援などやってくださった．
-
1430-1600 伊庭さん「アルゴリズムとしての MCMC」．
内容的には
岩波「統計科学のフロンティア12
計算統計 II
マルコフ連鎖モンテカルロ法の基礎」
の p.28-98 あたりの範囲のものが多かった．
ベイズの補足，
マルコフ連鎖の収束定理 (エルゴート性)，
レプリカ交換法とその背景，
MCMC のベイズじゃない応用 (魔方陣の探索とか)
……
-
16 時すぎ，
公開講座ぶじに終了．
本日きいた伊庭さんコメントいくつか
……
例によって久保の聞き書きなので，
まちがってたりする部分もあります．
-
統数研ではベイズや MCMC に関する長い研究の歴史があって，
そもそも赤池さん自身「AIC の次の研究」
としてベイズにとりくんでいた．
AIC によるモデル選択は離散的というか，
「モデル A かモデル B か，どちらか一方」
といった「かくかく」したものなので，
ベイズによってそのあたりなめらかにならないか，
といった方向性で．
-
このように「モデル選択はやらなくてもよいように」
といったことがベイズを研究する理由のひとつだった．
しかしモデル選択のようなことができれば便利な場合もある．
モデル選択については「ベイズ統計と統計物理」
でもふれている．
Bayes factor が使えるかも，
といったハナシはあるのだが，
事前分布に無情報事前分布があると
……
-
はずれ値は t 分布などスソの長い確率分布で
ベイズモデル化できるけれど，
その MCMC 計算は難しいものになりがち．
変化点の検出だの，
グループ数を指定しないグループわけだのも同様．
-
しかし予備検定とかそういうことをやらずに，
まとめてひとつの統計モデルであつかいたい，
というのがベイズな方向性のひとつ．
そういえば，
伊庭さんの講義の中で，
まだ MCMC わざが統計物理なヒトたちだけが使うものだった時代，
その計算技法はマンツーマンで伝授されていくものだったんだけど，
こんにちのように MCMC 使いが急増してしまってだいじょうぶかなぁ，
といったハナシも．
-
終了後，
統数研-ドイツ大使館前の坂をくだった
スープ屋
さんで伊庭さんとスープ (なぜ?) とコーヒー．
ふーむ，
岩波書店ってそんなに「攻撃的」なところだったんですか
……
-
ふと，
南麻布の統数研にくるのはこれが最後かも，
と思ったり．
今年
移転
するそーで．
-
1738 日比谷線広尾発．
1742 恵比寿着．
1748 JR 山手線に乗って発．
羽田空港まで 550 円．
1757 品川着．
1801 京急にのりかえて発
……
のはずが，
乗りそこなってしまった．
8 両編成なのにホームの一番前方で待ってたため．
そして人ごみなので移動できず．
1805 久里浜いき快特にのって発．
1815 蒲田着．
1822 京急羽田線に乗って発．
1832 羽田空港着．
-
1915 ADO の B767-300 に搭乗．
15 分おくれで
1930 taxing 開始．
なかなか離陸しなかった．
機内で晩飯，
あとはひたすら休む．
2110 新千歳空港着．
2130 空港発．
2210 JR 札幌着．
気温 -9°C．
2223 帰宅．
ぐったり．
-
帰宅したら不幸のメイルが．
A 棟 7-8F 内の PC 一台がウィルスのたぐいにやられて暴走，
あちこちの IP アドレス
445 番ポートに対する無差別大量アクセスやってるので，
北大ネットから切り離されてしまった．
やれやれ，
明日の朝一番はこれに対処しないといけないよーで
……
-
[今日の運動]
-
[今日の食卓]
- 朝 (0700):
こんびに握り飯．
- 昼 (1240):
統数研の公開講座会場で
こんびに握り飯．
- 晩 (2010):
AIRDO 機内で羽田空港サンドイッチ．
機内さーびすのゴボウスープ．
2009 年 02 月 24 日 (火)
-
0730 起床．
まだねむい．
朝飯．
コーヒー．
0855 自宅発．
曇．
0910 研究室着．
-
昨日も書いたように，
私が東京にいるあいだに A 棟 7-8F LAN (通称「闇ネット」)
が北大ネットから切り離されて音信不通状態．
原因は A 棟 7-8F のどこかで
ウィルスとかにやられた暴走 PC が
学外にむかって無差別大量の不正アクセスをやっていたため．
(2009-02-23 13:13)
ネットワークチームです。
久保先生が管理者として登録されています端末
IPアドレス：133.50.217.249
MACアドレス：00:A0:B0:3C:0F:6A
機種名：router box
設置場所：A802
より、学外の不特定多数IPアドレス445番ポートに対して、
短時間に多量の不正アクセスを行っているのを検知しました。
意図されていないアクセスの場合、当該端末がウィルス・ワーム・
マルウェアなどに感染、またはrootkitが組み込まれた可能性があります。
お手数ですが端末のご確認をお願い致します。
そしてなぜかこういうときだけは迅速に止められた．
(2009-02-23 14:04)
ネットワークチームです。
現在、学外への445番ポートに対する通信が検出され続けていますので、
影響を考えまして通信を停止しましたのでご了解ください。
-
とりあえずの対処療法を．
A 棟 7-8F ルータで TCP/IP フィルタリングを．
今までやってなかった私が LAN 管理者として無能，
という側面あったわけだが．
昨夜のうちにネット上で調べて見つけておいた
「 LAN→WANで必ず遮断しなければならないパケット一覧」:
プロトコル ポート番号 サービス名 コメント
TCP,UDP 135 epmap RPCに使用
TCP,UDP 137 netbios-ns コンピュータ名の名前解決に使用
TCP,UDP 138 netbios-dgm コンピュータの一覧表示に使用
TCP 139 netbios-ssn ファイル、プリンタ共有に使用
TCP,UDP 445 microsoft-ds ダイレクトホスティングSMBサービス
なるほど 445 番もとうぜんのごとくに含まれています．
えーい，
こういうゐんどーづ関係のポートはふさいでおく必要があったんだな．
ともかくルータの設定．
-
北海道大学情報基盤センター
ネットワークチーム
(内線 3456, 2948)
に電話かけて (何しろネットが使えないのでメイルがだせない)，
再び北大につないでもらう．
-
とりあえずうまくいった
……
と A 棟 8F 内をうろうろ点検してたら，
同時多発的にともうしますか，
別の事故が．
数年前に買った安物の
switching hub
(I-O data ET-FSWH8S2)
がまたまた突然死．
「またまた」というのは 8 台ぐらい同時に
まとめて買ったんだけどここ 1 年ぐらいでそのほとんどが
ナゾの「突然死」をとげる，
という事故がつづいていまして
……
とりあえず研究室のかたすみにころがってたばか hub に交換．
-
1030 より
研究室セミナー，
今日は平田さんで紫外線がサロベツのミカヅキグサ・ヌマガヤ
に何か影響あるか，
というハナシ．
まあ実験室内の実験とはちがうので，
いくらがんばってもそんなに顕著な「処理差 (紫外線カットフィルム差)」
なんてでないのだけど，
この二種の植物差についてはいろいろと考えさせられる．
ヌマガヤのほうが紫外線吸収物質の量が多い (らしく)，
同じ期間内のバイオマスも多い (?)，
いっぽうでミカヅキグサのほうが地上部への分配重視で，
繁殖開始齢が早そう (江川さんのデータ)
……
なんか関連ありませんかね?
-
研究室内に院生の人数がそろってきたので
(つまりこの中に「犯人」がいるだろうということで)，
セミナー終了後に「感染者」さがしを始める．
院生部屋の回線上にばか hub をかませて
(このハコはこういう状況では必要不可欠な道具)，
そこに自分の ThinkPad をつないで
tcpdump
という TCP/IP パケット盗聴ツールとでもいうべきもので
アヤしげなパケットをさがす．
-
tcpdump -nn
して他の host の 445 番ポートに大量不正アクセスしている流れをさがす
……
A802 室で あっさり
みつかった．
LAN 内 IP アドレス
192.168.1.149
だ．
やはり院生の PC だった．
下の
log
を見ればわかるように，
乱数をつかって「近傍」らしき IP アドレスを生成して，
なんとかもぐりこもうと努力している
……
効率わるいというべきなのだが，
これであっさり感染してしまう院生もいるわけで．
11:59:35.015672 IP 192.168.1.149.27181 > 192.168.59.208.445: S 1688427587:1688427587(0) win 65535 <mss 1460,nop,wscale 2,nop,nop,sackOK>
11:59:35.015876 IP 192.168.1.149.27180 > 192.168.59.207.445: S 1927696528:1927696528(0) win 65535 <mss 1460,nop,wscale 2,nop,nop,sackOK>
11:59:35.016765 IP 192.168.1.149.27182 > 192.168.59.209.445: S 3167003295:3167003295(0) win 65535 <mss 1460,nop,wscale 2,nop,nop,sackOK>
11:59:35.017949 IP 192.168.1.149.27183 > 192.168.59.210.445: S 1402010468:1402010468(0) win 65535 <mss 1460,nop,wscale 2,nop,nop,sackOK>
11:59:35.018161 IP 192.168.1.149.27184 > 192.168.59.211.445: S 1871805862:1871805862(0) win 65535 <mss 1460,nop,wscale 2,nop,nop,sackOK>
11:59:35.018369 IP 192.168.1.149.27185 > 192.168.59.212.445: S 2931914747:2931914747(0) win 65535 <mss 1460,nop,wscale 2,nop,nop,sackOK>
11:59:35.018562 IP 192.168.1.149.27186 > 192.168.59.213.445: S 3740762442:3740762442(0) win 65535 <mss 1460,nop,wscale 2,nop,nop,sackOK>
11:59:35.003877 IP 192.168.1.149.27177 > 65.54.245.72.25: . ack 89 win 32725
11:59:35.004052 IP 192.168.1.149.27177 > 65.54.245.72.25: F 7:7(0) ack 89 win 32725
11:59:35.004989 IP 65.54.245.72.25 > 192.168.1.149.27177: . ack 8 win 1460
11:59:35.015352 IP 192.168.1.149.27178 > 208.36.123.235.25: P 17:49(32) ack 34 win 32749
11:59:35.016949 IP 208.36.123.235.25 > 192.168.1.149.27178: P 34:42(8) ack 49 win 1460
11:59:35.017688 IP 192.168.1.149.27178 > 208.36.123.235.25: P 49:81(32) ack 42 win 32747
[地環研前の雪だるま]
[統数研ぐっず]
昨日，
統数研でおみやげにいただいた統数研ぐっず．
エコバッグとかポストイットとか．
秀逸なのは
「
日本人の国民性調査」
トランプ．
カードのオモテには国民性調査アンケートの質問が書いてある．
-
[今日の運動]
-
こういう日にこそうんどうすべきなんだろうけど
……
-
[今日の食卓]
- 朝 (0820):
米麦 0.6 合．
ネギ・卵の炒めもの．
- 昼 (1330):
研究室お茶部屋．
米麦 0.8 合．
レトルトパウチドのカレー．
- 晩 (2200):
米麦 0.8 合．
ダイコン・ニンジン・ネギ・豚肉の味噌汁．
2009 年 02 月 25 日 (水)
-
0730 起床．
自宅のコメとかがなくなったので，
0740 自宅発．
曇．
気温 +2°C
なのですごく暖かく感じる．
0755 研究室着．
朝飯．
コーヒー．
-
昨日のウィルスどたばたのぎょーむ日誌かき．
こういうバカげた事件に時間をとられて，
さらにその記録に時間をとられるのはどうよ，
とも思うわけだが，
こういう事件は今後も発生するのであり，
そのときにこういう記録が役にたつ．
昨日も，
以前の
ぎょーむ日誌
のおかげで作業時間はかなり短縮できたと思う．
-
メイルかき．
一昨日の公開講座の参加者のかたたちから
質問メイルなどいただく．
公開講座に参加した人たちは分野も経験もかなりさまざまなかんじで，
統計モデルとか初めてといったかんじのヒトから，
私などよりよほどベイズとかに詳しそうなヒトたちまで
……
-
で，
たとえば共役な事前分布を使えば事後分布は解析的に導出できるので，
Gibbs sampler は簡単に作れる．
しかし久保例題のように二項分布 * 正規分布といった統計モデル
の事後分布が WinBUGS でサンプリングできるのはなぜ?
といった疑問をもつヒトもいたり
(当日もそういう質問があった)．
とりあえずの回答例:
ご指摘のとおり，単純な Gibbs sampling では事後分布からサンプルは得られ
ません．そこで (例の質問の回答でもちょっとだけふれましたが) 数値的な
Gibbs sampling といった方法が使われます．要するに提案分布を使わずになん
とか事後分布からのサンプルを得たい，という発想ですね．アルゴリズムの名
前としては adaptive rejection sampler (ARS) とか slice sampler などあり
ます．
日本語で書かれた本でこのあたりにふれているのは
丹後俊郎. 2000. 統計モデル入門. 朝倉書店.
です．ARS などについて説明されています．英語本だといろいろあるのですが，
たとえば
Gilks, Richardson and Spiegelhalter. 1996. Markov chain Monte Carlo in
Practice. Chapman & Hall/CRC.
Gelman et al. 2004. Bayesian Data Analysis. Chapman & Hall/CRC.
は丁寧に説明されています．
当日も説明したんだけど，
WinBUGS は Gibbs sampler と呼ばれているが，
解析的に導出できる事後分布からの Gibbs sampling だけでなく，
じつは内部では上述の ARS や slice sampler による
いわば「数値的な」 Gibbs sampling さらに M-H 法も使われているそーで
(WinBUGS
user manual によると)．
伊庭さんのうまい説明を借用すると，
WinBUGS は「外からみて」
Gibbs sampler であるかのように見えているだけ，
ということ．
-
窓の外はまたいつのまにやら雪．
-
一昨日の公開講座の配布資料の修正にとりくむ．
皆さんにお配りしたもの，
けっこうヒドいもので
……
なぜ，
ヒドいかといういいわけのたぐいはぎょーむ日誌
2009-02-05
あたりに．
-
昼飯後も修正作業のつづき．
-
最低限の修正はひととおり終わったように思えるので，
公開講座のペイジ
からダウンロードできるようにする．
時刻は 1540．
ふう．
-
それをまちかまえていたかのように「蟻」解説の修正ファイルが
……
うう．
-
最初は技術的な難問にとりくんでみるか
……
引用文献書式を
jplain
ぢゃないやつ (ふつーのともうしますか apalike 的
ともうしますか)
にしなさい，
という指示．
-
それが簡単にできるのであれば最初からやってるわけだが
……
とりあえず，
makebst
(対話的な bst ファイルの
生成)
してみたんだけど，
やはり日本語文献データ混在では動作しない．
今まで，
ヘンてこなるぢゃーなるの悪夢のごとき reference
書式という難敵の数々を撃破してきた
makebst
も文字コードの絶壁はこえられないのか?
えーい，
どうしたらいいんだ
……
-
しかしながら，
捨てる神あらば拾う神あり，
というべきか
……
武田史郎さんが公開しておられる
jecon.bst: 経済学用BibTeXスタイルファイル，
これが今回の「正解」だった!
-
jecon
には詳しいチュートリアルがあり，
それをまず読んでから，
たいへんわかりやすく整理されている
bst ファイル内の関数の設定を変えていくと
(あと少しばかりそれ以外も改造すると)，
要求されている書式とほとんど同じ日英混在言語版ものができた．
助かった．
-
おっと，
ぢゃーなる名は省略名ではなくもとにもどさんといけませんな．
これは以前に作ったジャーナル名 → 省略名
変換スクリプト (in Perl)
をちょっと改造したら一瞬でできた．
-
しかしさっき
makebst
いじってたときに気づいたんだけど，
どうも BibTeX にはぢゃーなる省略名を生成する機能がある
(?)
らしく，
そういうオプション指定が可能だった．
つまり BibTeX の中でぢゃーなる省略名を
生成する，
というのがより正しい方法なのかもしれない．
-
技術的な難問はほとんど解決した．
あと一ヶ所ホントに小さい点で修正しないといけないかも
……
必要とあらば，
これは Perl を使えば (これまた一瞬で) 解決するんだけど，
まあ，
しばらく放置してみよう．
-
で，
次の問題は本文の修正というか
……
うーむ，
例によって私の作文がしょぼくてあちこち修正すべき点を
多々指摘していただきました．
-
しかしながら，
なおしのほうはほとんど進捗せぬまま撤退．
1910 研究室発．
買いもの．
1940 帰宅．
晩飯の準備．
晩飯．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0820):
研究室．
食パン．
- 昼 (1330):
研究室お茶部屋．
米麦 0.6 合．
ダイコン・ニンジン・ネギ・豚肉の味噌汁．
- 晩 (2120):
米麦 0.8 合．
ダイコンの味噌汁．
キャベツ・ネギ・ミズナ・ハムのサラダ．
2009 年 02 月 26 日 (木)
2009 年 02 月 27 日 (金)
-
0740 起床．
ねむい．
朝飯．
コーヒー．
洗濯．
0925 自宅発．
晴．
0940 研究室着．
-
「蟻」解説の修正作業のつづき
……
-
午前中いっぱいかかってひととーり修正．
「蟻」編集委員の恒岡さんにいろいろご指摘いただいたおかげで，
かなりマシになったような気がする．
原稿 PDF ファイルをダウンロードできるようにして，
連絡メイル送信．
ふう．
昼飯．
-
次．
とある論文原稿のチェック作業．
うう．
-
WinBUGS 問いあわせメイルが．
はい．
ゐんどーづ上で bugs code でベイズモデルを定義しているファイルへの
パス名 (
C:\Documents and settings\... なアレ)
中に日本語文字が混在してると WinBUGS は動きません．
ゐんばぐすのおかしなところ
参照．
-
「ばぐす」と書いてみて「ばぐず」かもと不安になったのだが
……
bugs は s の前が有声音だから「ばぐす」でいいのかな?
-
今ごろになって先日の統数研公開講座における兼業届をだせ，
とこちらの事務室から．
えーと，
私の時給は
……
と過去メイルを検索してみると 6200 円だった．
-
ふーむ，
世の中には
多重比較法国際会議
なるものがあるのか
……
-
うをー，
やはり「蟻」原稿は最終的には MS わーどファイル化しないといけないのか?
いやはやー
-
……
えーと，
本日午後のメインぎょーむである
論文原稿のチェック，
なんとかひととーり終わった．
勉強にはなった．
しかしコメントはたいへんそう．
窓の外は大雪．
時刻は 1645．
-
深谷君がやってきて 15 分間ほど
glm()
雑談．
-
論文原稿チェックのコメントかき，
けっこうたいへんそうになるかもしれない，
と危惧しつつとりくんでみる．
うう．
-
しかしあまりススまぬまま追いコンの時刻．
1750 研究室発．
雪．
札幌駅パセオ地下 1F の居酒屋「荒波商店」．
忘年会のときはむちゃくちゃに大量に料理がでたけど，
本日は適切な分量だった．
駅地下 + 酒重視の店なので料理すくなめ，
といったところか?
2030 帰宅．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0840):
米麦 0.8 合．
ネギ・卵の炒めもの．
- 昼 (1220):
研究室．
米麦 0.8 合．
味噌にぎりめし．
- 晩 (1810):
追いコン．
2009 年 02 月 28 日 (土)
-
0840 起床．
ちょっと寝すぎだったかな?
朝飯．
コーヒー．
0936 地震，
震度 3．
1000 自宅発．
曇．
1015 研究室着．
-
昨日よんだ論文原稿のコメントかき，
午前中におわらせなくては，
とがんばる．
-
書きあげてメイル送信．
一件落着．
時刻は 1128．
-
これから「
北海道若手生態学研究会」
に参加するべく札幌市南区の真駒内をめざすことに．
1145 研究室発．
1201 地下鉄南北線にのって札幌駅発．
1220 真駒内着．
ここまで 280 円．
-
会場にむかう途中，
国道 453 ぞいの蕎麦屋で昼飯．
おいしい蕎麦でした．
1230 蕎麦屋着の 1245 同発．
1255 研究会会場である北海道青少年会館着．
参加者は 20 名ぐらい．
[北海道青少年会館]
真駒内の住宅地の丘の上にあった．
[研究会会場]
司会の矢澤さんと最初の講演者の寺川さん．
参加者は 20 名ぐらい．
-
1315 より研究会はじまる．
最初は寺川眞理さん
「野外観察と遺伝解析にもとづくニホンザルによる
ヤマモモの種子散布特性の解明」．
屋久島のサル
(フンの中にヤマモモ種子が入っていて種子散布をする)
をずーっとおいかけたり，
種子島のヤマモモの樹木を 12 時間みはったり，
たくさん遺伝子を調べたりといった労力をかけてえられた
データにもとづくハナシ．
内容は多岐にわたるんだけど，
屋久島のヤマモモ好きなサルをおいかけると季節変化がわかり，
5 月後半は早めに熟した特定の個体にみな集まって交代で
木のぼりして食べ，
6 月前半あちこちでヤマモモが熟すと散らばって食べるといった変化
があるそーで．
サルはかなり優秀な散布者で結果的に数百メートルはこんでいた．
種子島は 50 年以上前にサルが絶滅したんだけど，
ヒヨドリ (これも優秀な散布者)
の努力のおかげか，
サル empty forest 化の影響はまだわからず
……
熱帯ではサルは重要な食物だそうで (サル肉屋さんの写真なんかも)，
狩猟による empty forest 化がすすんでいて，
どうなるんでしょうね．
-
次は竹内やよいさん
「東南アジア熱帯樹林 Shore 属の繁殖・更新過程とその空間スケール」．
ランビル 52 ha 調査地の
Shore laxa
(比較的低密度 1.4 本/ha)
15 本の結実とその樹冠下での食害の様子を調べられたもの．
GLMM で解析されているのだが，
NNI なる 2 パラメーター (range と power)
の指数 (あ，樹木サイズが入っていない)
の AIC 最小モデルを探すため range-power
二次元格子上の各点で AIC を評価していく，
という
(「AIC 地図」，圧倒されました)
……
莫大な個数の AIC 評価 → 最小 AIC
をさがすというのはけっこうアヤういんだけど
……
と思ってたけど，
考えてみればあの AIC 地図上ではモデル選択をしているわけではないので，
むしろ deviance 地図と呼ぶべきで，
ならばぼこぼこした対数尤度地形
(おもに range parameter の影響)
上の最尤推定ということか．
推定結果はなるほどなところ
(ハムシのような送粉者であっても熱帯林冠上部なら
風にふかれてけっこう飛ぶかも)
も難しそうなところも
(食害に関して power = 0, range = 200 m ぐらいの「円筒」なカーネル)
あるんだけど
……
やっぱり空間ものはおもしろくて難しく，
代替案がすぐには思いつかない．
-
本日最後は中村さん「ブナを利用する食植性昆虫の緯度勾配」．
主な調査地は黒松内 (ブナ北限)・椎葉 (南限)・秩父 (ブナほ場)．
秩父の結果をみると，
同一環境にもかかわらず産地によって葉の厚さは異なっていそう，
つまり遺伝的にきまっている部分がありそう
(ただし産地が同じでも厚さはさまざま)．
この南厚北薄ならびに
「とにかく厚ければ咀嚼性昆虫にはやられない」
傾向は椎葉・黒松内での観察とも一致．
つまりブナは環境傾度に対応して進化してそうだけど，
咀嚼性昆虫はそうではないかも，
南の厚い葉なんか食べなくても
……
そのあたりの理由は咀嚼性は generalist が多い
(同時に調べられた潜葉性やゴール形成の食植者は specialist が多い)
から?
私の勝手な憶測なんだけど，
ブナアオシャチホコのような咀嚼性な specialist
が卓越するためには，
(小林君いうところの) ブナ分布のフロントではだめ，
とか?
-
このあと，
晩飯たべたり，
夜中まで皆さんといろいろ議論したり
……
-
[今日の運動]
-
[今日の食卓]
- 朝 (0920):
シリアル．
- 昼 (1240):
真駒内の国道 453 ぞいの蕎麦屋
盛朗庵
の大もり 730 円．
おいしかった．
- 晩 (1800):
会館のレストランの洋風定食．
チキンカツ．
