小弟做了写了一个winbugs程序做跳跃利率模型的参数估计，哪位大哥能指点下是否正确？（程序可以运行）程序如下：

model
{
  for( i in 2 : n ) {

r[i] ~ dnorm(rmean[i],v)
rmean[i] <-  k1*mu1+(1-k1)*r[i-1]+j[i]*J[i]
j[i]~dbern(h)
J[i]~dnorm(muj,tauj)


}


k1 ~ dnorm( 0, 100)

mu1~ dnorm( 2,10)

v~dgamma(25,1)
h ~ dbeta(1,10)
tauj ~ dgamma(1,1)
sigma<-pow(v,(-0.5))
sigmaj<-pow(tauj,(-0.5))
muj~ dnorm( 0,100)


}