程序如下，红色加粗部分是需要修改的，感谢百忙之中帮我解答疑惑，或您发现了其他问题可私聊。

model
{
#likelihood
   for(i in 1:n){
          x ~ dnorm( mu, tau0 )
         mu<- mu1*t*I[0]+(mu1*tau-mu2*tau+mu2*t[i,j+1]-mu1*t[i,j])*I[1]+mu2*t*I[2]
         if tau> t[i, j+1]
            I[0]=1
         if t[i,j]<tau< t[i, j+1]
            I[1]=1
         if tau< t[i, j]
            I[2]=1
         sigma <-1/sqrt(tau0)
#prior
   tau~dnorm(mutau,tautau)
   mutau~dorm(0,  1.0E-3)
   tautau~dgamma(0.01, 0.01)
  sigmatau<-1/sqrt(tautau)
   }

   mu1~dnorm(mumu1, tau1)
   mumu1~dorm(0, 1.0E-3)
   tau1~dgamma(0.01,0.01)
   sigmatau1<-1/sqrt(tau1)

   mu2~dnorm(mumu2, tau2)
   mumu2~dorm(0, 1.0E-3)
   tau2~dgamma(0.01,0.01)
   sigmatau2<-1/sqrt(tau2)

   tau0~dgamma(alpha,beta)
   sigmatau0<-1/tau0

}