-
0750 起床．
朝飯．
コーヒー．
0920 自宅発．
曇．
0935 研究室着．
-
また ECOAS 八甲田湿原植物データ解析のつづき．
さすがに昨日あれだけいろいろ調べまわると，
何をどうやればいいのか，
少しは見当がつくようになってきた．
とはいえ，
やはり試行錯誤をつみかさねていくことでしか，
この統計モデルは改善されないわけだが．
-
計算まち時間に ECOAS 動物プランクトンに関する解析の結果をみなおす．
「脱」 ECOAS のためには，
-
得られた事後分布の表示が改善できないか検討する
-
この推定結果をもちいて，
何か ECOAS 的な予測というか憶測計算を示してみる
といったことが必要なのだろう．
-
八甲田湿原植物の試行錯誤．
競争力に上限・下限をもつような関数型にするとずいぶんとマシにはなるけど，
(昨日かいたような)
「談合で実現する種間競争力」
になってしまうところがあるねえ．
ここでは種間競争が完全に相対評価になってしまってるからなぁ
……
-
とはいえ，
絶対評価のがちんこ勝負にすると
「種間競争キビしすぎ」
な状況になる，
と (種多様性が高いほど現存量がどんどん減る
……
実際にそうなってる群落があっても不思議ではないけど)．
とすると，
バイオマス獲得の限界を種間競争だけではなく種内競争にも原因が，
と何らかの形で表現できればいいのか?
……
これは意味不明ぎみだな．
-
奇策は回避して，
相対的な種間競争に power を導入してみることに
……
これはよけーに事態を悪化させるとわかった．
-
いろいろ試行錯誤しつつ昼飯．
まずいなあ，
計算試行錯誤に時間をとられすぎてるよ
……
-
群落バイオマス分割，
ディリクレ分布
p[i, 1:N.sp] ~ ddirch(alpha[i,])
つかえばいいぢゃん，
と今さらながら気づいたんだが
……
unable to choose updater for node p[1,1]
なぜだー
……
dbern(...)
なんかとの組み合わせで無ければ動いてくれない，
ってコトかしらん?
共役事前分布でなきゃイヤイヤってのも横着なハナシだなぁ．
-
だったら自分で
ぢりくれ分布
つくればいいぢゃん，
と気づいたので，
alpha[i, j] ~ dgamma(mu.alpha[i, j], 1)
あれこれなどなどとしてみたんだが
cannot bracket slice for node alpha[12,11]
ゐんばぐす内の slice sampler がコケてるみたいだねえ
……
-
alpha[,]
に初期値 (1 とか) を与えたら計算するようにはなった
……
しかし 4000 MCMC step に 1165 秒も費やして
(slice sampler とかが苦闘してたんだろう)，
結果はやはり同じ，
と．
パラメーターとかの収束はいいのにねえ．
うーむ，
とはいえ，
やはり，
統計モデルがかなり良くない，
ってことかしらん?
-
で，
またまた今さらながら，
このコドラート内バイオマス「分割」モデルのよけーな計算に気づいた
……
「分割」だからって割算とかする必要ないぢゃん!
で，
憎むべき割算を取りさると計算は数倍速くなったんだが
(割算のせいで sampler によけーな心労・負担をかけていた，
ということだ)
……
うん，
列強談合体質はまだ完全には取り去れないね．
ちょー長時間計算すればいつかは一ヶ所にまとまるのかしらん?
-
どーしよーかなーとお茶部屋に逃げると，
そこにも湿原植物問題が．
といっても八甲田ではなくサロベツの湿原なんだけど．
江川さんから今シーズンの調査とその後について教えていただく
……
あいかわらずサロベツ的な状況
(A 棟 7-8F まわりでは「過酷・熾烈な重労働」を意味する)
のようで．
しかし昨年の春に播種した 3 種のうちヨシがめでたく全滅して，
残りはヌマガヤ・ミカヅキグサだけになったみたいだから，
解析したうけ業者としてはラクになったな．
-
1700 から研究院アワーに参加する．
気候力学分野の PD 奥西さんによる
「ニューラルネットワークを利用したマイワシの産卵回遊の再現」，
人工ニューラルネット (ANN) の脳をもつ
マイワシが北太平洋から土佐湾に帰ってくるように
人為選択 (genetic algorithm による)
かけてみました，
というもの．
わりとあらっぽい ANN でも行動を自動制御できたりするもんだな
……
-
1840 研究室発．
1855 帰宅．
晩飯の準備．
晩飯．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0820):
ヨーグルト．
リンゴ．
- 昼 (1330):
研究室お茶部屋．
食パン．
春雨スープ．
リンゴ．
- 晩 (2030):
米麦 0.7 合．
ジャガイモ・タマネギ・ニンジン・トマト・牛スネ肉のシチュー．
キャベツ・ニンジン・ネギのサラダ．