Question

我需要使用Winbugs估算四个参数：tau[i], mu1[i],mu2[i]and tau0[i]。

当我单击检查模型时，它在第四行的expected right parenthesis和tau0*之间显示(t[j+1]-t[j])。当我删除(t[j+1]-t[j])时，它在expected right parenthesis和c

之间显示相同的<-(1-step( tau[i]-t[j]))

model
{
#likelihood
   for(i in 1:n){   
    for(j in 1:m){      
 x[i,j]~ dnorm( mu[i,j],  tau0*(t[j+1]-t[j]))      
 mu[i,j]<- mu1*(t[j+1]-t[j])*a+mu1*(tau[i]-t[j])*b+mu2*(t[j+1]-tau[i])*b+mu2*(t[j+1]-t[j])*c       
 sigma <-1/sqrt(tau0)                                                                                       
 a<-step( tau[i]-t[ j+1] )
 b<-(step(tau[i]- t[j])-(step(tau[i]- t[j+1]))
 c<-(1-step( tau[i]-t[j]))
}
#prior
   tau[i]~dnorm(mutau,tautau)
   mutau~dorm(0,  0.001)
   tautau~dgamma(0.01, 0.01)
   sigmatau<-1/sqrt(tautau)
   mu1[i]~dnorm(mumu1, tau1)
   mumu1~dorm(0, 0.001)
   tau1~dgamma(0.01,0.01)
   sigmatau1<-1/sqrt(tau1)
   mu2[i]~dnorm(mumu2, tau2)
   mumu2~dorm(0, 0.001)
   tau2~dgamma(0.01,0.01)
   sigmatau2<-1/sqrt(tau2)

   tau0[i]~dgamma(alpha,beta)
   sigmatau0<-1/tau0
  }
}
DATA
list(n=5, m=16,   mutau=0, tautau=1,  mumu1=0, tau1=1,   mumu2=0, tau2=1, alpha=0.01, beta=0.01 )

我期望正确的模型语法

Answer 1

b <-（step（tau [i]-t [j + 1]）-（step（tau [i]-t [j + 1]））

(
    step(
            tau[i]- t[j]
        )
    -
    (
    step(
            tau[i]- t[j+1]
        )
    )

打开4个括号，关闭3个括号。第一括号未关闭。

如何解决在第四行中出现的“预期右括号”：winbugs中的tau0 *（t [j + 1] -t [j]）

1 个答案: