「ぎょーむ日誌」目次に戻る
|
KuboWeb top に戻る
|
twilog
|
atom
ぎょーむ日誌 2008-02-22
苦情・お叱りは, たいへんお手数かけて恐縮ですが, 久保 (
kubo@ees.hokudai.ac.jp
) までお知らせください.
本日
(
kubolog20080222
) |
次の日
|
1 日前
|
7 日前
|
31 日前
|
365 日前
|
top
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): 二種類の鍋物, カレイの煮魚, イカ刺身などごちそうになりました.
本日
(
kubolog20080222
) |
次の日
|
1 日前
|
7 日前
|
31 日前
|
365 日前
|
top
KuboLog
|
KuboWeb