能不能具体你的问题，相关代码随便参考一本论坛上的书就可以了

您好，我已经在原帖中增加了模型的相关内容，请问哪本书中有相关的软件代码？

模型好复杂。我觉得你可以试试WINBUGS（或OPENBUGS，两者只有很小的差别，几乎一样使用，但好像WINBUGS在高版本的WINDOWS上装不了，那就用OPENBUGS），均为免费下载（GOOGLE就能找到）。专门对付MCMC，好像在R也可调用上述软件（需要装好上述软件后），可能可以帮到你。论坛也有WINBUGS板块。
不过，依我看：模型太复杂（时变的东东太多），可能程序会运行不出，或估计的结果无多大的实际意义。

就两个时变参数，一个是常数项，是个是π前面的参数么。