[quote]
引用第42楼baifei1106于2007-03-30 11:01发表的“”:
cran兄:
这是我编的GARCH(1,1)-N模型的WinBUGS程序,请你看看问题出在哪里
model
{
.......[/quote]
比较难明白,我没有做过这种模型.有一个可能,就是你的起始数不好,比如sigma不能是负数.
还有
<br />
y[t]~dnorm(0,sigma[t])<br />
要注意,BUGS里,dnorm的argument不是这么用的.第二个agument是1/var
所以你要写
<br />
prec[t] <- pow(sigma[t],2)<br />
y[t]~dnorm(0,prec[t])<br />