苫小牧直径成長モデル,
推定計算プログラム in R
の改造.
光合成気象値・代謝気象値を別 list にする推定プログラムで,
これらの組み合わせを同時に選抜するルーチンを追加した.
試験運転もかねて,
「ねんのため」
計算開始.
……
# AIC = 2210.97
# min AIC = 2071.37
# --- wv.mb set = 7/27 ---
# --- wv.ps set = 39/54 ---
temp.ps: best=16.0/width=6.0
ppfd: factor=0.02
vpd: factor=1.00
1998 1999 2000 2001 2002 2003 2004
26.65381 23.59560 22.93142 24.22508 25.75016 24.64961 25.81366
# AIC = 2218.24, np = 4+8 (scale.log const ps.prev mb6)
scale.log const ps.prev mb6
-0.53700289 1.22804230 0.01813444 0.15193820
# AIC = 2218.24
# min AIC = 2071.37
# --- wv.mb set = 7/27 ---
# --- wv.ps set = 40/54 ---
temp.ps: best=17.0/width=6.0
ppfd: factor=0.02
vpd: factor=1.00
1998 1999 2000 2001 2002 2003 2004
26.64290 23.88847 23.50196 25.19602 25.54211 23.64804 24.48687
……
てなかんぢで進行中.
27 × 54 = 1458 気象値セットしか調べてないのに,
たぶん 7 時間ぐらい必要.
次はこの「同時選抜の高速化」
をやらんといかんね
……