ぎょーむ日誌 2008-02-(21-29)
2008 年 02 月 21 日 (木)
-
0740 起床．
朝飯．
コーヒー．
0845 自宅発．
雪．
ときどき突発的な吹雪．
0900 研究室着．
-
さてさて，
データ解析こんさるの山がちっとも減らない日々だ
……
とりあえず送っていただいた花タイプのデータに関して，
個体ごとの作図をいつものごとく
library(lattice)
の
xyplot()
してみる．
-
「とりあえず図を作りました．作図 R コードはこんなかんじです．
統計モデリングのほうはしばらくおまちください」
といった時間かせぎメイル送信
……
とはいえ，
上のような図を作ってみると統計モデリングしてみるまでもなく，
花のタイプの出現パターンなんかはおおよそわかってしまうのである
(ということで，
上の図はわざと見にくくしてある)．
-
窓の外は吹雪．
今日は気温は比較的たかいんだけど．
-
昨日は放置してしまったヒヨコ問題にとりくむ．
大学院生の書いた R & BUGS コードを読むのは楽しいけど，
なんといいますか，
「あー，こらこらいたずらしちゃダメだよ」
保父さん的寛容を発揮せねばならぬ．
そして経験の長短にかかわらず hacker 的な資質のあるヒトってのは
良くも悪くも暴れん坊みたいなもので
……
-
本日最初の件に関して，
また質問メイルいただく
……
そう，
まだ慣れてないヒトたちにとって呪われソフトウェア
ゑくせるに幽閉されてるデータを CSV 出力する，
というのはなかなか難問だったりするわけですよね
……
-
ヒヨコ作図プログラムとか整理できたように思えたので，
昼飯調達の旅にでる．
時刻はすでに 1355．
吹雪はやんでるのでいまのうちに．
今日は北大生協ではなく北 11 西 4 すなわち地環研前のサンクスに行ってみる．
なんだか A 棟 8F お茶部屋は今宵の「追いだしコンパ」
料理準備のため戦場のごときありさまのようなので，
こわくなって研究室で昼飯たべる．
-
ヒヨコなメイル書こうとすると
……
地環研 server 雑用．
あらー，
先日 (15 日) の
chown -R topeditor
をしくじっていたか．
物質環境科学全体のディレクトリがこんなところに隠されていたとは
……
-
直せるところは直して連絡．
とりあえずは一件落着．
しかし来月もしくは再来月にめんどうなことになりそうな予感．
-
で，
ひよこメイルを送信した瞬間に，
水中栄養塩類メイルがくる，
というのは偶然なのだろうか?
この部屋のどこかに監視カメラがついてるような気がしてきた
……
-
どんどん横からわりこみが．
A802 院生部屋のプリンター不調とのことで修理にきてもらっていたんだが
……
原因は CANON LBP2810 video 基板 (プリンターのメイン基板をこう呼ぶらしい)
とゆーことのようで，
修理みつもりは最低でも 10 万円．
買いかえたほうがマシかも，
というかんじなんだけど，
雪野さん発案で大学事務から CANON と直接交渉
(そもそも不良品では? 疑惑の追及)
するようなんで
……
まあ，
とりあえずその方向で．
-
現状ではモノクロ印刷はとくに問題なし
(やや紙づまりぎみ?)．
カラー印刷に問題あり
(とくに Mac や Linux 経由)．
そもそも LBP2810 は同社でもとらぶる続きのマシン，
ということになってるらしい．
-
だんだんばててきた
……
-
で，
天塩川の水中問題にもどる．
ああ，
われながらうかつなことに降水量の影響とかどわすれしてましたよ．
この要因はめんどうきわまりないなぁ
……
国土交通省の
水文水質データベイス
なんかだと
流量の時間変化はわかるから間接的に降水の影響はわかるんではないかな
……
ってそんなややこしいモデリングはやりたくないなぁ．
-
どうでもいいケド，
上記データベイス，
ひょっとして MSIE ぢゃないと動作しない，
とか?
嗚呼，
おやくしょでーたべいす
……
-
まあ，
降水量とかつまるところ短期的な流量変化は
(土壌栄養塩の溶脱速度は地下流量に比例するとか仮定することにして)
random effects としてあつかうしかないでしょう，
といったお返事メイルかき．
I hope some effects of land cover/use can be extracted by using
such simplified model from the massive flood of random effects.
-
1800 より A804 室の追いだしコンパに参加．
2105 研究室発．
2120 帰宅．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0800):
ダイコン・豆腐のコンソメスープ．
ヨーグルト．
リンゴ．
- 昼 (1415):
研究室．
玄米パン．
- 晩 (1830):
[追いコン]
今回も大学院生尽力の料理が並んだ．
ただし揚げものは最近 (?)
北 24 条にできたというスーパー
アークス
で買ってきたもの．
2008 年 02 月 22 日 (金)
-
0740 起床．
朝飯．
コーヒー．
ちょっと旅じたく．
0905 自宅発．
曇．
0920 研究室着．
-
またデータ解析こんさるメイルかき．
-
そしてトドマツ母子距離モデルの修正原稿のチェック．
いくつかのポイントチェックして 1055 ごろとりあえず返信．
-
同じトドマツ修正に関して追加できた問題への対応．
この観測データでは
-
たくさん子供がいるならそれらは近距離にいる
-
遠距離に子供がいるなら子供数は少ない
(注: 子供数が少ないからといって遠距離とはかぎらない)
といったパターンがあるんだけど，
査読者 2 が「それってホントか」と疑っており，
subject editor は「子までの平均距離の分散に関する検定やったら?」
という提案．
-
先日，
後藤さんがこちらにこられたときはめんどうだったので，
このあたりまるごと削除したら，
といったハナシになったんだけど，
「まるごと削除」
すると他にも影響がでるのでかえってめんどう．
-
そこでここで「けんてー」を持ちだすのは不本意であり，
えーい，
本来ならば子供数と散布距離の共通事前分布を二変量正規分布にして，
そこからこの二形質の相関の事後分布を推定すべきであったのだが
……
しかしこのトドマツ計算やってたころにはそのやりかたがわかんなかったんだよねえ．
-
そして現状で計算やりなおすと結果がそれなりに変わるだろうから
……
やむをえず
subject editor 殿の提案にしたがって
vmd <- var(mean(distance))
について調べてみるとことに．
考えればわかるようにこれは randomization (permutation)
test するほかない．
-
R
だとこのていどのプログラムは試行錯誤しつつ 15 分ほどで書ける．
そして 1000 回の sampling は数秒で終了．
-
で，
結果は観察された
vmd は 1963.9
だったんだけど，
母樹を無作為化した場合の vmd
95% 区間は (610.8, 1855.9) となり観察された値は vmd が大きすぎる
……
-
つまり全 82 母樹 (ただしここでいう母樹は一個体以上の子供をもつもの)
の中で 42 母樹が複数の log (倒木) 上に子供をもつんだけど，
その複数の log がランダムに分布しているのではなく，
母樹に近いところにあつまっているのでこのように
母樹 - log 間距離の分散が大きくなった，
と結論すればよいわけだ．
-
といった方向にて一件落着! メイルを書くとすでに出発すべき時間．
-
ということで宮城県の山の中にある
東北大学大学院農学研究科附属複合生態フィールド教育研究センターにむけて，
1143 研究室発．
むこうでは福岡大会のベイズ統計企画集会の事前勉強会，
みたいなのが開催されるんだよね．
1155 JR 札幌駅から快速エアポートにのって出発．
1232 新千歳空港着．
雪．
[B-767-300]
離陸準備中．
氷ついた主翼に湯をふきかけてるところ．
-
機内でくう昼飯調達して 1250 ANA 724 便 (また) B-767-300 に搭乗．
1259 taxing 開始．
離陸．
[機内]
やはりこの時刻は乗客すくなく，
しかも事前座席指定でいちばんうしろの席をとったので
「ひろびろー」といったかんぢで，
とてものんびりできた．
-
機中では Bayes 本などながめる．
昼飯．
1400 仙台空港に着陸．
ああ，
先週のつくばと同じくここも二ヶ月ぐらい季節がススんでいる．
1424 JR 空港アクセス線にのって同発．
[仙台空港アクセス線]
初めて乗った．
二両編成．
-
1443 仙台駅着．
なんとここで一時間ぐらい待たねばならぬようだ．
いつまでたっても覚えられぬ東北地方内の新幹線名，
本日は「やまびこ」が東北新幹線の各駅停車列車だと理解できた
(「こだま」にかけてるわけだな)．
そしてこれは 2 時間に一本ぐらいしかない，
ということか
……
-
一駅北の古川までいくだけなんだけどね．
えーい，
いっそのこと在来線で古川まで移動するか，
と思ったのだがネットに接続されてない私は無力というほかなく，
時刻表を買ったり立ち読みする気力もなく，
自由席特急券を買ってしまう．
840 円．
どうせ最終到達駅・川渡温泉駅における集合時刻 1647
より早く行っても意味はないわけで
……
-
ということで，
新幹線まちあいフロアでのんびりすることに．
その前に乗車券の清算しといた．
空港の切符自販機では古川まで 1350 円，
しか買えなかったのでここで川渡温泉までの追加料金 500 円支払う．
合計 1850 円．
-
1542 仙台発．
16 分後，
1556 古川着．
ここで JR 陸羽東線・鳴子温泉行普通にのって
1611 同発．
車内で関東方面からきたメンバー，
すなわち角谷さん・深澤さん・石濱さん・山道さん・島谷さん・竹中さんと合流．
[1647 川渡温泉着]
うーん，
一度きたことあるんだけど記憶がない
……
駅には富田君たちが車でむかえにきてくれていた．
仙台とかぜんぜん積雪なかったけど，
宮城・山形県境にちかいここでは 30-50 cm ぐらいつもっている．
[1700 東北大「川渡農場」着]
このへんはなんとなくおぼえているなぁ．
[院生部屋で]
なぜか竹中さんが
今晩の料理づくり．
手すきのものは鳴子温泉へ．
[鳴子温泉内の早稲田温泉]
こういう奇妙な建物で．
湯の温度がやや高いこともあって，
私は 30 分ほどで出たんだけど，
他の人々は一時間ほど湯につかりながら議論してた．
[晩ごはん]
今晩は院生部屋で．
こちらの院生の皆さんといっしょに．
2400 ごろまでハナシがつづく．
とあるエラい先生が
「ブナの北限はどのように決まっているか?
それは恐竜絶滅によって林内撹乱が……」
といったトンでも話を一般むけ講演でしたとかしなかったとか．
いやはや．
[寝室]
学生実習なんかで使われる寝室．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0750):
タマネギ・セロリ・ニンジン・ソーセイジのスープ．
シリアル．
- 昼 (1330):
B-767-300 機中．
空港でかったもりもとサンドイッチ．
ANA コンソメスープ．
- 晩 (1900):
二種類の鍋物，
カレイの煮魚，
イカ刺身などごちそうになりました．
2008 年 02 月 23 日 (土)
-
0650 起床．
東北大川渡農場ベイズ合宿二日目．
[川渡農場]
なぜか目のまえにシカの飼育檻が．
本日は朝は曇で昼すぎに大雨．
それが吹雪にかわったりときどき晴れたり，
最後はまた吹雪，
という山地らしいかわりやすい天気．
[川渡農場の模型]
山の上にも放牧地をもつ広大な農場．
低地には実験用水田なんかもあるそうで．
-
皆さんは朝風呂へ．
朝から風呂にいく元気のない私は院生部屋で朝飯の準備でも
……
と思ったら，
カギがかかってて入れない!
院生部屋のカギ，
というものはえてしてそのへんに隠してあるモノで
……
とヒトけのない庁舎内をうろうろと探しまわったけどみつからず．
とうとう朝風呂から帰ってきた富田君にうかがってみると
……
そこに隠してあったのか，
とがっくり．
-
ともかく院生部屋で朝飯の準備．
コーヒー．
朝飯．
-
予定よりややおくれて本日のセッション開始．
[ベイズ合宿ベイズ勉強会]
ベイズ企画集会で発表される内容に関して，
一件ずつ一時間以上かけての検討あれこれ．
いやー，
なんだかまにあっくな世界で
……
-
昼前に農場をでて全員車にのせてもらって JR 川渡駅へ．
本日合流の鈴木牧さん・高田まゆらさんをおむかえする．
そのまま国道ぞいのドライブインで昼飯．
-
また農場にもどって日没までベイズ議論．
-
1750 また車にのせていただいて農場発．
当地では移動のたびにこちらの院生のみなさんのおせわになっております．
昼と同じ場所で晩飯．
温泉．
-
2100 ごろから 2230 ごろまで，
本日最後のベイズ議論．
-
で，
そこからさらに夜おそくまで飲み会が継続する．
合気道弐段の高田さん・山道君による実技披露
(演武)
があったり
……
-
[今日の運動]
-
[今日の食卓]
- 朝 (0830):
米麦 0.6 合．
昨晩の鍋ののこりもので作った味噌汁．
豆腐に当地では「だし」とよばれるソース
(漬けものを細かくきざんでペイスト状にしたもの)
をかけたもの．
- 昼 (1230):
国道ぞいの「ドライブインおーとり」．
混んでたので観光バス乗務員用の隠し座敷でいただく．
カツカレー 780 円．
- 晩 (1840):
晩飯も国道ぞいの「ドライブインおーとり」．
[カキ鍋定食]
980 円．
味噌鍋じたてでした．
2008 年 02 月 24 日 (日)
-
0720 起床．
川渡ベイズ合宿 3 日目．
早朝まで飲んでた皆さん，
まだ寝ておられる．
ということで図書室へ．
昨日朝はうまくいかなかった川渡農場ネットへの侵入に再びとりくむ
……
なんとかうまくいった．
うーむ，
メイルがたまっている
……
[川渡農場庁舎中庭]
左側のバラックが院生研究室，
中庭をはさんで右側がベイズ合宿議論に使われてる会議室．
必ずしも若いとは言えないヒトたちがおきてきたので，
朝の作業．
竹中さんとこのへんの除雪作業を楽しんでみたり．
-
院生部屋でコーヒー．
朝飯．
寝室で荷物の撤収．
みなさん，
まだおきてこない．
[川渡農場ネコ]
この他にも二匹ほどこのへんをうろうろしてる．
寒そう．
[廊下のナゾのハコ]
なんなのだろう?
と思っていたのだけど室内暖房用の灯油が備蓄されてた．
-
会議室でうだうだする．
皆さん，
なかなかおきてこない．
-
1000 すぎから本日のべいづ議論．
ここまでのまとめ &
企画集会全体の構成なんかについて
……
-
私はこっそり某雄蕊数問題のデータ解析にとりくんでいたり
……
なンと!
花位置に依存してるのか．
そしてこの植物形態に関するデータ，
ある部分においてまったく個体差というものがなく
(データ数がそれなりであるにもかかわらず)
完全な二項分布になっている!
じつに教科書的な少数遺伝子によるメンデル分離の実例か?
……
とすごくびっくりした次の瞬間，
室内の議論にもどって，
「えー，
そもそも野外調査・生物あいての実験では
random effects がない，
ということはありえないわけで
……」
とかコメントしてる私はやはりいいかげんな研究者なのだろうか
……
-
1300 ごろまでクタびれたヒトたち (もちろん私も)
による熱心な議論がだらだらとつづき，
まあなんとなく今回はおひらき，
とゆーことで．
-
昼飯．
そのまま竹中さん・牧さん・高田まゆらさんは川渡温泉駅まで
送っていただく．
とりあえず仙台まで 1450 円，
の切符を買う．
-
ここから誰もが予想しなかった「先の見えぬ苦難の旅」が始まる．
そう，
川渡温泉のような外界から隔離された世界でわれわれがべいづしてるあいだに，
日本列島はとんでもない暴風雨と降雪 (これは北海道のみ?)
にみまわれていて
……
-
まず陸羽東線が遅れる．
1411 発のはずが 1419 に．
JR 古川での新幹線のりつぎがまにあうのかしらん
……
などと考えていたけど心配無用のことだった．
東北新幹線の運航はもっとめちゃくちゃになっていたからである．
-
1454 古川着．
1457 発の「やまびこ」は
……
1520 発．
しかも乗車率 200% ぐらいのすごいこみかたに．
[15:37 やまびこが JR 仙台駅着]
東京方面にもどられる皆さんはこのこみこみ新幹線であと
2 時間も，
と心配しつつ見送る
……
後続の列車はもっと遅れてるようなので，
このまま乗りつづける他ないのである．
とはいえ，
当方にもこのあとこれまた予想外の苦難がまちかまえていたのであった
……
-
JR の仙台空港アクセス線 1538 発にはまにあわないなあ，
と思ってたら発車を遅らせてまっていてくれた．
1542 仙台発．
各駅停車．
1605 仙台空港発．
川渡温泉からここまで乗車券 1850 円．
-
出発まで一時間あるのでのんびりと買いもの．
保安ゲイト通過すると，
何やらぶきみな放送が
……
ANA 727 便，
札幌千歳空港が雪のため，
着陸不能の場合は旭川空港，
または仙台空港へのひきかえし
……
定刻より 25 分おくれの出発?!
-
で，
実際には定刻より 35 分おくれの 1740 taxing 開始．
ANA 727 便 A-320．
搭乗率 100% のちょー満員．
なんで日曜日夕方の仙台 → 札幌便がこんなにこむんだろうね．
-
そして離陸したのは良いものの津軽海峡南側に終わりなき旋回飛行を命じられる．
新千歳空港の滑走路の除雪がどーのこーので
……
-
蛇足ながらなぜ新千歳空港 (もとは海軍航空隊の飛行場)
が人口密集地の札幌ちかくではなく
ずっと南の太平洋側にもうけられているかといえば，
日本海側の豪雪をさけるためである．
その千歳で雪とはねえ
……
春一番の南風と関係あるのかしらん?
-
1920 着陸．
予定では 75 分間のはずが (このあとの飛行場内移動もふくめて)
120 分ぐらいになってしまったなぁ
……
と思ってたら空港内の駐機場も満員御礼状態で，
乗客をおろせる場所をもとめて乗機はこれまた終わりなき空港内放浪を命じられる．
[19:46 新千歳空港で飛行機からおりる]
で，
飛行機からおりても空港内バスの旅になる，
と．
こんなかんじで積雪してます．
[19:52 新千歳空港内のバスの旅]
たしかに千歳ではありえない積雪量だなぁ
……
[19:56 新千歳空港内のバスの旅]
空港内バスの旅はなかなか時間がかかる．
というのも，
どこかで U ターンしないと乗客をおろせないんだけど，
そのターンできる場所がかなり限定されているので，
無駄な移動経路が長いのである．
ということで，
着陸から 40 分かかって到着ターミナルについた．
[20:05 新千歳空港，ありえない人ごみ]
ANA 到着口ふきん．
地下にある JR 駅から地上にとどく長蛇の列が形成されてて
(推定列長は 300 m 以上)，
こりゃダメだ状態．
高速バスのカウンター前にも
……
[20:28 新千歳空港 JAL 到着口]
ふと思いついて ANA から JAL の搭乗口のほうに移動してみると，
電光掲示板に示されてるとーりのきなみ欠航・遅延で
到着した乗客数がぜんぜん少なく状況が把握しやすい．
しかし最悪なのは最大の輸送力をもつ JR 北海道で，
(ふだんは何時発どこ行きと掲示される)
巨大な電光掲示板があるにもかかわらず「調整中」の三文字しか
表示してないんだよね．
[20:56 新千歳空港，高速バスカウンター前行列]
周りの人々のハナシによると，
やはり地下の JR 新千歳空港には数千人の人ごみ状態のようで．
JAL 口の中央バスカウンター前 のほうが短いのでならんでみること．
中央バスも臨時バスを増発して対応してるようだが，
やはり情報提示がダメなので状況がわからない．
-
一時間ほどならんでいたら札幌駅いきのバス切符 1000 円が買えて，
バスにのることができた．
たった一時間
ならんだだけで買えたのはまったくの幸運というべきだろう．
なぜなら私の後ろにはさらに待つべき人々が並んでいて，
行列の終端まで到達するにはあと数時間を要する，
と予測されるのだから
……
-
2121 新千歳空港発．
札幌につながる高速道路にのってはみたものの，
時速 50 km しかだせない状況
……
-
今回えられた個人的な教訓として，
かくのごとくどこに行っても情報は不明，
どこにいってもヒトが行列を作っている，
という状況において，
自分でもわりと意外だったんだけど
-
けっこうあわてる
……
あわててるつもりはないけど，
うっかりミスが増える
-
けっこういらいらする
……
「えーい，こんなとこでいらいらしても無駄無駄……」
と自制できるていどではあるけれど
やっぱり混乱した状況で情報がないと，
人間は挙動不審になるもんだねえ．
今回ばかりはけーたい電話をけーたいしたいと思った．
-
で，
高速バスはゆるゆると北進して
……
あらら，
北広島のインターチェンジでおりてしまった．
とすると，
ここから行程 20 km の札幌市内雪上移動が始まるのか
……
-
平均時速 20 km ぐらいのバス旅がイヤになってしまったので，
札幌駅までいくのはあきらめて，
札幌ドームちかくの福住バスターミナルで降りる．
時刻は 2237．
-
ここは地下鉄東豊線の駅があって
……
2250 地下鉄東豊線福住発．
さすがに地下鉄は地上の状況に影響をうけずに運航されていて，
15 km ぐらいの札幌市内南北移動は 15 分で終了．
2305 自宅すぐちかくの北 13 東着．
2310 帰宅．
ぐったり．
空腹．
晩飯．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0900):
院生部屋でこんびにおにぎり．
- 昼 (1220):
国道ぞいの「ドライブインおーとり」．
中華そば 500 円．
- 晩 (2350):
スパゲッティー．
レトルトパウチドソース．
2008 年 02 月 25 日 (月)
-
0855 起床．
うー
……
昨日のぐったりはかなり解消．
まあ，
たしか本日は川渡出張 (出張ということにしてました)
の代休，
すなわち書類上の休日のはずだよな
……
朝飯．
コーヒー．
1030 自宅発．
晴．
札幌市内もけっこう積雪してるなぁ．
1045 研究室着．
-
昨日の交通混乱のせいで，
北大の入試は一日おくれることに．
なんだか
たいへん
なことになっているなぁ
……
私の印象だけど，
北大はなけなしの事務処理能力を入試課に集中しているので，
かくのごとく迅速な対応ができている，
というような．
しかしネットにアクセスできない受験生はなかなかたいへんなことに
……
-
そしてなぜかしらこの変更が，
私の今週末の出張日程に影響するかもしれない．
そのあたりちょっと事務員の三浦さんに教えていただく．
なるほど，
事後的に「出張いけませんでした」
と返金も可能か
……
-
いろいろとメイルかき，
とか．
-
A 棟 8F の平尾君と大雪山 OTC データ解析の方法の相談．
まあ，
平尾君なら概要なんかを説明すればあとは何でもやってくれるわけで
……
-
昼飯．
-
宿題あれこれがあまりにもあるので何から手をつけたものやら，
と途方にくれる．
ああ，
生態学会福岡大会初日 (3/14 金)
まであと 18 日か．
えーと
……
最低限のやることリスト:
-
一般講演: アリポスターの準備
-
R プログラミング自由集会: 前口上的な部分
-
ベイズ企画集会: コメント (まあ 1-2 分で終わるような)
-
山地生態学企画集会: 占部さんの計算したうけ，PP-ZP-Fish 系の再計算
-
旅費関係というか，
航空券関係を旅費だしてくれる九大理学部等財務企画係に送る作業．
-
またこんさるメイル書き．
さすがに古い順から対応してみる
……
-
1600 ごろから次のデータ解析こんさる．
Nina さんと天塩川の
……
本日は栄養塩類濃度の土地利用のモデリング．
-
しかしデータが MS Office 2007 の
*.xlsx
なやつだったので当方では解毒できず．
けっきょく，
Nina さんはいったん研究室にもどる．
あとで露崎さんに
「Word/Excel/PowerPoint 2007 ファイル形式用 Microsoft Office 互換機能パック」
(その
アヤしげな URL)
なるモノを教えていただいた．
まち時間に A801 院生部屋の共用 PC
で実験してみるとたしかにこれで open
できるので CSV 変換とかで解毒できる．
ただし MS の説明そのままでは動かず，
ゐんどーづ 2000 + おひす 2003 という組みあわせでは，
例によって例のごとく理不尽なる「OS 再起動」なしでは機能しないようだ．
-
Nina さんもどってくる．
データあれこれながめるのに使う．
私は水文学・水質学とかまったく知らないんだけど，
常識的に考えて水量 (これはある河川断面を
一秒間に通過する水の体積で表現される)
が重要，
ということで水量
×
栄養塩類濃度，
つまり栄養塩類質量が土地利用によってどう変わるか，
といったモデリングを検討する．
-
そのためにデータのみかたなどあれこれとうかがう．
国土交通省は
水文水質データベイス
とやらはゐんどーづ & MSIE の組みあわせでも
データをうまく取得できないらしい，
といった「何のためのデータ公開なんだか」
といったハナシをうかがう．
-
いずれにせよ雪どけ期間のモデリングはやりたくない
……
他の月とぜんぜんちがうんだよねえ．
などと「とりあえず版」モデルにむけてモデリング対象のしぼりこみ，
とかで本日は終了．
今日・明日中ぐらいになんとかしないといけないような
……
-
先日のぎょーむ日誌に合気道「二段」とか書いてたんだけど，
これは「弐段」が正しいんですよと高田まゆら合気道弐段からご指摘いただき，
あわてて修正．
ふーむ，
合気会審査要項
なるものも公開されてるな
……
-
こんさるメイルかき，
というか夜は一般社会とちょっと接点のあるメイルかき．
なかなか仕事がススまなくて，
といった時間かせぎメイルにすぎないんだが．
-
2255 研究室発．
2310 帰宅．
今日もおそい晩飯．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0910):
シリアル．
- 昼 (1430):
研究室お茶部屋．
北 12 条生協の納豆巻．
リンゴ．
- 晩 (2400):
米麦 0.8 合．
ダイコン・ニンジン・コンブ・ワカメ・すりみのスープ．
リンゴ．
2008 年 02 月 26 日 (火)
-
0910 起床．
ねむい．
生活周期がずれてる
……
朝飯．
コーヒー．
1010 自宅発．
曇．
1025 研究室着．
-
1030 より
研究室セミナー，
本日は吉田君で葉寿命の数理モデル．
なぜかしら生涯光合成速度の最大化ではなく生涯純生産の最大化を目的にした
計算になっているので，
挙動がいろいろとヘンだ．
本人もよく理解してなかったようなので，
そのあたりを整理してみる．
-
昼飯．
-
1330 より 1.5 時間ほど M1 木村君のデータ解析こんさる．
火入れのある仏沼の種数が何で決まるかとか，
ヨシ・アオスゲ・イネ科草本の被度だの草高だの．
問題は例によって被度 (coverage)
なる「目でみて決める量」なんだが
……
ハナシを聞いてみると (目視で) 25 個の sub quadrat
に分割して，
ということなので 25 個中のいくつといったデータに換算できそうだ．
-
窓の外はまた吹雪．
-
金曜日の帯広出張は残念ながら取り消さざるをえない状況のようだ．
紺野さんにおわびメイルだす．
JR 北海道の予約キャンセル．
これはなぜか 530 円もとられてしまった．
帯広宿泊とりけしは無料だった
あとは大学まわりの出張とりけし，
なんだけどこれはまた後日に．
-
院生になんだかいきなり作文推敲を命じられる．
「月末に出さないといけない国際学会の要旨です」
「え」
「久保さん共著者ですから」
「いつのまに」
「要旨をなおしてください」
ということで解読・修正に苦闘することに．
-
日没時刻ごろから次第にばててきたので
……
撤退．
1850 研究室発．
1915 帰宅．
晩飯．
-
生態学会福岡大会初日 (3/14 金) まであと 17 日．
なにも準備がススんでいない．
-
いろいろと下うけデータ解析が堆積してるなぁ
……
ととりくんでみようとしたんだけど，
やはりススまず．
ばててるな．
-
さっさと寝なくては，
と思ったらメイルがきて仕事遅延ができない状況とわかる．
ということで，
夜中に解析結果の報告かき．
-
そして明日から三日間，
北大某所に隔離監禁されてなぞの下うけぎょーむ．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0920):
ヨーグルト．
リンゴ．
- 昼 (1300):
研究室お茶部屋．
米麦 0.6 合．
ダイコン・ニンジン・コンブ・ワカメ・すりみのスープ．
- 晩 (2000):
米麦 0.8 合．
ダイコン・ニンジン・コンブ・ワカメ・すりみのスープ．
リンゴ．
2008 年 02 月 27 日 (水)
-
0730 起床．
朝飯．
コーヒー．
0840 自宅発．
雪．
0855 北大内某所着．
-
1700 までナゾの下うけ労働．
-
1710 研究室着．
けっこう緊急性の高い
「お願い!」
メイルが二通も
……
-
とりあえず一件目に対処．
-
空腹になってきたので撤退．
1850 研究室発．
雪．
1900 帰宅．
晩飯の準備．
晩飯．
-
なかなか二件目に着手できず．
-
生態学会福岡大会初日 (3/14 金) まであと 16 日．
-
2500 すぎまでネバってはみたものの，
やはりなかなかよい結果はえられないわけで
……
-
そして明日もナゾぎょーむつづく．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0750):
ヨーグルト．
リンゴ．
- 昼 (1230):
北大生協北部食堂．
ライス M，
豚汁 S，
モロヘイヤオクラ，
豆腐ビビンバ．
647 キロカロリー．
367 円．
- 晩 (2120):
米麦 0.8 合．
ホッケのさつまあげ．
2008 年 02 月 28 日 (木)
-
0740 起床．
ねむい
……
朝飯．
コーヒー．
洗濯．
0845 自宅発．
また雪．
そしてまた昨晩もかなり積雪したようで．
0900 ナゾの下うけ労働会場着．
-
午前中でめでたく終了．
終了したので書くけど，
北大の学部の二次試験の採点作業なのであった．
-
のんびり昼飯．
1335 研究室着．
-
昨晩真夜中すぎのデータ解析でいささかばてぎみ．
その原因となった天塩川問題を少しだけススめる．
-
マイアミ大学の，
なんといいますか私からしますと
「甲山さん悪友」
とでもいうべき
(ただし私より若い)
Matt がくる．
今日・明日は以前と同じくこの部屋を拠点にそのへんふらふらとする予定．
なぜかしら，
さっそく彼の熱帯林データを階層ベイズモデル化する工夫に関する議論．
いわゆる種特性なるよくわからぬ値を fixed effects
的な説明変数としてもちいるのではなく，
多変量正規分布な事前分布をつかってハイパーパラメーターたる
分散共分散行列と同時に target 形質と種特性の同時分布みたいなもの
(ってふつーの事後分布なんだが)
を推定すればよろしかろう，
といったらいやはやすぐさま理解してしまって，
参考文献おしえてくれー状態に．
-
動物生態の松田君のとびこみデータ解析こんさる．
お茶部屋で．
うーむ，
テングザルはすごく面白いんだけど，
かくのごとく長期間にわたって群に密着して得られた時間・空間・社会構造のあるデータの
統計モデリングはすごーくたいへんなんだってば
……
たとえば，
ハーレムの重層性なるものをむりやりというか現時点での定義にそって
統計モデル的に考えると複数群の位置の時系列データ，
その相互作用つまりつるみ相手・つるみぐあいを推定せよってことになるんだよねえ
……
-
経済学分野の統計もでらーなヒトたちも
人間の購買行動のベイズモデリングばっかりやってないで，
こういうおさるさんモデリングをやってくれると助かるのだが
……
どっちも霊長類ぢゃありませんか．
とはいえ，
じつは上の重層性うんぬんはミクロなハナシというより，
むしろ日本円と他の通貨の相場変動のデータ解析
(相関のある複数時系列データ解析)
なんかに近いんだろうけど．
-
ともかく，
松田君の莫大なるデータを小分けにして論文として出版していく計画を拝聴し，
最初の「寝場所えらび」論文の投稿準備ができたら，
おさるさん寝場所 randomization な解析でもやろうか
……
とまた当方の時間かせぎを含む方向への誘導に成功する．
-
A802 「カラー印刷不能」状態におちいった彩色印刷機
CANON LBP-2810 問題
(なおしたけりゃマザーボード換装が必要)
に関する検討．
-
まずは (最高権力者?) 雪野さんからいろいろ教えていただく:
2002 年リリースの LBP-2810 の修理に 10 万円以上を投入するのは
わりがあわなさそう，
代替新型機の入手の予算規模・購入可能時期なんかについて
-
現在，
A 棟 7-8 F 闇ネット内で彩色 & 両面印刷が
唯一可能になってしまった EPSON ELP S6500 (7F 院生部屋)
に関してユーザーたちの意見をうかがう
……
なぜかしら Mac よりむしろゐんどーづからの印刷出力が遅い?
-
で，
EPSON サイトを調べてみたんだが
……
S6500 とトナーカートリッジを共有できる廉価マシンはナシ，
と確認
-
だったらまた CANON でもいいかなー
……
というあたりで本日はチカラつきて終了．
-
このプリンター調査のついでに，
小林君とたちばなし的ブナ林食植者昆虫の統計モデリングこんさる．
さらにさらについでながら，
なぜかしら某「恐竜絶滅とブナの北限」問題に関する，
よりフカい事情なんかも．
とりあえず，
ブナの起源は恐竜時代の末期といった
ホントかどうかよくわからぬ説があって
……
-
生態学会の和文誌に次の
自由集会
について何か書きませんかねメイルに関する質問メイルかき．
-
そして時間かせぎのきかぬ天塩川問題のつづき
……
だがすでに Nina さんの deadline にまにあわないよーなんで，
現状のまとめと今後の展望に関するやや長いメイルかき．
-
といったことで，
くたびれきったので
2015 研究室発．
2030 帰宅．
晩飯．
-
今晩もまた別のデータ解析
……
と思ったのだけど，
残念ながらこれは明日にまわすことに．
-
生態学会福岡大会初日 (3/14 金) まであと 15 日．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0800):
ヨーグルト．
リンゴ．
- 昼 (1200):
北大生協北部食堂．
豚汁．
サバの味噌煮．
コマツナゴマあえ．
1048 kcal．
514 円．
- 晩 (2110):
米麦 0.8 合．
ネギ・シイタケ・卵炒飯．
2008 年 02 月 29 日 (金)
-
0740 起床．
朝飯．
コーヒー．
0850 自宅発．
晴．
今日は北低南高の気圧配置で，
気温が高い．
+3°C
ぐらい．
0905 研究室着．
-
高田まゆらさんのクモ vs カメムシの水田内闘争の階層ベイズモデル，
WinBUGS + R2WinBUGS
による実装例の構築，
という宿題にとりくむ．
全データのごくごく一部しか使ってないのに，
BUGS コードを書いてみるとけっこうややこしいモデルになってしまった．
-
で 1100 からその作業から引きはがされて，
甲山さん・Matt・矢澤さんの熱帯会議になぜかまきこまれてしまう
……
うーむ，
よくもまあ次から次へとアヤしげな
……
-
甲山さん部屋の横の A702 で昼飯などくいつつ 1400 ごろまで継続．
うう．
-
クモ vs カメムシ宿題のほうはだいたいできたかな
……
(BUGS code,
事後分布表)
-
で，
あれこれ一式をダウンロードできるところに置いて，
WinBUGS のインストールからクモ vs カメムシモデルの注釈にいたる説明メイルかき．
この「説明メイルかき」って何であれ時間がかかるんだよねえ
……
-
で 1600 ごろに何とかかんとか一件落着して，
また別のメイルかきあれこれ．
-
1700 から
Trendy セミナー，
本日は Matt で熱帯林伐採方法に関する経済学よりの理論的なハナシ．
けっきょく人間活動の意志決定ってのは
こういうおおざっぱなモデルによることになるのかなぁ，
熱帯林の詳細なデータもないわけだし
……
というかんぢだ．
そうそう，
3 月から北大に着任する佐竹暁子さんも参加してましたよ．
-
1850 研究室発．
Matt 懇親会，
Asty の居酒屋で Matt 懇親会．
2105 帰宅．
-
生態学会福岡大会初日 (3/14 金) まであと 14 日．
-
えーと，
明日から大雪山ふもとで開催される
進化植物学研究会
に参加するので，
日曜日の夜まで連絡とれない予定です．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0750):
米麦 0.5 合．
海藻スープ．
- 昼 (1230):
研究室．
北大生協にぎり飯．
- 晩 (1920):
地鶏と旬魚 手づくり豆冨 ととと．
豆腐がうまかった．
