//このページが表示される方は、URLから「action=SOURCE&」を削除してみてください
[[状態空間モデリング - 記事一覧]]
!!!sspirサンプルコード
*投稿者： みゅ
*カテゴリ： なし
*優先度： 普通
*状態： 完了
*日時： 2008年06月12日 14時56分55秒

[sspirサンプルコード|http://r-project.matrix.jp/R/wiki.cgi?page=BugTrack%2D%BE%F5%C2%D6%B6%F5%B4%D6%A5%E2%A5%C7%A5%EA%A5%F3%A5%B0%2F1]

!!内容
RのsspirパッケージのSSオブジェクトのサンプルコード
サンプルコードだからいいんだけど、これで何がわかるんだろうという例・・・

 library(sspir)
 data(kurit)  ## West & Harrison, page 40
 m1 <- SS(y=kurit,
          Fmat=function(tt,x,phi) return(matrix(1)),
          Gmat=function(tt,x,phi) return(matrix(1)),
          Wmat=function(tt,x,phi) return(matrix(5)),
          Vmat=function(tt,x,phi) return(matrix(100)),
          m0=matrix(130),C0=matrix(400)
          )
 
 plot(m1$y)
 m1.f <- kfilter(m1)
 m1.s <- smoother(m1.f)
 lines(m1.f$m,lty=2,col=2)
{{ref_image sspir1.png}}
↑フィルターをおこなった系列のグラフ
 lines(m1.s$m,lty=2,col=2)
{{ref_image sspir2.png}}
↑平滑化までおこなった系列を描画

!!コメント
*あまり良いサンプルとはいえないような・・・ - みゅ (2008年06月12日 15時05分19秒)
//{{comment}}