状態遷移の式 α(t+1)=T*αt+ε (記号はKFASから） でTの中に2種類のQが入るだけなので、何も問題がないはずです。 nn=3 a=rep(c(1,2,3),5*nn) b=rep(c(5,5,0,0,0),3*nn) d=10+a+b ts.plot(d[1:30]) scode=" data{ int N; real y[N]; } parameters{ real sigma;…