ぎょーむ日誌 2006-07-26
2006 年 07 月 26 日 (水)
-
0800 起床.
朝飯.
コーヒー.
0920 自宅発.
曇.
0935 研究室着.
-
また九大事務とメイルやりとり.
先方はご立腹かも,
とおびえてみたり.
こわい.
-
昨日,
「等エネルギー面上でかぎられた領域しか蛇行できない状態」
図をさんざん見せられたせいか,
だんだん不安になってきたので,
母子里 Gibbs 林冠モデルを
レプリカ交換モンテカルロ法
(パラレルテンパリング法,
Metropolis-coupled MCMC)
に改造する方法を検討してみる.
意外と簡単にできるかも
……
-
プログラム修正とばぐとり試験運転に 3 時間ほど要した.
レプリカ数 8.
じゃ,
MCMC 計算で林冠状態サンプリングの試験運転やらせてみますか.
-
……
と放置しとこうとしたんだけど,
こんなにたくさんの林冠モデルを同時に動かすのははじめてなんで,
マクロなパラメーター設定に苦闘する.
具体的にはレプリカ間の「温度差」だ.
たとえば,
「温度差」をあけすぎると,
すぐにこういうふうに順位がかたまってしまう.
#CANOPY_ID: 1 3 4 5 2 6 7 8
#CANOPY_ID: 3 1 4 2 5 7 8 6
#CANOPY_ID: 3 4 2 1 7 8 5 6
#CANOPY_ID: 3 4 1 7 2 5 6 8
#CANOPY_ID: 3 1 7 2 4 6 8 5
#CANOPY_ID: 3 7 1 4 2 6 5 8
#CANOPY_ID: 3 7 1 2 6 5 4 8
#CANOPY_ID: 7 3 2 6 5 4 1 8
#CANOPY_ID: 7 3 2 5 4 1 8 6
#CANOPY_ID: 7 3 5 2 4 8 6 1
#CANOPY_ID: 7 3 2 5 8 4 1 6
#CANOPY_ID: 7 3 5 8 4 1 6 2
#CANOPY_ID: 7 3 8 4 5 1 2 6
#CANOPY_ID: 7 3 8 4 1 5 6 2
#CANOPY_ID: 7 3 8 4 1 6 5 2
#CANOPY_ID: 7 8 3 4 1 6 2 5
#CANOPY_ID: 7 8 4 3 6 1 2 5
#CANOPY_ID: 7 8 4 3 6 2 1 5
#CANOPY_ID: 7 8 4 6 3 2 1 5
... (以下略)
断続する計算まち時間に雑用あれこれと.
アカマツ原稿はラスヴェガスではなくインド送りになった.
かなーり「温度差」をせばめないといかん,
ということがわかってきた
……
するとレプリカ間でどんどん順番がいれかわる.
みょーなところにトラップされた林冠が他のものに追い抜かれ,
温度の高い領域にほうりこまれてかき混ぜられてまた敗者復活をとげる.
葉群配置に山あり谷あり栄枯盛衰ドラマあり.
これが MCMC 計算術の初目録 (?)
拡張アンサンブル系混合わざ
parallel tempering (並列焼きもどし) 法だ!
-
しかしこういう図を見てると,
統計的な推定にとってエルゴート仮説はどうでもよい,
というか (chain が一本でない場合は) ほとんど関係もないかも,
といまさらながら気づいた.
わりとあたりまえのことで,
ここで MCMC 計算に要求されていることは,
ただ単に
「Gibbs 分布からランダムサンプルせよ,
手段はとわない」
ってコトなんだから.
-
ぱられるの欠点としては時間かかること.
この場合 chain 数は 8 なので単純に言って従来の 8 倍以上の時間がかる.
のこり計算時間はあと 3 時間ちかく必要か,
ということで撤退.
1920 研究室発.
1940 帰宅.
-
2100 北大構内走.
2155 帰宅.
体重 74.2kg.
晩飯.
-
[今日の運動]
-
[今日の食卓]
- 朝 (0820):
米麦 0.6 合.
タマネギ・ニンジン・コンニャク・ワカメの味噌汁.
- 昼 :
食ってない
- 晩 (2210):
米麦 0.8 合.
ニラ卵.
タマネギ・ニンジン・コンニャク・ワカメの味噌汁.