rt，刚才贴的程序不对，改了下：
model{
    yisigma2[1]<-1/exp(theta[1])
    ymean[1]<-beta0+d*exp(theta[1])
    y[1]~dnorm(ymean[1],yisigma2[1])
for(t in 2:n){
    yisigma2[t]<-1/exp(theta[t])
    ymean[t]<-beta0+beta1*(y[t-1]-beta0)+d*exp(theta[t])
    y[t]~dnorm(ymean[t],yisigma2[t])
   }
mu~dnorm(0,0.1)
phi~dbeta(20,1.5)
itau2 ~dgamma(2.5,0.025)
beta0~dnorm(0,1.0E-6)
beta1~dnorm(0,1.0E-6)
d~dnorm(0,1.0E-6)
#betal1<- 2*beta1star-1
tau<-sqrt(1/itau2)
theta0~dnorm(mu,itau2)
thmean[1]<-mu+phi*theta0
theta[1]~dnorm(thmean[1],itau2)
for(t in 2:n){
    thmean[t]<-mu+phi*(theta[t-1]-mu)
    theta[t]~dnorm(thmean[t],itau2)
    }
}
1.关于初始值：用winbugs做SV-M模型，总是提示“this chain contains unitialized varibles”。
程序如下：
我设置的初始值是：list(itau2=0.5, mu=-10, phi=0.975, d=6,beta0=0,beta1=0,theta0=0.1)

2.我的DIC值出不来，翻了之前的帖子，看来是程序没有运行掉进了trap，但程序验证说是语法没有问题。所以不知道了。
3.关于收敛性问题，因为我做的是单链，G_R统计量使用多链问题，那单链的收敛问题如何解决呢。

谢谢大家了。