ときどき原著(2007)もチェックした。充実した解説を公開されている方がいる。私も何箇所か参考にさせていただきました。私とは再現の範囲が違うし、ところどころ間違えているところもあると思うので(もちろん私の方も)、見比べていただければと思います。RプログラムRのプログラムはこちら。実行にあたっては、著者らが配布しているデータとプログラムに含まれている、以下のデータファイルを作業ディレクトリに置く必要がある。わかりやすくするために、すべての例題について同じ変数名・関数名をつかうように心がけた。以下、このように私が命名した変数名なり関数名を、ブラケットをつけて[oFiltered.DLM]というように表記する。著者らが配布しているssfpackプログラムも実行し、その結果と比べてみた。結局、各例題について、dlmの出力、KFASの出力、ssfpackの出力、本の記載の4つを比べることができるわけである。(←ここ、すごく悩んだところなので強調しておきます)。どのようにずれるかを調べるため、すべての例題について次の5つを比較した。CK本の(8.7)式にあてはめて得た最大対数尤度KFASが出力する最大対数尤度CK本に記載されている最大対数尤度。これは著者らが配布しているssfpackプログラムで出力される最大対数尤度と同じであることを、全例題について確認した。KFASによる推定を比較した。特記ない限り、ほぼ同じ結果が得られていた。(p.94のe_t)は、基本的にdlmとKFASの間で一致しているようだ。干渉変数のあるモデル(6章,7章)ではズレが生じる。