0)'、と考えるのだ。固定してしまう場合は気にしなくて良いが(7.1節とか)、そうでない限り、Qのサイズには十分に気を配った方が良さそうだ。nrow=2)と指定している。matrix(0)でもかまわないけど、こうしておいたほうが、Q行列のどの要素がどの状態変数なのかがわかりやすいと思ったからである。KFASが出力する最大対数尤度[gLLbyFun.KFAS]は、2,3章ではCK本の(8.7)式とぴったり一致するが、4章以降は少しずれる。理由はよくわからない。なにか事情があるんでしょうね。パラメータ推定について。fitSSM()に与えるinitsはoptim()にそのまま渡される。KFASのマニュアルを眺めていると、どうやら意図している初期値の対数を渡すべきらしいのだが、はっきりした説明が見当たらないので自信がない(KFASってこういうところが不親切なんですよね...)。与えるべき順番は、先頭が不規則要素、残りが状態撹乱要素、らしい。いくつかの例題では、fitSSM()がデフォルトのままだとエラーになる。method="BFGS"をつけるとうまくいく。7.1節なので、要するに状態変数がすべて確定的なときにこの問題が生じるようだ。なにか事情があるんでしょうね。coef()とrstandard()の存在に気づきました。4章以降ではこれらの関数を使っています。7章)、KFS()の返し値[oEstimated.KFAS]に含まれる予測誤差系列やその分散の系列は、それらしい値を持っている。